Create geometry if not exists

demos/simple_pvloop.py

@@ -4,12 +4,13 @@
 # Local imports
 import logging
 import dolfin
+from pathlib import Path
 
 # from pulse2.problem import Problem
 from pulse2.itertarget import itertarget
 import pulse2
 from pulse2.geometry import LVGeometry
-from cardiac_geometries.geometry import Geometry
+import create_geometries
 
 # from pulse2.material import HolzapfelOgden, Guccione
 
@@ -22,8 +23,14 @@
 comm = dolfin.MPI.comm_world
 
 # Read geometry from file. If the file is not present we regenerate it.
+geofolder = Path("lv")
+if not geofolder.is_dir():
+    create_geometries.create_lv_ellipsoid(  # type: ignore
+        geofolder,
+        create_fibers=True,
+    )
 
-geo = Geometry.from_folder("lv")
+geo = create_geometries.geometry.Geometry.from_folder(geofolder)  # type: ignore
 geo = LVGeometry(
     mesh=geo.mesh,
     markers=geo.markers,

pyproject.toml

@@ -41,17 +41,19 @@ pypi = [
 ]
 demo = [
     "cardiac-geometries",
+    "ldrb"
 ]
 docs = [
    "jupyter-book",
    "jupytext",
-   "cardiac-geometries"
+   "pulse2[demo]"
 ]
 all = [
-   "pulse[test]",
-   "pulse[docs]",
-   "pulse[pypi]",
-   "pulse[dev]"
+   "pulse2[test]",
+   "pulse2[docs]",
+   "pulse2[pypi]",
+   "pulse2[demo]",
+   "pulse2[dev]"
 ]