svMultiPhysics
Loading...
Searching...
No Matches
Public Member Functions | Public Attributes | List of all members
CepModTtp Class Reference

Public Member Functions

void actv_strn (const double c_Ca, const double I4f, const double dt, double &gf)
 Compute macroscopic fiber strain based on sacromere force-length relationship and calcium concentration.
 
void actv_strs (const double c_Ca, const double dt, double &Tact, double &epsX)
 
void getf (const int nX, const int nG, const Vector< double > &X, const Vector< double > &Xg, Vector< double > &dX, const double I_stim, const double K_sac, Vector< double > &RPAR)
 Compute currents and time derivatives of state variables.
 
void getj (const int nX, const int nG, const Vector< double > &X, const Vector< double > &Xg, Array< double > &JAC, const double Ksac)
 
void init (const int nX, const int nG, Vector< double > &X, Vector< double > &Xg)
 
void set_initial_conditions (const TTPInitialConditionsParameters &params)
 
void integ_cn2 (const int imyo, const int nX, const int nG, Vector< double > &X, Vector< double > &Xg, const double Ts, const double dt, const double Istim, const double Ksac, Vector< int > &IPAR, Vector< double > &RPAR)
 Time integration performed using Crank-Nicholson method.
 
void integ_fe (const int imyo, const int nX, const int nG, Vector< double > &X, Vector< double > &Xg, const double Ts, const double dt, const double Istim, const double Ksac, Vector< double > &RPAR)
 
void integ_rk (const int imyo, const int nX, const int nG, Vector< double > &X, Vector< double > &Xg, const double Ts, const double dt, const double Istim, const double Ksac, Vector< double > &RPAR)
 
void update_g (const int imyo, const double dt, const int n, const int nG, const Vector< double > &X, Vector< double > &Xg)
 Update all the gating variables.
 
void copy_state_to_vectors (Vector< double > &X, Vector< double > &Xg) const
 
void distribute_conductance (const CmMod &cm_mod, const cmType &cm)
 
void distribute_initial_state (const CmMod &cm_mod, const cmType &cm)
 

Public Attributes

double Rc = 8314.472
 Gas constant [J/mol/K].
 
double Tc = 310.0
 Temperature [K].
 
double Fc = 96485.3415
 Faraday constant [C/mmol].
 
double Cm = 0.185
 Cell capacitance per unit surface area [uF/cm^{2}].
 
double sV = 0.2
 Surface to volume ratio [um^{-1}].
 
double rho = 162.0
 Cellular resistivity [ $\Omega$-cm].
 
double V_c = 16.404E-3
 Cytoplasmic volume [um^{3}].
 
double V_sr = 1.094E-3
 Sacroplasmic reticulum volume [um^{3}].
 
double V_ss = 5.468E-5
 Subspace volume [um^{3}].
 
double K_o = 5.4
 Extracellular K concentration [mM].
 
double Na_o = 140.0
 Extracellular Na concentration [mM].
 
double Ca_o = 2.0
 Extracellular Ca concentration [mM].
 
double G_Na = 14.838
 Maximal I_Na conductance [nS/pF].
 
double G_K1 = 5.405
 Maximal I_K1 conductance [nS/pF].
 
double G_to = 0.294
 Maximal I_to conductance [nS/pF].
 
double G_Kr = 0.153
 Maximal I_Kr conductance [nS/pF].
 
double G_Ks = 0.392
 Maximal I_Ks conductance [nS/pF].
 
double p_KNa = 3.E-2
 Relative I_Ks permeability to Na [-].
 
double G_CaL = 3.98E-5
 Maximal I_CaL conductance [cm^{3}/uF/ms].
 
double K_NaCa = 1000.
 Maximal I_NaCa [pA/pF].
 
double gamma = 0.35
 Voltage dependent parameter of I_NaCa [-].
 
double K_mCa = 1.38
 Ca_i half-saturation constant for I_NaCa [mM].
 
double K_mNai = 87.5
 Na_i half-saturation constant for I_NaCa [mM].
 
double K_sat = 0.1
 Saturation factor for I_NaCa [-].
 
double alpha = 2.5
 Factor enhancing outward nature of I_NaCa [-].
 
double p_NaK = 2.724
 Maximal I_NaK [pA/pF].
 
double K_mK = 1.
 K_o half-saturation constant of I_NaK [mM].
 
double K_mNa = 40.
 Na_i half-saturation constant of I_NaK [mM].
 
double G_pK = 1.46E-2
 Maximal I_pK conductance [nS/pF].
 
double G_pCa = 0.1238
 Maximal I_pCa conductance [pA/pF].
 
double K_pCa = 5.E-4
 Half-saturation constant of I_pCa [mM].
 
double G_bNa = 2.9E-4
 Maximal I_bNa conductance [nS/pF].
 
double G_bCa = 5.92E-4
 Maximal I_bCa conductance [nS/pF].
 
double Vmax_up = 6.375E-3
 Maximal I_up conductance [mM/ms].
 
double K_up = 2.5E-4
 Half-saturation constant of I_up [mM].
 
double V_rel = 0.102
 Maximal I_rel conductance [mM/ms].
 
double k1p = 0.15
 R to O and RI to I, I_rel transition rate [mM^{-2}/ms].
 
double k2p = 4.5E-2
 O to I and R to RI, I_rel transition rate [mM^{-1}/ms].
 
double k3 = 6.E-2
 O to R and I to RI, I_rel transition rate [ms^{-1}].
 
double k4 = 5.E-3
 I to O and Ri to I, I_rel transition rate [ms^{-1}].
 
double EC = 1.5
 Ca_sr half-saturation constant of k_casr [mM].
 
double max_sr = 2.5
 Maximum value of k_casr [-].
 
double min_sr = 1.
 Minimum value of k_casr [-].
 
double V_leak = 3.6E-4
 Maximal I_leak conductance [mM/ms].
 
double V_xfer = 3.8E-3
 Maximal I_xfer conductance [mM/ms].
 
double Buf_c = 0.2
 Total cytoplasmic buffer concentration [mM].
 
double K_bufc = 1.E-3
 Ca_i half-saturation constant for cytplasmic buffer [mM].
 
double Buf_sr = 10.
 Total sacroplasmic buffer concentration [mM].
 
double K_bufsr = 0.3
 Ca_sr half-saturation constant for subspace buffer [mM].
 
double Buf_ss = 0.4
 Total subspace buffer concentration [mM].
 
double K_bufss = 2.5E-4
 Ca_ss half-saturation constant for subspace buffer [mM].
 
double Vrest = -85.23
 Resting potential [mV].
 
double Ca_rest = 5.E-5
 Resting Ca concentration [mM].
 
double Ca_crit = 8.E-4
 Critical Ca concentration [mM].
 
double eta_T = 12.5
 Saturation of concentration [MPa/mM].
 
double eps_0 = 0.1
 Minimum activation [ms^{-1}].
 
double eps_i = 1.
 Maximum activation [ms^{-1}].
 
double xi_T = 4.E3
 Transition rate [mM^{-1}].
 
double alFa = -4.E6
 Active force of sacromere [-mM^{-2}].
 
double c_Ca0 = 2.155E-4
 Resting Ca concentration [mM].
 
double mu_Ca = 5.E6
 Viscous-type constant [ms-mM^{-2}].
 
double SL0 = 1.95
 Initial length of sacromeres [um].
 
double SLmin = 1.7
 Min. length of sacromeres [um].
 
double SLmax = 2.6
 Max. length of sacromeres [um].
 
double f0 = -4333.618335582119
 Fourier coefficients.
 
double fc1 = 2570.395355352195
 
double fs1 = -2051.827278991976
 
double fc2 = 1329.53611689133
 
double fs2 = 302.216784558222
 
double fc3 = 104.943770305116
 
double fs3 = 218.375174229422
 
double Vscale = 1.
 Voltage scaling.
 
double Tscale = 1.
 Time scaling.
 
double Voffset = 0.
 Voltage offset parameter.
 
double E_Na
 Reverse potentials for Na, K, Ca.
 
double E_K
 
double E_Ca
 
double E_Ks
 
double I_Na
 Fast sodium current.
 
double I_K1
 inward rectifier outward current
 
double I_to
 transient outward current
 
double I_Kr
 rapid delayed rectifier current
 
double I_Ks
 slow delayed rectifier current
 
double I_CaL
 L-type Ca current.
 
double I_NaCa
 Na-Ca exchanger current.
 
double I_NaK
 Na-K pump current.
 
double I_pCa
 plateau Ca current
 
double I_pK
 plateau K current
 
double I_bCa
 background Ca current
 
double I_bNa
 background Na current
 
double I_leak
 sacroplasmic reticulum Ca leak current
 
double I_up
 sacroplasmic reticulum Ca pump current
 
double I_rel
 Ca induced Ca release current.
 
double I_xfer
 diffusive Ca current
 
double V
 
double K_i
 
double Na_i
 
double Ca_i
 
double Ca_ss
 
double Ca_sr
 
double R_bar
 
double xr1
 
double xr1i
 
double xr2
 
double xr2i
 
double xs
 
double xsi
 
double m
 
double mi
 
double h
 
double hi
 
double j
 
double ji
 
double d
 
double di
 
double f
 
double fi
 
double f2
 
double f2i
 
double fcass
 
double fcassi
 
double s
 
double si
 
double r
 
double ri
 
double k1
 
double k2
 
double k_casr
 
double O
 
double E_Na_Nai
 
double E_K_Ki
 
double E_Ca_Cai
 
double E_Ks_Ki
 
double E_Ks_Nai
 
double I_Na_V
 
double I_Na_Nai
 
double I_to_V
 
double I_to_Ki
 
double I_K1_V
 
double I_K1_Ki
 
double I_Kr_V
 
double I_Kr_Ki
 
double I_Ks_V
 
double I_Ks_Ki
 
double I_Ks_Nai
 
double I_CaL_V
 
double I_CaL_Cass
 
double I_NaCa_V
 
double I_NaCa_Nai
 
double I_NaCa_Cai
 
double I_NaK_V
 
double I_NaK_Nai
 
double I_pCa_Cai
 
double I_pK_V
 
double I_pK_Ki
 
double I_bCa_V
 
double I_bCa_Cai
 
double I_bNa_V
 
double I_bNa_Nai
 
double I_leak_Cai
 
double I_leak_Casr
 
double I_up_Cai
 
double I_rel_Cass
 
double I_rel_Casr
 
double I_rel_Rbar
 
double I_xfer_Cai
 
double I_xfer_Cass
 
double k_casr_sr
 
double k1_casr
 
double O_Casr
 
double O_Cass
 
double O_Rbar
 
bool user_initial_state = false
 

Member Function Documentation

◆ actv_strn()

void CepModTtp::actv_strn ( const double  c_Ca,
const double  I4f,
const double  dt,
double &  gf 
)

Compute macroscopic fiber strain based on sacromere force-length relationship and calcium concentration.

◆ getf()

void CepModTtp::getf ( const int  nX,
const int  nG,
const Vector< double > &  X,
const Vector< double > &  Xg,
Vector< double > &  dX,
const double  I_stim,
const double  K_sac,
Vector< double > &  RPAR 
)

Compute currents and time derivatives of state variables.

Reproduces Fortran 'GETF()'.

◆ integ_cn2()

void CepModTtp::integ_cn2 ( const int  imyo,
const int  nX,
const int  nG,
Vector< double > &  X,
Vector< double > &  Xg,
const double  Ts,
const double  dt,
const double  Istim,
const double  Ksac,
Vector< int > &  IPAR,
Vector< double > &  RPAR 
)

Time integration performed using Crank-Nicholson method.

◆ update_g()

void CepModTtp::update_g ( const int  imyo,
const double  dt,
const int  n,
const int  nG,
const Vector< double > &  X,
Vector< double > &  Xg 
)

Update all the gating variables.

Member Data Documentation

◆ alFa

double CepModTtp::alFa = -4.E6

Active force of sacromere [-mM^{-2}].

◆ alpha

double CepModTtp::alpha = 2.5

Factor enhancing outward nature of I_NaCa [-].

◆ Buf_c

double CepModTtp::Buf_c = 0.2

Total cytoplasmic buffer concentration [mM].

◆ Buf_sr

double CepModTtp::Buf_sr = 10.

Total sacroplasmic buffer concentration [mM].

◆ Buf_ss

double CepModTtp::Buf_ss = 0.4

Total subspace buffer concentration [mM].

◆ c_Ca0

double CepModTtp::c_Ca0 = 2.155E-4

Resting Ca concentration [mM].

◆ Ca_crit

double CepModTtp::Ca_crit = 8.E-4

Critical Ca concentration [mM].

◆ Ca_o

double CepModTtp::Ca_o = 2.0

Extracellular Ca concentration [mM].

◆ Ca_rest

double CepModTtp::Ca_rest = 5.E-5

Resting Ca concentration [mM].

◆ Cm

double CepModTtp::Cm = 0.185

Cell capacitance per unit surface area [uF/cm^{2}].

◆ E_Na

double CepModTtp::E_Na

Reverse potentials for Na, K, Ca.

◆ EC

double CepModTtp::EC = 1.5

Ca_sr half-saturation constant of k_casr [mM].

◆ eps_0

double CepModTtp::eps_0 = 0.1

Minimum activation [ms^{-1}].

◆ eps_i

double CepModTtp::eps_i = 1.

Maximum activation [ms^{-1}].

◆ eta_T

double CepModTtp::eta_T = 12.5

Saturation of concentration [MPa/mM].

◆ f0

double CepModTtp::f0 = -4333.618335582119

Fourier coefficients.

◆ Fc

double CepModTtp::Fc = 96485.3415

Faraday constant [C/mmol].

◆ G_bCa

double CepModTtp::G_bCa = 5.92E-4

Maximal I_bCa conductance [nS/pF].

◆ G_bNa

double CepModTtp::G_bNa = 2.9E-4

Maximal I_bNa conductance [nS/pF].

◆ G_CaL

double CepModTtp::G_CaL = 3.98E-5

Maximal I_CaL conductance [cm^{3}/uF/ms].

◆ G_K1

double CepModTtp::G_K1 = 5.405

Maximal I_K1 conductance [nS/pF].

◆ G_Kr

double CepModTtp::G_Kr = 0.153

Maximal I_Kr conductance [nS/pF].

◆ G_Ks

double CepModTtp::G_Ks = 0.392

Maximal I_Ks conductance [nS/pF].

◆ G_Na

double CepModTtp::G_Na = 14.838

Maximal I_Na conductance [nS/pF].

◆ G_pCa

double CepModTtp::G_pCa = 0.1238

Maximal I_pCa conductance [pA/pF].

◆ G_pK

double CepModTtp::G_pK = 1.46E-2

Maximal I_pK conductance [nS/pF].

◆ G_to

double CepModTtp::G_to = 0.294

Maximal I_to conductance [nS/pF].

◆ gamma

double CepModTtp::gamma = 0.35

Voltage dependent parameter of I_NaCa [-].

◆ I_bCa

double CepModTtp::I_bCa

background Ca current

◆ I_bNa

double CepModTtp::I_bNa

background Na current

◆ I_CaL

double CepModTtp::I_CaL

L-type Ca current.

◆ I_K1

double CepModTtp::I_K1

inward rectifier outward current

◆ I_Kr

double CepModTtp::I_Kr

rapid delayed rectifier current

◆ I_Ks

double CepModTtp::I_Ks

slow delayed rectifier current

◆ I_leak

double CepModTtp::I_leak

sacroplasmic reticulum Ca leak current

◆ I_Na

double CepModTtp::I_Na

Fast sodium current.

◆ I_NaCa

double CepModTtp::I_NaCa

Na-Ca exchanger current.

◆ I_NaK

double CepModTtp::I_NaK

Na-K pump current.

◆ I_pCa

double CepModTtp::I_pCa

plateau Ca current

◆ I_pK

double CepModTtp::I_pK

plateau K current

◆ I_rel

double CepModTtp::I_rel

Ca induced Ca release current.

◆ I_to

double CepModTtp::I_to

transient outward current

◆ I_up

double CepModTtp::I_up

sacroplasmic reticulum Ca pump current

◆ I_xfer

double CepModTtp::I_xfer

diffusive Ca current

◆ k1p

double CepModTtp::k1p = 0.15

R to O and RI to I, I_rel transition rate [mM^{-2}/ms].

◆ k2p

double CepModTtp::k2p = 4.5E-2

O to I and R to RI, I_rel transition rate [mM^{-1}/ms].

◆ k3

double CepModTtp::k3 = 6.E-2

O to R and I to RI, I_rel transition rate [ms^{-1}].

◆ k4

double CepModTtp::k4 = 5.E-3

I to O and Ri to I, I_rel transition rate [ms^{-1}].

◆ K_bufc

double CepModTtp::K_bufc = 1.E-3

Ca_i half-saturation constant for cytplasmic buffer [mM].

◆ K_bufsr

double CepModTtp::K_bufsr = 0.3

Ca_sr half-saturation constant for subspace buffer [mM].

◆ K_bufss

double CepModTtp::K_bufss = 2.5E-4

Ca_ss half-saturation constant for subspace buffer [mM].

◆ K_mCa

double CepModTtp::K_mCa = 1.38

Ca_i half-saturation constant for I_NaCa [mM].

◆ K_mK

double CepModTtp::K_mK = 1.

K_o half-saturation constant of I_NaK [mM].

◆ K_mNa

double CepModTtp::K_mNa = 40.

Na_i half-saturation constant of I_NaK [mM].

◆ K_mNai

double CepModTtp::K_mNai = 87.5

Na_i half-saturation constant for I_NaCa [mM].

◆ K_NaCa

double CepModTtp::K_NaCa = 1000.

Maximal I_NaCa [pA/pF].

◆ K_o

double CepModTtp::K_o = 5.4

Extracellular K concentration [mM].

◆ K_pCa

double CepModTtp::K_pCa = 5.E-4

Half-saturation constant of I_pCa [mM].

◆ K_sat

double CepModTtp::K_sat = 0.1

Saturation factor for I_NaCa [-].

◆ K_up

double CepModTtp::K_up = 2.5E-4

Half-saturation constant of I_up [mM].

◆ max_sr

double CepModTtp::max_sr = 2.5

Maximum value of k_casr [-].

◆ min_sr

double CepModTtp::min_sr = 1.

Minimum value of k_casr [-].

◆ mu_Ca

double CepModTtp::mu_Ca = 5.E6

Viscous-type constant [ms-mM^{-2}].

◆ Na_o

double CepModTtp::Na_o = 140.0

Extracellular Na concentration [mM].

◆ p_KNa

double CepModTtp::p_KNa = 3.E-2

Relative I_Ks permeability to Na [-].

◆ p_NaK

double CepModTtp::p_NaK = 2.724

Maximal I_NaK [pA/pF].

◆ Rc

double CepModTtp::Rc = 8314.472

Gas constant [J/mol/K].

◆ rho

double CepModTtp::rho = 162.0

Cellular resistivity [ $\Omega$-cm].

◆ SL0

double CepModTtp::SL0 = 1.95

Initial length of sacromeres [um].

◆ SLmax

double CepModTtp::SLmax = 2.6

Max. length of sacromeres [um].

◆ SLmin

double CepModTtp::SLmin = 1.7

Min. length of sacromeres [um].

◆ sV

double CepModTtp::sV = 0.2

Surface to volume ratio [um^{-1}].

◆ Tc

double CepModTtp::Tc = 310.0

Temperature [K].

◆ Tscale

double CepModTtp::Tscale = 1.

Time scaling.

◆ V_c

double CepModTtp::V_c = 16.404E-3

Cytoplasmic volume [um^{3}].

◆ V_leak

double CepModTtp::V_leak = 3.6E-4

Maximal I_leak conductance [mM/ms].

◆ V_rel

double CepModTtp::V_rel = 0.102

Maximal I_rel conductance [mM/ms].

◆ V_sr

double CepModTtp::V_sr = 1.094E-3

Sacroplasmic reticulum volume [um^{3}].

◆ V_ss

double CepModTtp::V_ss = 5.468E-5

Subspace volume [um^{3}].

◆ V_xfer

double CepModTtp::V_xfer = 3.8E-3

Maximal I_xfer conductance [mM/ms].

◆ Vmax_up

double CepModTtp::Vmax_up = 6.375E-3

Maximal I_up conductance [mM/ms].

◆ Voffset

double CepModTtp::Voffset = 0.

Voltage offset parameter.

◆ Vrest

double CepModTtp::Vrest = -85.23

Resting potential [mV].

◆ Vscale

double CepModTtp::Vscale = 1.

Voltage scaling.

◆ xi_T

double CepModTtp::xi_T = 4.E3

Transition rate [mM^{-1}].

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.8