adding fiber generation codes

javijv4 · javijv4 · commit b010f79c0701 · 2025-12-05T11:29:56.000-08:00
diff --git a/utilities/fiber_generation/README.md b/utilities/fiber_generation/README.md
@@ -0,0 +1,35 @@
+# SV-fibergen
+Python + SVmultiphysics codes for fiber generation. Two methods are implemented:
+* Bayer et al. (2012). [link](https://doi.org/10.1007/s10439-012-0593-5)
+* Doste et al. (2018). [link](https://doi.org/10.1002/cnm.3185)
+
+### Examples
+The `main_bayer.py` and `main_doste.py` are scripts to run both methods in the geometry described in the `example/truncated` and `example/ot` folders respectively.
+
+<img src="example/truncated/bayer.png" alt="Results for truncated BiV (Bayer)" width="640" />
+<img src="example/ot/doste.png" alt="Results for BiV w/ outflow tracts (Doste)" width="640" />
+
+Note that the Doste methods needs a geometry with outflow tracts to be run (each valve needs to be defined as a separated surface). Bayer can be run in any biventricular geometry.
+
+
+### Updates to the old code
+* All operations are vectorized now.
+* The SVmultiphysics solver now solves a Laplace equation.
+* In Bayer: For the bislerp interpolation, instead of using the correction described in Bayer et al. (that returns a discontinuity), the basis are flipped to maintain a coherent fiber direction (see function `generate_fibers_BiV_Bayer_cells` in `FibGen.py`).
+* In Bayer: The beta angles were not being included correctly. The second rotation was being applied respect the first vector (circumferential) when it should be respect to the second vector (longitudinal) (see function `generate_fibers_BiV_Bayer_cells` in `FibGen.py`).
+
+### Notes on SVmultiphysics solver
+
+To solve a Laplace equation directly from the transient HEAT solver in SVmultiphysics, in `<GeneralSimulationParameters>` we need to set,
+```
+<Number_of_time_steps> 1 </Number_of_time_steps>
+<Time_step_size> 1 </Time_step_size>
+<Spectral_radius_of_infinite_time_step> 0. </Spectral_radius_of_infinite_time_step>
+```
+and in `<Add_equation type="heatS" >`,
+```
+<Conductivity> 1.0 </Conductivity>
+<Source_term> 0.0 </Source_term>
+<Density> 0.0 </Density>
+```
+This will allow us to solve the Laplace equation directly in 1 timestep and 1 iteration.
diff --git a/utilities/fiber_generation/example/ot/doste.png b/utilities/fiber_generation/example/ot/doste.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4da5c28acd2e2b9745ba972ec783b7341afa2156a48244059684ca28f45652ac
+size 494459
diff --git a/utilities/fiber_generation/example/ot/mesh-complete.mesh.vtu b/utilities/fiber_generation/example/ot/mesh-complete.mesh.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:9347d6e7907816749cc16471130397f264c7721d029626fb1b9f7fdd6ac1e02c
+size 28209696
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/EPI_APEX.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/EPI_APEX.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:1431a90e42af727af648dacde59095e9c3fd8e37d4a7dfef3fef57b633eaf139
+size 2867
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/av.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/av.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:ee7d4d8eb31b2639f75599c188cb2af9fa4fcc466d27e34b65dc3ec4f4f3633b
+size 14489
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/endo_lv.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/endo_lv.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:fa09ea16a80da5813b133a22f42b94b7aae48cd414e9c8ce174d46a4d7eb3a19
+size 563261
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/endo_rv.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/endo_rv.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b0c07b7ca1ad1b8c38d826d363f5572e1e6359595f5701a2b3f8f25606e30073
+size 734382
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/epi.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/epi.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:2168042f6f99283d747e5895df8bb30d443631b5be7e31b51c55e1a7be83fa29
+size 1293881
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/mv.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/mv.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:712904ededfdd630a736abfc213758155f8dd68a7f1512c2128c2b0305eeecb6
+size 90751
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/pv.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/pv.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:1b99698dc66260cd1c63d4612402f07e98b0a547549e42d896048bf7bdbe2afb
+size 12559
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/top.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/top.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:edec0ae0c98ef7db7add744db0340aeba1c925613403ec4bc1e580d0c1454767
+size 87143
diff --git a/utilities/fiber_generation/example/ot/mesh-surfaces/tv.vtp b/utilities/fiber_generation/example/ot/mesh-surfaces/tv.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b33291fd2b529845d68d4584c2f4f5241c288dda1eb6a5e76f5bc839e21e2f62
+size 19198
diff --git a/utilities/fiber_generation/example/truncated/VOLUME.vtu b/utilities/fiber_generation/example/truncated/VOLUME.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:bcd011125435852298177578498ebcf2b5aac1080b824163d1facdefc48c2a27
+size 24524913
diff --git a/utilities/fiber_generation/example/truncated/bayer.png b/utilities/fiber_generation/example/truncated/bayer.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c01ec2a92c75b6047fa9b6602924a619d1e20b577cc26de7852a91cc06d19b6a
+size 313768
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/BASE.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/BASE.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b56eb592388c9cc2226c6372ffc8e0199952698f84ffa81330672d272503d7aa
+size 191256
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:2c5ecc05750d086d04dfe9902e7a7e200831ed5cb11ffe6339186112d8bac0ce
+size 526632
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI_APEX.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI_APEX.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c6e68548080a65b62e5ea62af5d685cafa9181604f43ae45fa9680b0f6e09938
+size 2833
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI_APEX2.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI_APEX2.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c6e68548080a65b62e5ea62af5d685cafa9181604f43ae45fa9680b0f6e09938
+size 2833
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI_MID.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/EPI_MID.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:cc6a73358eff58939c32e73dfd9326b9a53cf25afa409e5c27b7a87ef3ec08ad
+size 525056
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/LV.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/LV.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e16a121071598c4696ea25623b52595c349bca27823c6d9cfb6e1bdbc7b77767
+size 145544
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/RV.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/RV.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c6028ab47b1af4947b9ee47cd79dfc58105d8e2bf86f1167aed877609023935c
+size 128852
diff --git a/utilities/fiber_generation/example/truncated/mesh-surfaces/exterior.vtp b/utilities/fiber_generation/example/truncated/mesh-surfaces/exterior.vtp
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:889f15170752b6310d35a52f9e1b207ab84afa1e3ce7c4be24c114403015fe3f
+size 978206
diff --git a/utilities/fiber_generation/main_bayer.py b/utilities/fiber_generation/main_bayer.py
@@ -0,0 +1,64 @@
+#!/usr/bin/env python
+# -*-coding:utf-8 -*-
+'''
+Created on 2025/11/21 20:38:14
+
+@author: Javiera Jilberto Vallejos 
+'''
+
+import os
+import src.FibGen as fg
+from time import time
+
+###########################################################
+############  USER INPUTS  ################################
+###########################################################
+
+run_flag = True
+svfsi_exec = "svmultiphysics "
+
+mesh_path = "example/truncated/VOLUME.vtu"
+surfaces_dir = f"example/truncated/mesh-surfaces"
+outdir = "example/truncated/output_b"
+
+surface_names = {'epi': 'EPI.vtp',
+                 'epi_apex': 'EPI_APEX.vtp',    # New surface
+                 'base': 'BASE.vtp',
+                 'endo_lv': 'LV.vtp',
+                 'endo_rv': 'RV.vtp'}
+
+# Parameters for the Bayer et al. method https://doi.org/10.1007/s10439-012-0593-5. 
+params = {
+    "ALFA_END": 60.0,
+    "ALFA_EPI": -60.0,
+    "BETA_END": 20.0,
+    "BETA_EPI": -20.0,
+}
+
+
+###########################################################
+############  FIBER GENERATION  ###########################
+###########################################################
+
+# Make sure the paths are full paths
+mesh_path = os.path.abspath(mesh_path)
+surfaces_dir = os.path.abspath(surfaces_dir)
+outdir = os.path.abspath(outdir)
+
+start = time()
+fg.generate_epi_apex(mesh_path, surfaces_dir, surface_names)
+
+# Run the Laplace solver
+if run_flag:
+    template_file = os.path.join(os.path.dirname(os.path.abspath(__file__)), "src", "templates", "solver_bayer.xml")
+    laplace_results_file = fg.runLaplaceSolver(mesh_path, surfaces_dir, mesh_path, svfsi_exec, template_file, outdir, surface_names)
+laplace_results_file = outdir + '/result_001.vtu'
+
+# Generate the fiber directions
+result_mesh = fg.generate_fibers_BiV_Bayer_cells(outdir, laplace_results_file, params, return_angles=True, return_intermediate=True)
+
+print(f"generate fibers (Bayer method) elapsed time: {time() - start:.3f} s")
+
+# Optional, save the result mesh with intermediate field and angles for checking
+result_mesh_path = os.path.join(outdir, "results_bayer.vtu")
+result_mesh.save(result_mesh_path)
diff --git a/utilities/fiber_generation/main_doste.py b/utilities/fiber_generation/main_doste.py
@@ -0,0 +1,85 @@
+#!/usr/bin/env python
+# -*-coding:utf-8 -*-
+'''
+Created on 2025/11/21 20:38:14
+
+@author: Javiera Jilberto Vallejos 
+'''
+
+import os
+import src.FibGen as fg
+from time import time
+
+###########################################################
+############  USER INPUTS  ################################
+###########################################################
+
+run_flag = True
+method = 'doste'
+svfsi_exec = "svmultiphysics "
+
+mesh_path = "example/ot/mesh-complete.mesh.vtu"
+surfaces_dir = f"example/ot/mesh-surfaces"
+outdir = "example/ot/output_d"
+
+surface_names = {'epi': 'epi.vtp',
+                 'epi_apex': 'epi_apex.vtp',    # New surface
+                 'av': 'av.vtp',
+                 'mv': 'mv.vtp',
+                 'tv': 'tv.vtp',
+                 'pv': 'pv.vtp',
+                 'base': 'top.vtp',             # This is all the valves together, it is used to find the apex.
+                 'endo_lv': 'endo_lv.vtp',
+                 'endo_rv': 'endo_rv.vtp'}
+
+# Parameters from the Doste paper https://doi.org/10.1002/cnm.3185
+params = {
+    # A = alpha angle
+    'AENDORV' : 90,
+    'AEPIRV' : -25,
+    'AENDOLV' : 60,
+    'AEPILV' : -60,
+
+    'AOTENDOLV' : 90, 
+    'AOTENDORV' : 90,
+    'AOTEPILV' : 0,
+    'AOTEPIRV' : 0,
+
+    # B = beta angle (this have an opposite sign to the Doste paper, 
+    # but it's because the longitudinal direction is opposite)
+    'BENDORV' : 20,
+    'BEPIRV' : -20,
+    'BENDOLV' : 20,
+    'BEPILV' : -20,
+}
+
+
+###########################################################
+############  FIBER GENERATION  ###########################
+###########################################################
+
+# Make sure the paths are full paths
+mesh_path = os.path.abspath(mesh_path)
+surfaces_dir = os.path.abspath(surfaces_dir)
+outdir = os.path.abspath(outdir)
+
+# Generate the apex surface
+start = time()
+
+start = time()
+fg.generate_epi_apex(mesh_path, surfaces_dir, surface_names)
+
+# Run the Laplace solver
+if run_flag:
+    template_file = os.path.join(os.path.dirname(os.path.abspath(__file__)), "src", "templates", "solver_doste.xml")
+    laplace_results_file = fg.runLaplaceSolver(mesh_path, surfaces_dir, mesh_path, svfsi_exec, template_file, outdir, surface_names)
+laplace_results_file = outdir + '/result_001.vtu'
+
+# Generate the fiber directions
+result_mesh = fg.generate_fibers_BiV_Doste_cells(outdir, laplace_results_file, params, return_angles=True, return_intermediate=False)
+
+print(f"generate fibers (Doste method) elapsed time: {time() - start:.3f} s")
+
+# Optional, save the result mesh with intermediate field and angles for checking
+result_mesh_path = os.path.join(outdir, "results_doste.vtu")
+result_mesh.save(result_mesh_path)
diff --git a/utilities/fiber_generation/src/FibGen.py b/utilities/fiber_generation/src/FibGen.py
diff --git a/utilities/fiber_generation/src/templates/solver_bayer.xml b/utilities/fiber_generation/src/templates/solver_bayer.xml
diff --git a/utilities/fiber_generation/src/templates/solver_doste.xml b/utilities/fiber_generation/src/templates/solver_doste.xml

Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:9347d6e7907816749cc16471130397f264c7721d029626fb1b9f7fdd6ac1e02c`
	`3`	`+size 28209696`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:1431a90e42af727af648dacde59095e9c3fd8e37d4a7dfef3fef57b633eaf139`
	`3`	`+size 2867`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:ee7d4d8eb31b2639f75599c188cb2af9fa4fcc466d27e34b65dc3ec4f4f3633b`
	`3`	`+size 14489`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:fa09ea16a80da5813b133a22f42b94b7aae48cd414e9c8ce174d46a4d7eb3a19`
	`3`	`+size 563261`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:b0c07b7ca1ad1b8c38d826d363f5572e1e6359595f5701a2b3f8f25606e30073`
	`3`	`+size 734382`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:2168042f6f99283d747e5895df8bb30d443631b5be7e31b51c55e1a7be83fa29`
	`3`	`+size 1293881`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:712904ededfdd630a736abfc213758155f8dd68a7f1512c2128c2b0305eeecb6`
	`3`	`+size 90751`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,3 @@`
	`1`	`+version https://git-lfs.github.com/spec/v1`
	`2`	`+oid sha256:1b99698dc66260cd1c63d4612402f07e98b0a547549e42d896048bf7bdbe2afb`
	`3`	`+size 12559`