#include <LuoRudy91.h>
Inheritance diagram for Nektar::LuoRudy91:
Public Member Functions | |
| LuoRudy91 (const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField) | |
| Constructor. More... | |
| virtual | ~LuoRudy91 () |
| 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 | |
| virtual 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... | |
| virtual void | v_GenerateSummary (SummaryList &s) override |
| Prints a summary of the model parameters. More... | |
| virtual 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()
|
inlinevirtual |
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:168
References Nektar::MemoryManager< DataType >::AllocateSharedPtr().
◆ v_GenerateSummary()
|
overrideprotectedvirtual |
Prints a summary of the model parameters.
Implements Nektar::CellModel.
Definition at line 543 of file LuoRudy91.cpp.
544{
546}
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 |
Set initial conditions for the cell model.
Implements Nektar::CellModel.
Definition at line 551 of file LuoRudy91.cpp.
552{
561}
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 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.
80{
81 boost::ignore_unused(time);
82
84 {
85 // Inputs:
86 // Time units: millisecond
87 NekDouble var_chaste_interface__membrane__V = inarray[0][i];
88 // Units: millivolt; Initial value: -84.3801107371
89 NekDouble var_chaste_interface__fast_sodium_current_m_gate__m =
90 inarray[1][i];
91 // Units: dimensionless; Initial value: 0.00171338077730188
92 NekDouble var_chaste_interface__fast_sodium_current_h_gate__h =
93 inarray[2][i];
94 // Units: dimensionless; Initial value: 0.982660523699656
95 NekDouble var_chaste_interface__fast_sodium_current_j_gate__j =
96 inarray[3][i];
97 // Units: dimensionless; Initial value: 0.989108212766685
98 NekDouble var_chaste_interface__slow_inward_current_d_gate__d =
99 inarray[4][i];
100 // Units: dimensionless; Initial value: 0.00302126301779861
101 NekDouble var_chaste_interface__slow_inward_current_f_gate__f =
102 inarray[5][i];
103 // Units: dimensionless; Initial value: 0.999967936476325
104 NekDouble
105 var_chaste_interface__time_dependent_potassium_current_X_gate__X =
106 inarray[6][i];
107 // Units: dimensionless; Initial value: 0.0417603108167287
108 NekDouble
109 var_chaste_interface__intracellular_calcium_concentration__Cai =
110 inarray[7][i];
111 // Units: millimolar; Initial value: 0.00017948816388306
112
113 // Mathematics
114 NekDouble d_dt_chaste_interface__membrane__V;
116 var_chaste_interface__fast_sodium_current_j_gate__j; // dimensionless
118 var_chaste_interface__fast_sodium_current_h_gate__h; // dimensionless
120 var_chaste_interface__fast_sodium_current_m_gate__m; // dimensionless
122 var_chaste_interface__membrane__V; // millivolt
124 var_chaste_interface__slow_inward_current_d_gate__d; // dimensionless
126 var_chaste_interface__slow_inward_current_f_gate__f; // dimensionless
128 var_chaste_interface__membrane__V; // millivolt
130 var_chaste_interface__intracellular_calcium_concentration__Cai; // millimolar
132 7.7 -
135 0.09 * var_slow_inward_current__d * var_slow_inward_current__f *
136 (var_slow_inward_current__V -
137 var_slow_inward_current__E_si); // microA_per_cm2
139 var_chaste_interface__membrane__V; // millivolt
141 var_chaste_interface__time_dependent_potassium_current_X_gate__X; // dimensionless
142#if 0
143 const NekDouble var_fast_sodium_current_m_gate__m = var_fast_sodium_current__m; // dimensionless
144#endif
146 var_fast_sodium_current__V; // millivolt
148 (0.32 * (var_fast_sodium_current_m_gate__V + 47.13)) /
149 (1.0 - exp((-0.1) * (var_fast_sodium_current_m_gate__V +
150 47.13))); // per_millisecond
152 0.08 *
153 exp((-var_fast_sodium_current_m_gate__V) / 11.0); // per_millisecond
154#if 0
155 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
156 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
157#endif
159 var_fast_sodium_current__V; // millivolt
161 (var_fast_sodium_current_h_gate__V < (-40.0))
162 ? ((3.56 * exp(0.079 * var_fast_sodium_current_h_gate__V)) +
163 (310000.0 * exp(0.35 * var_fast_sodium_current_h_gate__V)))
164 : (1.0 / (0.13 * (1.0 + exp((var_fast_sodium_current_h_gate__V +
165 10.66) /
166 (-11.1))))); // per_millisecond
168 (var_fast_sodium_current_h_gate__V < (-40.0))
169 ? (0.135 *
170 exp((80.0 + var_fast_sodium_current_h_gate__V) / (-6.8)))
171 : 0.0; // per_millisecond
172#if 0
173 const NekDouble var_fast_sodium_current_h_gate__h = var_fast_sodium_current__h; // dimensionless
174 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
175 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
176#endif
178 var_fast_sodium_current__V; // millivolt
180 (var_fast_sodium_current_j_gate__V < (-40.0))
181 ? (((((-127140.0) *
182 exp(0.2444 * var_fast_sodium_current_j_gate__V)) -
183 (3.474e-05 *
184 exp((-0.04391) * var_fast_sodium_current_j_gate__V))) *
185 (var_fast_sodium_current_j_gate__V + 37.78)) /
186 (1.0 +
187 exp(0.311 * (var_fast_sodium_current_j_gate__V + 79.23))))
188 : 0.0; // per_millisecond
190 (var_fast_sodium_current_j_gate__V < (-40.0))
191 ? ((0.1212 *
192 exp((-0.01052) * var_fast_sodium_current_j_gate__V)) /
193 (1.0 + exp((-0.1378) *
194 (var_fast_sodium_current_j_gate__V + 40.14))))
195 : ((0.3 *
196 exp((-2.535e-07) * var_fast_sodium_current_j_gate__V)) /
197 (1.0 + exp((-0.1) * (var_fast_sodium_current_j_gate__V +
198 32.0)))); // per_millisecond
199#if 0
200 const NekDouble var_fast_sodium_current_j_gate__j = var_fast_sodium_current__j; // dimensionless
201 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
202 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
203#endif
205 var_slow_inward_current__V; // millivolt
207 (0.095 * exp((-0.01) * (var_slow_inward_current_d_gate__V - 5.0))) /
208 (1.0 + exp((-0.072) * (var_slow_inward_current_d_gate__V -
209 5.0))); // per_millisecond
210#if 0
211 const NekDouble var_slow_inward_current_d_gate__d = var_slow_inward_current__d; // dimensionless
212#endif
214 (0.07 *
215 exp((-0.017) * (var_slow_inward_current_d_gate__V + 44.0))) /
216 (1.0 + exp(0.05 * (var_slow_inward_current_d_gate__V +
217 44.0))); // per_millisecond
218#if 0
219 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
220 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
221 const NekDouble var_slow_inward_current_f_gate__f = var_slow_inward_current__f; // dimensionless
222#endif
224 var_slow_inward_current__V; // millivolt
226 (0.012 *
227 exp((-0.008) * (var_slow_inward_current_f_gate__V + 28.0))) /
228 (1.0 + exp(0.15 * (var_slow_inward_current_f_gate__V +
229 28.0))); // per_millisecond
231 (0.0065 *
232 exp((-0.02) * (var_slow_inward_current_f_gate__V + 30.0))) /
233 (1.0 + exp((-0.2) * (var_slow_inward_current_f_gate__V +
234 30.0))); // per_millisecond
235#if 0
236 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
237 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
238 const NekDouble var_time_dependent_potassium_current_X_gate__X = var_time_dependent_potassium_current__X; // dimensionless
239#endif
241 var_time_dependent_potassium_current__V; // millivolt
243 (0.0013 *
244 exp((-0.06) *
245 (var_time_dependent_potassium_current_X_gate__V + 20.0))) /
246 (1.0 +
247 exp((-0.04) * (var_time_dependent_potassium_current_X_gate__V +
248 20.0))); // per_millisecond
250 (0.0005 *
251 exp(0.083 *
252 (var_time_dependent_potassium_current_X_gate__V + 50.0))) /
253 (1.0 + exp(0.057 * (var_time_dependent_potassium_current_X_gate__V +
254 50.0))); // per_millisecond
255#if 0
256 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
257 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
258#endif
260 var_chaste_interface__intracellular_calcium_concentration__Cai; // millimolar
262 var_slow_inward_current__i_si; // microA_per_cm2
264 var_intracellular_calcium_concentration__d_Cai_d_environment__time =
265 (((-0.0001) / 1.0) *
266 var_intracellular_calcium_concentration__i_si) +
267 (0.07 *
268 (0.0001 -
269 var_intracellular_calcium_concentration__Cai)); // 'millimole
270 // per litre
271 // per
272 // millisecond'
273#if 0
274 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
275 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
276 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
277 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
278 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
279 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
280#endif
282 var_chaste_interface__intracellular_calcium_concentration__d_Cai_d_environment__time =
283 var_intracellular_calcium_concentration__d_Cai_d_environment__time; // ___units_16
284#if 0
285 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
286 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
287 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
288 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
289 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
290 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
291#endif
293 d_dt_chaste_interface__intracellular_calcium_concentration__Cai =
294 var_chaste_interface__intracellular_calcium_concentration__d_Cai_d_environment__time; // 'millimole per litre per millisecond'
295
302 var_chaste_interface__membrane__I_stim; // microA_per_cm2
305 var_membrane__R; // joule_per_kilomole_kelvin
307 var_membrane__F; // coulomb_per_mole
310 var_ionic_concentrations__Nao; // millimolar
313 var_ionic_concentrations__Nai; // millimolar
316 ((var_fast_sodium_current__R * var_fast_sodium_current__T) /
317 var_fast_sodium_current__F) *
318 log(var_fast_sodium_current__Nao /
319 var_fast_sodium_current__Nai); // millivolt
321 var_fast_sodium_current__g_Na *
322 pow(var_fast_sodium_current__m, 3.0) * var_fast_sodium_current__h *
323 var_fast_sodium_current__j *
324 (var_fast_sodium_current__V -
325 var_fast_sodium_current__E_Na); // microA_per_cm2
327 var_fast_sodium_current__i_Na; // microA_per_cm2
329 var_slow_inward_current__i_si; // microA_per_cm2
331 var_time_dependent_potassium_current__V; // millivolt
333 (var_time_dependent_potassium_current_Xi_gate__V > (-100.0))
334 ? ((2.837 *
335 (exp(0.04 *
336 (var_time_dependent_potassium_current_Xi_gate__V +
337 77.0)) -
338 1.0)) /
339 ((var_time_dependent_potassium_current_Xi_gate__V + 77.0) *
340 exp(0.04 *
341 (var_time_dependent_potassium_current_Xi_gate__V +
342 35.0))))
343 : 1.0; // dimensionless
345 var_time_dependent_potassium_current_Xi_gate__Xi; // dimensionless
348 var_ionic_concentrations__Ko; // millimolar
350 0.282 * sqrt(var_time_dependent_potassium_current__Ko /
351 5.4); // milliS_per_cm2
353 0.01833; // dimensionless
355 var_membrane__F; // coulomb_per_mole
357 var_ionic_concentrations__Nao; // millimolar
360 var_ionic_concentrations__Ki; // millimolar
362 var_ionic_concentrations__Nai; // millimolar
364 var_membrane__T; // kelvin
366 var_membrane__R; // joule_per_kilomole_kelvin
368 ((var_time_dependent_potassium_current__R *
369 var_time_dependent_potassium_current__T) /
370 var_time_dependent_potassium_current__F) *
371 log((var_time_dependent_potassium_current__Ko +
372 (var_time_dependent_potassium_current__PR_NaK *
373 var_time_dependent_potassium_current__Nao)) /
374 (var_time_dependent_potassium_current__Ki +
375 (var_time_dependent_potassium_current__PR_NaK *
376 var_time_dependent_potassium_current__Nai))); // millivolt
378 var_time_dependent_potassium_current__g_K *
379 var_time_dependent_potassium_current__X *
380 var_time_dependent_potassium_current__Xi *
381 (var_time_dependent_potassium_current__V -
382 var_time_dependent_potassium_current__E_K); // microA_per_cm2
384 var_time_dependent_potassium_current__i_K; // microA_per_cm2
386 var_chaste_interface__membrane__V; // millivolt
388 var_time_independent_potassium_current__V; // millivolt
390 var_ionic_concentrations__Ki; // millimolar
392 var_membrane__R; // joule_per_kilomole_kelvin
394 var_membrane__F; // coulomb_per_mole
396 var_ionic_concentrations__Ko; // millimolar
398 var_membrane__T; // kelvin
400 ((var_time_independent_potassium_current__R *
401 var_time_independent_potassium_current__T) /
402 var_time_independent_potassium_current__F) *
403 log(var_time_independent_potassium_current__Ko /
404 var_time_independent_potassium_current__Ki); // millivolt
406 var_time_independent_potassium_current__E_K1; // millivolt
408 var_time_independent_potassium_current_K1_gate__beta_K1 =
409 ((0.49124 *
410 exp(0.08032 *
411 ((var_time_independent_potassium_current_K1_gate__V +
412 5.476) -
413 var_time_independent_potassium_current_K1_gate__E_K1))) +
414 (1.0 *
415 exp(0.06175 *
416 (var_time_independent_potassium_current_K1_gate__V -
417 (var_time_independent_potassium_current_K1_gate__E_K1 +
418 594.31))))) /
419 (1.0 +
420 exp((-0.5143) *
421 ((var_time_independent_potassium_current_K1_gate__V -
422 var_time_independent_potassium_current_K1_gate__E_K1) +
423 4.753))); // per_millisecond
425 var_time_independent_potassium_current_K1_gate__alpha_K1 =
426 1.02 /
427 (1.0 +
428 exp(0.2385 *
429 ((var_time_independent_potassium_current_K1_gate__V -
430 var_time_independent_potassium_current_K1_gate__E_K1) -
431 59.215))); // per_millisecond
433 var_time_independent_potassium_current_K1_gate__K1_infinity =
434 var_time_independent_potassium_current_K1_gate__alpha_K1 /
435 (var_time_independent_potassium_current_K1_gate__alpha_K1 +
436 var_time_independent_potassium_current_K1_gate__beta_K1); // dimensionless
438 var_time_independent_potassium_current_K1_gate__K1_infinity; // dimensionless
440 0.6047 * sqrt(var_time_independent_potassium_current__Ko /
441 5.4); // milliS_per_cm2
443 var_time_independent_potassium_current__g_K1 *
444 var_time_independent_potassium_current__K1_infinity *
445 (var_time_independent_potassium_current__V -
446 var_time_independent_potassium_current__E_K1); // microA_per_cm2
448 var_time_independent_potassium_current__i_K1; // microA_per_cm2
450 0.0183; // milliS_per_cm2
452 var_chaste_interface__membrane__V; // millivolt
454 1.0 / (1.0 + exp((7.488 - var_plateau_potassium_current__V) /
455 5.98)); // dimensionless
457 var_time_independent_potassium_current__E_K1; // millivolt
459 var_plateau_potassium_current__E_K1; // millivolt
461 var_plateau_potassium_current__g_Kp *
462 var_plateau_potassium_current__Kp *
463 (var_plateau_potassium_current__V -
464 var_plateau_potassium_current__E_Kp); // microA_per_cm2
466 var_plateau_potassium_current__i_Kp; // microA_per_cm2
470 var_chaste_interface__membrane__V; // millivolt
472 var_background_current__g_b *
473 (var_background_current__V -
474 var_background_current__E_b); // microA_per_cm2
476 var_background_current__i_b; // microA_per_cm2
478 ((-1.0) / var_membrane__C) *
479 (var_membrane__I_stim + var_membrane__i_Na + var_membrane__i_si +
480 var_membrane__i_K + var_membrane__i_K1 + var_membrane__i_Kp +
481 var_membrane__i_b); // 'millivolt per millisecond'
483 var_chaste_interface__membrane__d_V_d_environment__time =
484 var_membrane__d_V_d_environment__time; // ___units_1
485 d_dt_chaste_interface__membrane__V =
486 var_chaste_interface__membrane__d_V_d_environment__time; // 'millivolt
487 // per
488 // millisecond'
490 (var_fast_sodium_current_m_gate__alpha_m +
491 var_fast_sodium_current_m_gate__beta_m);
493 var_fast_sodium_current_m_gate__beta_m);
495 (var_fast_sodium_current_h_gate__alpha_h +
496 var_fast_sodium_current_h_gate__beta_h);
498 var_fast_sodium_current_h_gate__beta_h);
500 (var_fast_sodium_current_j_gate__alpha_j +
501 var_fast_sodium_current_j_gate__beta_j);
503 var_fast_sodium_current_j_gate__beta_j);
505 (var_slow_inward_current_d_gate__alpha_d +
506 var_slow_inward_current_d_gate__beta_d);
508 var_slow_inward_current_d_gate__beta_d);
510 (var_slow_inward_current_f_gate__alpha_f +
511 var_slow_inward_current_f_gate__beta_f);
513 var_slow_inward_current_f_gate__beta_f);
515 var_time_dependent_potassium_current_X_gate__alpha_X /
516 (var_time_dependent_potassium_current_X_gate__alpha_X +
517 var_time_dependent_potassium_current_X_gate__beta_X);
519 1.0 / (var_time_dependent_potassium_current_X_gate__alpha_X +
520 var_time_dependent_potassium_current_X_gate__beta_X);
521
522 outarray[0][i] = d_dt_chaste_interface__membrane__V;
523 outarray[1][i] = m_inf;
524 m_gates_tau[0][i] = m_tau;
525 outarray[2][i] = h_inf;
526 m_gates_tau[1][i] = h_tau;
527 outarray[3][i] = j_inf;
528 m_gates_tau[2][i] = j_tau;
529 outarray[4][i] = d_inf;
530 m_gates_tau[3][i] = d_tau;
531 outarray[5][i] = f_inf;
532 m_gates_tau[4][i] = f_tau;
533 outarray[6][i] = X_inf;
534 m_gates_tau[5][i] = X_tau;
535 outarray[7][i] =
536 d_dt_chaste_interface__intracellular_calcium_concentration__Cai;
537 }
538}
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:198
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.
