search

finsberg / ldrb

A software for assigning myocardial fiber orientations based on the Laplace Dirichlet Ruled-Based algorithm
http://finsberg.github.io/ldrb
GNU Lesser General Public License v3.0
22 stars 8 forks source link

Fiber generation problem #88

Closed wang-hairui closed 1 year ago

wang-hairui commented 1 year ago

Hi，Dr. Finsberg!

​ Thank you very much for providing the Python implementation of LDRB. I have a question and would like to consult you.

​ When running the example "Life X - idealized_LV" (https://finsberg.github.io/ldrb/demo_lifex_idealized_LV.html), I noticed that the number of generated fiber vectors corresponds to the number of points in the mesh (read from the stored fiber.xdmf using meshio, where fibers are also stored in ‘point_data’). Can I understand that the algorithm you provided assigns a fiber direction to each point in the mesh?

​ If I want to assign a specific fiber direction to each tetrahedron in the mesh, how should I proceed? Do you have any suggestions or examples I can refer to?

finsberg commented 1 year ago

If you use the fiber space "P_1" it will assign a fiber direction to each point in the mesh. To assign one fiber direction to each tetrahedron, you need to use "DG_0" space

wang-hairui commented 1 year ago

Thank you very much for your response！！This answer solved my problem. May you provide some links or documents for me to further understand the other available options for the "fiber_space"?

finsberg commented 1 year ago

The "fiber_space" is basically a string on the form "{family}_{degree}", where family is the finite element space and degree is the degree. I guess good resources for this are https://defelement.com and https://femtable.org.

The relevant code is found here: https://github.com/finsberg/ldrb/blob/f12e5a6cb769bc51d4aa4ead68c32b8a3db61777/ldrb/utils.py#L188C2-L188C2

wang-hairui commented 1 year ago

I will carefully read these references. Thanks again for your response !!!