00009 :
00010 tdim = cell2dim[celltype]
00011 gdim = tdim
00012 mesh = Mesh()
00013 editor = MeshEditor()
00014 editor.open(mesh, celltype, tdim, gdim)
00015 if celltype == "interval":
00016 vertices = [(0.0,),
00017 (1.0,)]
00018 if celltype == "triangle":
00019 vertices = [(0.0, 0.0),
00020 (1.0, 0.0),
00021 (0.0, 1.0)]
00022 if celltype == "tetrahedron":
00023 vertices = [(0.0, 0.0, 0.0),
00024 (1.0, 0.0, 0.0),
00025 (0.0, 1.0, 0.0),
00026 (0.0, 0.0, 1.0)]
00027 if celltype == "quadrilateral":
00028 vertices = [(0.0, 0.0),
00029 (1.0, 0.0),
00030 (1.0, 1.0),
00031 (0.0, 1.0)]
00032 if celltype == "hexahedron":
00033 vertices = [(0.0, 0.0, 0.0),
00034 (1.0, 0.0, 0.0),
00035 (1.0, 1.0, 0.0),
00036 (0.0, 1.0, 0.0),
00037 (0.0, 0.0, 1.0),
00038 (1.0, 0.0, 1.0),
00039 (1.0, 1.0, 1.0),
00040 (0.0, 1.0, 1.0)]
00041 editor.initVertices(len(vertices))
00042 editor.initCells(1)
00043 for i, p in enumerate(vertices):
00044 editor.addVertex(i, *p)
00045 editor.addCell(0, *range(len(vertices)))
00046 editor.close()
00047 return mesh
00048