cardiac-geometriesx 0.4.1__tar.gz → 0.4.2__tar.gz

{cardiac_geometriesx-0.4.1/src/cardiac_geometriesx.egg-info → cardiac_geometriesx-0.4.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: cardiac-geometriesx
-Version: 0.4.1
+Version: 0.4.2
 Summary: A python library for cardiac geometries
 Author-email: Henrik Finsberg <henriknf@simula.no>
 License: MIT

{cardiac_geometriesx-0.4.1 → cardiac_geometriesx-0.4.2}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "cardiac-geometriesx"
-version = "0.4.1"
+version = "0.4.2"
 description = "A python library for cardiac geometries"
 authors = [{name = "Henrik Finsberg", email = "henriknf@simula.no"}]
 license = {text = "MIT"}
@@ -174,7 +174,7 @@ tag = true
 sign_tags = false
 tag_name = "v{new_version}"
 tag_message = "Bump version: {current_version} → {new_version}"
-current_version = "0.4.1"
+current_version = "0.4.2"
 
 
 [[tool.bumpversion.files]]

{cardiac_geometriesx-0.4.1 → cardiac_geometriesx-0.4.2}/src/cardiac_geometries/fibers/lv_ellipsoid.py

@@ -8,7 +8,12 @@ from . import utils
 
 
 def mu_theta(
-    x: np.ndarray, y: np.ndarray, z: np.ndarray, long_axis: int = 0
+    x: np.ndarray,
+    y: np.ndarray,
+    z: np.ndarray,
+    rs: np.ndarray,
+    rl: np.ndarray,
+    long_axis: int = 0,
 ) -> tuple[np.ndarray, np.ndarray, list[int]]:
     """Get the angles mu and theta from the coordinates x, y, z
     given the long axis.
@@ -21,6 +26,10 @@ def mu_theta(
         The y-coordinates
     z : np.ndarray
         The z-coordinates
+    rs : np.ndarray
+        The short radius
+    rl : np.ndarray
+        The long radius
     long_axis : int, optional
         The long axis, by default 0 (x-axis)
 
@@ -35,18 +44,18 @@ def mu_theta(
         If the long axis is not 0, 1 or 2
     """
     if long_axis == 0:
-        a = np.sqrt(y**2 + z**2)
-        b = x
+        a = np.sqrt(y**2 + z**2) / rs
+        b = x / rl
         theta = np.pi - np.arctan2(z, -y)
         perm = [0, 1, 2]
     elif long_axis == 1:
-        a = np.sqrt(x**2 + z**2)
-        b = y
+        a = np.sqrt(x**2 + z**2) / rs
+        b = y / rl
         theta = np.pi - np.arctan2(z, -x)
         perm = [1, 0, 2]
     elif long_axis == 2:
-        a = np.sqrt(x**2 + y**2)
-        b = z
+        a = np.sqrt(x**2 + y**2) / rs
+        b = z / rl
         theta = np.pi - np.arctan2(x, -y)
         perm = [2, 1, 0]
     else:
@@ -120,7 +129,7 @@ def compute_system(
     y = dof_coordinates[:, 1]
     z = dof_coordinates[:, 2]
 
-    mu, theta, perm = mu_theta(x, y, z, long_axis=long_axis)
+    mu, theta, perm = mu_theta(x, y, z, rs, rl, long_axis=long_axis)
 
     e_t = np.array(
         [

{cardiac_geometriesx-0.4.1 → cardiac_geometriesx-0.4.2}/src/cardiac_geometries/fibers/utils.py

@@ -25,6 +25,7 @@ def save_microstructure(
     if len(functions) == 0:
         return
     # Save for paraview visualization
+
     if functions[0].function_space.ufl_element().family_name == "quadrature":
         from scifem.xdmf import create_pointcloud
 
@@ -35,8 +36,8 @@ def save_microstructure(
                 mesh.comm, Path(outdir) / "microstructure-viz.bp", functions, engine="BP4"
             ) as file:
                 file.write(0.0)
-        except RuntimeError:
-            pass
+        except RuntimeError as ex:
+            print(f"Failed to write microstructure: {ex}")
 
     # Save with proper function space
     filename = Path(outdir) / "microstructure.bp"

{cardiac_geometriesx-0.4.1 → cardiac_geometriesx-0.4.2/src/cardiac_geometriesx.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: cardiac-geometriesx
-Version: 0.4.1
+Version: 0.4.2
 Summary: A python library for cardiac geometries
 Author-email: Henrik Finsberg <henriknf@simula.no>
 License: MIT