Electrophysiology benchmark proposed in Niederer et al. 2011. More...

#include <lifev/electrophysiology/solver/ElectroETAMonodomainSolver.hpp>
#include <lifev/electrophysiology/testsuite/test_benchmark/benchmarkUtility.hpp>
#include <sys/stat.h>

Include dependency graph for lifev/electrophysiology/testsuite/test_benchmark/main.cpp:

Go to the source code of this file.

Macros
#define	finalActivationTime 53.1

Functions
Int	main (Int argc, char **argv)

Detailed Description

Electrophysiology benchmark proposed in Niederer et al. 2011.

This test is the basic test to start using the electrophysiology module. We consider a parallelepiped excited in a small region in one corner. The excitation wave propagates and excites the whole domain. We check the time it takes for the wave to activate the whole domain. We provide two coarse meshes with this test. Be aware that with such mesh sizes we are far from convergence.

The parameter list is an xml file. To work properly it should be called MonodomainSolverParamList.xml as the provided one, or change the name of the file (see below in the code.)

Run it using e.g.: mpirun -n 2 Electrophysiology_benchmark -o OutputFolder

WARNING: The provided mesh is coarse!!! You should not expect nice results using such mesh. A finer mesh can be obtained just using gmsh and using the refine by splitting function. The required mesh would be too large and we decided to provide only a coarse one.

NOTICE: The value of the membrane capacitance has a big influence on the propagation. For example the TenTusscher model has Cm = 2 and final activation time (in this benchmark) about twice the one computed with Luo-Rudy I or the Minimal Model. Imposing Cm = 1, then, these ionic models give similar results in terms of activation times. Be aware of what you are doing! Don't use the default parameters as black boxes!!!