43 T& make_ref(T&& x) {
return x; }
71 double Fc = 96485.3415;
83 double V_c = 16.404E-3;
274 double f0 = -4333.618335582119;
275 double fc1 = 2570.395355352195;
276 double fs1 = -2051.827278991976;
277 double fc2 = 1329.53611689133;
278 double fs2 = 302.216784558222;
279 double fc3 = 104.943770305116;
280 double fs3 = 218.375174229422;
368 double fcass, fcassi;
379 double E_Na_Nai, E_K_Ki, E_Ca_Cai, E_Ks_Ki, E_Ks_Nai;
380 double I_Na_V, I_Na_Nai;
381 double I_to_V, I_to_Ki;
382 double I_K1_V, I_K1_Ki;
383 double I_Kr_V, I_Kr_Ki;
384 double I_Ks_V, I_Ks_Ki, I_Ks_Nai;
385 double I_CaL_V, I_CaL_Cass;
386 double I_NaCa_V, I_NaCa_Nai, I_NaCa_Cai;
387 double I_NaK_V, I_NaK_Nai;
389 double I_pK_V, I_pK_Ki;
390 double I_bCa_V, I_bCa_Cai;
391 double I_bNa_V, I_bNa_Nai;
392 double I_leak_Cai, I_leak_Casr;
394 double I_rel_Cass, I_rel_Casr, I_rel_Rbar;
395 double I_xfer_Cai, I_xfer_Cass;
396 double k_casr_sr, k1_casr, O_Casr, O_Cass, O_Rbar;
398 void actv_strn(
const double c_Ca,
const double I4f,
const double dt,
double& gf);
399 void actv_strs(
const double c_Ca,
const double dt,
double& Tact,
double& epsX);
405 Array<double>& JAC,
const double Ksac);
413 const double Ts,
const double dt,
const double Istim,
const double Ksac,
417 const double Ts,
const double dt,
const double Istim,
const double Ksac,
Vector<double>& RPAR);
420 const double Ts,
const double dt,
const double Istim,
const double Ksac,
Vector<double>& RPAR);
This module defines data structures for ten Tusscher-Panfilov epicardial cellular activation model fo...
Definition: CepModTtp.h:51
double I_NaK
Na-K pump current.
Definition: CepModTtp.h:323
double K_o
Extracellular K concentration [mM].
Definition: CepModTtp.h:92
double mu_Ca
Viscous-type constant [ms-mM^{-2}].
Definition: CepModTtp.h:261
double SL0
Initial length of sacromeres [um].
Definition: CepModTtp.h:265
double Cm
Cell capacitance per unit surface area [uF/cm^{2}].
Definition: CepModTtp.h:74
double K_bufsr
Ca_sr half-saturation constant for subspace buffer [mM].
Definition: CepModTtp.h:221
double alpha
Factor enhancing outward nature of I_NaCa [-].
Definition: CepModTtp.h:143
double eps_i
Maximum activation [ms^{-1}].
Definition: CepModTtp.h:246
double V_rel
Maximal I_rel conductance [mM/ms].
Definition: CepModTtp.h:182
void getf(const int i, 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.
Definition: CepModTtp.cpp:77
double SLmax
Max. length of sacromeres [um].
Definition: CepModTtp.h:271
double k1p
R to O and RI to I, I_rel transition rate [mM^{-2}/ms].
Definition: CepModTtp.h:185
double I_up
sacroplasmic reticulum Ca pump current
Definition: CepModTtp.h:341
double K_mCa
Ca_i half-saturation constant for I_NaCa [mM].
Definition: CepModTtp.h:134
double I_bNa
background Na current
Definition: CepModTtp.h:335
double K_NaCa
Maximal I_NaCa [pA/pF].
Definition: CepModTtp.h:128
double K_bufc
Ca_i half-saturation constant for cytplasmic buffer [mM].
Definition: CepModTtp.h:215
double c_Ca0
Resting Ca concentration [mM].
Definition: CepModTtp.h:258
double I_bCa
background Ca current
Definition: CepModTtp.h:332
double max_sr
Maximum value of k_casr [-].
Definition: CepModTtp.h:200
double Tscale
Time scaling.
Definition: CepModTtp.h:288
double I_to
transient outward current
Definition: CepModTtp.h:308
double eta_T
Saturation of concentration [MPa/mM].
Definition: CepModTtp.h:240
double alFa
Active force of sacromere [-mM^{-2}].
Definition: CepModTtp.h:255
double sV
Surface to volume ratio [um^{-1}].
Definition: CepModTtp.h:77
double V_leak
Maximal I_leak conductance [mM/ms].
Definition: CepModTtp.h:206
double V_xfer
Maximal I_xfer conductance [mM/ms].
Definition: CepModTtp.h:209
double Voffset
Voltage offset parameter.
Definition: CepModTtp.h:291
double K_mK
K_o half-saturation constant of I_NaK [mM].
Definition: CepModTtp.h:149
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.
Definition: CepModTtp.cpp:575
double rho
Cellular resistivity [ -cm].
Definition: CepModTtp.h:80
double I_rel
Ca induced Ca release current.
Definition: CepModTtp.h:344
double G_pCa
Maximal I_pCa conductance [pA/pF].
Definition: CepModTtp.h:161
double V_sr
Sacroplasmic reticulum volume [um^{3}].
Definition: CepModTtp.h:86
double G_pK
Maximal I_pK conductance [nS/pF].
Definition: CepModTtp.h:155
double G_Na
Maximal I_Na conductance [nS/pF].
Definition: CepModTtp.h:101
double G_bCa
Maximal I_bCa conductance [nS/pF].
Definition: CepModTtp.h:173
CepModTtp()
Definition: CepModTtp.cpp:37
Vector< double > G_to
Maximal epicardial I_to conductance [nS/pF].
Definition: CepModTtp.h:107
double xi_T
Transition rate [mM^{-1}].
Definition: CepModTtp.h:249
double K_mNai
Na_i half-saturation constant for I_NaCa [mM].
Definition: CepModTtp.h:137
double Buf_sr
Total sacroplasmic buffer concentration [mM].
Definition: CepModTtp.h:218
double Vscale
Voltage scaling.
Definition: CepModTtp.h:285
double k3
O to R and I to RI, I_rel transition rate [ms^{-1}].
Definition: CepModTtp.h:191
double I_Kr
rapid delayed rectifier current
Definition: CepModTtp.h:311
double I_NaCa
Na-Ca exchanger current.
Definition: CepModTtp.h:320
double p_NaK
Maximal I_NaK [pA/pF].
Definition: CepModTtp.h:146
double I_Na
Fast sodium current.
Definition: CepModTtp.h:302
double K_mNa
Na_i half-saturation constant of I_NaK [mM].
Definition: CepModTtp.h:152
double SLmin
Min. length of sacromeres [um].
Definition: CepModTtp.h:268
double V_ss
Subspace volume [um^{3}].
Definition: CepModTtp.h:89
double I_K1
inward rectifier outward current
Definition: CepModTtp.h:305
double p_KNa
Relative I_Ks permeability to Na [-].
Definition: CepModTtp.h:122
double Rc
Gas constant [J/mol/K].
Definition: CepModTtp.h:65
double min_sr
Minimum value of k_casr [-].
Definition: CepModTtp.h:203
double I_xfer
diffusive Ca current
Definition: CepModTtp.h:347
double G_bNa
Maximal I_bNa conductance [nS/pF].
Definition: CepModTtp.h:170
double Buf_ss
Total subspace buffer concentration [mM].
Definition: CepModTtp.h:224
double K_up
Half-saturation constant of I_up [mM].
Definition: CepModTtp.h:179
double Tc
Temperature [K].
Definition: CepModTtp.h:68
double Buf_c
Total cytoplasmic buffer concentration [mM].
Definition: CepModTtp.h:212
double Na_o
Extracellular Na concentration [mM].
Definition: CepModTtp.h:95
double I_pK
plateau K current
Definition: CepModTtp.h:329
double f0
Fourier coefficients.
Definition: CepModTtp.h:274
double Vrest
Resting potential [mV].
Definition: CepModTtp.h:230
double Ca_crit
Critical Ca concentration [mM].
Definition: CepModTtp.h:237
double E_Na
Reverse potentials for Na, K, Ca.
Definition: CepModTtp.h:296
void update_g(const int i, const double dt, const int n, const int nG, const Vector< double > &X, Vector< double > &Xg)
Update all the gating variables.
Definition: CepModTtp.cpp:691
double G_Kr
Maximal I_Kr conductance [nS/pF].
Definition: CepModTtp.h:110
double Ca_rest
Resting Ca concentration [mM].
Definition: CepModTtp.h:234
double Fc
Faraday constant [C/mmol].
Definition: CepModTtp.h:71
double k4
I to O and Ri to I, I_rel transition rate [ms^{-1}].
Definition: CepModTtp.h:194
double eps_0
Minimum activation [ms^{-1}].
Definition: CepModTtp.h:243
double V_c
Cytoplasmic volume [um^{3}].
Definition: CepModTtp.h:83
double k2p
O to I and R to RI, I_rel transition rate [mM^{-1}/ms].
Definition: CepModTtp.h:188
double I_Ks
slow delayed rectifier current
Definition: CepModTtp.h:314
double I_CaL
L-type Ca current.
Definition: CepModTtp.h:317
double G_CaL
Maximal I_CaL conductance [cm^{3}/uF/ms].
Definition: CepModTtp.h:125
double K_sat
Saturation factor for I_NaCa [-].
Definition: CepModTtp.h:140
double K_pCa
Half-saturation constant of I_pCa [mM].
Definition: CepModTtp.h:167
double gamma
Voltage dependent parameter of I_NaCa [-].
Definition: CepModTtp.h:131
Vector< double > G_Ks
Maximal epicardial I_Ks conductance [nS/pF].
Definition: CepModTtp.h:116
double EC
Ca_sr half-saturation constant of k_casr [mM].
Definition: CepModTtp.h:197
double Ca_o
Extracellular Ca concentration [mM].
Definition: CepModTtp.h:98
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 concentrati...
Definition: CepModTtp.cpp:46
double I_pCa
plateau Ca current
Definition: CepModTtp.h:326
double K_bufss
Ca_ss half-saturation constant for subspace buffer [mM].
Definition: CepModTtp.h:227
double G_K1
Maximal I_K1 conductance [nS/pF].
Definition: CepModTtp.h:104
double Vmax_up
Maximal I_up conductance [mM/ms].
Definition: CepModTtp.h:176
double I_leak
sacroplasmic reticulum Ca leak current
Definition: CepModTtp.h:338