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
|
import pathlib
import tempfile
import helpers
import meshio
import pygalmesh
def test_remesh_surface():
this_dir = pathlib.Path(__file__).resolve().parent
# mesh = pygalmesh.remesh_surface(
# this_dir / "meshes" / "lion-head.off",
# max_edge_size_at_feature_edges=0.025,
# min_facet_angle=25,
# max_radius_surface_delaunay_ball=0.1,
# max_facet_distance=0.001,
# verbose=False,
# )
with tempfile.TemporaryDirectory() as tmp:
out_filename = str(pathlib.Path(tmp) / "out.vtk")
pygalmesh._cli.remesh_surface(
[
str(this_dir / "meshes" / "lion-head.off"),
out_filename,
"--max-edge-size-at-feature-edges",
"0.025",
"--min-facet-angle",
"25",
"--max-radius-surface-delaunay-ball",
"0.1",
"--max-facet-distance",
"0.001",
"--quiet",
]
)
mesh = meshio.read(out_filename)
tol = 1.0e-3
ref = [
3.705640e-01,
-3.711630e-01,
4.754940e-01,
-4.728260e-01,
4.998690e-01,
-4.998350e-01,
]
assert abs(max(mesh.points[:, 0]) - ref[0]) < tol * abs(ref[0])
assert abs(min(mesh.points[:, 0]) - ref[1]) < tol * abs(ref[1])
assert abs(max(mesh.points[:, 1]) - ref[2]) < tol * abs(ref[2])
assert abs(min(mesh.points[:, 1]) - ref[3]) < tol * abs(ref[3])
assert abs(max(mesh.points[:, 2]) - ref[4]) < tol * abs(ref[4])
assert abs(min(mesh.points[:, 2]) - ref[5]) < tol * abs(ref[5])
print(helpers.compute_triangle_areas(mesh.points, mesh.get_cells_type("triangle")))
vol = sum(
helpers.compute_triangle_areas(mesh.points, mesh.get_cells_type("triangle"))
)
ref = 1.917942005786831
assert abs(vol - ref) < ref * 1.0e-3
|
