36#ifndef LUORUDYSPIRALWAVECELLFACTORY_HPP_
37#define LUORUDYSPIRALWAVECELLFACTORY_HPP_
39#include "LuoRudy1991BackwardEuler.hpp"
40#include "AbstractCardiacCellFactory.hpp"
41#include "MultiStimulus.hpp"
42#include "SimpleStimulus.hpp"
95 double x_threshold_for_S1 = 0.1 + 1e-6;
96 double x_threshold_for_S2 =
mXExtent*0.6;
97 double y_threshold_for_S2 =
mYExtent*0.5;
100 if (x < x_threshold_for_S1)
102 if (y < y_threshold_for_S2)
111 else if ((x < x_threshold_for_S2) && (y < y_threshold_for_S2))
122 p_cell->
SetParameter(
"membrane_L_type_calcium_current_conductance",0);
124 unsigned node_index = pNode->
GetIndex();
131 std::cout <<
"G_si = " << p_cell->
GetAnyVariable(
"membrane_L_type_calcium_current_conductance") <<
"\n";
boost::shared_ptr< AbstractIvpOdeSolver > mpSolver
boost::shared_ptr< ZeroStimulus > mpZeroStimulus
void SetParameter(const std::string &rParameterName, double value)
double GetAnyVariable(const std::string &rName, double time=0.0)
LuoRudySpiralWaveCellFactory(double xExtent, double yExtent)
boost::shared_ptr< SimpleStimulus > mpS1Stimulus
boost::shared_ptr< SimpleStimulus > mpS2Stimulus
boost::shared_ptr< MultiStimulus > mpBothStimulus
AbstractCardiacCell * CreateCardiacCellForTissueNode(Node< 2 > *pNode)
const c_vector< double, SPACE_DIM > & rGetLocation() const
unsigned GetIndex() const