def test_triangle_mesh(): mesh = Mesh( MPI.comm_world, CellType.triangle, numpy.array([[0.0, 0.0], [1.0, 0.0], [0.0, 1.0]], dtype=numpy.float64), numpy.array([[0, 1, 2]], dtype=numpy.int32), [], cpp.mesh.GhostMode.none) assert mesh.num_entities_global(0) == 3 assert mesh.num_entities_global(2) == 1
def test_read_mesh_data(tempdir, tdim, n): filename = os.path.join(tempdir, "mesh.xdmf") mesh = mesh_factory(tdim, n) encoding = XDMFFile.Encoding.HDF5 ghost_mode = cpp.mesh.GhostMode.none with XDMFFile(mesh.mpi_comm(), filename, encoding) as file: file.write(mesh) with XDMFFile(MPI.comm_world, filename) as file: cell_type, points, cells, indices = file.read_mesh_data(MPI.comm_world) mesh2 = Mesh(MPI.comm_world, cell_type, points, cells, indices, ghost_mode) assert (mesh.topology.cell_type == mesh2.topology.cell_type) assert mesh.num_entities_global(0) == mesh2.num_entities_global(0) dim = mesh.topology.dim assert mesh.num_entities_global(dim) == mesh2.num_entities_global(dim)