search

rsachetto / MonoAlg3D_C

Extensible Monodomain simulator
MIT License
16 stars 13 forks source link

Updating Purkinje-library - Package 8 + Bug Fix #19

Closed bergolho closed 4 years ago

bergolho commented 4 years ago

Hello Sachetto,

In this pull request, a multiple include error that was happening in the "grid.c" was fixed. This issue was happening when you recompile the whole project from scrath or when you try to compile for the first time.

In addition, I added some simulations in order to get the steady-state vector of any given cellular model and a script that can plot all the state variables. Using these new features I was able to update the initial conditions with the steady-state vector for the ToRORd and Trovato models. There is also a new function that can calculate the propagation velocity across a cable Purkinje fiber. I still need to improve this function to work on any given Purkinje network, but for that, I will need to implement a queue data structure, so it can be used on the Dijkstra algorithm.

The next steps are to allow the Purkinje-coupling examples to be executed using batch simulations and implement the Dijkstra algorithm to calculate the propagation velocity.

rsachetto commented 4 years ago

I dont think that the print_purkinje_propagation_velocity function should be in the mondomain_solver.c file. Maybe in the save_library?

bergolho commented 4 years ago

You are right. The 'print_purkinje_propagation_velocity' only works when you have the activation times from the cells and it makes more sense to put the function alongside the functions that write the LAT and APD map in the 'save_mesh_helper.c'. This way, it will be possible to turn on/off the propagation velocity calculus for the Purkinje in the configuration file.