Merge pull request #110 from nschloe/test-coverage

Test coverage

Author: nschloe

pygalmesh/__init__.py

@@ -22,6 +22,7 @@
     Union,
 )
 
+from . import _cli
 from .__about__ import __version__
 from .main import (
     generate_from_array,
@@ -37,6 +38,8 @@
 __all__ = [
     "__version__",
     #
+    "_cli",
+    #
     "DomainBase",
     "SizingFieldBase",
     "Translate",

test/test_inr.py

@@ -1,15 +1,29 @@
-import os
+import pathlib
+import tempfile
 
 import helpers
+import meshio
 
 import pygalmesh
 
 
 def test_inr():
-    this_dir = os.path.dirname(os.path.abspath(__file__))
-    mesh = pygalmesh.generate_from_inr(
-        os.path.join(this_dir, "meshes", "skull_2.9.inr"), cell_size=5.0, verbose=False
-    )
+    this_dir = pathlib.Path(__file__).resolve().parent
+    # mesh = pygalmesh.generate_from_inr(
+    #     this_dir / "meshes" / "skull_2.9.inr", cell_size=5.0, verbose=False
+    # )
+    with tempfile.TemporaryDirectory() as tmp:
+        out_filename = str(pathlib.Path(tmp) / "out.vtk")
+        pygalmesh._cli.inr(
+            [
+                str(this_dir / "meshes" / "skull_2.9.inr"),
+                out_filename,
+                "--cell-size",
+                "5.0",
+                "--quiet",
+            ]
+        )
+        mesh = meshio.read(out_filename)
 
     tol = 2.0e-3
     ref = [2.031053e02, 3.739508e01, 2.425594e02, 2.558910e01, 2.300883e02, 1.775010e00]

test/test_remesh_surface.py

@@ -1,20 +1,40 @@
-import os
+import pathlib
+import tempfile
 
 import helpers
+import meshio
 
 import pygalmesh
 
 
 def test_inr():
-    this_dir = os.path.dirname(os.path.abspath(__file__))
-    mesh = pygalmesh.remesh_surface(
-        os.path.join(this_dir, "meshes", "lion-head.off"),
-        edge_size=0.025,
-        facet_angle=25,
-        facet_size=0.1,
-        facet_distance=0.001,
-        verbose=False,
-    )
+    this_dir = pathlib.Path(__file__).resolve().parent
+    # mesh = pygalmesh.remesh_surface(
+    #     this_dir / "meshes" / "lion-head.off",
+    #     edge_size=0.025,
+    #     facet_angle=25,
+    #     facet_size=0.1,
+    #     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,
+                "--edge-size",
+                "0.025",
+                "--facet-angle",
+                "25",
+                "--facet-size",
+                "0.1",
+                "--facet-distance",
+                "0.001",
+                "--quiet",
+            ]
+        )
+        mesh = meshio.read(out_filename)
 
     tol = 1.0e-3
     ref = [

test/test_volume_from_surface.py

@@ -1,20 +1,40 @@
-import os
+import pathlib
+import tempfile
 
 import helpers
+import meshio
 
 import pygalmesh
 
 
 def test_volume_from_surface():
-    this_dir = os.path.dirname(os.path.abspath(__file__))
-    mesh = pygalmesh.generate_volume_mesh_from_surface_mesh(
-        os.path.join(this_dir, "meshes", "elephant.vtu"),
-        facet_angle=25.0,
-        facet_size=0.15,
-        facet_distance=0.008,
-        cell_radius_edge_ratio=3.0,
-        verbose=False,
-    )
+    this_dir = pathlib.Path(__file__).resolve().parent
+    # mesh = pygalmesh.generate_volume_mesh_from_surface_mesh(
+    #     this_dir / "meshes" / "elephant.vtu",
+    #     facet_angle=25.0,
+    #     facet_size=0.15,
+    #     facet_distance=0.008,
+    #     cell_radius_edge_ratio=3.0,
+    #     verbose=False,
+    # )
+    with tempfile.TemporaryDirectory() as tmp:
+        out_filename = str(pathlib.Path(tmp) / "out.vtk")
+        pygalmesh._cli.volume_from_surface(
+            [
+                str(this_dir / "meshes" / "elephant.vtu"),
+                out_filename,
+                "--facet-angle",
+                "0.5",
+                "--facet-size",
+                "0.15",
+                "--facet-distance",
+                "0.008",
+                "--cell-radius-edge-ratio",
+                "3.0",
+                "--quiet",
+            ]
+        )
+        mesh = meshio.read(out_filename)
 
     tol = 2.0e-2
     assert abs(max(mesh.points[:, 0]) - 0.357612477657) < tol