#include <TenTusscher06.h>
Inheritance diagram for Nektar::TenTusscher06:
Public Member Functions | |
| TenTusscher06 (const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField) | |
| Constructor. More... | |
| virtual | ~TenTusscher06 () |
| Desctructor. More... | |
 Public Member Functions inherited from Nektar::CellModel | |
| CellModel (const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField) | |
| virtual | ~CellModel () |
| void | Initialise () |
| Initialise the cell model storage and set initial conditions. More... | |
| void | TimeIntegrate (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time) |
| Time integrate the cell model by one PDE timestep. More... | |
| void | Update (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time) |
| Compute the derivatives of cell model variables. More... | |
| void | GenerateSummary (SummaryList &s) |
| Print a summary of the cell model. More... | |
| size_t | GetNumCellVariables () |
| std::string | GetCellVarName (size_t idx) |
| Array< OneD, NekDouble > | GetCellSolutionCoeffs (size_t idx) |
| Array< OneD, NekDouble > | GetCellSolution (size_t idx) |
Static Public Member Functions | |
| static CellModelSharedPtr | create (const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField) |
| Creates an instance of this class. More... | |
Static Public Attributes | |
| static std::string | className |
| Name of class. More... | |
Protected Types | |
| enum | Variants { eEpicardium , eEndocardium , eMid , eIschemia } |
Protected Member Functions | |
| virtual void | v_Update (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time) override |
| virtual void | v_GenerateSummary (SummaryList &s) override |
| Prints a summary of the model parameters. More... | |
| virtual void | v_SetInitialConditions () override |
| Protected Member Functions inherited from Nektar::CellModel | |
| virtual void | v_Update (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)=0 |
| virtual void | v_GenerateSummary (SummaryList &s)=0 |
| virtual std::string | v_GetCellVarName (size_t idx) |
| virtual void | v_SetInitialConditions ()=0 |
| void | LoadCellModel () |
Protected Attributes | |
| NekDouble | g_to |
| NekDouble | g_Ks |
| NekDouble | s_inf_factor |
| NekDouble | s_tau_f1 |
| NekDouble | s_tau_f2 |
| NekDouble | s_tau_f3 |
| NekDouble | s_tau_f4 |
| NekDouble | s_tau_f5 |
| NekDouble | k_0 |
| enum Variants | model_variant |
| Protected Attributes inherited from Nektar::CellModel | |
| LibUtilities::SessionReaderSharedPtr | m_session |
| Session. More... | |
| MultiRegions::ExpListSharedPtr | m_field |
| Transmembrane potential field from PDE system. More... | |
| size_t | m_nq |
| Number of physical points. More... | |
| size_t | m_nvar |
| Number of variables in cell model (inc. transmembrane voltage) More... | |
| NekDouble | m_lastTime |
| Timestep for pde model. More... | |
| size_t | m_substeps |
| Number of substeps to take. More... | |
| Array< OneD, Array< OneD, NekDouble > > | m_cellSol |
| Cell model solution variables. More... | |
| Array< OneD, Array< OneD, NekDouble > > | m_wsp |
| Cell model integration workspace. More... | |
| bool | m_useNodal |
| Flag indicating whether nodal projection in use. More... | |
| StdRegions::StdNodalTriExpSharedPtr | m_nodalTri |
| StdNodalTri for cell model calculations. More... | |
| StdRegions::StdNodalTetExpSharedPtr | m_nodalTet |
| Array< OneD, Array< OneD, NekDouble > > | m_nodalTmp |
| Temporary array for nodal projection. More... | |
| std::vector< int > | m_concentrations |
| Indices of cell model variables which are concentrations. More... | |
| std::vector< int > | m_gates |
| Indices of cell model variables which are gates. More... | |
| Array< OneD, Array< OneD, NekDouble > > | m_gates_tau |
| Storage for gate tau values. More... | |
Static Protected Attributes | |
| static std::string | lookupIds [] |
| static std::string | def |
Detailed Description
Definition at line 42 of file TenTusscher06.h.
Member Enumeration Documentation
◆ Variants
|
protected |
| Enumerator | |
|---|---|
| eEpicardium | |
| eEndocardium | |
| eMid | |
| eIschemia | |
Definition at line 88 of file TenTusscher06.h.
Constructor & Destructor Documentation
◆ TenTusscher06()
| Nektar::TenTusscher06::TenTusscher06 | ( | const LibUtilities::SessionReaderSharedPtr & | pSession, |
| const MultiRegions::ExpListSharedPtr & | pField | ||
| ) |
Constructor.
Definition at line 64 of file TenTusscher06.cpp.
67 : CellModel(pSession, pField)
68{
70 "CellModelVariant");
71
73 {
75 g_to = 0.294;
76 g_Ks = 0.392;
77 s_inf_factor = 20.0;
78 s_tau_f1 = 85.0;
79 s_tau_f2 = 45.0;
80 s_tau_f3 = 320.0;
81 s_tau_f4 = 0.0;
82 s_tau_f5 = 1.0;
83 k_0 = 5.4;
84 break;
86 g_to = 0.073;
87 g_Ks = 0.392;
88 s_inf_factor = 28.0;
89 s_tau_f1 = 1000.0;
90 s_tau_f2 = 67.0;
91 s_tau_f3 = 1000.0;
92 s_tau_f4 = 8.0;
93 s_tau_f5 = 0.0;
94 k_0 = 5.4;
95 break;
97 g_to = 0.294;
98 g_Ks = 0.098;
99 s_inf_factor = 20.0;
100 s_tau_f1 = 85.0;
101 s_tau_f2 = 45.0;
102 s_tau_f3 = 320.0;
103 s_tau_f4 = 0.0;
104 s_tau_f5 = 1.0;
105 k_0 = 5.4;
106 break;
108 g_to = 0.294;
109 g_Ks = 0.392;
110 s_inf_factor = 20.0;
111 s_tau_f1 = 85.0;
112 s_tau_f2 = 45.0;
113 s_tau_f3 = 320.0;
114 s_tau_f4 = 0.0;
115 s_tau_f5 = 1.0;
116 k_0 = 9.0;
117 break;
118 }
119 m_nvar = 19;
120
121 m_gates.push_back(1);
122 m_gates.push_back(2);
123 m_gates.push_back(3);
124 m_gates.push_back(4);
125 m_gates.push_back(5);
126 m_gates.push_back(6);
127 m_gates.push_back(7);
128 m_gates.push_back(8);
129 m_gates.push_back(9);
130 m_gates.push_back(10);
131 m_gates.push_back(11);
132 m_gates.push_back(12);
133 m_concentrations.push_back(13);
134 m_concentrations.push_back(14);
135 m_concentrations.push_back(15);
136 m_concentrations.push_back(16);
137 m_concentrations.push_back(17);
138 m_concentrations.push_back(18);
139}
std::vector< int > m_concentrations
Indices of cell model variables which are concentrations.
Definition: CellModel.h:139
std::vector< int > m_gates
Indices of cell model variables which are gates.
Definition: CellModel.h:141
size_t m_nvar
Number of variables in cell model (inc. transmembrane voltage)
Definition: CellModel.h:119
CellModel(const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField)
Definition: CellModel.cpp:66
References eEndocardium, eEpicardium, eIschemia, eMid, g_Ks, g_to, k_0, Nektar::CellModel::m_concentrations, Nektar::CellModel::m_gates, Nektar::CellModel::m_nvar, model_variant, s_inf_factor, s_tau_f1, s_tau_f2, s_tau_f3, s_tau_f4, and s_tau_f5.
◆ ~TenTusscher06()
|
inlinevirtual |
Member Function Documentation
◆ create()
|
inlinestatic |
Creates an instance of this class.
Definition at line 47 of file TenTusscher06.h.
50 {
52 pField);
53 }
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:168
References Nektar::MemoryManager< DataType >::AllocateSharedPtr().
◆ v_GenerateSummary()
|
overrideprotectedvirtual |
Prints a summary of the model parameters.
Implements Nektar::CellModel.
Definition at line 1079 of file TenTusscher06.cpp.
1080{
1082}
void AddSummaryItem(SummaryList &l, const std::string &name, const std::string &value)
Adds a summary item to the summary info list.
Definition: Misc.cpp:49
References Nektar::SolverUtils::AddSummaryItem().
◆ v_SetInitialConditions()
|
overrideprotectedvirtual |
Implements Nektar::CellModel.
Definition at line 1087 of file TenTusscher06.cpp.
1088{
1089 // Values taken from CellML website
1091 {
1112 break;
1133 break;
1154 break;
1175 break;
1176 }
1177}
Array< OneD, Array< OneD, NekDouble > > m_cellSol
Cell model solution variables.
Definition: CellModel.h:126
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:43
References eEndocardium, eEpicardium, eIschemia, eMid, Vmath::Fill(), Nektar::CellModel::m_cellSol, Nektar::CellModel::m_nq, and model_variant.
◆ v_Update()
|
overrideprotectedvirtual |
Implements Nektar::CellModel.
Definition at line 144 of file TenTusscher06.cpp.
147{
148 boost::ignore_unused(time);
149
151 {
152 // Inputs:
153 // Time units: millisecond
154 NekDouble var_chaste_interface__membrane__V = inarray[0][i];
155 // Units: millivolt; Initial value: -85.23
156 NekDouble
157 var_chaste_interface__rapid_time_dependent_potassium_current_Xr1_gate__Xr1 =
158 inarray[1][i];
159 // Units: dimensionless; Initial value: 0.00621
160 NekDouble
161 var_chaste_interface__rapid_time_dependent_potassium_current_Xr2_gate__Xr2 =
162 inarray[2][i];
163 // Units: dimensionless; Initial value: 0.4712
164 NekDouble
165 var_chaste_interface__slow_time_dependent_potassium_current_Xs_gate__Xs =
166 inarray[3][i];
167 // Units: dimensionless; Initial value: 0.0095
168 NekDouble var_chaste_interface__fast_sodium_current_m_gate__m =
169 inarray[4][i];
170 // Units: dimensionless; Initial value: 0.00172
171 NekDouble var_chaste_interface__fast_sodium_current_h_gate__h =
172 inarray[5][i];
173 // Units: dimensionless; Initial value: 0.7444
174 NekDouble var_chaste_interface__fast_sodium_current_j_gate__j =
175 inarray[6][i];
176 // Units: dimensionless; Initial value: 0.7045
177 NekDouble var_chaste_interface__L_type_Ca_current_d_gate__d =
178 inarray[7][i];
179 // Units: dimensionless; Initial value: 3.373e-5
180 NekDouble var_chaste_interface__L_type_Ca_current_f_gate__f =
181 inarray[8][i];
182 // Units: dimensionless; Initial value: 0.7888
183 NekDouble var_chaste_interface__L_type_Ca_current_f2_gate__f2 =
184 inarray[9][i];
185 // Units: dimensionless; Initial value: 0.9755
186 NekDouble var_chaste_interface__L_type_Ca_current_fCass_gate__fCass =
187 inarray[10][i];
188 // Units: dimensionless; Initial value: 0.9953
189 NekDouble var_chaste_interface__transient_outward_current_s_gate__s =
190 inarray[11][i];
191 // Units: dimensionless; Initial value: 0.999998
192 NekDouble var_chaste_interface__transient_outward_current_r_gate__r =
193 inarray[12][i];
194 // Units: dimensionless; Initial value: 2.42e-8
195 NekDouble var_chaste_interface__calcium_dynamics__Ca_i = inarray[13][i];
196 // Units: millimolar; Initial value: 0.000126
197 NekDouble var_chaste_interface__calcium_dynamics__Ca_SR =
198 inarray[14][i];
199 // Units: millimolar; Initial value: 3.64
200 NekDouble var_chaste_interface__calcium_dynamics__Ca_ss =
201 inarray[15][i];
202 // Units: millimolar; Initial value: 0.00036
203 NekDouble var_chaste_interface__calcium_dynamics__R_prime =
204 inarray[16][i];
205 // Units: dimensionless; Initial value: 0.9073
206 NekDouble var_chaste_interface__sodium_dynamics__Na_i = inarray[17][i];
207 // Units: millimolar; Initial value: 8.604
208 NekDouble var_chaste_interface__potassium_dynamics__K_i =
209 inarray[18][i];
210 // Units: millimolar; Initial value: 136.89
211
212 // Mathematics
213 NekDouble d_dt_chaste_interface__membrane__V;
220 var_chaste_interface__membrane__V; // millivolt
222 var_chaste_interface__potassium_dynamics__K_i; // millimolar
224 var_membrane__R; // joule_per_mole_kelvin
227 var_membrane__F; // coulomb_per_millimole
230 var_potassium_dynamics__K_o; // millimolar
232 ((var_reversal_potentials__R * var_reversal_potentials__T) /
233 var_reversal_potentials__F) *
234 log(var_reversal_potentials__K_o /
235 var_reversal_potentials__K_i); // millivolt
237 var_reversal_potentials__E_K; // millivolt
239 ((3.0 *
240 exp(0.0002 * ((var_inward_rectifier_potassium_current__V -
241 var_inward_rectifier_potassium_current__E_K) +
242 100.0))) +
243 exp(0.1 * ((var_inward_rectifier_potassium_current__V -
244 var_inward_rectifier_potassium_current__E_K) -
245 10.0))) /
246 (1.0 +
247 exp((-0.5) *
248 (var_inward_rectifier_potassium_current__V -
249 var_inward_rectifier_potassium_current__E_K))); // dimensionless
251 0.1 /
252 (1.0 + exp(0.06 * ((var_inward_rectifier_potassium_current__V -
253 var_inward_rectifier_potassium_current__E_K) -
254 200.0))); // dimensionless
256 var_inward_rectifier_potassium_current__alpha_K1 /
257 (var_inward_rectifier_potassium_current__alpha_K1 +
258 var_inward_rectifier_potassium_current__beta_K1); // dimensionless
260 var_potassium_dynamics__K_o; // millimolar
262 5.405; // nanoS_per_picoF
264 var_inward_rectifier_potassium_current__g_K1 *
265 var_inward_rectifier_potassium_current__xK1_inf *
266 sqrt(var_inward_rectifier_potassium_current__K_o / 5.4) *
267 (var_inward_rectifier_potassium_current__V -
268 var_inward_rectifier_potassium_current__E_K); // picoA_per_picoF
270 var_chaste_interface__transient_outward_current_s_gate__s; // dimensionless
272 var_chaste_interface__transient_outward_current_r_gate__r; // dimensionless
276 var_chaste_interface__membrane__V; // millivolt
278 var_reversal_potentials__E_K; // millivolt
280 var_transient_outward_current__g_to *
281 var_transient_outward_current__r *
282 var_transient_outward_current__s *
283 (var_transient_outward_current__V -
284 var_transient_outward_current__E_K); // picoA_per_picoF
286 var_chaste_interface__rapid_time_dependent_potassium_current_Xr1_gate__Xr1; // dimensionless
288 var_chaste_interface__rapid_time_dependent_potassium_current_Xr2_gate__Xr2; // dimensionless
290 var_chaste_interface__membrane__V; // millivolt
292 var_potassium_dynamics__K_o; // millimolar
294 var_reversal_potentials__E_K; // millivolt
296 0.153; // nanoS_per_picoF
298 var_rapid_time_dependent_potassium_current__g_Kr *
299 sqrt(var_rapid_time_dependent_potassium_current__K_o / 5.4) *
300 var_rapid_time_dependent_potassium_current__Xr1 *
301 var_rapid_time_dependent_potassium_current__Xr2 *
302 (var_rapid_time_dependent_potassium_current__V -
303 var_rapid_time_dependent_potassium_current__E_K); // picoA_per_picoF
308 var_sodium_dynamics__Na_o; // millimolar
310 var_chaste_interface__sodium_dynamics__Na_i; // millimolar
313 ((var_reversal_potentials__R * var_reversal_potentials__T) /
314 var_reversal_potentials__F) *
315 log((var_reversal_potentials__K_o +
316 (var_reversal_potentials__P_kna *
317 var_reversal_potentials__Na_o)) /
318 (var_reversal_potentials__K_i +
319 (var_reversal_potentials__P_kna *
320 var_reversal_potentials__Na_i))); // millivolt
322 var_reversal_potentials__E_Ks; // millivolt
324 var_chaste_interface__membrane__V; // millivolt
326 var_chaste_interface__slow_time_dependent_potassium_current_Xs_gate__Xs; // dimensionless
328 var_slow_time_dependent_potassium_current__g_Ks *
329 pow(var_slow_time_dependent_potassium_current__Xs, 2.0) *
330 (var_slow_time_dependent_potassium_current__V -
331 var_slow_time_dependent_potassium_current__E_Ks); // picoA_per_picoF
333 var_chaste_interface__calcium_dynamics__Ca_ss; // millimolar
335 3.98e-05; // litre_per_farad_second
337 var_chaste_interface__L_type_Ca_current_f_gate__f; // dimensionless
339 var_chaste_interface__L_type_Ca_current_d_gate__d; // dimensionless
341 var_membrane__F; // coulomb_per_millimole
343 var_chaste_interface__L_type_Ca_current_f2_gate__f2; // dimensionless
345 var_chaste_interface__L_type_Ca_current_fCass_gate__fCass; // dimensionless
347 var_chaste_interface__membrane__V; // millivolt
351 var_calcium_dynamics__Ca_o; // millimolar
353 var_membrane__R; // joule_per_mole_kelvin
355 (((var_L_type_Ca_current__g_CaL * var_L_type_Ca_current__d *
356 var_L_type_Ca_current__f * var_L_type_Ca_current__f2 *
357 var_L_type_Ca_current__fCass * 4.0 *
358 (var_L_type_Ca_current__V - 15.0) *
359 pow(var_L_type_Ca_current__F, 2.0)) /
360 (var_L_type_Ca_current__R * var_L_type_Ca_current__T)) *
361 ((0.25 * var_L_type_Ca_current__Ca_ss *
362 exp((2.0 * (var_L_type_Ca_current__V - 15.0) *
363 var_L_type_Ca_current__F) /
364 (var_L_type_Ca_current__R * var_L_type_Ca_current__T))) -
365 var_L_type_Ca_current__Ca_o)) /
366 (exp((2.0 * (var_L_type_Ca_current__V - 15.0) *
367 var_L_type_Ca_current__F) /
368 (var_L_type_Ca_current__R * var_L_type_Ca_current__T)) -
369 1.0); // picoA_per_picoF
371 var_chaste_interface__sodium_dynamics__Na_i; // millimolar
373 var_membrane__R; // joule_per_mole_kelvin
375 var_membrane__T; // kelvin
377 var_chaste_interface__membrane__V; // millivolt
379 1.0; // millimolar
381 2.724; // picoA_per_picoF
383 40.0; // millimolar
385 var_membrane__F; // coulomb_per_millimole
387 var_potassium_dynamics__K_o; // millimolar
389 ((((var_sodium_potassium_pump_current__P_NaK *
390 var_sodium_potassium_pump_current__K_o) /
391 (var_sodium_potassium_pump_current__K_o +
392 var_sodium_potassium_pump_current__K_mk)) *
393 var_sodium_potassium_pump_current__Na_i) /
394 (var_sodium_potassium_pump_current__Na_i +
395 var_sodium_potassium_pump_current__K_mNa)) /
396 (1.0 +
397 (0.1245 * exp(((-0.1) * var_sodium_potassium_pump_current__V *
398 var_sodium_potassium_pump_current__F) /
399 (var_sodium_potassium_pump_current__R *
400 var_sodium_potassium_pump_current__T))) +
401 (0.0353 *
402 exp(((-var_sodium_potassium_pump_current__V) *
403 var_sodium_potassium_pump_current__F) /
404 (var_sodium_potassium_pump_current__R *
405 var_sodium_potassium_pump_current__T)))); // picoA_per_picoF
407 var_chaste_interface__fast_sodium_current_j_gate__j; // dimensionless
409 var_chaste_interface__fast_sodium_current_h_gate__h; // dimensionless
411 14.838; // nanoS_per_picoF
413 var_chaste_interface__fast_sodium_current_m_gate__m; // dimensionless
415 var_chaste_interface__membrane__V; // millivolt
417 ((var_reversal_potentials__R * var_reversal_potentials__T) /
418 var_reversal_potentials__F) *
419 log(var_reversal_potentials__Na_o /
420 var_reversal_potentials__Na_i); // millivolt
422 var_reversal_potentials__E_Na; // millivolt
424 var_fast_sodium_current__g_Na *
425 pow(var_fast_sodium_current__m, 3.0) * var_fast_sodium_current__h *
426 var_fast_sodium_current__j *
427 (var_fast_sodium_current__V -
428 var_fast_sodium_current__E_Na); // picoA_per_picoF
430 var_chaste_interface__membrane__V; // millivolt
432 var_reversal_potentials__E_Na; // millivolt
434 0.00029; // nanoS_per_picoF
436 var_sodium_background_current__g_bna *
437 (var_sodium_background_current__V -
438 var_sodium_background_current__E_Na); // picoA_per_picoF
440 2.5; // dimensionless
442 0.35; // dimensionless
444 0.1; // dimensionless
446 1.38; // millimolar
448 1000.0; // picoA_per_picoF
450 var_membrane__F; // coulomb_per_millimole
452 var_chaste_interface__calcium_dynamics__Ca_i; // millimolar
454 var_calcium_dynamics__Ca_o; // millimolar
456 var_chaste_interface__membrane__V; // millivolt
458 var_membrane__R; // joule_per_mole_kelvin
460 87.5; // millimolar
462 var_sodium_dynamics__Na_o; // millimolar
464 var_chaste_interface__sodium_dynamics__Na_i; // millimolar
466 var_membrane__T; // kelvin
468 (var_sodium_calcium_exchanger_current__K_NaCa *
469 ((exp((var_sodium_calcium_exchanger_current__gamma *
470 var_sodium_calcium_exchanger_current__V *
471 var_sodium_calcium_exchanger_current__F) /
472 (var_sodium_calcium_exchanger_current__R *
473 var_sodium_calcium_exchanger_current__T)) *
474 pow(var_sodium_calcium_exchanger_current__Na_i, 3.0) *
475 var_sodium_calcium_exchanger_current__Ca_o) -
476 (exp(((var_sodium_calcium_exchanger_current__gamma - 1.0) *
477 var_sodium_calcium_exchanger_current__V *
478 var_sodium_calcium_exchanger_current__F) /
479 (var_sodium_calcium_exchanger_current__R *
480 var_sodium_calcium_exchanger_current__T)) *
481 pow(var_sodium_calcium_exchanger_current__Na_o, 3.0) *
482 var_sodium_calcium_exchanger_current__Ca_i *
483 var_sodium_calcium_exchanger_current__alpha))) /
484 ((pow(var_sodium_calcium_exchanger_current__Km_Nai, 3.0) +
485 pow(var_sodium_calcium_exchanger_current__Na_o, 3.0)) *
486 (var_sodium_calcium_exchanger_current__Km_Ca +
487 var_sodium_calcium_exchanger_current__Ca_o) *
488 (1.0 +
489 (var_sodium_calcium_exchanger_current__K_sat *
490 exp(((var_sodium_calcium_exchanger_current__gamma - 1.0) *
491 var_sodium_calcium_exchanger_current__V *
492 var_sodium_calcium_exchanger_current__F) /
493 (var_sodium_calcium_exchanger_current__R *
494 var_sodium_calcium_exchanger_current__T))))); // picoA_per_picoF
496 var_calcium_dynamics__Ca_o; // millimolar
498 var_chaste_interface__calcium_dynamics__Ca_i; // millimolar
500 ((0.5 * var_reversal_potentials__R * var_reversal_potentials__T) /
501 var_reversal_potentials__F) *
502 log(var_reversal_potentials__Ca_o /
503 var_reversal_potentials__Ca_i); // millivolt
505 var_reversal_potentials__E_Ca; // millivolt
507 0.000592; // nanoS_per_picoF
509 var_chaste_interface__membrane__V; // millivolt
511 var_calcium_background_current__g_bca *
512 (var_calcium_background_current__V -
513 var_calcium_background_current__E_Ca); // picoA_per_picoF
515 0.0146; // nanoS_per_picoF
517 var_chaste_interface__membrane__V; // millivolt
519 var_reversal_potentials__E_K; // millivolt
521 (var_potassium_pump_current__g_pK *
522 (var_potassium_pump_current__V -
523 var_potassium_pump_current__E_K)) /
524 (1.0 + exp((25.0 - var_potassium_pump_current__V) /
525 5.98)); // picoA_per_picoF
528 var_chaste_interface__calcium_dynamics__Ca_i; // millimolar
530 0.1238; // picoA_per_picoF
532 (var_calcium_pump_current__g_pCa * var_calcium_pump_current__Ca_i) /
533 (var_calcium_pump_current__Ca_i +
534 var_calcium_pump_current__K_pCa); // picoA_per_picoF
537 var_rapid_time_dependent_potassium_current__V; // millivolt
539 var_rapid_time_dependent_potassium_current_Xr1_gate__alpha_xr1 =
540 450.0 /
541 (1.0 +
542 exp(((-45.0) -
543 var_rapid_time_dependent_potassium_current_Xr1_gate__V) /
544 10.0)); // dimensionless
546 var_rapid_time_dependent_potassium_current_Xr1_gate__beta_xr1 =
547 6.0 /
548 (1.0 +
549 exp((var_rapid_time_dependent_potassium_current_Xr1_gate__V +
550 30.0) /
551 11.5)); // dimensionless
553 var_rapid_time_dependent_potassium_current_Xr1_gate__tau_xr1 =
554 1.0 *
555 var_rapid_time_dependent_potassium_current_Xr1_gate__alpha_xr1 *
556 var_rapid_time_dependent_potassium_current_Xr1_gate__beta_xr1; // millisecond
558 var_rapid_time_dependent_potassium_current_Xr1_gate__xr1_inf =
559 1.0 /
560 (1.0 +
561 exp(((-26.0) -
562 var_rapid_time_dependent_potassium_current_Xr1_gate__V) /
563 7.0)); // dimensionless
565 var_rapid_time_dependent_potassium_current__V; // millivolt
567 var_rapid_time_dependent_potassium_current_Xr2_gate__alpha_xr2 =
568 3.0 /
569 (1.0 +
570 exp(((-60.0) -
571 var_rapid_time_dependent_potassium_current_Xr2_gate__V) /
572 20.0)); // dimensionless
574 var_rapid_time_dependent_potassium_current_Xr2_gate__beta_xr2 =
575 1.12 /
576 (1.0 +
577 exp((var_rapid_time_dependent_potassium_current_Xr2_gate__V -
578 60.0) /
579 20.0)); // dimensionless
581 var_rapid_time_dependent_potassium_current_Xr2_gate__tau_xr2 =
582 1.0 *
583 var_rapid_time_dependent_potassium_current_Xr2_gate__alpha_xr2 *
584 var_rapid_time_dependent_potassium_current_Xr2_gate__beta_xr2; // millisecond
586 var_rapid_time_dependent_potassium_current_Xr2_gate__xr2_inf =
587 1.0 /
588 (1.0 +
589 exp((var_rapid_time_dependent_potassium_current_Xr2_gate__V +
590 88.0) /
591 24.0)); // dimensionless
593 var_slow_time_dependent_potassium_current__V; // millivolt
595 var_slow_time_dependent_potassium_current_Xs_gate__beta_xs =
596 1.0 /
597 (1.0 +
598 exp((var_slow_time_dependent_potassium_current_Xs_gate__V -
599 35.0) /
600 15.0)); // dimensionless
602 var_slow_time_dependent_potassium_current_Xs_gate__alpha_xs =
603 1400.0 /
604 sqrt(
605 1.0 +
606 exp((5.0 -
607 var_slow_time_dependent_potassium_current_Xs_gate__V) /
608 6.0)); // dimensionless
610 var_slow_time_dependent_potassium_current_Xs_gate__tau_xs =
611 (1.0 *
612 var_slow_time_dependent_potassium_current_Xs_gate__alpha_xs *
613 var_slow_time_dependent_potassium_current_Xs_gate__beta_xs) +
614 80.0; // millisecond
616 var_slow_time_dependent_potassium_current_Xs_gate__xs_inf =
617 1.0 /
618 (1.0 +
619 exp(((-5.0) -
620 var_slow_time_dependent_potassium_current_Xs_gate__V) /
621 14.0)); // dimensionless
623 var_fast_sodium_current__V; // millivolt
625 1.0 / (1.0 + exp(((-60.0) - var_fast_sodium_current_m_gate__V) /
626 5.0)); // dimensionless
628 (0.1 /
629 (1.0 + exp((var_fast_sodium_current_m_gate__V + 35.0) / 5.0))) +
630 (0.1 / (1.0 + exp((var_fast_sodium_current_m_gate__V - 50.0) /
631 200.0))); // dimensionless
633 1.0 * var_fast_sodium_current_m_gate__alpha_m *
634 var_fast_sodium_current_m_gate__beta_m; // millisecond
636 1.0 / pow(1.0 + exp(((-56.86) - var_fast_sodium_current_m_gate__V) /
637 9.03),
638 2.0); // dimensionless
640 var_fast_sodium_current__V; // millivolt
642 1.0 /
643 pow(1.0 + exp((var_fast_sodium_current_h_gate__V + 71.55) / 7.43),
644 2.0); // dimensionless
646 (var_fast_sodium_current_h_gate__V < (-40.0))
647 ? ((2.7 * exp(0.079 * var_fast_sodium_current_h_gate__V)) +
648 (310000.0 * exp(0.3485 * var_fast_sodium_current_h_gate__V)))
649 : (0.77 /
650 (0.13 *
651 (1.0 + exp((var_fast_sodium_current_h_gate__V + 10.66) /
652 (-11.1))))); // per_millisecond
654 (var_fast_sodium_current_h_gate__V < (-40.0))
655 ? (0.057 *
656 exp((-(var_fast_sodium_current_h_gate__V + 80.0)) / 6.8))
657 : 0.0; // per_millisecond
659 1.0 / (var_fast_sodium_current_h_gate__alpha_h +
660 var_fast_sodium_current_h_gate__beta_h); // millisecond
662 var_fast_sodium_current__V; // millivolt
664 1.0 /
665 pow(1.0 + exp((var_fast_sodium_current_j_gate__V + 71.55) / 7.43),
666 2.0); // dimensionless
668 (var_fast_sodium_current_j_gate__V < (-40.0))
669 ? ((((((-25428.0) *
670 exp(0.2444 * var_fast_sodium_current_j_gate__V)) -
671 (6.948e-06 *
672 exp((-0.04391) * var_fast_sodium_current_j_gate__V))) *
673 (var_fast_sodium_current_j_gate__V + 37.78)) /
674 1.0) /
675 (1.0 +
676 exp(0.311 * (var_fast_sodium_current_j_gate__V + 79.23))))
677 : 0.0; // per_millisecond
679 (var_fast_sodium_current_j_gate__V < (-40.0))
680 ? ((0.02424 *
681 exp((-0.01052) * var_fast_sodium_current_j_gate__V)) /
682 (1.0 + exp((-0.1378) *
683 (var_fast_sodium_current_j_gate__V + 40.14))))
684 : ((0.6 * exp(0.057 * var_fast_sodium_current_j_gate__V)) /
685 (1.0 + exp((-0.1) * (var_fast_sodium_current_j_gate__V +
686 32.0)))); // per_millisecond
688 1.0 / (var_fast_sodium_current_j_gate__alpha_j +
689 var_fast_sodium_current_j_gate__beta_j); // millisecond
691 var_L_type_Ca_current__V; // millivolt
693 (1.4 /
694 (1.0 + exp(((-35.0) - var_L_type_Ca_current_d_gate__V) / 13.0))) +
695 0.25; // dimensionless
697 1.0 / (1.0 + exp((50.0 - var_L_type_Ca_current_d_gate__V) /
698 20.0)); // millisecond
700 1.4 / (1.0 + exp((var_L_type_Ca_current_d_gate__V + 5.0) /
701 5.0)); // dimensionless
703 (1.0 * var_L_type_Ca_current_d_gate__alpha_d *
704 var_L_type_Ca_current_d_gate__beta_d) +
705 var_L_type_Ca_current_d_gate__gamma_d; // millisecond
707 1.0 / (1.0 + exp(((-8.0) - var_L_type_Ca_current_d_gate__V) /
708 7.5)); // dimensionless
710 var_L_type_Ca_current__V; // millivolt
712 (1102.5 *
713 exp((-pow(var_L_type_Ca_current_f_gate__V + 27.0, 2.0)) / 225.0)) +
714 (200.0 /
715 (1.0 + exp((13.0 - var_L_type_Ca_current_f_gate__V) / 10.0))) +
716 (180.0 /
717 (1.0 + exp((var_L_type_Ca_current_f_gate__V + 30.0) / 10.0))) +
718 20.0; // millisecond
720 1.0 / (1.0 + exp((var_L_type_Ca_current_f_gate__V + 20.0) /
721 7.0)); // dimensionless
723 var_L_type_Ca_current__V; // millivolt
725 (0.67 /
726 (1.0 + exp((var_L_type_Ca_current_f2_gate__V + 35.0) / 7.0))) +
727 0.33; // dimensionless
729 (562.0 * exp((-pow(var_L_type_Ca_current_f2_gate__V + 27.0, 2.0)) /
730 240.0)) +
731 (31.0 /
732 (1.0 + exp((25.0 - var_L_type_Ca_current_f2_gate__V) / 10.0))) +
733 (80.0 / (1.0 + exp((var_L_type_Ca_current_f2_gate__V + 30.0) /
734 10.0))); // millisecond
736 var_L_type_Ca_current__Ca_ss; // millimolar
738 (80.0 /
739 (1.0 + pow(var_L_type_Ca_current_fCass_gate__Ca_ss / 0.05, 2.0))) +
740 2.0; // millisecond
742 (0.6 /
743 (1.0 + pow(var_L_type_Ca_current_fCass_gate__Ca_ss / 0.05, 2.0))) +
744 0.4; // dimensionless
746 var_transient_outward_current__V; // millivolt
748 1.0 / (1.0 + exp((var_transient_outward_current_s_gate__V +
749 s_inf_factor) /
750 5.0)); // dimensionless
752 (s_tau_f1 *
753 exp((-pow(var_transient_outward_current_s_gate__V + s_tau_f2,
754 2.0)) /
755 s_tau_f3)) +
756 s_tau_f4 +
757 s_tau_f5 *
758 ((5.0 /
759 (1.0 + exp((var_transient_outward_current_s_gate__V - 20.0) /
760 5.0))) +
761 3.0); // millisecond
763 var_transient_outward_current__V; // millivolt
765 1.0 / (1.0 + exp((20.0 - var_transient_outward_current_r_gate__V) /
766 6.0)); // dimensionless
768 (9.5 *
769 exp((-pow(var_transient_outward_current_r_gate__V + 40.0, 2.0)) /
770 1800.0)) +
771 0.8; // millisecond
773 var_chaste_interface__calcium_dynamics__Ca_i; // millimolar
775 var_chaste_interface__calcium_dynamics__Ca_SR; // millimolar
777 var_chaste_interface__calcium_dynamics__Ca_ss; // millimolar
780 var_chaste_interface__calcium_dynamics__R_prime; // dimensionless
782 0.15; // per_millimolar2_per_millisecond
787 var_calcium_dynamics__max_sr -
788 ((var_calcium_dynamics__max_sr - var_calcium_dynamics__min_sr) /
789 (1.0 + pow(var_calcium_dynamics__EC / var_calcium_dynamics__Ca_SR,
790 2.0))); // dimensionless
792 var_calcium_dynamics__k1_prime /
793 var_calcium_dynamics__kcasr; // per_millimolar2_per_millisecond
796 (var_calcium_dynamics__k1 * pow(var_calcium_dynamics__Ca_ss, 2.0) *
797 var_calcium_dynamics__R_prime) /
798 (var_calcium_dynamics__k3 +
799 (var_calcium_dynamics__k1 *
800 pow(var_calcium_dynamics__Ca_ss, 2.0))); // dimensionless
802 var_calcium_dynamics__V_rel * var_calcium_dynamics__O *
803 (var_calcium_dynamics__Ca_SR -
804 var_calcium_dynamics__Ca_ss); // millimolar_per_millisecond
806 0.006375; // millimolar_per_millisecond
809 var_calcium_dynamics__Vmax_up /
810 (1.0 + (pow(var_calcium_dynamics__K_up, 2.0) /
811 pow(var_calcium_dynamics__Ca_i,
812 2.0))); // millimolar_per_millisecond
814 0.00036; // per_millisecond
816 var_calcium_dynamics__V_leak *
817 (var_calcium_dynamics__Ca_SR -
818 var_calcium_dynamics__Ca_i); // millimolar_per_millisecond
820 0.0038; // per_millisecond
822 var_calcium_dynamics__V_xfer *
823 (var_calcium_dynamics__Ca_ss -
824 var_calcium_dynamics__Ca_i); // millimolar_per_millisecond
826 0.045; // per_millimolar_per_millisecond
828 var_calcium_dynamics__k2_prime *
829 var_calcium_dynamics__kcasr; // per_millimolar_per_millisecond
834 1.0 /
835 (1.0 +
836 ((var_calcium_dynamics__Buf_c * var_calcium_dynamics__K_buf_c) /
837 pow(var_calcium_dynamics__Ca_i + var_calcium_dynamics__K_buf_c,
838 2.0))); // dimensionless
842 1.0 /
843 (1.0 +
844 ((var_calcium_dynamics__Buf_sr * var_calcium_dynamics__K_buf_sr) /
845 pow(var_calcium_dynamics__Ca_SR + var_calcium_dynamics__K_buf_sr,
846 2.0))); // dimensionless
850 1.0 /
851 (1.0 +
852 ((var_calcium_dynamics__Buf_ss * var_calcium_dynamics__K_buf_ss) /
853 pow(var_calcium_dynamics__Ca_ss + var_calcium_dynamics__K_buf_ss,
854 2.0))); // dimensionless
858 var_membrane__V_c; // micrometre3
860 var_membrane__F; // coulomb_per_millimole
863 var_L_type_Ca_current__i_CaL; // picoA_per_picoF
865 var_sodium_calcium_exchanger_current__i_NaCa; // picoA_per_picoF
867 var_calcium_pump_current__i_p_Ca; // picoA_per_picoF
869 var_calcium_background_current__i_b_Ca; // picoA_per_picoF
871 var_calcium_dynamics__Ca_i_bufc *
872 (((((var_calcium_dynamics__i_leak - var_calcium_dynamics__i_up) *
873 var_calcium_dynamics__V_sr) /
874 var_calcium_dynamics__V_c) +
875 var_calcium_dynamics__i_xfer) -
876 ((1.0 *
877 ((var_calcium_dynamics__i_b_Ca + var_calcium_dynamics__i_p_Ca) -
878 (2.0 * var_calcium_dynamics__i_NaCa)) *
879 var_calcium_dynamics__Cm) /
880 (2.0 * 1.0 * var_calcium_dynamics__V_c *
881 var_calcium_dynamics__F))); // 'millimole per litre per
882 // millisecond'
884 var_calcium_dynamics__Ca_sr_bufsr *
885 (var_calcium_dynamics__i_up -
886 (var_calcium_dynamics__i_rel +
887 var_calcium_dynamics__i_leak)); // 'millimole per litre per
888 // millisecond'
890 var_calcium_dynamics__Ca_ss_bufss *
891 (((((-1.0) * var_calcium_dynamics__i_CaL *
892 var_calcium_dynamics__Cm) /
893 (2.0 * 1.0 * var_calcium_dynamics__V_ss *
894 var_calcium_dynamics__F)) +
895 ((var_calcium_dynamics__i_rel * var_calcium_dynamics__V_sr) /
896 var_calcium_dynamics__V_ss)) -
897 ((var_calcium_dynamics__i_xfer * var_calcium_dynamics__V_c) /
898 var_calcium_dynamics__V_ss)); // 'millimole per litre per
899 // millisecond'
901 ((-var_calcium_dynamics__k2) * var_calcium_dynamics__Ca_ss *
902 var_calcium_dynamics__R_prime) +
903 (var_calcium_dynamics__k4 *
904 (1.0 - var_calcium_dynamics__R_prime)); // per_millisecond
906 var_membrane__F; // coulomb_per_millimole
909 var_membrane__V_c; // micrometre3
911 var_fast_sodium_current__i_Na; // picoA_per_picoF
913 var_sodium_calcium_exchanger_current__i_NaCa; // picoA_per_picoF
915 var_sodium_potassium_pump_current__i_NaK; // picoA_per_picoF
917 var_sodium_background_current__i_b_Na; // picoA_per_picoF
919 (((-1.0) *
920 (var_sodium_dynamics__i_Na + var_sodium_dynamics__i_b_Na +
921 (3.0 * var_sodium_dynamics__i_NaK) +
922 (3.0 * var_sodium_dynamics__i_NaCa))) /
923 (1.0 * var_sodium_dynamics__V_c * var_sodium_dynamics__F)) *
924 var_sodium_dynamics__Cm; // 'millimole per litre per millisecond'
926 var_membrane__F; // coulomb_per_millimole
929 var_membrane__V_c; // micrometre3
931 var_inward_rectifier_potassium_current__i_K1; // picoA_per_picoF
933 var_transient_outward_current__i_to; // picoA_per_picoF
935 var_sodium_potassium_pump_current__i_NaK; // picoA_per_picoF
937 var_rapid_time_dependent_potassium_current__i_Kr; // picoA_per_picoF
939 var_slow_time_dependent_potassium_current__i_Ks; // picoA_per_picoF
941 var_potassium_pump_current__i_p_K; // picoA_per_picoF
943 var_potassium_dynamics__chaste_interface__chaste_membrane_capacitance =
944 1.0; // uF_per_cm2
946 var_chaste_interface__membrane__i_Stim; // uA_per_cm2
948 var_potassium_dynamics__i_Stim_converter /
949 var_potassium_dynamics__chaste_interface__chaste_membrane_capacitance; // picoA_per_picoF
951 (((-1.0) *
952 ((var_potassium_dynamics__i_K1 + var_potassium_dynamics__i_to +
953 var_potassium_dynamics__i_Kr + var_potassium_dynamics__i_Ks +
954 var_potassium_dynamics__i_p_K +
955 var_potassium_dynamics__i_Stim) -
956 (2.0 * var_potassium_dynamics__i_NaK))) /
957 (1.0 * var_potassium_dynamics__V_c * var_potassium_dynamics__F)) *
958 var_potassium_dynamics__Cm; // 'millimole per litre per millisecond'
960 var_chaste_interface__calcium_dynamics__d_Ca_i_d_environment__time =
961 var_calcium_dynamics__d_Ca_i_d_environment__time; // millimolar_per_millisecond
963 var_chaste_interface__calcium_dynamics__d_Ca_SR_d_environment__time =
964 var_calcium_dynamics__d_Ca_SR_d_environment__time; // millimolar_per_millisecond
966 var_chaste_interface__calcium_dynamics__d_Ca_ss_d_environment__time =
967 var_calcium_dynamics__d_Ca_ss_d_environment__time; // millimolar_per_millisecond
969 var_chaste_interface__calcium_dynamics__d_R_prime_d_environment__time =
970 var_calcium_dynamics__d_R_prime_d_environment__time; // per_millisecond
972 var_chaste_interface__sodium_dynamics__d_Na_i_d_environment__time =
973 var_sodium_dynamics__d_Na_i_d_environment__time; // millimolar_per_millisecond
975 var_chaste_interface__potassium_dynamics__d_K_i_d_environment__time =
976 var_potassium_dynamics__d_K_i_d_environment__time; // millimolar_per_millisecond
978 var_chaste_interface__calcium_dynamics__d_Ca_i_d_environment__time; // 'millimole per litre per millisecond'
980 var_chaste_interface__calcium_dynamics__d_Ca_SR_d_environment__time; // 'millimole per litre per millisecond'
982 var_chaste_interface__calcium_dynamics__d_Ca_ss_d_environment__time; // 'millimole per litre per millisecond'
984 var_chaste_interface__calcium_dynamics__d_R_prime_d_environment__time; // per_millisecond
986 var_chaste_interface__sodium_dynamics__d_Na_i_d_environment__time; // 'millimole per litre per millisecond'
988 var_chaste_interface__potassium_dynamics__d_K_i_d_environment__time; // 'millimole per litre per millisecond'
989
991 var_inward_rectifier_potassium_current__i_K1; // picoA_per_picoF
993 var_transient_outward_current__i_to; // picoA_per_picoF
995 var_rapid_time_dependent_potassium_current__i_Kr; // picoA_per_picoF
997 var_slow_time_dependent_potassium_current__i_Ks; // picoA_per_picoF
999 var_L_type_Ca_current__i_CaL; // picoA_per_picoF
1001 var_sodium_potassium_pump_current__i_NaK; // picoA_per_picoF
1003 var_fast_sodium_current__i_Na; // picoA_per_picoF
1005 var_sodium_background_current__i_b_Na; // picoA_per_picoF
1007 var_sodium_calcium_exchanger_current__i_NaCa; // picoA_per_picoF
1009 var_calcium_background_current__i_b_Ca; // picoA_per_picoF
1011 var_potassium_pump_current__i_p_K; // picoA_per_picoF
1013 var_calcium_pump_current__i_p_Ca; // picoA_per_picoF
1015 var_chaste_interface__membrane__i_Stim; // uA_per_cm2
1017 var_membrane__chaste_interface__chaste_membrane_capacitance =
1018 1.0; // uF_per_cm2
1020 var_membrane__i_Stim_converter /
1021 var_membrane__chaste_interface__chaste_membrane_capacitance; // picoA_per_picoF
1023 ((-1.0) / 1.0) *
1024 (var_membrane__i_K1 + var_membrane__i_to + var_membrane__i_Kr +
1025 var_membrane__i_Ks + var_membrane__i_CaL + var_membrane__i_NaK +
1026 var_membrane__i_Na + var_membrane__i_b_Na + var_membrane__i_NaCa +
1027 var_membrane__i_b_Ca + var_membrane__i_p_K + var_membrane__i_p_Ca +
1028 var_membrane__i_Stim); // 'millivolt per millisecond'
1030 var_chaste_interface__membrane__d_V_d_environment__time =
1031 var_membrane__d_V_d_environment__time; // ___units_1
1032 d_dt_chaste_interface__membrane__V =
1033 var_chaste_interface__membrane__d_V_d_environment__time; // 'millivolt
1034 // per
1035 // millisecond'
1036 outarray[0][i] = d_dt_chaste_interface__membrane__V;
1037 outarray[1][i] =
1038 var_rapid_time_dependent_potassium_current_Xr1_gate__xr1_inf;
1039 m_gates_tau[0][i] =
1040 var_rapid_time_dependent_potassium_current_Xr1_gate__tau_xr1;
1041 outarray[2][i] =
1042 var_rapid_time_dependent_potassium_current_Xr2_gate__xr2_inf;
1043 m_gates_tau[1][i] =
1044 var_rapid_time_dependent_potassium_current_Xr2_gate__tau_xr2;
1045 outarray[3][i] =
1046 var_slow_time_dependent_potassium_current_Xs_gate__xs_inf;
1047 m_gates_tau[2][i] =
1048 var_slow_time_dependent_potassium_current_Xs_gate__tau_xs;
1049 outarray[4][i] = var_fast_sodium_current_m_gate__m_inf;
1050 m_gates_tau[3][i] = var_fast_sodium_current_m_gate__tau_m;
1051 outarray[5][i] = var_fast_sodium_current_h_gate__h_inf;
1052 m_gates_tau[4][i] = var_fast_sodium_current_h_gate__tau_h;
1053 outarray[6][i] = var_fast_sodium_current_j_gate__j_inf;
1054 m_gates_tau[5][i] = var_fast_sodium_current_j_gate__tau_j;
1055 outarray[7][i] = var_L_type_Ca_current_d_gate__d_inf;
1056 m_gates_tau[6][i] = var_L_type_Ca_current_d_gate__tau_d;
1057 outarray[8][i] = var_L_type_Ca_current_f_gate__f_inf;
1058 m_gates_tau[7][i] = var_L_type_Ca_current_f_gate__tau_f;
1059 outarray[9][i] = var_L_type_Ca_current_f2_gate__f2_inf;
1060 m_gates_tau[8][i] = var_L_type_Ca_current_f2_gate__tau_f2;
1061 outarray[10][i] = var_L_type_Ca_current_fCass_gate__fCass_inf;
1062 m_gates_tau[9][i] = var_L_type_Ca_current_fCass_gate__tau_fCass;
1063 outarray[11][i] = var_transient_outward_current_s_gate__s_inf;
1064 m_gates_tau[10][i] = var_transient_outward_current_s_gate__tau_s;
1065 outarray[12][i] = var_transient_outward_current_r_gate__r_inf;
1066 m_gates_tau[11][i] = var_transient_outward_current_r_gate__tau_r;
1067 outarray[13][i] = d_dt_chaste_interface__calcium_dynamics__Ca_i;
1068 outarray[14][i] = d_dt_chaste_interface__calcium_dynamics__Ca_SR;
1069 outarray[15][i] = d_dt_chaste_interface__calcium_dynamics__Ca_ss;
1070 outarray[16][i] = d_dt_chaste_interface__calcium_dynamics__R_prime;
1071 outarray[17][i] = d_dt_chaste_interface__sodium_dynamics__Na_i;
1072 outarray[18][i] = d_dt_chaste_interface__potassium_dynamics__K_i;
1073 }
1074}
Array< OneD, Array< OneD, NekDouble > > m_gates_tau
Storage for gate tau values.
Definition: CellModel.h:143
References g_Ks, g_to, k_0, tinysimd::log(), Nektar::CellModel::m_gates_tau, Nektar::CellModel::m_nq, s_inf_factor, s_tau_f1, s_tau_f2, s_tau_f3, s_tau_f4, s_tau_f5, and tinysimd::sqrt().
Member Data Documentation
◆ className
|
static |
Initial value:
=
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Definition: NekFactory.hpp:198
static CellModelSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField)
Creates an instance of this class.
Definition: TenTusscher06.h:47
Name of class.
Definition at line 56 of file TenTusscher06.h.
◆ def
|
staticprotected |
Initial value:
=
LibUtilities::SessionReader::RegisterDefaultSolverInfo("CellModelVariant",
"Epicardium")
static std::string RegisterDefaultSolverInfo(const std::string &pName, const std::string &pValue)
Registers the default string value of a solver info property.
Definition: SessionReader.h:617
Definition at line 98 of file TenTusscher06.h.
◆ g_Ks
|
protected |
Definition at line 79 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ g_to
|
protected |
Definition at line 78 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ k_0
|
protected |
Definition at line 86 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ lookupIds
|
staticprotected |
Initial value:
= {
static std::string RegisterEnumValue(std::string pEnum, std::string pString, int pEnumValue)
Registers an enumeration value.
Definition: SessionReader.h:583
Definition at line 97 of file TenTusscher06.h.
◆ model_variant
|
protected |
Definition at line 95 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_SetInitialConditions().
◆ s_inf_factor
|
protected |
Definition at line 80 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ s_tau_f1
|
protected |
Definition at line 81 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ s_tau_f2
|
protected |
Definition at line 82 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ s_tau_f3
|
protected |
Definition at line 83 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ s_tau_f4
|
protected |
Definition at line 84 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
◆ s_tau_f5
|
protected |
Definition at line 85 of file TenTusscher06.h.
Referenced by TenTusscher06(), and v_Update().
