Add a new version of the 'calc_apd_full_grid' script using Elnaz's al…

…gorithm
bergolho · Mar 31, 2019 · 2c70498 · 2c70498
1 parent c499bb4
commit 2c70498
Show file tree

Hide file tree

Showing 11 changed files with 11,133 additions and 10,603 deletions.
diff --git a/ToDo b/ToDo
@@ -7,12 +7,10 @@ Sachetto
 - Rewrite some library APIs (maybe we need only the config to be passed to the library function)
 
 Berg
-- Rewrite the 'calc_apd.py' script so that we could calculate multiple APDs
-	- Use Elnaz idea for this one ...
-
 - Implement 'save_vtk_polydata_grid_as_vtu_compressed' function
 - Implement coupling between Purkinje and tissue
 - Implement multi-celular-model solver (use the 'extra_data' data structure as Sachetto suggested)
+- Try to reduce the size of the files generated by the "calc_apd_full_grid.py" script
 
 Pedro
 - Upgrade and refactor the DDM code.

diff --git a/example_configs/plain_mesh_ddm_2d.ini b/example_configs/plain_mesh_ddm_2d.ini
@@ -2,7 +2,7 @@
 [main]
 num_threads=2
 dt_pde=0.01
-simulation_time=200.0
+simulation_time=1000.0
 abort_on_no_activity=true
 use_adaptivity=false
 

diff --git a/example_configs/plain_mesh_no_fibrosis_example.ini b/example_configs/plain_mesh_no_fibrosis_example.ini
@@ -1,7 +1,7 @@
 [main]
 num_threads=2
 dt_pde=0.02
-simulation_time=500.0
+simulation_time=1000.0
 abort_on_no_activity=false
 use_adaptivity=false
 
@@ -11,7 +11,7 @@ function=update_monodomain_default
 [save_result]
 ;/////mandatory/////////
 print_rate=100
-output_dir=./outputs/plain_100_100_100_tentusscher
+output_dir=./outputs/plain_100_100_100_fhn_4aps
 function=save_as_vtu
 save_pvd=true
 ;//////////////////
@@ -52,22 +52,24 @@ dt_ode=0.02
 use_gpu=yes
 gpu_id=0
 ;////////////////////////////////////////////////
-library_file=shared_libs/libten_tusscher_2006.so
+;library_file=shared_libs/libten_tusscher_2006.so
 ;////////////////////////////////////////////////
-;library_file=shared_libs/libfhn_mod.so
+library_file=shared_libs/libfhn_mod.so
 
 ;////////////////////////////////////////////////
-[stim_plain]
-start = 0.0
-duration = 2.0
-current = -50.0f
-x_limit = 100.0
-function=stim_if_x_less_than
-;////////////////////////////////////////////////
-
 ;[stim_plain]
 ;start = 0.0
 ;duration = 2.0
-;current = 0.5f
+;period = 250.0
+;current = -50.0f
 ;x_limit = 100.0
 ;function=stim_if_x_less_than
+;////////////////////////////////////////////////
+
+[stim_plain]
+start = 0.0
+duration = 2.0
+period = 250.0
+current = 1.0
+x_limit = 100.0
+function=stim_if_x_less_than