Merge pull request #131 from tataratat/no-splines

move spline functionalities to splinepy

examples/.gitignore

@@ -0,0 +1 @@
+samples

examples/load_sample_file.py

@@ -0,0 +1,34 @@
+"""Helper file for examples. Loads sample files from the tataratat/samples."""
+import pathlib
+
+import requests
+
+base_url = "https://raw.githubusercontent.com/tataratat/samples/main/"
+local_path = pathlib.Path(__file__).parent / "samples/"
+
+
+def load_sample_file(filename: str, force_reload: bool = False) -> bool:
+    """Loads a sample file from the tataratat/samples repository.
+
+    Args:
+        filename (str): Path of the file to load.
+        force_reload (bool, optional): Loads the file even if it already
+        exists. Defaults to False.
+
+    Returns:
+        bool: File could be loaded.
+    """
+    local_file = local_path / filename
+    if local_file.exists() and not force_reload:
+        return True
+    pathlib.Path.mkdir(local_file.parent, parents=True, exist_ok=True)
+
+    url = base_url + filename
+    response = requests.get(url)
+    if response.status_code == 200:
+        with open(local_file, "wb") as f:
+            f.write(response.content)
+        return True
+    else:
+        print(f"Cant load {url}, status code: {response.status_code}.")
+        return False

examples/show_gmsh.py

@@ -5,31 +5,31 @@
 """
 import pathlib
 
+import load_sample_file
+
 from gustaf import io
 
 if __name__ == "__main__":
-    base_samples_path = pathlib.Path("samples/faces")
+    mesh_file_tri = pathlib.Path("faces/tri/2DChannelTria.msh")
+    mesh_file_quad = pathlib.Path("faces/quad/2DChannelQuad.msh")
+
+    base_samples_path = pathlib.Path("samples")
     if not base_samples_path.exists():
-        raise RuntimeError(
-            "The geometries could not be found. Please initialize the "
-            "samples submodule, instructions can be found in the "
-            "README.md."
-        )
-    mesh_file_tri = pathlib.Path("samples/faces/tri/2DChannelTria.msh")
-    mesh_file_quad = pathlib.Path("samples/faces/quad/2DChannelQuad.msh")
+        load_sample_file.load_sample_file(str(mesh_file_tri))
+        load_sample_file.load_sample_file(str(mesh_file_quad))
 
     # load the .msh file directly with the correct io module (meshio)
-    loaded_mesh_tri = io.meshio.load(mesh_file_tri)
+    loaded_mesh_tri = io.meshio.load(base_samples_path / mesh_file_tri)
 
     loaded_mesh_tri.show()
 
     # load the .msh file directly with the correct io module (meshio)
-    loaded_mesh_quad = io.meshio.load(mesh_file_quad)
+    loaded_mesh_quad = io.meshio.load(base_samples_path / mesh_file_quad)
 
     loaded_mesh_quad.show()
 
     # load the .msh file with the default load function which needs to find out
     # it self which module is the correct one.
-    loaded_mesh_default = io.load(mesh_file_tri)
+    loaded_mesh_default = io.load(base_samples_path / mesh_file_tri)
 
     loaded_mesh_default.show()