46 "Pandit-Giles-Demir 2003 cell model.");
65 int nvariables = inarray.num_elements();
67 for (
unsigned int i = 0; i < nq; ++i)
71 NekDouble var_chaste_interface__membrane__V = inarray[0][i];
73 NekDouble var_chaste_interface__sodium_current_m_gate__m = inarray[2][i];
75 NekDouble var_chaste_interface__sodium_current_h_gate__h = inarray[3][i];
77 NekDouble var_chaste_interface__sodium_current_j_gate__j = inarray[4][i];
79 NekDouble var_chaste_interface__L_type_Ca_channel_d_gate__d = inarray[5][i];
81 NekDouble var_chaste_interface__L_type_Ca_channel_f_11_gate__f_11 = inarray[6][i];
83 NekDouble var_chaste_interface__L_type_Ca_channel_f_12_gate__f_12 = inarray[7][i];
85 NekDouble var_chaste_interface__L_type_Ca_channel_Ca_inact_gate__Ca_inact = inarray[8][i];
87 NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__r = inarray[9][i];
89 NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__s = inarray[10][i];
91 NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__s_slow = inarray[11][i];
93 NekDouble var_chaste_interface__steady_state_outward_K_current_r_ss_gate__r_ss = inarray[12][i];
95 NekDouble var_chaste_interface__steady_state_outward_K_current_s_ss_gate__s_ss = inarray[13][i];
97 NekDouble var_chaste_interface__hyperpolarisation_activated_current_y_gate__y = inarray[14][i];
99 NekDouble var_chaste_interface__SR_Ca_release_channel__P_O1 = inarray[15][i];
101 NekDouble var_chaste_interface__SR_Ca_release_channel__P_O2 = inarray[16][i];
103 NekDouble var_chaste_interface__SR_Ca_release_channel__P_C1 = inarray[17][i];
105 NekDouble var_chaste_interface__SR_Ca_release_channel__P_C2 = inarray[18][i];
107 NekDouble var_chaste_interface__intracellular_and_SR_Ca_fluxes__HTRPNCa = inarray[19][i];
109 NekDouble var_chaste_interface__intracellular_and_SR_Ca_fluxes__LTRPNCa = inarray[20][i];
111 NekDouble var_chaste_interface__intracellular_ion_concentrations__Na_i = inarray[21][i];
113 NekDouble var_chaste_interface__intracellular_ion_concentrations__Ca_i = inarray[22][i];
115 NekDouble var_chaste_interface__intracellular_ion_concentrations__K_i = inarray[23][i];
117 NekDouble var_chaste_interface__intracellular_ion_concentrations__Ca_ss = inarray[24][i];
119 NekDouble var_chaste_interface__intracellular_ion_concentrations__Ca_JSR = inarray[25][i];
121 NekDouble var_chaste_interface__intracellular_ion_concentrations__Ca_NSR = inarray[26][i];
126 NekDouble d_dt_chaste_interface__membrane__V;
127 const NekDouble var_membrane__R = 8314.5;
129 const NekDouble var_membrane__F = 96487.0;
130 const NekDouble var_sodium_current__F = var_membrane__F;
131 const NekDouble var_standard_ionic_concentrations__Na_o = 140.0;
132 const NekDouble var_sodium_current__Na_o = var_standard_ionic_concentrations__Na_o;
133 const NekDouble var_sodium_current__T = var_membrane__T;
134 const NekDouble var_sodium_current__R = var_membrane__R;
135 const NekDouble var_sodium_current__Na_i = var_chaste_interface__intracellular_ion_concentrations__Na_i;
136 const NekDouble var_sodium_current__E_Na = ((var_sodium_current__R * var_sodium_current__T) / var_sodium_current__F) * log(var_sodium_current__Na_o / var_sodium_current__Na_i);
137 const NekDouble var_sodium_current__m = var_chaste_interface__sodium_current_m_gate__m;
138 const NekDouble var_sodium_current__j = var_chaste_interface__sodium_current_j_gate__j;
139 const NekDouble var_sodium_current__h = var_chaste_interface__sodium_current_h_gate__h;
140 const NekDouble var_sodium_current__V = var_chaste_interface__membrane__V;
141 const NekDouble var_sodium_current__g_Na = 1.064;
142 const NekDouble var_sodium_current__i_Na = var_sodium_current__g_Na * pow(var_sodium_current__m, 3.0) * var_sodium_current__h * var_sodium_current__j * (var_sodium_current__V - var_sodium_current__E_Na);
143 const NekDouble var_L_type_Ca_channel__d = var_chaste_interface__L_type_Ca_channel_d_gate__d;
144 const NekDouble var_L_type_Ca_channel__E_Ca_L = 65.0;
145 const NekDouble var_L_type_Ca_channel__g_Ca_L_normal = 0.0341;
146 const NekDouble var_membrane__Diabetes = 0.0;
147 const NekDouble var_L_type_Ca_channel__Diabetes = var_membrane__Diabetes;
148 const NekDouble var_L_type_Ca_channel__g_Ca_L = (var_L_type_Ca_channel__Diabetes == 0.0) ? var_L_type_Ca_channel__g_Ca_L_normal : (0.76 * var_L_type_Ca_channel__g_Ca_L_normal);
149 const NekDouble var_L_type_Ca_channel__f_12 = var_chaste_interface__L_type_Ca_channel_f_12_gate__f_12;
150 const NekDouble var_L_type_Ca_channel__f_11 = var_chaste_interface__L_type_Ca_channel_f_11_gate__f_11;
151 const NekDouble var_L_type_Ca_channel__V = var_chaste_interface__membrane__V;
152 const NekDouble var_L_type_Ca_channel__Ca_inact = var_chaste_interface__L_type_Ca_channel_Ca_inact_gate__Ca_inact;
153 const NekDouble var_L_type_Ca_channel__i_Ca_L = var_L_type_Ca_channel__g_Ca_L * var_L_type_Ca_channel__d * (((0.9 + (var_L_type_Ca_channel__Ca_inact / 10.0)) * var_L_type_Ca_channel__f_11) + ((0.1 - (var_L_type_Ca_channel__Ca_inact / 10.0)) * var_L_type_Ca_channel__f_12)) * (var_L_type_Ca_channel__V - var_L_type_Ca_channel__E_Ca_L);
154 const NekDouble var_Ca_independent_transient_outward_K_current__Diabetes = var_membrane__Diabetes;
155 const NekDouble var_Ca_independent_transient_outward_K_current__b = (var_Ca_independent_transient_outward_K_current__Diabetes == 0.0) ? 0.114 : 0.31;
156 const NekDouble var_Ca_independent_transient_outward_K_current__a = 1.0 - var_Ca_independent_transient_outward_K_current__b;
157 const NekDouble var_Ca_independent_transient_outward_K_current__g_t_normal = 0.04375;
158 const NekDouble var_Ca_independent_transient_outward_K_current__g_t = (var_Ca_independent_transient_outward_K_current__Diabetes == 0.0) ? var_Ca_independent_transient_outward_K_current__g_t_normal : (0.68 * var_Ca_independent_transient_outward_K_current__g_t_normal);
159 const NekDouble var_Ca_independent_transient_outward_K_current__V = var_chaste_interface__membrane__V;
160 const NekDouble var_Ca_independent_transient_outward_K_current__r = var_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__r;
161 const NekDouble var_Ca_independent_transient_outward_K_current__s = var_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__s;
162 const NekDouble var_Ca_independent_transient_outward_K_current__F = var_membrane__F;
163 const NekDouble var_Ca_independent_transient_outward_K_current__T = var_membrane__T;
164 const NekDouble var_Ca_independent_transient_outward_K_current__K_i = var_chaste_interface__intracellular_ion_concentrations__K_i;
165 const NekDouble var_Ca_independent_transient_outward_K_current__R = var_membrane__R;
166 const NekDouble var_standard_ionic_concentrations__K_o = 5.4;
167 const NekDouble var_Ca_independent_transient_outward_K_current__K_o = var_standard_ionic_concentrations__K_o;
168 const NekDouble var_Ca_independent_transient_outward_K_current__E_K = ((var_Ca_independent_transient_outward_K_current__R * var_Ca_independent_transient_outward_K_current__T) / var_Ca_independent_transient_outward_K_current__F) * log(var_Ca_independent_transient_outward_K_current__K_o / var_Ca_independent_transient_outward_K_current__K_i);
169 const NekDouble var_Ca_independent_transient_outward_K_current__s_slow = var_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__s_slow;
170 const NekDouble var_Ca_independent_transient_outward_K_current__i_t = var_Ca_independent_transient_outward_K_current__g_t * var_Ca_independent_transient_outward_K_current__r * ((var_Ca_independent_transient_outward_K_current__a * var_Ca_independent_transient_outward_K_current__s) + (var_Ca_independent_transient_outward_K_current__b * var_Ca_independent_transient_outward_K_current__s_slow)) * (var_Ca_independent_transient_outward_K_current__V - var_Ca_independent_transient_outward_K_current__E_K);
171 const NekDouble var_steady_state_outward_K_current__V = var_chaste_interface__membrane__V;
172 const NekDouble var_steady_state_outward_K_current__s_ss = var_chaste_interface__steady_state_outward_K_current_s_ss_gate__s_ss;
173 const NekDouble var_steady_state_outward_K_current__E_K = var_Ca_independent_transient_outward_K_current__E_K;
174 const NekDouble var_steady_state_outward_K_current__g_ss_normal = 0.0077;
175 const NekDouble var_steady_state_outward_K_current__Diabetes = var_membrane__Diabetes;
176 const NekDouble var_steady_state_outward_K_current__g_ss = (var_steady_state_outward_K_current__Diabetes == 0.0) ? var_steady_state_outward_K_current__g_ss_normal : (0.77 * var_steady_state_outward_K_current__g_ss_normal);
177 const NekDouble var_steady_state_outward_K_current__r_ss = var_chaste_interface__steady_state_outward_K_current_r_ss_gate__r_ss;
178 const NekDouble var_steady_state_outward_K_current__i_ss = var_steady_state_outward_K_current__g_ss * var_steady_state_outward_K_current__r_ss * var_steady_state_outward_K_current__s_ss * (var_steady_state_outward_K_current__V - var_steady_state_outward_K_current__E_K);
179 const NekDouble var_hyperpolarisation_activated_current__f_Na = 0.2;
180 const NekDouble var_hyperpolarisation_activated_current__E_Na = var_sodium_current__E_Na;
181 const NekDouble var_hyperpolarisation_activated_current__y = var_chaste_interface__hyperpolarisation_activated_current_y_gate__y;
182 const NekDouble var_hyperpolarisation_activated_current__V = var_chaste_interface__membrane__V;
183 const NekDouble var_hyperpolarisation_activated_current__g_f = 0.00145;
184 const NekDouble var_hyperpolarisation_activated_current__i_f_Na = var_hyperpolarisation_activated_current__g_f * var_hyperpolarisation_activated_current__y * var_hyperpolarisation_activated_current__f_Na * (var_hyperpolarisation_activated_current__V - var_hyperpolarisation_activated_current__E_Na);
185 const NekDouble var_hyperpolarisation_activated_current__f_K = 1.0 - var_hyperpolarisation_activated_current__f_Na;
186 const NekDouble var_hyperpolarisation_activated_current__E_K = var_Ca_independent_transient_outward_K_current__E_K;
187 const NekDouble var_hyperpolarisation_activated_current__i_f_K = var_hyperpolarisation_activated_current__g_f * var_hyperpolarisation_activated_current__y * var_hyperpolarisation_activated_current__f_K * (var_hyperpolarisation_activated_current__V - var_hyperpolarisation_activated_current__E_K);
188 const NekDouble var_inward_rectifier__T = var_membrane__T;
189 const NekDouble var_inward_rectifier__V = var_chaste_interface__membrane__V;
190 const NekDouble var_inward_rectifier__g_K1 = 0.024;
191 const NekDouble var_inward_rectifier__R = var_membrane__R;
192 const NekDouble var_inward_rectifier__K_o = var_standard_ionic_concentrations__K_o;
193 const NekDouble var_inward_rectifier__F = var_membrane__F;
194 const NekDouble var_inward_rectifier__E_K = var_Ca_independent_transient_outward_K_current__E_K;
195 const NekDouble var_inward_rectifier__i_K1 = ((((48.0 / (exp((var_inward_rectifier__V + 37.0) / 25.0) + exp((var_inward_rectifier__V + 37.0) / (-25.0)))) + 10.0) * 0.0001) / (1.0 + exp((var_inward_rectifier__V - (var_inward_rectifier__E_K + 76.77)) / (-17.0)))) + ((var_inward_rectifier__g_K1 * (var_inward_rectifier__V - (var_inward_rectifier__E_K + 1.73))) / ((1.0 + exp((1.613 * var_inward_rectifier__F * (var_inward_rectifier__V - (var_inward_rectifier__E_K + 1.73))) / (var_inward_rectifier__R * var_inward_rectifier__T))) * (1.0 + exp((var_inward_rectifier__K_o - 0.9988) / (-0.124)))));
196 const NekDouble var_background_currents__R = var_membrane__R;
197 const NekDouble var_background_currents__Ca_i = var_chaste_interface__intracellular_ion_concentrations__Ca_i;
198 const NekDouble var_standard_ionic_concentrations__Ca_o = 1.2;
199 const NekDouble var_background_currents__Ca_o = var_standard_ionic_concentrations__Ca_o;
200 const NekDouble var_background_currents__F = var_membrane__F;
201 const NekDouble var_background_currents__T = var_membrane__T;
202 const NekDouble var_background_currents__E_Ca = ((0.5 * var_background_currents__R * var_background_currents__T) / var_background_currents__F) * log(var_background_currents__Ca_o / var_background_currents__Ca_i);
203 const NekDouble var_background_currents__V = var_chaste_interface__membrane__V;
204 const NekDouble var_background_currents__g_B_Ca_normal = 3.24e-05;
205 const NekDouble var_background_currents__Diabetes = var_membrane__Diabetes;
206 const NekDouble var_background_currents__g_B_Ca = (var_background_currents__Diabetes == 0.0) ? var_background_currents__g_B_Ca_normal : (0.5 * var_background_currents__g_B_Ca_normal);
207 const NekDouble var_background_currents__i_B_Ca = var_background_currents__g_B_Ca * (var_background_currents__V - var_background_currents__E_Ca);
208 const NekDouble var_background_currents__g_B_K = 0.000138;
209 const NekDouble var_background_currents__E_K = var_Ca_independent_transient_outward_K_current__E_K;
210 const NekDouble var_background_currents__i_B_K = var_background_currents__g_B_K * (var_background_currents__V - var_background_currents__E_K);
211 const NekDouble var_background_currents__E_Na = var_sodium_current__E_Na;
212 const NekDouble var_background_currents__g_B_Na_normal = 8.015e-05;
213 const NekDouble var_background_currents__g_B_Na = (var_background_currents__Diabetes == 0.0) ? var_background_currents__g_B_Na_normal : (1.25 * var_background_currents__g_B_Na_normal);
214 const NekDouble var_background_currents__i_B_Na = var_background_currents__g_B_Na * (var_background_currents__V - var_background_currents__E_Na);
215 const NekDouble var_sodium_potassium_pump__Na_o = var_standard_ionic_concentrations__Na_o;
216 const NekDouble var_sodium_potassium_pump__sigma = (exp(var_sodium_potassium_pump__Na_o / 67.3) - 1.0) / 7.0;
217 const NekDouble var_sodium_potassium_pump__K_m_Na = 10.0;
218 const NekDouble var_sodium_potassium_pump__K_m_K = 1.5;
219 const NekDouble var_sodium_potassium_pump__V = var_chaste_interface__membrane__V;
220 const NekDouble var_sodium_potassium_pump__T = var_membrane__T;
221 const NekDouble var_sodium_potassium_pump__R = var_membrane__R;
222 const NekDouble var_sodium_potassium_pump__Na_i = var_chaste_interface__intracellular_ion_concentrations__Na_i;
223 const NekDouble var_sodium_potassium_pump__F = var_membrane__F;
224 const NekDouble var_sodium_potassium_pump__i_NaK_max_normal = 0.08;
225 const NekDouble var_sodium_potassium_pump__Diabetes = var_membrane__Diabetes;
226 const NekDouble var_sodium_potassium_pump__i_NaK_max = (var_sodium_potassium_pump__Diabetes == 0.0) ? var_sodium_potassium_pump__i_NaK_max_normal : (0.63 * var_sodium_potassium_pump__i_NaK_max_normal);
227 const NekDouble var_sodium_potassium_pump__K_o = var_standard_ionic_concentrations__K_o;
228 const NekDouble var_sodium_potassium_pump__i_NaK = (((var_sodium_potassium_pump__i_NaK_max / (1.0 + (0.1245 * exp(((-0.1) * var_sodium_potassium_pump__V * var_sodium_potassium_pump__F) / (var_sodium_potassium_pump__R * var_sodium_potassium_pump__T))) + (0.0365 * var_sodium_potassium_pump__sigma * exp(((-var_sodium_potassium_pump__V) * var_sodium_potassium_pump__F) / (var_sodium_potassium_pump__R * var_sodium_potassium_pump__T))))) * var_sodium_potassium_pump__K_o) / (var_sodium_potassium_pump__K_o + var_sodium_potassium_pump__K_m_K)) / (1.0 + pow(var_sodium_potassium_pump__K_m_Na / var_sodium_potassium_pump__Na_i, 1.5));
229 const NekDouble var_Na_Ca_ion_exchanger_current__Na_o = var_standard_ionic_concentrations__Na_o;
230 const NekDouble var_Na_Ca_ion_exchanger_current__Na_i = var_chaste_interface__intracellular_ion_concentrations__Na_i;
231 const NekDouble var_Na_Ca_ion_exchanger_current__gamma_NaCa = 0.5;
232 const NekDouble var_Na_Ca_ion_exchanger_current__Ca_i = var_chaste_interface__intracellular_ion_concentrations__Ca_i;
233 const NekDouble var_Na_Ca_ion_exchanger_current__Ca_o = var_standard_ionic_concentrations__Ca_o;
234 const NekDouble var_Na_Ca_ion_exchanger_current__K_NaCa = 9.984e-06;
235 const NekDouble var_Na_Ca_ion_exchanger_current__d_NaCa = 0.0001;
236 const NekDouble var_Na_Ca_ion_exchanger_current__V = var_chaste_interface__membrane__V;
237 const NekDouble var_Na_Ca_ion_exchanger_current__i_NaCa = (var_Na_Ca_ion_exchanger_current__K_NaCa * ((pow(var_Na_Ca_ion_exchanger_current__Na_i, 3.0) * var_Na_Ca_ion_exchanger_current__Ca_o * exp(0.03743 * var_Na_Ca_ion_exchanger_current__V * var_Na_Ca_ion_exchanger_current__gamma_NaCa)) - (pow(var_Na_Ca_ion_exchanger_current__Na_o, 3.0) * var_Na_Ca_ion_exchanger_current__Ca_i * exp(0.03743 * var_Na_Ca_ion_exchanger_current__V * (var_Na_Ca_ion_exchanger_current__gamma_NaCa - 1.0))))) / (1.0 + (var_Na_Ca_ion_exchanger_current__d_NaCa * ((var_Na_Ca_ion_exchanger_current__Ca_i * pow(var_Na_Ca_ion_exchanger_current__Na_o, 3.0)) + (var_Na_Ca_ion_exchanger_current__Ca_o * pow(var_Na_Ca_ion_exchanger_current__Na_i, 3.0)))));
238 const NekDouble var_sarcolemmal_calcium_pump_current__Ca_i = var_chaste_interface__intracellular_ion_concentrations__Ca_i;
239 const NekDouble var_sarcolemmal_calcium_pump_current__i_Ca_P_max = 0.004;
240 const NekDouble var_sarcolemmal_calcium_pump_current__i_Ca_P = (var_sarcolemmal_calcium_pump_current__i_Ca_P_max * var_sarcolemmal_calcium_pump_current__Ca_i) / (var_sarcolemmal_calcium_pump_current__Ca_i + 0.0004);
241 const NekDouble var_sodium_current_m_gate__V = var_sodium_current__V;
242 const NekDouble var_sodium_current_m_gate__m_infinity = 1.0 / (1.0 + exp((var_sodium_current_m_gate__V + 45.0) / (-6.5)));
243 const NekDouble var_sodium_current_m_gate__tau_m = 0.00136 / (((0.32 * (var_sodium_current_m_gate__V + 47.13)) / (1.0 - exp((-0.1) * (var_sodium_current_m_gate__V + 47.13)))) + (0.08 * exp((-var_sodium_current_m_gate__V) / 11.0)));
244 const NekDouble var_sodium_current_m_gate__m = var_sodium_current__m;
245 const NekDouble var_sodium_current_m_gate__d_m_d_environment__time = (var_sodium_current_m_gate__m_infinity - var_sodium_current_m_gate__m) / var_sodium_current_m_gate__tau_m;
246 const NekDouble var_sodium_current__sodium_current_m_gate__d_m_d_environment__time = var_sodium_current_m_gate__d_m_d_environment__time;
247 const NekDouble var_sodium_current_h_gate__V = var_sodium_current__V;
248 const NekDouble var_sodium_current_h_gate__h_infinity = 1.0 / (1.0 + exp((var_sodium_current_h_gate__V + 76.1) / 6.07));
249 const NekDouble var_sodium_current_h_gate__h = var_sodium_current__h;
250 const NekDouble var_sodium_current_h_gate__tau_h = (var_sodium_current_h_gate__V >= (-40.0)) ? (0.0004537 * (1.0 + exp((-(var_sodium_current_h_gate__V + 10.66)) / 11.1))) : (0.00349 / ((0.135 * exp((-(var_sodium_current_h_gate__V + 80.0)) / 6.8)) + (3.56 * exp(0.079 * var_sodium_current_h_gate__V)) + (310000.0 * exp(0.35 * var_sodium_current_h_gate__V))));
251 const NekDouble var_sodium_current_h_gate__d_h_d_environment__time = (var_sodium_current_h_gate__h_infinity - var_sodium_current_h_gate__h) / var_sodium_current_h_gate__tau_h;
252 const NekDouble var_sodium_current__sodium_current_h_gate__d_h_d_environment__time = var_sodium_current_h_gate__d_h_d_environment__time;
253 const NekDouble var_sodium_current_j_gate__V = var_sodium_current__V;
254 const NekDouble var_sodium_current_j_gate__j_infinity = 1.0 / (1.0 + exp((var_sodium_current_j_gate__V + 76.1) / 6.07));
255 const NekDouble var_sodium_current_j_gate__j = var_sodium_current__j;
256 const NekDouble var_sodium_current_j_gate__tau_j = (var_sodium_current_j_gate__V >= (-40.0)) ? ((0.01163 * (1.0 + exp((-0.1) * (var_sodium_current_j_gate__V + 32.0)))) / exp((-2.535e-07) * var_sodium_current_j_gate__V)) : (0.00349 / ((((var_sodium_current_j_gate__V + 37.78) / (1.0 + exp(0.311 * (var_sodium_current_j_gate__V + 79.23)))) * (((-127140.0) * exp(0.2444 * var_sodium_current_j_gate__V)) - (3.474e-05 * exp((-0.04391) * var_sodium_current_j_gate__V)))) + ((0.1212 * exp((-0.01052) * var_sodium_current_j_gate__V)) / (1.0 + exp((-0.1378) * (var_sodium_current_j_gate__V + 40.14))))));
257 const NekDouble var_sodium_current_j_gate__d_j_d_environment__time = (var_sodium_current_j_gate__j_infinity - var_sodium_current_j_gate__j) / var_sodium_current_j_gate__tau_j;
258 const NekDouble var_sodium_current__sodium_current_j_gate__d_j_d_environment__time = var_sodium_current_j_gate__d_j_d_environment__time;
259 const NekDouble var_L_type_Ca_channel__Ca_ss = var_chaste_interface__intracellular_ion_concentrations__Ca_ss;
260 const NekDouble var_L_type_Ca_channel_d_gate__d = var_L_type_Ca_channel__d;
261 const NekDouble var_L_type_Ca_channel_d_gate__V = var_L_type_Ca_channel__V;
262 const NekDouble var_L_type_Ca_channel_d_gate__d_infinity = 1.0 / (1.0 + exp((var_L_type_Ca_channel_d_gate__V + 15.3) / (-5.0)));
263 const NekDouble var_L_type_Ca_channel_d_gate__tau_d = (0.00305 * exp((-0.0045) * pow(var_L_type_Ca_channel_d_gate__V + 7.0, 2.0))) + (0.00105 * exp((-0.002) * pow(var_L_type_Ca_channel_d_gate__V - 18.0, 2.0))) + 0.00025;
264 const NekDouble var_L_type_Ca_channel_d_gate__d_d_d_environment__time = (var_L_type_Ca_channel_d_gate__d_infinity - var_L_type_Ca_channel_d_gate__d) / var_L_type_Ca_channel_d_gate__tau_d;
265 const NekDouble var_L_type_Ca_channel__L_type_Ca_channel_d_gate__d_d_d_environment__time = var_L_type_Ca_channel_d_gate__d_d_d_environment__time;
266 const NekDouble var_L_type_Ca_channel_f_11_gate__f_11 = var_L_type_Ca_channel__f_11;
267 const NekDouble var_L_type_Ca_channel_f_11_gate__Diabetes = var_L_type_Ca_channel__Diabetes;
268 const NekDouble var_L_type_Ca_channel_f_11_gate__V = var_L_type_Ca_channel__V;
269 const NekDouble var_L_type_Ca_channel_f_11_gate__tau_f_11 = (var_L_type_Ca_channel_f_11_gate__Diabetes == 0.0) ? ((0.105 * exp(-pow((var_L_type_Ca_channel_f_11_gate__V + 45.0) / 12.0, 2.0))) + (0.04 / (1.0 + exp(((-var_L_type_Ca_channel_f_11_gate__V) + 25.0) / 25.0))) + (0.015 / (1.0 + exp((var_L_type_Ca_channel_f_11_gate__V + 75.0) / 25.0))) + 0.0017) : ((0.9 * 0.105 * exp(-pow((var_L_type_Ca_channel_f_11_gate__V + 45.0) / 12.0, 2.0))) + (0.04 / (1.0 + exp(((-var_L_type_Ca_channel_f_11_gate__V) + 25.0) / 25.0))) + (0.015 / (1.0 + exp((var_L_type_Ca_channel_f_11_gate__V + 75.0) / 25.0))) + 0.0017);
270 const NekDouble var_L_type_Ca_channel_f_11_gate__f_11_infinity = 1.0 / (1.0 + exp((var_L_type_Ca_channel_f_11_gate__V + 26.7) / 5.4));
271 const NekDouble var_L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time = (var_L_type_Ca_channel_f_11_gate__f_11_infinity - var_L_type_Ca_channel_f_11_gate__f_11) / var_L_type_Ca_channel_f_11_gate__tau_f_11;
272 const NekDouble var_L_type_Ca_channel__L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time = var_L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time;
273 const NekDouble var_L_type_Ca_channel_f_12_gate__f_12 = var_L_type_Ca_channel__f_12;
274 const NekDouble var_L_type_Ca_channel_f_12_gate__V = var_L_type_Ca_channel__V;
275 const NekDouble var_L_type_Ca_channel_f_12_gate__f_12_infinity = 1.0 / (1.0 + exp((var_L_type_Ca_channel_f_12_gate__V + 26.7) / 5.4));
276 const NekDouble var_L_type_Ca_channel_f_12_gate__tau_f_12 = (0.041 * exp(-pow((var_L_type_Ca_channel_f_12_gate__V + 47.0) / 12.0, 2.0))) + (0.08 / (1.0 + exp((var_L_type_Ca_channel_f_12_gate__V + 55.0) / (-5.0)))) + (0.015 / (1.0 + exp((var_L_type_Ca_channel_f_12_gate__V + 75.0) / 25.0))) + 0.0017;
277 const NekDouble var_L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time = (var_L_type_Ca_channel_f_12_gate__f_12_infinity - var_L_type_Ca_channel_f_12_gate__f_12) / var_L_type_Ca_channel_f_12_gate__tau_f_12;
278 const NekDouble var_L_type_Ca_channel__L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time = var_L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time;
279 const NekDouble var_L_type_Ca_channel_Ca_inact_gate__Ca_ss = var_L_type_Ca_channel__Ca_ss;
280 const NekDouble var_L_type_Ca_channel_Ca_inact_gate__Ca_inact_infinity = 1.0 / (1.0 + (var_L_type_Ca_channel_Ca_inact_gate__Ca_ss / 0.01));
281 const NekDouble var_L_type_Ca_channel_Ca_inact_gate__tau_Ca_inact = 0.009;
282 const NekDouble var_L_type_Ca_channel_Ca_inact_gate__Ca_inact = var_L_type_Ca_channel__Ca_inact;
283 const NekDouble var_L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time = (var_L_type_Ca_channel_Ca_inact_gate__Ca_inact_infinity - var_L_type_Ca_channel_Ca_inact_gate__Ca_inact) / var_L_type_Ca_channel_Ca_inact_gate__tau_Ca_inact;
284 const NekDouble var_L_type_Ca_channel__L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time = var_L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time;
285 const NekDouble var_Ca_independent_transient_outward_K_current_r_gate__V = var_Ca_independent_transient_outward_K_current__V;
286 const NekDouble var_Ca_independent_transient_outward_K_current_r_gate__r_infinity = 1.0 / (1.0 + exp((var_Ca_independent_transient_outward_K_current_r_gate__V + 10.6) / (-11.42)));
287 const NekDouble var_Ca_independent_transient_outward_K_current_r_gate__r = var_Ca_independent_transient_outward_K_current__r;
288 const NekDouble var_Ca_independent_transient_outward_K_current_r_gate__tau_r = 1.0 / ((45.16 * exp(0.03577 * (var_Ca_independent_transient_outward_K_current_r_gate__V + 50.0))) + (98.9 * exp((-0.1) * (var_Ca_independent_transient_outward_K_current_r_gate__V + 38.0))));
289 const NekDouble var_Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time = (var_Ca_independent_transient_outward_K_current_r_gate__r_infinity - var_Ca_independent_transient_outward_K_current_r_gate__r) / var_Ca_independent_transient_outward_K_current_r_gate__tau_r;
290 const NekDouble var_Ca_independent_transient_outward_K_current__Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time = var_Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time;
291 const NekDouble var_Ca_independent_transient_outward_K_current_s_gate__V = var_Ca_independent_transient_outward_K_current__V;
292 const NekDouble var_Ca_independent_transient_outward_K_current_s_gate__s_infinity = 1.0 / (1.0 + exp((var_Ca_independent_transient_outward_K_current_s_gate__V + 45.3) / 6.8841));
293 const NekDouble var_Ca_independent_transient_outward_K_current_s_gate__tau_s = (0.35 * exp(-pow((var_Ca_independent_transient_outward_K_current_s_gate__V + 70.0) / 15.0, 2.0))) + 0.035;
294 const NekDouble var_Ca_independent_transient_outward_K_current_s_gate__s = var_Ca_independent_transient_outward_K_current__s;
295 const NekDouble var_Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time = (var_Ca_independent_transient_outward_K_current_s_gate__s_infinity - var_Ca_independent_transient_outward_K_current_s_gate__s) / var_Ca_independent_transient_outward_K_current_s_gate__tau_s;
296 const NekDouble var_Ca_independent_transient_outward_K_current__Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time = var_Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time;
297 const NekDouble var_Ca_independent_transient_outward_K_current_s_slow_gate__V = var_Ca_independent_transient_outward_K_current__V;
298 const NekDouble var_Ca_independent_transient_outward_K_current_s_slow_gate__s_slow_infinity = 1.0 / (1.0 + exp((var_Ca_independent_transient_outward_K_current_s_slow_gate__V + 45.3) / 6.8841));
299 const NekDouble var_Ca_independent_transient_outward_K_current_s_slow_gate__s_slow = var_Ca_independent_transient_outward_K_current__s_slow;
300 const NekDouble var_Ca_independent_transient_outward_K_current_s_slow_gate__tau_s_slow = (3.7 * exp(-pow((var_Ca_independent_transient_outward_K_current_s_slow_gate__V + 70.0) / 30.0, 2.0))) + 0.035;
301 const NekDouble var_Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time = (var_Ca_independent_transient_outward_K_current_s_slow_gate__s_slow_infinity - var_Ca_independent_transient_outward_K_current_s_slow_gate__s_slow) / var_Ca_independent_transient_outward_K_current_s_slow_gate__tau_s_slow;
302 const NekDouble var_Ca_independent_transient_outward_K_current__Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time = var_Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time;
303 const NekDouble var_steady_state_outward_K_current_r_ss_gate__V = var_steady_state_outward_K_current__V;
304 const NekDouble var_steady_state_outward_K_current_r_ss_gate__r_ss_infinity = 1.0 / (1.0 + exp((var_steady_state_outward_K_current_r_ss_gate__V + 11.5) / (-11.82)));
305 const NekDouble var_steady_state_outward_K_current_r_ss_gate__r_ss = var_steady_state_outward_K_current__r_ss;
306 const NekDouble var_steady_state_outward_K_current_r_ss_gate__tau_r_ss = 10.0 / ((45.16 * exp(0.03577 * (var_steady_state_outward_K_current_r_ss_gate__V + 50.0))) + (98.9 * exp((-0.1) * (var_steady_state_outward_K_current_r_ss_gate__V + 38.0))));
307 const NekDouble var_steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time = (var_steady_state_outward_K_current_r_ss_gate__r_ss_infinity - var_steady_state_outward_K_current_r_ss_gate__r_ss) / var_steady_state_outward_K_current_r_ss_gate__tau_r_ss;
308 const NekDouble var_steady_state_outward_K_current__steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time = var_steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time;
309 const NekDouble var_steady_state_outward_K_current_s_ss_gate__s_ss = var_steady_state_outward_K_current__s_ss;
310 const NekDouble var_steady_state_outward_K_current_s_ss_gate__tau_s_ss = 2.1;
311 const NekDouble var_steady_state_outward_K_current_s_ss_gate__V = var_steady_state_outward_K_current__V;
312 const NekDouble var_steady_state_outward_K_current_s_ss_gate__s_ss_infinity = 1.0 / (1.0 + exp((var_steady_state_outward_K_current_s_ss_gate__V + 87.5) / 10.3));
313 const NekDouble var_steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time = (var_steady_state_outward_K_current_s_ss_gate__s_ss_infinity - var_steady_state_outward_K_current_s_ss_gate__s_ss) / var_steady_state_outward_K_current_s_ss_gate__tau_s_ss;
314 const NekDouble var_steady_state_outward_K_current__steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time = var_steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time;
315 const NekDouble var_hyperpolarisation_activated_current_y_gate__V = var_hyperpolarisation_activated_current__V;
316 const NekDouble var_hyperpolarisation_activated_current_y_gate__tau_y = 1.0 / ((0.11885 * exp((var_hyperpolarisation_activated_current_y_gate__V + 80.0) / 28.37)) + (0.5623 * exp((var_hyperpolarisation_activated_current_y_gate__V + 80.0) / (-14.19))));
317 const NekDouble var_hyperpolarisation_activated_current_y_gate__y = var_hyperpolarisation_activated_current__y;
318 const NekDouble var_hyperpolarisation_activated_current_y_gate__y_infinity = 1.0 / (1.0 + exp((var_hyperpolarisation_activated_current_y_gate__V + 138.6) / 10.48));
319 const NekDouble var_hyperpolarisation_activated_current_y_gate__d_y_d_environment__time = (var_hyperpolarisation_activated_current_y_gate__y_infinity - var_hyperpolarisation_activated_current_y_gate__y) / var_hyperpolarisation_activated_current_y_gate__tau_y;
320 const NekDouble var_hyperpolarisation_activated_current__hyperpolarisation_activated_current_y_gate__d_y_d_environment__time = var_hyperpolarisation_activated_current_y_gate__d_y_d_environment__time;
321 const NekDouble var_SR_Ca_release_channel__Ca_ss = var_chaste_interface__intracellular_ion_concentrations__Ca_ss;
322 const NekDouble var_SR_Ca_release_channel__v1 = 1800.0;
323 const NekDouble var_SR_Ca_release_channel__P_O2 = var_chaste_interface__SR_Ca_release_channel__P_O2;
324 const NekDouble var_SR_Ca_release_channel__P_O1 = var_chaste_interface__SR_Ca_release_channel__P_O1;
325 const NekDouble var_SR_Ca_release_channel__Ca_JSR = var_chaste_interface__intracellular_ion_concentrations__Ca_JSR;
326 const NekDouble var_SR_Ca_release_channel__J_rel = var_SR_Ca_release_channel__v1 * (var_SR_Ca_release_channel__P_O1 + var_SR_Ca_release_channel__P_O2) * (var_SR_Ca_release_channel__Ca_JSR - var_SR_Ca_release_channel__Ca_ss);
327 const NekDouble var_SR_Ca_release_channel__k_a_plus = 1.215e+13;
328 const NekDouble var_SR_Ca_release_channel__k_a_minus = 576.0;
329 const NekDouble var_SR_Ca_release_channel__k_b_plus = 4050000000.0;
330 const NekDouble var_SR_Ca_release_channel__k_b_minus = 1930.0;
331 const NekDouble var_SR_Ca_release_channel__k_c_plus = 100.0;
332 const NekDouble var_SR_Ca_release_channel__k_c_minus = 0.8;
333 const NekDouble var_SR_Ca_release_channel__P_C1 = var_chaste_interface__SR_Ca_release_channel__P_C1;
334 const NekDouble var_SR_Ca_release_channel__P_C2 = var_chaste_interface__SR_Ca_release_channel__P_C2;
335 const NekDouble var_SR_Ca_release_channel__n = 4.0;
336 const NekDouble var_SR_Ca_release_channel__m = 3.0;
337 const NekDouble var_SR_Ca_release_channel__d_P_O1_d_environment__time = ((var_SR_Ca_release_channel__k_a_plus * pow(var_SR_Ca_release_channel__Ca_ss / 1.0, var_SR_Ca_release_channel__n) * var_SR_Ca_release_channel__P_C1) - ((var_SR_Ca_release_channel__k_a_minus * var_SR_Ca_release_channel__P_O1) + (var_SR_Ca_release_channel__k_b_plus * pow(var_SR_Ca_release_channel__Ca_ss / 1.0, var_SR_Ca_release_channel__m) * var_SR_Ca_release_channel__P_O1) + (var_SR_Ca_release_channel__k_c_plus * var_SR_Ca_release_channel__P_O1))) + (var_SR_Ca_release_channel__k_b_minus * var_SR_Ca_release_channel__P_O2) + (var_SR_Ca_release_channel__k_c_minus * var_SR_Ca_release_channel__P_C2);
338 const NekDouble var_SR_Ca_release_channel__d_P_O2_d_environment__time = (var_SR_Ca_release_channel__k_b_plus * pow(var_SR_Ca_release_channel__Ca_ss / 1.0, var_SR_Ca_release_channel__m) * var_SR_Ca_release_channel__P_O1) - (var_SR_Ca_release_channel__k_b_minus * var_SR_Ca_release_channel__P_O2);
339 const NekDouble var_SR_Ca_release_channel__d_P_C1_d_environment__time = ((-var_SR_Ca_release_channel__k_a_plus) * pow(var_SR_Ca_release_channel__Ca_ss / 1.0, var_SR_Ca_release_channel__n) * var_SR_Ca_release_channel__P_C1) + (var_SR_Ca_release_channel__k_a_minus * var_SR_Ca_release_channel__P_O1);
340 const NekDouble var_SR_Ca_release_channel__d_P_C2_d_environment__time = (var_SR_Ca_release_channel__k_c_plus * var_SR_Ca_release_channel__P_O1) - (var_SR_Ca_release_channel__k_c_minus * var_SR_Ca_release_channel__P_C2);
341 const NekDouble var_SERCA2a_pump__N_fb = 1.2;
342 const NekDouble var_SERCA2a_pump__K_fb = 0.000168;
343 const NekDouble var_SERCA2a_pump__Ca_i = var_chaste_interface__intracellular_ion_concentrations__Ca_i;
344 const NekDouble var_SERCA2a_pump__fb = pow(var_SERCA2a_pump__Ca_i / var_SERCA2a_pump__K_fb, var_SERCA2a_pump__N_fb);
345 const NekDouble var_SERCA2a_pump__Diabetes = var_membrane__Diabetes;
346 const NekDouble var_SERCA2a_pump__Vmaxf_normal = 0.04;
347 const NekDouble var_SERCA2a_pump__Vmaxf = (var_SERCA2a_pump__Diabetes == 0.0) ? var_SERCA2a_pump__Vmaxf_normal : (0.8 * var_SERCA2a_pump__Vmaxf_normal);
348 const NekDouble var_SERCA2a_pump__K_SR_normal = 1.0;
349 const NekDouble var_SERCA2a_pump__K_SR = (var_SERCA2a_pump__Diabetes == 0.0) ? var_SERCA2a_pump__K_SR_normal : (0.55 * var_SERCA2a_pump__K_SR_normal);
350 const NekDouble var_SERCA2a_pump__Ca_NSR = var_chaste_interface__intracellular_ion_concentrations__Ca_NSR;
351 const NekDouble var_SERCA2a_pump__K_rb = 3.29;
352 const NekDouble var_SERCA2a_pump__N_rb = 1.0;
353 const NekDouble var_SERCA2a_pump__rb = pow(var_SERCA2a_pump__Ca_NSR / var_SERCA2a_pump__K_rb, var_SERCA2a_pump__N_rb);
354 const NekDouble var_SERCA2a_pump__Vmaxr = 0.9;
355 const NekDouble var_SERCA2a_pump__J_up = (var_SERCA2a_pump__K_SR * ((var_SERCA2a_pump__Vmaxf * var_SERCA2a_pump__fb) - (var_SERCA2a_pump__Vmaxr * var_SERCA2a_pump__rb))) / (1.0 + var_SERCA2a_pump__fb + var_SERCA2a_pump__rb);
356 const NekDouble var_intracellular_and_SR_Ca_fluxes__tau_tr = 0.0005747;
357 const NekDouble var_intracellular_and_SR_Ca_fluxes__Ca_JSR = var_chaste_interface__intracellular_ion_concentrations__Ca_JSR;
358 const NekDouble var_intracellular_and_SR_Ca_fluxes__Ca_NSR = var_chaste_interface__intracellular_ion_concentrations__Ca_NSR;
359 const NekDouble var_intracellular_and_SR_Ca_fluxes__J_tr = (var_intracellular_and_SR_Ca_fluxes__Ca_NSR - var_intracellular_and_SR_Ca_fluxes__Ca_JSR) / var_intracellular_and_SR_Ca_fluxes__tau_tr;
360 const NekDouble var_intracellular_and_SR_Ca_fluxes__Ca_ss = var_chaste_interface__intracellular_ion_concentrations__Ca_ss;
361 const NekDouble var_intracellular_and_SR_Ca_fluxes__Ca_i = var_chaste_interface__intracellular_ion_concentrations__Ca_i;
362 const NekDouble var_intracellular_and_SR_Ca_fluxes__tau_xfer = 0.0267;
363 const NekDouble var_intracellular_and_SR_Ca_fluxes__J_xfer = (var_intracellular_and_SR_Ca_fluxes__Ca_ss - var_intracellular_and_SR_Ca_fluxes__Ca_i) / var_intracellular_and_SR_Ca_fluxes__tau_xfer;
364 const NekDouble var_intracellular_and_SR_Ca_fluxes__k_htrpn_plus = 200000.0;
365 const NekDouble var_intracellular_and_SR_Ca_fluxes__k_htrpn_minus = 0.066;
366 const NekDouble var_intracellular_and_SR_Ca_fluxes__HTRPNCa = var_chaste_interface__intracellular_and_SR_Ca_fluxes__HTRPNCa;
367 const NekDouble var_intracellular_and_SR_Ca_fluxes__HTRPN_tot = 0.14;
368 const NekDouble var_intracellular_and_SR_Ca_fluxes__J_HTRPNCa = (var_intracellular_and_SR_Ca_fluxes__k_htrpn_plus * var_intracellular_and_SR_Ca_fluxes__Ca_i * (var_intracellular_and_SR_Ca_fluxes__HTRPN_tot - var_intracellular_and_SR_Ca_fluxes__HTRPNCa)) - (var_intracellular_and_SR_Ca_fluxes__k_htrpn_minus * var_intracellular_and_SR_Ca_fluxes__HTRPNCa);
369 const NekDouble var_intracellular_and_SR_Ca_fluxes__LTRPNCa = var_chaste_interface__intracellular_and_SR_Ca_fluxes__LTRPNCa;
370 const NekDouble var_intracellular_and_SR_Ca_fluxes__LTRPN_tot = 0.07;
371 const NekDouble var_intracellular_and_SR_Ca_fluxes__k_ltrpn_plus = 40000.0;
372 const NekDouble var_intracellular_and_SR_Ca_fluxes__k_ltrpn_minus = 40.0;
373 const NekDouble var_intracellular_and_SR_Ca_fluxes__J_LTRPNCa = (var_intracellular_and_SR_Ca_fluxes__k_ltrpn_plus * var_intracellular_and_SR_Ca_fluxes__Ca_i * (var_intracellular_and_SR_Ca_fluxes__LTRPN_tot - var_intracellular_and_SR_Ca_fluxes__LTRPNCa)) - (var_intracellular_and_SR_Ca_fluxes__k_ltrpn_minus * var_intracellular_and_SR_Ca_fluxes__LTRPNCa);
374 const NekDouble var_intracellular_and_SR_Ca_fluxes__J_trpn = var_intracellular_and_SR_Ca_fluxes__J_HTRPNCa + var_intracellular_and_SR_Ca_fluxes__J_LTRPNCa;
375 const NekDouble var_intracellular_and_SR_Ca_fluxes__d_HTRPNCa_d_environment__time = (var_intracellular_and_SR_Ca_fluxes__k_htrpn_plus * var_intracellular_and_SR_Ca_fluxes__Ca_i * (var_intracellular_and_SR_Ca_fluxes__HTRPN_tot - var_intracellular_and_SR_Ca_fluxes__HTRPNCa)) - (var_intracellular_and_SR_Ca_fluxes__k_htrpn_minus * var_intracellular_and_SR_Ca_fluxes__HTRPNCa);
376 const NekDouble var_intracellular_and_SR_Ca_fluxes__d_LTRPNCa_d_environment__time = (var_intracellular_and_SR_Ca_fluxes__k_ltrpn_plus * var_intracellular_and_SR_Ca_fluxes__Ca_i * (var_intracellular_and_SR_Ca_fluxes__LTRPN_tot - var_intracellular_and_SR_Ca_fluxes__LTRPNCa)) - (var_intracellular_and_SR_Ca_fluxes__k_ltrpn_minus * var_intracellular_and_SR_Ca_fluxes__LTRPNCa);
377 const NekDouble var_intracellular_ion_concentrations__Ca_i = var_chaste_interface__intracellular_ion_concentrations__Ca_i;
378 const NekDouble var_intracellular_ion_concentrations__Ca_ss = var_chaste_interface__intracellular_ion_concentrations__Ca_ss;
379 const NekDouble var_intracellular_ion_concentrations__Ca_JSR = var_chaste_interface__intracellular_ion_concentrations__Ca_JSR;
380 const NekDouble var_intracellular_ion_concentrations__V_myo = 9.36e-06;
381 const NekDouble var_intracellular_ion_concentrations__V_JSR = 5.6e-07;
382 const NekDouble var_intracellular_ion_concentrations__V_NSR = 5.04e-06;
383 const NekDouble var_intracellular_ion_concentrations__V_SS = 1.2e-08;
384 const NekDouble var_intracellular_ion_concentrations__K_mCMDN = 0.00238;
385 const NekDouble var_intracellular_ion_concentrations__K_mCSQN = 0.8;
386 const NekDouble var_intracellular_ion_concentrations__K_mEGTA = 0.00015;
387 const NekDouble var_intracellular_ion_concentrations__CMDN_tot = 0.05;
388 const NekDouble var_intracellular_ion_concentrations__CSQN_tot = 15.0;
389 const NekDouble var_intracellular_ion_concentrations__EGTA_tot = 10.0;
390 const NekDouble var_intracellular_ion_concentrations__beta_i = 1.0 / (1.0 + ((var_intracellular_ion_concentrations__CMDN_tot * var_intracellular_ion_concentrations__K_mCMDN) / pow(var_intracellular_ion_concentrations__K_mCMDN + var_intracellular_ion_concentrations__Ca_i, 2.0)) + ((var_intracellular_ion_concentrations__EGTA_tot * var_intracellular_ion_concentrations__K_mEGTA) / pow(var_intracellular_ion_concentrations__K_mEGTA + var_intracellular_ion_concentrations__Ca_i, 2.0)));
391 const NekDouble var_intracellular_ion_concentrations__beta_SS = 1.0 / (1.0 + ((var_intracellular_ion_concentrations__CMDN_tot * var_intracellular_ion_concentrations__K_mCMDN) / pow(var_intracellular_ion_concentrations__K_mCMDN + var_intracellular_ion_concentrations__Ca_ss, 2.0)));
392 const NekDouble var_intracellular_ion_concentrations__beta_JSR = 1.0 / (1.0 + ((var_intracellular_ion_concentrations__CSQN_tot * var_intracellular_ion_concentrations__K_mCSQN) / pow(var_intracellular_ion_concentrations__K_mCSQN + var_intracellular_ion_concentrations__Ca_JSR, 2.0)));
393 const NekDouble var_intracellular_ion_concentrations__F = var_membrane__F;
394 const NekDouble var_intracellular_ion_concentrations__i_Na = var_sodium_current__i_Na;
395 const NekDouble var_intracellular_ion_concentrations__i_Ca_L = var_L_type_Ca_channel__i_Ca_L;
396 const NekDouble var_intracellular_ion_concentrations__i_B_Na = var_background_currents__i_B_Na;
397 const NekDouble var_intracellular_ion_concentrations__i_NaCa = var_Na_Ca_ion_exchanger_current__i_NaCa;
398 const NekDouble var_intracellular_ion_concentrations__i_NaK = var_sodium_potassium_pump__i_NaK;
399 const NekDouble var_intracellular_ion_concentrations__i_f_Na = var_hyperpolarisation_activated_current__i_f_Na;
400 const NekDouble var_intracellular_ion_concentrations__i_f_K = var_hyperpolarisation_activated_current__i_f_K;
401 const NekDouble var_intracellular_ion_concentrations__i_B_K = var_background_currents__i_B_K;
402 const NekDouble var_intracellular_ion_concentrations__i_K1 = var_inward_rectifier__i_K1;
403 const NekDouble var_intracellular_ion_concentrations__i_t = var_Ca_independent_transient_outward_K_current__i_t;
404 const NekDouble var_intracellular_ion_concentrations__i_ss = var_steady_state_outward_K_current__i_ss;
405 const NekDouble var_intracellular_ion_concentrations__i_Ca_P = var_sarcolemmal_calcium_pump_current__i_Ca_P;
406 const NekDouble var_intracellular_ion_concentrations__i_B_Ca = var_background_currents__i_B_Ca;
407 const NekDouble var_intracellular_ion_concentrations__J_up = var_SERCA2a_pump__J_up;
408 const NekDouble var_intracellular_ion_concentrations__J_rel = var_SR_Ca_release_channel__J_rel;
409 const NekDouble var_intracellular_ion_concentrations__J_xfer = var_intracellular_and_SR_Ca_fluxes__J_xfer;
410 const NekDouble var_intracellular_ion_concentrations__J_trpn = var_intracellular_and_SR_Ca_fluxes__J_trpn;
411 const NekDouble var_intracellular_ion_concentrations__J_tr = var_intracellular_and_SR_Ca_fluxes__J_tr;
412 const NekDouble var_intracellular_ion_concentrations__d_Na_i_d_environment__time = (-(var_intracellular_ion_concentrations__i_Na + var_intracellular_ion_concentrations__i_B_Na + (var_intracellular_ion_concentrations__i_NaCa * 3.0) + (var_intracellular_ion_concentrations__i_NaK * 3.0) + var_intracellular_ion_concentrations__i_f_Na)) / (var_intracellular_ion_concentrations__V_myo * var_intracellular_ion_concentrations__F);
413 const NekDouble var_intracellular_ion_concentrations__d_Ca_i_d_environment__time = var_intracellular_ion_concentrations__beta_i * (var_intracellular_ion_concentrations__J_xfer - (var_intracellular_ion_concentrations__J_up + var_intracellular_ion_concentrations__J_trpn + (((var_intracellular_ion_concentrations__i_B_Ca - (2.0 * var_intracellular_ion_concentrations__i_NaCa)) + var_intracellular_ion_concentrations__i_Ca_P) / (2.0 * var_intracellular_ion_concentrations__V_myo * var_intracellular_ion_concentrations__F))));
414 const NekDouble var_intracellular_ion_concentrations__d_K_i_d_environment__time = (-(var_intracellular_ion_concentrations__i_ss + var_intracellular_ion_concentrations__i_B_K + var_intracellular_ion_concentrations__i_t + var_intracellular_ion_concentrations__i_K1 + var_intracellular_ion_concentrations__i_f_K + (var_intracellular_ion_concentrations__i_NaK * (-2.0)))) / (var_intracellular_ion_concentrations__V_myo * var_intracellular_ion_concentrations__F);
415 const NekDouble var_intracellular_ion_concentrations__d_Ca_ss_d_environment__time = var_intracellular_ion_concentrations__beta_SS * ((((var_intracellular_ion_concentrations__J_rel * var_intracellular_ion_concentrations__V_JSR) / var_intracellular_ion_concentrations__V_SS) - ((var_intracellular_ion_concentrations__J_xfer * var_intracellular_ion_concentrations__V_myo) / var_intracellular_ion_concentrations__V_SS)) - (var_intracellular_ion_concentrations__i_Ca_L / (2.0 * var_intracellular_ion_concentrations__V_SS * var_intracellular_ion_concentrations__F)));
416 const NekDouble var_intracellular_ion_concentrations__d_Ca_JSR_d_environment__time = var_intracellular_ion_concentrations__beta_JSR * (var_intracellular_ion_concentrations__J_tr - var_intracellular_ion_concentrations__J_rel);
417 const NekDouble var_intracellular_ion_concentrations__d_Ca_NSR_d_environment__time = ((var_intracellular_ion_concentrations__J_up * var_intracellular_ion_concentrations__V_myo) / var_intracellular_ion_concentrations__V_NSR) - ((var_intracellular_ion_concentrations__J_tr * var_intracellular_ion_concentrations__V_JSR) / var_intracellular_ion_concentrations__V_NSR);
418 const NekDouble var_chaste_interface__sodium_current_m_gate__d_m_d_environment__time_converter = var_sodium_current__sodium_current_m_gate__d_m_d_environment__time;
419 const NekDouble var_chaste_interface__sodium_current_m_gate__d_m_d_environment__time = 0.001 * var_chaste_interface__sodium_current_m_gate__d_m_d_environment__time_converter;
420 const NekDouble var_chaste_interface__sodium_current_h_gate__d_h_d_environment__time_converter = var_sodium_current__sodium_current_h_gate__d_h_d_environment__time;
421 const NekDouble var_chaste_interface__sodium_current_h_gate__d_h_d_environment__time = 0.001 * var_chaste_interface__sodium_current_h_gate__d_h_d_environment__time_converter;
422 const NekDouble var_chaste_interface__sodium_current_j_gate__d_j_d_environment__time_converter = var_sodium_current__sodium_current_j_gate__d_j_d_environment__time;
423 const NekDouble var_chaste_interface__sodium_current_j_gate__d_j_d_environment__time = 0.001 * var_chaste_interface__sodium_current_j_gate__d_j_d_environment__time_converter;
424 const NekDouble var_chaste_interface__L_type_Ca_channel_d_gate__d_d_d_environment__time_converter = var_L_type_Ca_channel__L_type_Ca_channel_d_gate__d_d_d_environment__time;
425 const NekDouble var_chaste_interface__L_type_Ca_channel_d_gate__d_d_d_environment__time = 0.001 * var_chaste_interface__L_type_Ca_channel_d_gate__d_d_d_environment__time_converter;
426 const NekDouble var_chaste_interface__L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time_converter = var_L_type_Ca_channel__L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time;
427 const NekDouble var_chaste_interface__L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time = 0.001 * var_chaste_interface__L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time_converter;
428 const NekDouble var_chaste_interface__L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time_converter = var_L_type_Ca_channel__L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time;
429 const NekDouble var_chaste_interface__L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time = 0.001 * var_chaste_interface__L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time_converter;
430 const NekDouble var_chaste_interface__L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time_converter = var_L_type_Ca_channel__L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time;
431 const NekDouble var_chaste_interface__L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time = 0.001 * var_chaste_interface__L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time_converter;
432 const NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time_converter = var_Ca_independent_transient_outward_K_current__Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time;
433 const NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time = 0.001 * var_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time_converter;
434 const NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time_converter = var_Ca_independent_transient_outward_K_current__Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time;
435 const NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time = 0.001 * var_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time_converter;
436 const NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time_converter = var_Ca_independent_transient_outward_K_current__Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time;
437 const NekDouble var_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time = 0.001 * var_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time_converter;
438 const NekDouble var_chaste_interface__steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time_converter = var_steady_state_outward_K_current__steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time;
439 const NekDouble var_chaste_interface__steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time = 0.001 * var_chaste_interface__steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time_converter;
440 const NekDouble var_chaste_interface__steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time_converter = var_steady_state_outward_K_current__steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time;
441 const NekDouble var_chaste_interface__steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time = 0.001 * var_chaste_interface__steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time_converter;
442 const NekDouble var_chaste_interface__hyperpolarisation_activated_current_y_gate__d_y_d_environment__time_converter = var_hyperpolarisation_activated_current__hyperpolarisation_activated_current_y_gate__d_y_d_environment__time;
443 const NekDouble var_chaste_interface__hyperpolarisation_activated_current_y_gate__d_y_d_environment__time = 0.001 * var_chaste_interface__hyperpolarisation_activated_current_y_gate__d_y_d_environment__time_converter;
444 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_O1_d_environment__time_converter = var_SR_Ca_release_channel__d_P_O1_d_environment__time;
445 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_O1_d_environment__time = 0.001 * var_chaste_interface__SR_Ca_release_channel__d_P_O1_d_environment__time_converter;
446 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_O2_d_environment__time_converter = var_SR_Ca_release_channel__d_P_O2_d_environment__time;
447 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_O2_d_environment__time = 0.001 * var_chaste_interface__SR_Ca_release_channel__d_P_O2_d_environment__time_converter;
448 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_C1_d_environment__time_converter = var_SR_Ca_release_channel__d_P_C1_d_environment__time;
449 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_C1_d_environment__time = 0.001 * var_chaste_interface__SR_Ca_release_channel__d_P_C1_d_environment__time_converter;
450 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_C2_d_environment__time_converter = var_SR_Ca_release_channel__d_P_C2_d_environment__time;
451 const NekDouble var_chaste_interface__SR_Ca_release_channel__d_P_C2_d_environment__time = 0.001 * var_chaste_interface__SR_Ca_release_channel__d_P_C2_d_environment__time_converter;
452 const NekDouble var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_HTRPNCa_d_environment__time_converter = var_intracellular_and_SR_Ca_fluxes__d_HTRPNCa_d_environment__time;
453 const NekDouble var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_HTRPNCa_d_environment__time = 0.001 * var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_HTRPNCa_d_environment__time_converter;
454 const NekDouble var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_LTRPNCa_d_environment__time_converter = var_intracellular_and_SR_Ca_fluxes__d_LTRPNCa_d_environment__time;
455 const NekDouble var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_LTRPNCa_d_environment__time = 0.001 * var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_LTRPNCa_d_environment__time_converter;
456 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Na_i_d_environment__time_converter = var_intracellular_ion_concentrations__d_Na_i_d_environment__time;
457 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Na_i_d_environment__time = 0.001 * var_chaste_interface__intracellular_ion_concentrations__d_Na_i_d_environment__time_converter;
458 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_i_d_environment__time_converter = var_intracellular_ion_concentrations__d_Ca_i_d_environment__time;
459 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_i_d_environment__time = 0.001 * var_chaste_interface__intracellular_ion_concentrations__d_Ca_i_d_environment__time_converter;
460 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_K_i_d_environment__time_converter = var_intracellular_ion_concentrations__d_K_i_d_environment__time;
461 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_K_i_d_environment__time = 0.001 * var_chaste_interface__intracellular_ion_concentrations__d_K_i_d_environment__time_converter;
462 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_ss_d_environment__time_converter = var_intracellular_ion_concentrations__d_Ca_ss_d_environment__time;
463 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_ss_d_environment__time = 0.001 * var_chaste_interface__intracellular_ion_concentrations__d_Ca_ss_d_environment__time_converter;
464 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_JSR_d_environment__time_converter = var_intracellular_ion_concentrations__d_Ca_JSR_d_environment__time;
465 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_JSR_d_environment__time = 0.001 * var_chaste_interface__intracellular_ion_concentrations__d_Ca_JSR_d_environment__time_converter;
466 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_NSR_d_environment__time_converter = var_intracellular_ion_concentrations__d_Ca_NSR_d_environment__time;
467 const NekDouble var_chaste_interface__intracellular_ion_concentrations__d_Ca_NSR_d_environment__time = 0.001 * var_chaste_interface__intracellular_ion_concentrations__d_Ca_NSR_d_environment__time_converter;
468 const NekDouble d_dt_chaste_interface__sodium_current_m_gate__m = var_chaste_interface__sodium_current_m_gate__d_m_d_environment__time;
469 const NekDouble d_dt_chaste_interface__sodium_current_h_gate__h = var_chaste_interface__sodium_current_h_gate__d_h_d_environment__time;
470 const NekDouble d_dt_chaste_interface__sodium_current_j_gate__j = var_chaste_interface__sodium_current_j_gate__d_j_d_environment__time;
471 const NekDouble d_dt_chaste_interface__L_type_Ca_channel_d_gate__d = var_chaste_interface__L_type_Ca_channel_d_gate__d_d_d_environment__time;
472 const NekDouble d_dt_chaste_interface__L_type_Ca_channel_f_11_gate__f_11 = var_chaste_interface__L_type_Ca_channel_f_11_gate__d_f_11_d_environment__time;
473 const NekDouble d_dt_chaste_interface__L_type_Ca_channel_f_12_gate__f_12 = var_chaste_interface__L_type_Ca_channel_f_12_gate__d_f_12_d_environment__time;
474 const NekDouble d_dt_chaste_interface__L_type_Ca_channel_Ca_inact_gate__Ca_inact = var_chaste_interface__L_type_Ca_channel_Ca_inact_gate__d_Ca_inact_d_environment__time;
475 const NekDouble d_dt_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__r = var_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__d_r_d_environment__time;
476 const NekDouble d_dt_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__s = var_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__d_s_d_environment__time;
477 const NekDouble d_dt_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__s_slow = var_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__d_s_slow_d_environment__time;
478 const NekDouble d_dt_chaste_interface__steady_state_outward_K_current_r_ss_gate__r_ss = var_chaste_interface__steady_state_outward_K_current_r_ss_gate__d_r_ss_d_environment__time;
479 const NekDouble d_dt_chaste_interface__steady_state_outward_K_current_s_ss_gate__s_ss = var_chaste_interface__steady_state_outward_K_current_s_ss_gate__d_s_ss_d_environment__time;
480 const NekDouble d_dt_chaste_interface__hyperpolarisation_activated_current_y_gate__y = var_chaste_interface__hyperpolarisation_activated_current_y_gate__d_y_d_environment__time;
481 const NekDouble d_dt_chaste_interface__SR_Ca_release_channel__P_O1 = var_chaste_interface__SR_Ca_release_channel__d_P_O1_d_environment__time;
482 const NekDouble d_dt_chaste_interface__SR_Ca_release_channel__P_O2 = var_chaste_interface__SR_Ca_release_channel__d_P_O2_d_environment__time;
483 const NekDouble d_dt_chaste_interface__SR_Ca_release_channel__P_C1 = var_chaste_interface__SR_Ca_release_channel__d_P_C1_d_environment__time;
484 const NekDouble d_dt_chaste_interface__SR_Ca_release_channel__P_C2 = var_chaste_interface__SR_Ca_release_channel__d_P_C2_d_environment__time;
485 const NekDouble d_dt_chaste_interface__intracellular_and_SR_Ca_fluxes__HTRPNCa = var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_HTRPNCa_d_environment__time;
486 const NekDouble d_dt_chaste_interface__intracellular_and_SR_Ca_fluxes__LTRPNCa = var_chaste_interface__intracellular_and_SR_Ca_fluxes__d_LTRPNCa_d_environment__time;
487 const NekDouble d_dt_chaste_interface__intracellular_ion_concentrations__Na_i = var_chaste_interface__intracellular_ion_concentrations__d_Na_i_d_environment__time;
488 const NekDouble d_dt_chaste_interface__intracellular_ion_concentrations__Ca_i = var_chaste_interface__intracellular_ion_concentrations__d_Ca_i_d_environment__time;
489 const NekDouble d_dt_chaste_interface__intracellular_ion_concentrations__K_i = var_chaste_interface__intracellular_ion_concentrations__d_K_i_d_environment__time;
490 const NekDouble d_dt_chaste_interface__intracellular_ion_concentrations__Ca_ss = var_chaste_interface__intracellular_ion_concentrations__d_Ca_ss_d_environment__time;
491 const NekDouble d_dt_chaste_interface__intracellular_ion_concentrations__Ca_JSR = var_chaste_interface__intracellular_ion_concentrations__d_Ca_JSR_d_environment__time;
492 const NekDouble d_dt_chaste_interface__intracellular_ion_concentrations__Ca_NSR = var_chaste_interface__intracellular_ion_concentrations__d_Ca_NSR_d_environment__time;
494 const NekDouble var_membrane__Cm = 0.0001;
495 const NekDouble var_membrane__i_Na = var_sodium_current__i_Na;
496 const NekDouble var_membrane__i_Ca_L = var_L_type_Ca_channel__i_Ca_L;
497 const NekDouble var_membrane__i_t = var_Ca_independent_transient_outward_K_current__i_t;
498 const NekDouble var_membrane__i_ss = var_steady_state_outward_K_current__i_ss;
499 const NekDouble var_hyperpolarisation_activated_current__i_f = var_hyperpolarisation_activated_current__i_f_Na + var_hyperpolarisation_activated_current__i_f_K;
500 const NekDouble var_membrane__i_f = var_hyperpolarisation_activated_current__i_f;
501 const NekDouble var_membrane__i_K1 = var_inward_rectifier__i_K1;
502 const NekDouble var_background_currents__i_B = var_background_currents__i_B_Na + var_background_currents__i_B_Ca + var_background_currents__i_B_K;
503 const NekDouble var_membrane__i_B = var_background_currents__i_B;
504 const NekDouble var_membrane__i_NaK = var_sodium_potassium_pump__i_NaK;
505 const NekDouble var_membrane__i_NaCa = var_Na_Ca_ion_exchanger_current__i_NaCa;
506 const NekDouble var_membrane__i_Ca_P = var_sarcolemmal_calcium_pump_current__i_Ca_P;
507 const NekDouble var_chaste_interface__membrane__i_Stim = 0.0;
508 const NekDouble var_membrane__i_Stim_converter = var_chaste_interface__membrane__i_Stim;
509 const NekDouble var_membrane__chaste_interface__chaste_membrane_capacitance = 1.0;
510 const NekDouble var_membrane__i_Stim = 1000.0 * ((var_membrane__i_Stim_converter * var_membrane__Cm) / var_membrane__chaste_interface__chaste_membrane_capacitance);
511 const NekDouble var_membrane__d_V_d_environment__time = (-(var_membrane__i_Na + var_membrane__i_Ca_L + var_membrane__i_t + var_membrane__i_ss + var_membrane__i_f + var_membrane__i_K1 + var_membrane__i_B + var_membrane__i_NaK + var_membrane__i_NaCa + var_membrane__i_Ca_P + var_membrane__i_Stim)) / var_membrane__Cm;
512 const NekDouble var_chaste_interface__membrane__d_V_d_environment__time_converter = var_membrane__d_V_d_environment__time;
513 const NekDouble var_chaste_interface__membrane__d_V_d_environment__time = 0.001 * var_chaste_interface__membrane__d_V_d_environment__time_converter;
514 d_dt_chaste_interface__membrane__V = var_chaste_interface__membrane__d_V_d_environment__time;
515 outarray[0][i] = d_dt_chaste_interface__membrane__V;
516 outarray[1][i] = 0.0;
517 outarray[2][i] = d_dt_chaste_interface__sodium_current_m_gate__m;
518 outarray[3][i] = d_dt_chaste_interface__sodium_current_h_gate__h;
519 outarray[4][i] = d_dt_chaste_interface__sodium_current_j_gate__j;
520 outarray[5][i] = d_dt_chaste_interface__L_type_Ca_channel_d_gate__d;
521 outarray[6][i] = d_dt_chaste_interface__L_type_Ca_channel_f_11_gate__f_11;
522 outarray[7][i] = d_dt_chaste_interface__L_type_Ca_channel_f_12_gate__f_12;
523 outarray[8][i] = d_dt_chaste_interface__L_type_Ca_channel_Ca_inact_gate__Ca_inact;
524 outarray[9][i] = d_dt_chaste_interface__Ca_independent_transient_outward_K_current_r_gate__r;
525 outarray[10][i] = d_dt_chaste_interface__Ca_independent_transient_outward_K_current_s_gate__s;
526 outarray[11][i] = d_dt_chaste_interface__Ca_independent_transient_outward_K_current_s_slow_gate__s_slow;
527 outarray[12][i] = d_dt_chaste_interface__steady_state_outward_K_current_r_ss_gate__r_ss;
528 outarray[13][i] = d_dt_chaste_interface__steady_state_outward_K_current_s_ss_gate__s_ss;
529 outarray[14][i] = d_dt_chaste_interface__hyperpolarisation_activated_current_y_gate__y;
530 outarray[15][i] = d_dt_chaste_interface__SR_Ca_release_channel__P_O1;
531 outarray[16][i] = d_dt_chaste_interface__SR_Ca_release_channel__P_O2;
532 outarray[17][i] = d_dt_chaste_interface__SR_Ca_release_channel__P_C1;
533 outarray[18][i] = d_dt_chaste_interface__SR_Ca_release_channel__P_C2;
534 outarray[19][i] = d_dt_chaste_interface__intracellular_and_SR_Ca_fluxes__HTRPNCa;
535 outarray[20][i] = d_dt_chaste_interface__intracellular_and_SR_Ca_fluxes__LTRPNCa;
536 outarray[21][i] = d_dt_chaste_interface__intracellular_ion_concentrations__Na_i;
537 outarray[22][i] = d_dt_chaste_interface__intracellular_ion_concentrations__Ca_i;
538 outarray[23][i] = d_dt_chaste_interface__intracellular_ion_concentrations__K_i;
539 outarray[24][i] = d_dt_chaste_interface__intracellular_ion_concentrations__Ca_ss;
540 outarray[25][i] = d_dt_chaste_interface__intracellular_ion_concentrations__Ca_JSR;
541 outarray[26][i] = d_dt_chaste_interface__intracellular_ion_concentrations__Ca_NSR;
virtual void v_GenerateSummary(SummaryList &s)
Prints a summary of the model parameters.
static std::string className
Name of class.
int m_nq
Number of physical points.
std::vector< std::pair< std::string, std::string > > SummaryList
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
PanditGilesDemir03(const LibUtilities::SessionReaderSharedPtr &pSession, const int nq)
Constructor.
void AddSummaryItem(SummaryList &l, const std::string &name, const std::string &value)
Adds a summary item to the summary info list.
static CellModelSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const int nq)
Creates an instance of this class.
CellModelFactory & GetCellModelFactory()
virtual void v_Update(const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
Computes the reaction terms $f(u,v)$ and $g(u,v)$.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.