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)