How to specify cell type in MeshEditor().open()?

Question

When I try to open a MeshEditor, using

editor.open( mesh, mesh.type(), mesh.topology().dim(), mesh.geometry().dim() )

I get the following TypeError:

in method 'MeshEditor_open', argument 3 of type 'std::string'

From inspecting the documentation, I would expect this to work. It works when I replace the second argument with "triangle" for 2d meshes. Am I missing something? If not, what's the best way of generically specifying the cell type for both 2d and 3d meshes?

Answer 1

Hi, in cpp there are three versions of MeshEditor.open. On python side, the function taking 5 arguments whose signature involves CellType is ignored. You get your type error because the remaining function taking five arguments expects string. If you insist on using that signature consider

from dolfin import *

mesh0 = UnitSquareMesh(2, 2)

mesh = Mesh()
editor = MeshEditor()

gdim = mesh0.geometry().dim()
tdim = mesh0.topology().dim()
c_type = mesh0.type()
c_str = c_type.type2string(c_type.cell_type())

editor.open(mesh, c_str, tdim, gdim)
editor.close()