#include <LuoRudy91.h>
Inheritance diagram for Nektar::LuoRudy91:
Public Member Functions | |
| LuoRudy91 (const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField) | |
| Constructor. More... | |
| ~LuoRudy91 () override | |
| Destructor. 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 Member Functions | |
| void | v_Update (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time) override |
| Computes the reaction terms $f(u,v)$ and $g(u,v)$. More... | |
| void | v_GenerateSummary (SummaryList &s) override |
| Prints a summary of the model parameters. More... | |
| void | v_SetInitialConditions () override |
| Set initial conditions for the cell model. More... | |
| 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 () |
Additional Inherited Members | |
| 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... | |
Detailed Description
Definition at line 41 of file LuoRudy91.h.
Constructor & Destructor Documentation
◆ LuoRudy91()
| Nektar::LuoRudy91::LuoRudy91 | ( | const LibUtilities::SessionReaderSharedPtr & | pSession, |
| const MultiRegions::ExpListSharedPtr & | pField | ||
| ) |
Constructor.
Definition at line 61 of file LuoRudy91.cpp.
63 : CellModel(pSession, pField)
64{
65 m_nq = pField->GetNpoints();
66
67 m_nvar = 8;
68 m_gates.push_back(1);
69 m_gates.push_back(2);
70 m_gates.push_back(3);
71 m_gates.push_back(4);
72 m_gates.push_back(5);
73 m_gates.push_back(6);
74 m_concentrations.push_back(7);
75}
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 Nektar::CellModel::m_concentrations, Nektar::CellModel::m_gates, Nektar::CellModel::m_nq, and Nektar::CellModel::m_nvar.
◆ ~LuoRudy91()
|
inlineoverride |
Member Function Documentation
◆ create()
|
inlinestatic |
Creates an instance of this class.
Definition at line 46 of file LuoRudy91.h.
49 {
51 }
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:166
References Nektar::MemoryManager< DataType >::AllocateSharedPtr().
◆ v_GenerateSummary()
|
overrideprotectedvirtual |
Prints a summary of the model parameters.
Implements Nektar::CellModel.
Definition at line 542 of file LuoRudy91.cpp.
543{
545}
void AddSummaryItem(SummaryList &l, const std::string &name, const std::string &value)
Adds a summary item to the summary info list.
Definition: Misc.cpp:47
References Nektar::SolverUtils::AddSummaryItem().
◆ v_SetInitialConditions()
|
overrideprotectedvirtual |
Set initial conditions for the cell model.
Implements Nektar::CellModel.
Definition at line 550 of file LuoRudy91.cpp.
551{
560}
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.hpp:54
References Vmath::Fill(), Nektar::CellModel::m_cellSol, and Nektar::CellModel::m_nq.
◆ v_Update()
|
overrideprotectedvirtual |
Computes the reaction terms $f(u,v)$ and $g(u,v)$.
Implements Nektar::CellModel.
Definition at line 77 of file LuoRudy91.cpp.
81{
83 {
84 // Inputs:
85 // Time units: millisecond
86 NekDouble var_chaste_interface__membrane__V = inarray[0][i];
87 // Units: millivolt; Initial value: -84.3801107371
88 NekDouble var_chaste_interface__fast_sodium_current_m_gate__m =
89 inarray[1][i];
90 // Units: dimensionless; Initial value: 0.00171338077730188
91 NekDouble var_chaste_interface__fast_sodium_current_h_gate__h =
92 inarray[2][i];
93 // Units: dimensionless; Initial value: 0.982660523699656
94 NekDouble var_chaste_interface__fast_sodium_current_j_gate__j =
95 inarray[3][i];
96 // Units: dimensionless; Initial value: 0.989108212766685
97 NekDouble var_chaste_interface__slow_inward_current_d_gate__d =
98 inarray[4][i];
99 // Units: dimensionless; Initial value: 0.00302126301779861
100 NekDouble var_chaste_interface__slow_inward_current_f_gate__f =
101 inarray[5][i];
102 // Units: dimensionless; Initial value: 0.999967936476325
103 NekDouble
104 var_chaste_interface__time_dependent_potassium_current_X_gate__X =
105 inarray[6][i];
106 // Units: dimensionless; Initial value: 0.0417603108167287
107 NekDouble
108 var_chaste_interface__intracellular_calcium_concentration__Cai =
109 inarray[7][i];
110 // Units: millimolar; Initial value: 0.00017948816388306
111
112 // Mathematics
113 NekDouble d_dt_chaste_interface__membrane__V;
115 var_chaste_interface__fast_sodium_current_j_gate__j; // dimensionless
117 var_chaste_interface__fast_sodium_current_h_gate__h; // dimensionless
119 var_chaste_interface__fast_sodium_current_m_gate__m; // dimensionless
121 var_chaste_interface__membrane__V; // millivolt
123 var_chaste_interface__slow_inward_current_d_gate__d; // dimensionless
125 var_chaste_interface__slow_inward_current_f_gate__f; // dimensionless
127 var_chaste_interface__membrane__V; // millivolt
129 var_chaste_interface__intracellular_calcium_concentration__Cai; // millimolar
131 7.7 -
134 0.09 * var_slow_inward_current__d * var_slow_inward_current__f *
135 (var_slow_inward_current__V -
136 var_slow_inward_current__E_si); // microA_per_cm2
138 var_chaste_interface__membrane__V; // millivolt
140 var_chaste_interface__time_dependent_potassium_current_X_gate__X; // dimensionless
141#if 0
142 const NekDouble var_fast_sodium_current_m_gate__m = var_fast_sodium_current__m; // dimensionless
143#endif
145 var_fast_sodium_current__V; // millivolt
147 (0.32 * (var_fast_sodium_current_m_gate__V + 47.13)) /
148 (1.0 - exp((-0.1) * (var_fast_sodium_current_m_gate__V +
149 47.13))); // per_millisecond
151 0.08 *
152 exp((-var_fast_sodium_current_m_gate__V) / 11.0); // per_millisecond
153#if 0
154 const NekDouble var_fast_sodium_current_m_gate__d_m_d_environment__time = (var_fast_sodium_current_m_gate__alpha_m * (1.0 - var_fast_sodium_current_m_gate__m)) - (var_fast_sodium_current_m_gate__beta_m * var_fast_sodium_current_m_gate__m); // per_millisecond
155 const NekDouble var_fast_sodium_current__fast_sodium_current_m_gate__d_m_d_environment__time = var_fast_sodium_current_m_gate__d_m_d_environment__time; // per_millisecond
156#endif
158 var_fast_sodium_current__V; // millivolt
160 (var_fast_sodium_current_h_gate__V < (-40.0))
161 ? ((3.56 * exp(0.079 * var_fast_sodium_current_h_gate__V)) +
162 (310000.0 * exp(0.35 * var_fast_sodium_current_h_gate__V)))
163 : (1.0 / (0.13 * (1.0 + exp((var_fast_sodium_current_h_gate__V +
164 10.66) /
165 (-11.1))))); // per_millisecond
167 (var_fast_sodium_current_h_gate__V < (-40.0))
168 ? (0.135 *
169 exp((80.0 + var_fast_sodium_current_h_gate__V) / (-6.8)))
170 : 0.0; // per_millisecond
171#if 0
172 const NekDouble var_fast_sodium_current_h_gate__h = var_fast_sodium_current__h; // dimensionless
173 const NekDouble var_fast_sodium_current_h_gate__d_h_d_environment__time = (var_fast_sodium_current_h_gate__alpha_h * (1.0 - var_fast_sodium_current_h_gate__h)) - (var_fast_sodium_current_h_gate__beta_h * var_fast_sodium_current_h_gate__h); // per_millisecond
174 const NekDouble var_fast_sodium_current__fast_sodium_current_h_gate__d_h_d_environment__time = var_fast_sodium_current_h_gate__d_h_d_environment__time; // per_millisecond
175#endif
177 var_fast_sodium_current__V; // millivolt
179 (var_fast_sodium_current_j_gate__V < (-40.0))
180 ? (((((-127140.0) *
181 exp(0.2444 * var_fast_sodium_current_j_gate__V)) -
182 (3.474e-05 *
183 exp((-0.04391) * var_fast_sodium_current_j_gate__V))) *
184 (var_fast_sodium_current_j_gate__V + 37.78)) /
185 (1.0 +
186 exp(0.311 * (var_fast_sodium_current_j_gate__V + 79.23))))
187 : 0.0; // per_millisecond
189 (var_fast_sodium_current_j_gate__V < (-40.0))
190 ? ((0.1212 *
191 exp((-0.01052) * var_fast_sodium_current_j_gate__V)) /
192 (1.0 + exp((-0.1378) *
193 (var_fast_sodium_current_j_gate__V + 40.14))))
194 : ((0.3 *
195 exp((-2.535e-07) * var_fast_sodium_current_j_gate__V)) /
196 (1.0 + exp((-0.1) * (var_fast_sodium_current_j_gate__V +
197 32.0)))); // per_millisecond
198#if 0
199 const NekDouble var_fast_sodium_current_j_gate__j = var_fast_sodium_current__j; // dimensionless
200 const NekDouble var_fast_sodium_current_j_gate__d_j_d_environment__time = (var_fast_sodium_current_j_gate__alpha_j * (1.0 - var_fast_sodium_current_j_gate__j)) - (var_fast_sodium_current_j_gate__beta_j * var_fast_sodium_current_j_gate__j); // per_millisecond
201 const NekDouble var_fast_sodium_current__fast_sodium_current_j_gate__d_j_d_environment__time = var_fast_sodium_current_j_gate__d_j_d_environment__time; // per_millisecond
202#endif
204 var_slow_inward_current__V; // millivolt
206 (0.095 * exp((-0.01) * (var_slow_inward_current_d_gate__V - 5.0))) /
207 (1.0 + exp((-0.072) * (var_slow_inward_current_d_gate__V -
208 5.0))); // per_millisecond
209#if 0
210 const NekDouble var_slow_inward_current_d_gate__d = var_slow_inward_current__d; // dimensionless
211#endif
213 (0.07 *
214 exp((-0.017) * (var_slow_inward_current_d_gate__V + 44.0))) /
215 (1.0 + exp(0.05 * (var_slow_inward_current_d_gate__V +
216 44.0))); // per_millisecond
217#if 0
218 const NekDouble var_slow_inward_current_d_gate__d_d_d_environment__time = (var_slow_inward_current_d_gate__alpha_d * (1.0 - var_slow_inward_current_d_gate__d)) - (var_slow_inward_current_d_gate__beta_d * var_slow_inward_current_d_gate__d); // per_millisecond
219 const NekDouble var_slow_inward_current__slow_inward_current_d_gate__d_d_d_environment__time = var_slow_inward_current_d_gate__d_d_d_environment__time; // per_millisecond
220 const NekDouble var_slow_inward_current_f_gate__f = var_slow_inward_current__f; // dimensionless
221#endif
223 var_slow_inward_current__V; // millivolt
225 (0.012 *
226 exp((-0.008) * (var_slow_inward_current_f_gate__V + 28.0))) /
227 (1.0 + exp(0.15 * (var_slow_inward_current_f_gate__V +
228 28.0))); // per_millisecond
230 (0.0065 *
231 exp((-0.02) * (var_slow_inward_current_f_gate__V + 30.0))) /
232 (1.0 + exp((-0.2) * (var_slow_inward_current_f_gate__V +
233 30.0))); // per_millisecond
234#if 0
235 const NekDouble var_slow_inward_current_f_gate__d_f_d_environment__time = (var_slow_inward_current_f_gate__alpha_f * (1.0 - var_slow_inward_current_f_gate__f)) - (var_slow_inward_current_f_gate__beta_f * var_slow_inward_current_f_gate__f); // per_millisecond
236 const NekDouble var_slow_inward_current__slow_inward_current_f_gate__d_f_d_environment__time = var_slow_inward_current_f_gate__d_f_d_environment__time; // per_millisecond
237 const NekDouble var_time_dependent_potassium_current_X_gate__X = var_time_dependent_potassium_current__X; // dimensionless
238#endif
240 var_time_dependent_potassium_current__V; // millivolt
242 (0.0013 *
243 exp((-0.06) *
244 (var_time_dependent_potassium_current_X_gate__V + 20.0))) /
245 (1.0 +
246 exp((-0.04) * (var_time_dependent_potassium_current_X_gate__V +
247 20.0))); // per_millisecond
249 (0.0005 *
250 exp(0.083 *
251 (var_time_dependent_potassium_current_X_gate__V + 50.0))) /
252 (1.0 + exp(0.057 * (var_time_dependent_potassium_current_X_gate__V +
253 50.0))); // per_millisecond
254#if 0
255 const NekDouble var_time_dependent_potassium_current_X_gate__d_X_d_environment__time = (var_time_dependent_potassium_current_X_gate__alpha_X * (1.0 - var_time_dependent_potassium_current_X_gate__X)) - (var_time_dependent_potassium_current_X_gate__beta_X * var_time_dependent_potassium_current_X_gate__X); // per_millisecond
256 const NekDouble var_time_dependent_potassium_current__time_dependent_potassium_current_X_gate__d_X_d_environment__time = var_time_dependent_potassium_current_X_gate__d_X_d_environment__time; // per_millisecond
257#endif
259 var_chaste_interface__intracellular_calcium_concentration__Cai; // millimolar
261 var_slow_inward_current__i_si; // microA_per_cm2
263 var_intracellular_calcium_concentration__d_Cai_d_environment__time =
264 (((-0.0001) / 1.0) *
265 var_intracellular_calcium_concentration__i_si) +
266 (0.07 *
267 (0.0001 -
268 var_intracellular_calcium_concentration__Cai)); // 'millimole
269 // per litre
270 // per
271 // millisecond'
272#if 0
273 const NekDouble var_chaste_interface__fast_sodium_current_m_gate__d_m_d_environment__time = var_fast_sodium_current__fast_sodium_current_m_gate__d_m_d_environment__time; // per_millisecond
274 const NekDouble var_chaste_interface__fast_sodium_current_h_gate__d_h_d_environment__time = var_fast_sodium_current__fast_sodium_current_h_gate__d_h_d_environment__time; // per_millisecond
275 const NekDouble var_chaste_interface__fast_sodium_current_j_gate__d_j_d_environment__time = var_fast_sodium_current__fast_sodium_current_j_gate__d_j_d_environment__time; // per_millisecond
276 const NekDouble var_chaste_interface__slow_inward_current_d_gate__d_d_d_environment__time = var_slow_inward_current__slow_inward_current_d_gate__d_d_d_environment__time; // per_millisecond
277 const NekDouble var_chaste_interface__slow_inward_current_f_gate__d_f_d_environment__time = var_slow_inward_current__slow_inward_current_f_gate__d_f_d_environment__time; // per_millisecond
278 const NekDouble var_chaste_interface__time_dependent_potassium_current_X_gate__d_X_d_environment__time = var_time_dependent_potassium_current__time_dependent_potassium_current_X_gate__d_X_d_environment__time; // per_millisecond
279#endif
281 var_chaste_interface__intracellular_calcium_concentration__d_Cai_d_environment__time =
282 var_intracellular_calcium_concentration__d_Cai_d_environment__time; // ___units_16
283#if 0
284 const NekDouble d_dt_chaste_interface__fast_sodium_current_m_gate__m = var_chaste_interface__fast_sodium_current_m_gate__d_m_d_environment__time; // per_millisecond
285 const NekDouble d_dt_chaste_interface__fast_sodium_current_h_gate__h = var_chaste_interface__fast_sodium_current_h_gate__d_h_d_environment__time; // per_millisecond
286 const NekDouble d_dt_chaste_interface__fast_sodium_current_j_gate__j = var_chaste_interface__fast_sodium_current_j_gate__d_j_d_environment__time; // per_millisecond
287 const NekDouble d_dt_chaste_interface__slow_inward_current_d_gate__d = var_chaste_interface__slow_inward_current_d_gate__d_d_d_environment__time; // per_millisecond
288 const NekDouble d_dt_chaste_interface__slow_inward_current_f_gate__f = var_chaste_interface__slow_inward_current_f_gate__d_f_d_environment__time; // per_millisecond
289 const NekDouble d_dt_chaste_interface__time_dependent_potassium_current_X_gate__X = var_chaste_interface__time_dependent_potassium_current_X_gate__d_X_d_environment__time; // per_millisecond
290#endif
292 d_dt_chaste_interface__intracellular_calcium_concentration__Cai =
293 var_chaste_interface__intracellular_calcium_concentration__d_Cai_d_environment__time; // 'millimole per litre per millisecond'
294
301 var_chaste_interface__membrane__I_stim; // microA_per_cm2
304 var_membrane__R; // joule_per_kilomole_kelvin
306 var_membrane__F; // coulomb_per_mole
309 var_ionic_concentrations__Nao; // millimolar
312 var_ionic_concentrations__Nai; // millimolar
315 ((var_fast_sodium_current__R * var_fast_sodium_current__T) /
316 var_fast_sodium_current__F) *
317 log(var_fast_sodium_current__Nao /
318 var_fast_sodium_current__Nai); // millivolt
320 var_fast_sodium_current__g_Na *
321 pow(var_fast_sodium_current__m, 3.0) * var_fast_sodium_current__h *
322 var_fast_sodium_current__j *
323 (var_fast_sodium_current__V -
324 var_fast_sodium_current__E_Na); // microA_per_cm2
326 var_fast_sodium_current__i_Na; // microA_per_cm2
328 var_slow_inward_current__i_si; // microA_per_cm2
330 var_time_dependent_potassium_current__V; // millivolt
332 (var_time_dependent_potassium_current_Xi_gate__V > (-100.0))
333 ? ((2.837 *
334 (exp(0.04 *
335 (var_time_dependent_potassium_current_Xi_gate__V +
336 77.0)) -
337 1.0)) /
338 ((var_time_dependent_potassium_current_Xi_gate__V + 77.0) *
339 exp(0.04 *
340 (var_time_dependent_potassium_current_Xi_gate__V +
341 35.0))))
342 : 1.0; // dimensionless
344 var_time_dependent_potassium_current_Xi_gate__Xi; // dimensionless
347 var_ionic_concentrations__Ko; // millimolar
349 0.282 * sqrt(var_time_dependent_potassium_current__Ko /
350 5.4); // milliS_per_cm2
352 0.01833; // dimensionless
354 var_membrane__F; // coulomb_per_mole
356 var_ionic_concentrations__Nao; // millimolar
359 var_ionic_concentrations__Ki; // millimolar
361 var_ionic_concentrations__Nai; // millimolar
363 var_membrane__T; // kelvin
365 var_membrane__R; // joule_per_kilomole_kelvin
367 ((var_time_dependent_potassium_current__R *
368 var_time_dependent_potassium_current__T) /
369 var_time_dependent_potassium_current__F) *
370 log((var_time_dependent_potassium_current__Ko +
371 (var_time_dependent_potassium_current__PR_NaK *
372 var_time_dependent_potassium_current__Nao)) /
373 (var_time_dependent_potassium_current__Ki +
374 (var_time_dependent_potassium_current__PR_NaK *
375 var_time_dependent_potassium_current__Nai))); // millivolt
377 var_time_dependent_potassium_current__g_K *
378 var_time_dependent_potassium_current__X *
379 var_time_dependent_potassium_current__Xi *
380 (var_time_dependent_potassium_current__V -
381 var_time_dependent_potassium_current__E_K); // microA_per_cm2
383 var_time_dependent_potassium_current__i_K; // microA_per_cm2
385 var_chaste_interface__membrane__V; // millivolt
387 var_time_independent_potassium_current__V; // millivolt
389 var_ionic_concentrations__Ki; // millimolar
391 var_membrane__R; // joule_per_kilomole_kelvin
393 var_membrane__F; // coulomb_per_mole
395 var_ionic_concentrations__Ko; // millimolar
397 var_membrane__T; // kelvin
399 ((var_time_independent_potassium_current__R *
400 var_time_independent_potassium_current__T) /
401 var_time_independent_potassium_current__F) *
402 log(var_time_independent_potassium_current__Ko /
403 var_time_independent_potassium_current__Ki); // millivolt
405 var_time_independent_potassium_current__E_K1; // millivolt
407 var_time_independent_potassium_current_K1_gate__beta_K1 =
408 ((0.49124 *
409 exp(0.08032 *
410 ((var_time_independent_potassium_current_K1_gate__V +
411 5.476) -
412 var_time_independent_potassium_current_K1_gate__E_K1))) +
413 (1.0 *
414 exp(0.06175 *
415 (var_time_independent_potassium_current_K1_gate__V -
416 (var_time_independent_potassium_current_K1_gate__E_K1 +
417 594.31))))) /
418 (1.0 +
419 exp((-0.5143) *
420 ((var_time_independent_potassium_current_K1_gate__V -
421 var_time_independent_potassium_current_K1_gate__E_K1) +
422 4.753))); // per_millisecond
424 var_time_independent_potassium_current_K1_gate__alpha_K1 =
425 1.02 /
426 (1.0 +
427 exp(0.2385 *
428 ((var_time_independent_potassium_current_K1_gate__V -
429 var_time_independent_potassium_current_K1_gate__E_K1) -
430 59.215))); // per_millisecond
432 var_time_independent_potassium_current_K1_gate__K1_infinity =
433 var_time_independent_potassium_current_K1_gate__alpha_K1 /
434 (var_time_independent_potassium_current_K1_gate__alpha_K1 +
435 var_time_independent_potassium_current_K1_gate__beta_K1); // dimensionless
437 var_time_independent_potassium_current_K1_gate__K1_infinity; // dimensionless
439 0.6047 * sqrt(var_time_independent_potassium_current__Ko /
440 5.4); // milliS_per_cm2
442 var_time_independent_potassium_current__g_K1 *
443 var_time_independent_potassium_current__K1_infinity *
444 (var_time_independent_potassium_current__V -
445 var_time_independent_potassium_current__E_K1); // microA_per_cm2
447 var_time_independent_potassium_current__i_K1; // microA_per_cm2
449 0.0183; // milliS_per_cm2
451 var_chaste_interface__membrane__V; // millivolt
453 1.0 / (1.0 + exp((7.488 - var_plateau_potassium_current__V) /
454 5.98)); // dimensionless
456 var_time_independent_potassium_current__E_K1; // millivolt
458 var_plateau_potassium_current__E_K1; // millivolt
460 var_plateau_potassium_current__g_Kp *
461 var_plateau_potassium_current__Kp *
462 (var_plateau_potassium_current__V -
463 var_plateau_potassium_current__E_Kp); // microA_per_cm2
465 var_plateau_potassium_current__i_Kp; // microA_per_cm2
469 var_chaste_interface__membrane__V; // millivolt
471 var_background_current__g_b *
472 (var_background_current__V -
473 var_background_current__E_b); // microA_per_cm2
475 var_background_current__i_b; // microA_per_cm2
477 ((-1.0) / var_membrane__C) *
478 (var_membrane__I_stim + var_membrane__i_Na + var_membrane__i_si +
479 var_membrane__i_K + var_membrane__i_K1 + var_membrane__i_Kp +
480 var_membrane__i_b); // 'millivolt per millisecond'
482 var_chaste_interface__membrane__d_V_d_environment__time =
483 var_membrane__d_V_d_environment__time; // ___units_1
484 d_dt_chaste_interface__membrane__V =
485 var_chaste_interface__membrane__d_V_d_environment__time; // 'millivolt
486 // per
487 // millisecond'
489 (var_fast_sodium_current_m_gate__alpha_m +
490 var_fast_sodium_current_m_gate__beta_m);
492 var_fast_sodium_current_m_gate__beta_m);
494 (var_fast_sodium_current_h_gate__alpha_h +
495 var_fast_sodium_current_h_gate__beta_h);
497 var_fast_sodium_current_h_gate__beta_h);
499 (var_fast_sodium_current_j_gate__alpha_j +
500 var_fast_sodium_current_j_gate__beta_j);
502 var_fast_sodium_current_j_gate__beta_j);
504 (var_slow_inward_current_d_gate__alpha_d +
505 var_slow_inward_current_d_gate__beta_d);
507 var_slow_inward_current_d_gate__beta_d);
509 (var_slow_inward_current_f_gate__alpha_f +
510 var_slow_inward_current_f_gate__beta_f);
512 var_slow_inward_current_f_gate__beta_f);
514 var_time_dependent_potassium_current_X_gate__alpha_X /
515 (var_time_dependent_potassium_current_X_gate__alpha_X +
516 var_time_dependent_potassium_current_X_gate__beta_X);
518 1.0 / (var_time_dependent_potassium_current_X_gate__alpha_X +
519 var_time_dependent_potassium_current_X_gate__beta_X);
520
521 outarray[0][i] = d_dt_chaste_interface__membrane__V;
522 outarray[1][i] = m_inf;
523 m_gates_tau[0][i] = m_tau;
524 outarray[2][i] = h_inf;
525 m_gates_tau[1][i] = h_tau;
526 outarray[3][i] = j_inf;
527 m_gates_tau[2][i] = j_tau;
528 outarray[4][i] = d_inf;
529 m_gates_tau[3][i] = d_tau;
530 outarray[5][i] = f_inf;
531 m_gates_tau[4][i] = f_tau;
532 outarray[6][i] = X_inf;
533 m_gates_tau[5][i] = X_tau;
534 outarray[7][i] =
535 d_dt_chaste_interface__intracellular_calcium_concentration__Cai;
536 }
537}
Array< OneD, Array< OneD, NekDouble > > m_gates_tau
Storage for gate tau values.
Definition: CellModel.h:143
References tinysimd::log(), Nektar::CellModel::m_gates_tau, Nektar::CellModel::m_nq, 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:197
static CellModelSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField)
Creates an instance of this class.
Definition: LuoRudy91.h:46
Name of class.
Definition at line 54 of file LuoRudy91.h.
