def create_example_mesh(): points = array( [ [0, 0, 0], [1, 0, 0], [0, 1, 0], [0, 0, 1], [2, 0, 0], [3, 0, 0], [3, 1, 0], [2, 1, 0], [2, 0, 1], [3, 0, 1], [3, 1, 1], [2, 1, 1], ], "f", ) cells = [[0, 1, 2, 3], [4, 5, 6, 7, 8, 9, 10, 11]] # tetra # hex faces = [[2, 7, 11]] edges = [[1, 4], [3, 8]] mesh = Mesh("example") # add points uids = mesh.add_points( Point(coordinates=point, data=DataContainer(TEMPERATURE=index)) for index, point in enumerate(points) ) # add edges mesh.add_edges( Edge(points=[uids[index] for index in element], data=DataContainer(TEMPERATURE=index)) for index, element in enumerate(edges) ) # add faces mesh.add_faces( Face(points=[uids[index] for index in element], data=DataContainer(TEMPERATURE=index)) for index, element in enumerate(faces) ) # add cells mesh.add_cells( Cell(points=[uids[index] for index in element], data=DataContainer(TEMPERATURE=index)) for index, element in enumerate(cells) ) return mesh
[3, 1, 1], [2, 1, 1], ], "f", ) cells = [[0, 1, 2, 3], [4, 5, 6, 7, 8, 9, 10, 11]] # tetra # hex faces = [[2, 7, 11]] edges = [[1, 4], [3, 8]] mesh = Mesh("example") # add points uids = mesh.add_points( Point(coordinates=point, data=DataContainer(TEMPERATURE=index)) for index, point in enumerate(points) ) # add edges edge_uids = mesh.add_edges(Edge(points=[uids[index] for index in element]) for index, element in enumerate(edges)) # add faces face_uids = mesh.add_faces(Face(points=[uids[index] for index in element]) for index, element in enumerate(faces)) # add cells cell_uids = mesh.add_cells(Cell(points=[uids[index] for index in element]) for index, element in enumerate(cells)) if __name__ == "__main__": from simphony.visualisation import paraview_tools # Visualise the Mesh object