def test_msh(tmpdir): msh = str(tmpdir.join("magnsec.msh")) msh2 = str(tmpdir.join("magnsec2.msh")) convert.to_msh("tests/data/magnsec.ISA7", msh) isa = isa7.read("tests/data/magnsec.ISA7") convert.to_msh(isa, msh2) with open(msh) as f: file1 = f.read() with open(msh2) as f: file2 = f.read() assert file1 == file2 m = meshio.read(msh) assert len(m.points) == 2340 assert len(m.point_data["potential"]) == len(m.points) assert len(m.cells["triangle"]) == 1522 assert len(m.cells["quad"]) == 1506 assert len(m.field_data.keys()) == 28 num_cells = sum(map(len, m.cells.values())) for data in m.cell_data["triangle"]: len_data = sum([len(m.cell_data[ctype][data]) for ctype in m.cell_data]) assert len_data == num_cells
def test_nastran_quads(tmpdir): msh = str(tmpdir.join("quads.msh")) convert.to_msh("tests/data/quads.nas", msh) mesh = meshio.read(msh) assert len(mesh.points) == 6 assert mesh.cells[0][0] == "quad" and len(mesh.cells[0][1]) == 2
def test_nastran_triangles(tmpdir): msh = str(tmpdir.join("triangles.msh")) convert.to_msh("tests/data/triangles.nas", msh) mesh = meshio.read(msh) assert len(mesh.points) == 4 assert mesh.cells[0][0] == "triangle" and len(mesh.cells[0][1]) == 2
def test_nastran_superelements(tmpdir): msh = str(tmpdir.join("superelements.msh")) convert.to_msh("tests/data/superelements.nas", msh) mesh = meshio.read(msh) assert len(mesh.points) == 6 assert mesh.cells[0][0] == "triangle" and len(mesh.cells[0][1]) == 2 assert mesh.cells[1][0] == "quad" and len(mesh.cells[1][1]) == 1 assert mesh.cell_data["gmsh:geometrical"][0][0] == 2 assert mesh.cell_data["gmsh:geometrical"][0][1] == 1 assert mesh.cell_data["gmsh:geometrical"][1][0] == 2
def test_msh(tmpdir): msh = str(tmpdir.join("magnsec.msh")) msh2 = str(tmpdir.join("magnsec2.msh")) convert.to_msh("tests/data/magnsec.ISA7", msh) isa = isa7.read("tests/data/magnsec.ISA7") convert.to_msh(isa, msh2) with open(msh) as f: file1 = f.read() with open(msh2) as f: file2 = f.read() assert file1 == file2 m = meshio.read(msh) assert len(m.points) == 2340 assert len(m.point_data["potential"]) == len(m.points) assert len(m.cells[1][1]) == 1522 assert len(m.cells[2][1]) == 1506 assert len(m.field_data.keys()) == 28
def test_msh_quads(tmpdir): msh = str(tmpdir.join("quads.msh")) convert.to_msh("tests/data/quads.ISA7", msh)
def test_msh_triangles(tmpdir): msh = str(tmpdir.join("triangles.msh")) convert.to_msh("tests/data/triangles.ISA7", msh)