#include <Nash2004ContractionModel.hpp>

Inheritance diagram for Nash2004ContractionModel:

Collaboration diagram for Nash2004ContractionModel:

Public Member Functions
	Nash2004ContractionModel ()
void	EvaluateYDerivatives (double time, const std::vector< double > &rY, std::vector< double > &rDY)
void	SetInputParameters (ContractionModelInputParameters &rInputParameters)
void	SetStretchAndStretchRate (double stretch, double stretchRate)
double	GetActiveTension ()
double	GetNextActiveTension ()
bool	IsStretchDependent ()
bool	IsStretchRateDependent ()
Private Attributes
double	mScaledVoltage
Static Private Attributes
static const double	kTa = 47.9
static const double	e0ByT0 = 1.0/100

Detailed Description

Nash2004 contraction model - Nash & Panfilov "Electromechanical mode of excitable tissue to study reentrant cardiac arrthrymias", Progress in Biophysics and Molecular Biology, 2004.

A simple, stretch- and stretch-rate-independent contraction model, just dependent on the voltage. If v is the non-dimensionalised voltage, Ta is given by dTa/dt = eps(v) (kTa*v - Ta); where eps(v) = e0 if v < 0.05, = 10*e0 v >= 0.05.

We use the non-dimensionalisation: V in [-85,40] ---> v in [0,1], ie v=(V+85)/125

Not sure what the appropriate value of e0 is: as the paper uses non-dimensionalised time the e0 above corresponds to "eps0/t0", where eps0 is value used in the paper (1.0) and t0 is the characteristic time. The paper suggests using t0=25.9 (?), which gives Ta growing too quickly at the beginning (as rapidly as the voltage). t0 = 100 coded at the moment..

Definition at line 60 of file Nash2004ContractionModel.hpp.

Constructor & Destructor Documentation

Nash2004ContractionModel::Nash2004ContractionModel ( ) [inline]

Constructor

Definition at line 73 of file Nash2004ContractionModel.hpp.

References AbstractUntemplatedParameterisedSystem::mpSystemInfo, mScaledVoltage, and AbstractParameterisedSystem< std::vector< double > >::mStateVariables.

Member Function Documentation

void Nash2004ContractionModel::EvaluateYDerivatives	(	double	time,
		const std::vector< double > &	rY,
		std::vector< double > &	rDY
	)		`[inline, virtual]`

Calculate the derivative of the Ta (the state variable)

Parameters:

time	time
rY	1D vector containing Ta
rDY	1D vector in which dTa/dt is set

Implements AbstractOdeSystem.

Definition at line 87 of file Nash2004ContractionModel.hpp.

References e0ByT0, kTa, and mScaledVoltage.

double Nash2004ContractionModel::GetActiveTension ( ) [inline, virtual]

Get the current active tension

Implements AbstractContractionModel.

Definition at line 115 of file Nash2004ContractionModel.hpp.

References AbstractParameterisedSystem< std::vector< double > >::rGetStateVariables().

double Nash2004ContractionModel::GetNextActiveTension ( ) [inline, virtual]

Get the active tension corresponding to the temporary stored state variables produced by calling RunDoNotUpdate (and before calling UpdateStateVariables())

Implements AbstractContractionModel.

Definition at line 124 of file Nash2004ContractionModel.hpp.

References AbstractOdeBasedContractionModel::mTemporaryStateVariables.

bool Nash2004ContractionModel::IsStretchDependent ( ) [inline, virtual]

This model is stretch-independent

Implements AbstractContractionModel.

Definition at line 133 of file Nash2004ContractionModel.hpp.

bool Nash2004ContractionModel::IsStretchRateDependent ( ) [inline, virtual]

This model is stretch-rate-independent

Implements AbstractContractionModel.

Definition at line 141 of file Nash2004ContractionModel.hpp.

void Nash2004ContractionModel::SetInputParameters ( ContractionModelInputParameters & rInputParameters ) [inline, virtual]