std::set< _Type, _Compare, _Alloc > CepsSet
C++ set.
float CepsReal
Need single precision floating point.
int32_t CepsInt
Need 32 bit integer.
Represents a ionic model for a group of cells, i.e. multiple systems of ODEs.
Base class for creating PDEs to solve.
Beeler Reuter model with bounded iNa current to avoid diverging solution because of.
static constexpr const CepsInt _vB
Maximum voltage above which current is bounded.
virtual CepsReal getINa(CepsReal *y, CepsReal gNaCa, CepsReal gNa, CepsReal v) const override
Returns the bounded current I_Na BR current.
BR77Modified(Unknown *u, AbstractPdeProblem *pb, const CepsSet< CepsAttribute > &attrs={}, InputParameters *params=nullptr)
Constructor (sets constants)
void setupWithParameters(InputParameters *p, FunctionDictionary *dico, AbstractPdeProblem *pb)
Sets the constants and the space dependant parameters from text inputs.
static constexpr const CepsInt _vA
Minimum voltage below which current is bounded.
Beeler Reuter (1977) ionic model.
static constexpr const CepsInt _ENa
Index alias.
void computeRates(CepsReal t, CepsReal *y, CepsReal *v, CepsReal *dtyL, CepsReal *dtyNL, CepsReal *dtv, DegreeOfFreedom *dof) const override
Get the linear and non linear part of the evolution function f. Also computes the ionic current.
static constexpr const CepsInt _CaI
Index alias.
ScalarField< DegreeOfFreedom > * m_gNa
This parameter may vary with space.
static constexpr const CepsInt _f
Index alias.
void setupWithParameters(InputParameters *p, FunctionDictionary *dico, AbstractPdeProblem *pb)
Sets the constants and the space dependant parameters from text inputs.
CepsReal convertCmFromCepsUnit(const CepsReal &cm) const final
Convert capacitance from ceps units (uF/cm2) to ionic model units. Does nothing by default.
CepsReal convertCurrentFromCepsUnit(const CepsReal &v) const final
Convert from muA per cm2 to muA per mm2.
BR77(Unknown *u, AbstractPdeProblem *pb, const CepsSet< CepsAttribute > &attrs={}, InputParameters *params=nullptr)
Constructor (sets constants)
void getInitialCondition(CepsReal *v, CepsReal *y) const final
Sets initial values of state variables and transmembrane voltage for a single point....
virtual CepsReal getINa(CepsReal *y, CepsReal gNaCa, CepsReal gNa, CepsReal v) const
Returns the iNa current.
static constexpr const CepsInt _m
Index alias.
CepsReal convertCmToCepsUnit(const CepsReal &cm) const final
Convert capacitance from ionic model units to ceps units (uF/cm2). Does nothing by default.
static constexpr const CepsInt _h
Index alias.
ScalarField< DegreeOfFreedom > * m_gNaCa
This parameter may vary with space.
~BR77()=default
Destructor.
static constexpr const CepsInt _j
Index alias.
static constexpr const CepsInt _d
Index alias.
static constexpr const CepsInt _x1
Index alias.
static constexpr const CepsInt _gS
Index alias.
A degree of freedom for any kind of problem The dof can be associated to a geometrical element or not...
A Field is an object wrapped around a SAFunc functor, defined on at least one domain.
FunctionDictionary that holds functions which can be used to define source terms, boundary conditions...
A class used to defined an unknown of a PDE problem The unknown can be defined on a specific region,...