forked from nschloe/meshio
-
Notifications
You must be signed in to change notification settings - Fork 0
/
test_flac3d.py
68 lines (59 loc) · 1.7 KB
/
test_flac3d.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
import copy
import pathlib
import sys
import helpers
import numpy
import pytest
import meshio
@pytest.mark.parametrize(
"mesh, binary, data",
[
(helpers.tet_mesh, False, []),
(helpers.hex_mesh, False, []),
(helpers.tet_mesh, False, [1, 2]),
(helpers.tet_mesh, True, []),
(helpers.hex_mesh, True, []),
(helpers.tet_mesh, True, [1, 2]),
],
)
def test(mesh, binary, data):
if data:
mesh = copy.deepcopy(mesh)
mesh.cell_data["flac3d:group"] = [numpy.array(data)]
helpers.write_read(
lambda f, m: meshio.flac3d.write(f, m, binary=binary),
meshio.flac3d.read,
mesh,
1.0e-15,
)
# the failure perhaps has to do with dictionary ordering
@pytest.mark.skipif(sys.version_info < (3, 6), reason="Fails with 3.5")
@pytest.mark.parametrize(
"filename",
["flac3d_mesh_ex.f3grid", "flac3d_mesh_ex_bin.f3grid"],
)
def test_reference_file(filename):
this_dir = pathlib.Path(__file__).resolve().parent
filename = this_dir / "meshes" / "flac3d" / filename
mesh = meshio.read(filename)
# points
assert numpy.isclose(mesh.points.sum(), 307.0)
# cells
ref_num_cells = [
("hexahedron", 45),
("pyramid", 9),
("hexahedron", 18),
("wedge", 9),
("hexahedron", 6),
("wedge", 3),
("hexahedron", 6),
("wedge", 3),
("pyramid", 6),
("tetra", 3),
("quad", 15),
("triangle", 3),
]
assert [(k, len(v)) for k, v in mesh.cells] == ref_num_cells
# Cell data
ref_sum_cell_data = [num_cell[1] for num_cell in ref_num_cells]
assert [len(arr) for arr in mesh.cell_data["flac3d:group"]] == ref_sum_cell_data