search

rsachetto / MonoAlg3D_C

Extensible Monodomain simulator
MIT License
16 stars 12 forks source link

Oxford models fixes and extra features #54

Closed bergolho closed 1 year ago

bergolho commented 1 year ago

Hello Sachetto,

In this pull request, I made a few changes to the Oxford cellular models ("Trovato" and "ToRORd_dynCl") to avoid NaN errors that arise when a Purkinje is coupled to the tissue in some corner cases. Also, I updated the initial conditions from these models as well.

In addition, there is a new "matrix_assembly" function that considers a fast endocardium layer for only ventricular tissue simulations. Previously, there was only a Purkinje coupled version for this.

Furthermore, a new "save_mesh" function to save the transmurality of a given mesh in VTU format is available. I made a new function, "new_vtk_unstructured_grid_from_alg_grid_using_values" to allow any float array to be passed as the CELL_DATA to be written in the VTU file.

Finally, there are some new scripts for trace plots.

Best regards, Lucas

bergolho commented 1 year ago

Code duplication avoided using the suggested if statement in the "new_vtk_unstructured_grid_from_alg_grid" function.