Generated Code
The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
There are a total of 71 entries in the algebraic variable array.
There are a total of 18 entries in each of the rate and state variable arrays.
There are a total of 73 entries in the constant variable array.
*/
/*
* VOI is time in component environment (second).
* STATES[0] is Vm in component Membrane (volt).
* ALGEBRAIC[17] is i_CaL in component i_CaL (A_per_F).
* ALGEBRAIC[49] is i_K1 in component i_K1 (A_per_F).
* ALGEBRAIC[50] is i_f in component i_f (A_per_F).
* ALGEBRAIC[4] is i_Na in component i_Na (A_per_F).
* ALGEBRAIC[32] is i_Kr in component i_Kr (A_per_F).
* ALGEBRAIC[41] is i_Ks in component i_Ks (A_per_F).
* ALGEBRAIC[58] is i_to in component i_to (A_per_F).
* ALGEBRAIC[57] is i_PCa in component i_PCa (A_per_F).
* ALGEBRAIC[55] is i_NaK in component i_NaK (A_per_F).
* ALGEBRAIC[56] is i_NaCa in component i_NaCa (A_per_F).
* ALGEBRAIC[54] is i_b_Ca in component i_b_Ca (A_per_F).
* ALGEBRAIC[53] is i_b_Na in component i_b_Na (A_per_F).
* ALGEBRAIC[0] is i_stim in component stim_mode (A_per_F).
* CONSTANTS[0] is Cm in component model_parameters (farad).
* CONSTANTS[1] is stim_flag in component stim_mode (dimensionless).
* CONSTANTS[2] is i_stim_Start in component stim_mode (second).
* CONSTANTS[3] is i_stim_End in component stim_mode (second).
* CONSTANTS[4] is i_stim_Amplitude in component stim_mode (ampere).
* CONSTANTS[5] is i_stim_frequency in component stim_mode (per_second).
* CONSTANTS[65] is i_stim_Period in component stim_mode (second).
* CONSTANTS[6] is i_stim_PulseDuration in component stim_mode (second).
* CONSTANTS[7] is TTX_3uM in component current_blockers (dimensionless).
* CONSTANTS[8] is TTX_10uM in component current_blockers (dimensionless).
* CONSTANTS[9] is TTX_30uM in component current_blockers (dimensionless).
* CONSTANTS[10] is E4031_30nM in component current_blockers (dimensionless).
* CONSTANTS[11] is E4031_100nM in component current_blockers (dimensionless).
* CONSTANTS[12] is nifed_3nM in component current_blockers (dimensionless).
* CONSTANTS[13] is nifed_10nM in component current_blockers (dimensionless).
* CONSTANTS[14] is nifed_30nM in component current_blockers (dimensionless).
* CONSTANTS[15] is nifed_100nM in component current_blockers (dimensionless).
* CONSTANTS[16] is Chromanol_iKs30 in component current_blockers (dimensionless).
* CONSTANTS[17] is Chromanol_iKs50 in component current_blockers (dimensionless).
* CONSTANTS[18] is Chromanol_iKs70 in component current_blockers (dimensionless).
* CONSTANTS[19] is Chromanol_iKs90 in component current_blockers (dimensionless).
* ALGEBRAIC[1] is E_Na in component electric_potentials (volt).
* CONSTANTS[66] is E_K in component electric_potentials (volt).
* ALGEBRAIC[2] is E_Ks in component electric_potentials (volt).
* ALGEBRAIC[3] is E_Ca in component electric_potentials (volt).
* CONSTANTS[20] is R in component model_parameters (joule_per_mole_kelvin).
* CONSTANTS[21] is T in component model_parameters (kelvin).
* CONSTANTS[22] is F in component model_parameters (coulomb_per_mole).
* STATES[1] is Nai in component sodium_dynamics (millimolar).
* CONSTANTS[23] is Nao in component model_parameters (millimolar).
* STATES[2] is Cai in component calcium_dynamics (millimolar).
* CONSTANTS[24] is Cao in component model_parameters (millimolar).
* CONSTANTS[25] is Ki in component model_parameters (millimolar).
* CONSTANTS[26] is Ko in component model_parameters (millimolar).
* CONSTANTS[27] is PkNa in component electric_potentials (dimensionless).
* CONSTANTS[28] is g_Na in component i_Na (S_per_F).
* STATES[3] is m in component i_Na_m_gate (dimensionless).
* STATES[4] is h in component i_Na_h_gate (dimensionless).
* STATES[5] is j in component i_Na_j_gate (dimensionless).
* CONSTANTS[67] is TTX_coeff in component i_Na (dimensionless).
* ALGEBRAIC[5] is m_inf in component i_Na_m_gate (dimensionless).
* ALGEBRAIC[8] is tau_m in component i_Na_m_gate (second).
* ALGEBRAIC[6] is alpha_m in component i_Na_m_gate (dimensionless).
* ALGEBRAIC[7] is beta_m in component i_Na_m_gate (dimensionless).
* ALGEBRAIC[9] is h_inf in component i_Na_h_gate (dimensionless).
* ALGEBRAIC[10] is alpha_h in component i_Na_h_gate (dimensionless).
* ALGEBRAIC[11] is beta_h in component i_Na_h_gate (dimensionless).
* ALGEBRAIC[12] is tau_h in component i_Na_h_gate (second).
* ALGEBRAIC[13] is j_inf in component i_Na_j_gate (dimensionless).
* ALGEBRAIC[14] is alpha_j in component i_Na_j_gate (dimensionless).
* ALGEBRAIC[15] is beta_j in component i_Na_j_gate (dimensionless).
* ALGEBRAIC[16] is tau_j in component i_Na_j_gate (second).
* CONSTANTS[29] is g_CaL in component i_CaL (metre_cube_per_F_per_s).
* STATES[6] is d in component i_CaL_d_gate (dimensionless).
* STATES[7] is f1 in component i_CaL_f1_gate (dimensionless).
* STATES[8] is f2 in component i_CaL_f2_gate (dimensionless).
* STATES[9] is fCa in component i_CaL_fCa_gate (dimensionless).
* CONSTANTS[68] is nifed_coeff in component i_CaL (dimensionless).
* ALGEBRAIC[18] is d_infinity in component i_CaL_d_gate (dimensionless).
* ALGEBRAIC[22] is tau_d in component i_CaL_d_gate (second).
* ALGEBRAIC[19] is alpha_d in component i_CaL_d_gate (dimensionless).
* ALGEBRAIC[20] is beta_d in component i_CaL_d_gate (dimensionless).
* ALGEBRAIC[21] is gamma_d in component i_CaL_d_gate (dimensionless).
* ALGEBRAIC[23] is f1_inf in component i_CaL_f1_gate (dimensionless).
* ALGEBRAIC[70] is tau_f1 in component i_CaL_f1_gate (second).
* ALGEBRAIC[24] is constf1 in component i_CaL_f1_gate (dimensionless).
* ALGEBRAIC[25] is f2_inf in component i_CaL_f2_gate (dimensionless).
* ALGEBRAIC[26] is tau_f2 in component i_CaL_f2_gate (second).
* CONSTANTS[69] is constf2 in component i_CaL_f2_gate (dimensionless).
* ALGEBRAIC[31] is constfCa in component i_CaL_fCa_gate (dimensionless).
* ALGEBRAIC[27] is alpha_fCa in component i_CaL_fCa_gate (dimensionless).
* ALGEBRAIC[28] is beta_fCa in component i_CaL_fCa_gate (dimensionless).
* ALGEBRAIC[29] is gamma_fCa in component i_CaL_fCa_gate (dimensionless).
* ALGEBRAIC[30] is fCa_inf in component i_CaL_fCa_gate (dimensionless).
* CONSTANTS[30] is tau_fCa in component i_CaL_fCa_gate (second).
* CONSTANTS[31] is g_Kr in component i_Kr (S_per_F).
* STATES[10] is Xr1 in component i_Kr_Xr1_gate (dimensionless).
* STATES[11] is Xr2 in component i_Kr_Xr2_gate (dimensionless).
* CONSTANTS[70] is E4031_coeff in component i_Kr (dimensionless).
* ALGEBRAIC[33] is Xr1_inf in component i_Kr_Xr1_gate (dimensionless).
* ALGEBRAIC[34] is alpha_Xr1 in component i_Kr_Xr1_gate (dimensionless).
* ALGEBRAIC[35] is beta_Xr1 in component i_Kr_Xr1_gate (dimensionless).
* ALGEBRAIC[36] is tau_Xr1 in component i_Kr_Xr1_gate (second).
* CONSTANTS[32] is L0 in component i_Kr_Xr1_gate (dimensionless).
* CONSTANTS[71] is V_half in component i_Kr_Xr1_gate (millivolt).
* CONSTANTS[33] is Q in component i_Kr_Xr1_gate (dimensionless).
* ALGEBRAIC[37] is Xr2_infinity in component i_Kr_Xr2_gate (dimensionless).
* ALGEBRAIC[38] is alpha_Xr2 in component i_Kr_Xr2_gate (dimensionless).
* ALGEBRAIC[39] is beta_Xr2 in component i_Kr_Xr2_gate (dimensionless).
* ALGEBRAIC[40] is tau_Xr2 in component i_Kr_Xr2_gate (second).
* CONSTANTS[34] is g_Ks in component i_Ks (S_per_F).
* STATES[12] is Xs in component i_Ks_Xs_gate (dimensionless).
* CONSTANTS[72] is Chromanol_coeff in component i_Ks (dimensionless).
* ALGEBRAIC[42] is Xs_infinity in component i_Ks_Xs_gate (dimensionless).
* ALGEBRAIC[43] is alpha_Xs in component i_Ks_Xs_gate (dimensionless).
* ALGEBRAIC[44] is beta_Xs in component i_Ks_Xs_gate (dimensionless).
* ALGEBRAIC[45] is tau_Xs in component i_Ks_Xs_gate (second).
* CONSTANTS[35] is g_K1 in component i_K1 (S_per_F).
* ALGEBRAIC[48] is XK1_inf in component i_K1 (dimensionless).
* ALGEBRAIC[46] is alpha_K1 in component i_K1 (dimensionless).
* ALGEBRAIC[47] is beta_K1 in component i_K1 (dimensionless).
* CONSTANTS[36] is g_f in component i_f (S_per_F).
* CONSTANTS[37] is E_f in component i_f (volt).
* STATES[13] is Xf in component i_f_Xf_gate (dimensionless).
* ALGEBRAIC[51] is Xf_infinity in component i_f_Xf_gate (dimensionless).
* ALGEBRAIC[52] is tau_Xf in component i_f_Xf_gate (second).
* CONSTANTS[38] is g_b_Na in component i_b_Na (S_per_F).
* CONSTANTS[39] is g_b_Ca in component i_b_Ca (S_per_F).
* CONSTANTS[40] is Km_K in component i_NaK (millimolar).
* CONSTANTS[41] is Km_Na in component i_NaK (millimolar).
* CONSTANTS[42] is PNaK in component i_NaK (A_per_F).
* CONSTANTS[43] is kNaCa in component i_NaCa (A_per_F).
* CONSTANTS[44] is alpha in component i_NaCa (dimensionless).
* CONSTANTS[45] is gamma in component i_NaCa (dimensionless).
* CONSTANTS[46] is Ksat in component i_NaCa (dimensionless).
* CONSTANTS[47] is KmCa in component i_NaCa (millimolar).
* CONSTANTS[48] is KmNai in component i_NaCa (millimolar).
* CONSTANTS[49] is g_PCa in component i_PCa (A_per_F).
* CONSTANTS[50] is KPCa in component i_PCa (millimolar).
* CONSTANTS[51] is g_to in component i_to (S_per_F).
* STATES[14] is q in component i_to_q_gate (dimensionless).
* STATES[15] is r in component i_to_r_gate (dimensionless).
* ALGEBRAIC[59] is q_inf in component i_to_q_gate (dimensionless).
* ALGEBRAIC[60] is tau_q in component i_to_q_gate (second).
* ALGEBRAIC[61] is r_inf in component i_to_r_gate (dimensionless).
* ALGEBRAIC[62] is tau_r in component i_to_r_gate (second).
* CONSTANTS[52] is Vc in component model_parameters (micrometre_cube).
* CONSTANTS[53] is V_SR in component model_parameters (micrometre_cube).
* STATES[16] is Ca_SR in component calcium_dynamics (millimolar).
* CONSTANTS[54] is a_rel in component calcium_dynamics (millimolar_per_second).
* CONSTANTS[55] is b_rel in component calcium_dynamics (millimolar).
* CONSTANTS[56] is c_rel in component calcium_dynamics (millimolar_per_second).
* STATES[17] is g in component calcium_dynamics (dimensionless).
* CONSTANTS[57] is tau_g in component calcium_dynamics (second).
* ALGEBRAIC[66] is g_inf in component calcium_dynamics (dimensionless).
* CONSTANTS[58] is Kup in component calcium_dynamics (millimolar).
* CONSTANTS[59] is Buf_C in component calcium_dynamics (millimolar).
* CONSTANTS[60] is Buf_SR in component calcium_dynamics (millimolar).
* CONSTANTS[61] is Kbuf_C in component calcium_dynamics (millimolar).
* CONSTANTS[62] is Kbuf_SR in component calcium_dynamics (millimolar).
* ALGEBRAIC[68] is Cai_bufc in component calcium_dynamics (dimensionless).
* ALGEBRAIC[69] is Ca_SR_bufSR in component calcium_dynamics (dimensionless).
* CONSTANTS[63] is VmaxUp in component calcium_dynamics (millimolar_per_second).
* ALGEBRAIC[67] is const2 in component calcium_dynamics (dimensionless).
* CONSTANTS[64] is V_leak in component calcium_dynamics (per_second).
* ALGEBRAIC[63] is i_rel in component calcium_dynamics (millimolar_per_second).
* ALGEBRAIC[64] is i_up in component calcium_dynamics (millimolar_per_second).
* ALGEBRAIC[65] is i_leak in component calcium_dynamics (millimolar_per_second).
* RATES[0] is d/dt Vm in component Membrane (volt).
* RATES[3] is d/dt m in component i_Na_m_gate (dimensionless).
* RATES[4] is d/dt h in component i_Na_h_gate (dimensionless).
* RATES[5] is d/dt j in component i_Na_j_gate (dimensionless).
* RATES[6] is d/dt d in component i_CaL_d_gate (dimensionless).
* RATES[7] is d/dt f1 in component i_CaL_f1_gate (dimensionless).
* RATES[8] is d/dt f2 in component i_CaL_f2_gate (dimensionless).
* RATES[9] is d/dt fCa in component i_CaL_fCa_gate (dimensionless).
* RATES[10] is d/dt Xr1 in component i_Kr_Xr1_gate (dimensionless).
* RATES[11] is d/dt Xr2 in component i_Kr_Xr2_gate (dimensionless).
* RATES[12] is d/dt Xs in component i_Ks_Xs_gate (dimensionless).
* RATES[13] is d/dt Xf in component i_f_Xf_gate (dimensionless).
* RATES[14] is d/dt q in component i_to_q_gate (dimensionless).
* RATES[15] is d/dt r in component i_to_r_gate (dimensionless).
* RATES[1] is d/dt Nai in component sodium_dynamics (millimolar).
* RATES[17] is d/dt g in component calcium_dynamics (dimensionless).
* RATES[2] is d/dt Cai in component calcium_dynamics (millimolar).
* RATES[16] is d/dt Ca_SR in component calcium_dynamics (millimolar).
* There are a total of 14 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -0.0743340057623841;
CONSTANTS[0] = 9.87109e-11;
CONSTANTS[1] = 0;
CONSTANTS[2] = 0;
CONSTANTS[3] = 800;
CONSTANTS[4] = 5.5e-10;
CONSTANTS[5] = 60;
CONSTANTS[6] = 0.005;
CONSTANTS[7] = 0;
CONSTANTS[8] = 0;
CONSTANTS[9] = 0;
CONSTANTS[10] = 0;
CONSTANTS[11] = 0;
CONSTANTS[12] = 0;
CONSTANTS[13] = 0;
CONSTANTS[14] = 0;
CONSTANTS[15] = 0;
CONSTANTS[16] = 0;
CONSTANTS[17] = 0;
CONSTANTS[18] = 0;
CONSTANTS[19] = 0;
CONSTANTS[20] = 8.314472;
CONSTANTS[21] = 310;
CONSTANTS[22] = 96485.3415;
STATES[1] = 10.9248496211574;
CONSTANTS[23] = 151;
STATES[2] = 1.80773974140477e-5;
CONSTANTS[24] = 1.8;
CONSTANTS[25] = 150;
CONSTANTS[26] = 5.4;
CONSTANTS[27] = 0.03;
CONSTANTS[28] = 3671.2302;
STATES[3] = 0.102953468725004;
STATES[4] = 0.786926637881461;
STATES[5] = 0.253943221774722;
CONSTANTS[29] = 8.635702e-5;
STATES[6] = 8.96088425225182e-5;
STATES[7] = 0.970411811263976;
STATES[8] = 0.999965815466749;
STATES[9] = 0.998925296531804;
CONSTANTS[30] = 0.002;
CONSTANTS[31] = 29.8667;
STATES[10] = 0.00778547011240132;
STATES[11] = 0.432162576531617;
CONSTANTS[32] = 0.025;
CONSTANTS[33] = 2.3;
CONSTANTS[34] = 2.041;
STATES[12] = 0.0322944866983666;
CONSTANTS[35] = 28.1492;
CONSTANTS[36] = 30.10312;
CONSTANTS[37] = -0.017;
STATES[13] = 0.100615100568753;
CONSTANTS[38] = 0.9;
CONSTANTS[39] = 0.69264;
CONSTANTS[40] = 1;
CONSTANTS[41] = 40;
CONSTANTS[42] = 1.841424;
CONSTANTS[43] = 4900;
CONSTANTS[44] = 2.8571432;
CONSTANTS[45] = 0.35;
CONSTANTS[46] = 0.1;
CONSTANTS[47] = 1.38;
CONSTANTS[48] = 87.5;
CONSTANTS[49] = 0.4125;
CONSTANTS[50] = 0.0005;
CONSTANTS[51] = 29.9038;
STATES[14] = 0.839295925773219;
STATES[15] = 0.00573289893326379;
CONSTANTS[52] = 8800;
CONSTANTS[53] = 583.73;
STATES[16] = 0.2734234751931;
CONSTANTS[54] = 16.464;
CONSTANTS[55] = 0.25;
CONSTANTS[56] = 8.232;
STATES[17] = 0.999999981028517;
CONSTANTS[57] = 0.002;
CONSTANTS[58] = 0.00025;
CONSTANTS[59] = 0.25;
CONSTANTS[60] = 10;
CONSTANTS[61] = 0.001;
CONSTANTS[62] = 0.3;
CONSTANTS[63] = 0.56064;
CONSTANTS[64] = 0.00044444;
CONSTANTS[65] = 60.0000/CONSTANTS[5];
CONSTANTS[66] = (( CONSTANTS[20]*CONSTANTS[21])/CONSTANTS[22])*log(CONSTANTS[26]/CONSTANTS[25]);
CONSTANTS[67] = (CONSTANTS[7]==1.00000 ? 0.180000 : CONSTANTS[8]==1.00000 ? 0.0600000 : CONSTANTS[9]==1.00000 ? 0.0200000 : 1.00000);
CONSTANTS[68] = (CONSTANTS[12]==1.00000 ? 0.930000 : CONSTANTS[13]==1.00000 ? 0.790000 : CONSTANTS[14]==1.00000 ? 0.560000 : CONSTANTS[15]==1.00000 ? 0.280000 : 1.00000);
CONSTANTS[69] = 1.00000;
CONSTANTS[70] = (CONSTANTS[10]==1.00000 ? 0.770000 : CONSTANTS[11]==1.00000 ? 0.500000 : 1.00000);
CONSTANTS[71] = 1000.00*( (( - CONSTANTS[20]*CONSTANTS[21])/( CONSTANTS[22]*CONSTANTS[33]))*log(pow(1.00000+CONSTANTS[24]/2.60000, 4.00000)/( CONSTANTS[32]*pow(1.00000+CONSTANTS[24]/0.580000, 4.00000))) - 0.0190000);
CONSTANTS[72] = (CONSTANTS[16]==1.00000 ? 0.700000 : CONSTANTS[17]==1.00000 ? 0.500000 : CONSTANTS[18]==1.00000 ? 0.300000 : CONSTANTS[19]==1.00000 ? 0.100000 : 1.00000);
RATES[0] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[13] = 0.1001;
RATES[14] = 0.1001;
RATES[15] = 0.1001;
RATES[1] = 0.1001;
RATES[17] = 0.1001;
RATES[2] = 0.1001;
RATES[16] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - - ((ALGEBRAIC[49]+(ALGEBRAIC[58]+(ALGEBRAIC[32]+(ALGEBRAIC[41]+(ALGEBRAIC[17]+(ALGEBRAIC[55]+(ALGEBRAIC[4]+(ALGEBRAIC[56]+(ALGEBRAIC[57]+(ALGEBRAIC[50]+(ALGEBRAIC[53]+ALGEBRAIC[54]))))))))))) - ALGEBRAIC[0]);
resid[1] = RATES[3] - (ALGEBRAIC[5] - STATES[3])/ALGEBRAIC[8];
resid[2] = RATES[4] - (ALGEBRAIC[9] - STATES[4])/ALGEBRAIC[12];
resid[3] = RATES[5] - (ALGEBRAIC[13] - STATES[5])/ALGEBRAIC[16];
resid[4] = RATES[6] - (ALGEBRAIC[18] - STATES[6])/ALGEBRAIC[22];
resid[5] = RATES[7] - (ALGEBRAIC[23] - STATES[7])/ALGEBRAIC[70];
resid[6] = RATES[8] - (ALGEBRAIC[25] - STATES[8])/ALGEBRAIC[26];
resid[7] = RATES[9] - ( ALGEBRAIC[31]*(ALGEBRAIC[30] - STATES[9]))/CONSTANTS[30];
resid[8] = RATES[10] - (ALGEBRAIC[33] - STATES[10])/ALGEBRAIC[36];
resid[9] = RATES[11] - (ALGEBRAIC[37] - STATES[11])/ALGEBRAIC[40];
resid[10] = RATES[12] - (ALGEBRAIC[42] - STATES[12])/ALGEBRAIC[45];
resid[11] = RATES[13] - (ALGEBRAIC[51] - STATES[13])/ALGEBRAIC[52];
resid[12] = RATES[14] - (ALGEBRAIC[59] - STATES[14])/ALGEBRAIC[60];
resid[13] = RATES[15] - (ALGEBRAIC[61] - STATES[15])/ALGEBRAIC[62];
resid[14] = RATES[1] - ( - CONSTANTS[0]*(ALGEBRAIC[4]+(ALGEBRAIC[53]+( 3.00000*ALGEBRAIC[55]+ 3.00000*ALGEBRAIC[56]))))/( CONSTANTS[22]*( CONSTANTS[52]*1.00000e-18));
resid[15] = RATES[17] - ( ALGEBRAIC[67]*(ALGEBRAIC[66] - STATES[17]))/CONSTANTS[57];
resid[16] = RATES[2] - ALGEBRAIC[68]*(((ALGEBRAIC[65] - ALGEBRAIC[64])+ALGEBRAIC[63]) - ( ((ALGEBRAIC[17]+(ALGEBRAIC[54]+ALGEBRAIC[57])) - 2.00000*ALGEBRAIC[56])*CONSTANTS[0])/( 2.00000*( CONSTANTS[52]*( CONSTANTS[22]*1.00000e-18))));
resid[17] = RATES[16] - (( ALGEBRAIC[69]*CONSTANTS[52])/CONSTANTS[53])*(ALGEBRAIC[64] - (ALGEBRAIC[63]+ALGEBRAIC[65]));
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (CONDVAR[0]>=0.00000&&(CONDVAR[1]<=0.00000&&CONDVAR[2]<=0.00000) ? ( CONSTANTS[1]*CONSTANTS[4])/CONSTANTS[0] : 0.00000);
ALGEBRAIC[1] = (( CONSTANTS[20]*CONSTANTS[21])/CONSTANTS[22])*log(CONSTANTS[23]/STATES[1]);
ALGEBRAIC[4] = CONSTANTS[67]*( CONSTANTS[28]*( pow(STATES[3], 3.00000)*( STATES[4]*( STATES[5]*(STATES[0] - ALGEBRAIC[1])))));
ALGEBRAIC[5] = 1.00000/pow(1.00000+exp(( - STATES[0]*1000.00 - 34.1000)/5.90000), 1.00000/3.00000);
ALGEBRAIC[6] = 1.00000/(1.00000+exp(( - STATES[0]*1000.00 - 60.0000)/5.00000));
ALGEBRAIC[7] = 0.100000/(1.00000+exp(( STATES[0]*1000.00+35.0000)/5.00000))+0.100000/(1.00000+exp(( STATES[0]*1000.00 - 50.0000)/200.000));
ALGEBRAIC[8] = ( 1.00000*( ALGEBRAIC[6]*ALGEBRAIC[7]))/1000.00;
ALGEBRAIC[9] = 1.00000/ pow((1.00000+exp(( STATES[0]*1000.00+72.1000)/5.70000)), 1.0 / 2);
ALGEBRAIC[10] = (CONDVAR[3]<0.00000 ? 0.0570000*exp(- ( STATES[0]*1000.00+80.0000)/6.80000) : 0.00000);
ALGEBRAIC[11] = (CONDVAR[4]<0.00000 ? 2.70000*exp( 0.0790000*( STATES[0]*1000.00))+ 3.10000*( pow(10.0000, 5.00000)*exp( 0.348500*( STATES[0]*1000.00))) : 0.770000/( 0.130000*(1.00000+exp(( STATES[0]*1000.00+10.6600)/- 11.1000))));
ALGEBRAIC[12] = (CONDVAR[5]<0.00000 ? 1.50000/( (ALGEBRAIC[10]+ALGEBRAIC[11])*1000.00) : 2.54200/1000.00);
ALGEBRAIC[13] = 1.00000/ pow((1.00000+exp(( STATES[0]*1000.00+72.1000)/5.70000)), 1.0 / 2);
ALGEBRAIC[14] = (CONDVAR[6]<0.00000 ? ( ( - 25428.0*exp( 0.244400*( STATES[0]*1000.00)) - 6.94800*( pow(10.0000, - 6.00000)*exp( - 0.0439100*( STATES[0]*1000.00))))*( STATES[0]*1000.00+37.7800))/(1.00000+exp( 0.311000*( STATES[0]*1000.00+79.2300))) : 0.00000);
ALGEBRAIC[15] = (CONDVAR[7]<0.00000 ? ( 0.0242400*exp( - 0.0105200*( STATES[0]*1000.00)))/(1.00000+exp( - 0.137800*( STATES[0]*1000.00+40.1400))) : ( 0.600000*exp( 0.0570000*( STATES[0]*1000.00)))/(1.00000+exp( - 0.100000*( STATES[0]*1000.00+32.0000))));
ALGEBRAIC[16] = 7.00000/( (ALGEBRAIC[14]+ALGEBRAIC[15])*1000.00);
ALGEBRAIC[17] = (( (( CONSTANTS[29]*( 4.00000*( STATES[0]*pow(CONSTANTS[22], 2.00000))))/( CONSTANTS[20]*CONSTANTS[21]))*( STATES[2]*exp(( 2.00000*( STATES[0]*CONSTANTS[22]))/( CONSTANTS[20]*CONSTANTS[21])) - 0.341000*CONSTANTS[24]))/(exp(( 2.00000*( STATES[0]*CONSTANTS[22]))/( CONSTANTS[20]*CONSTANTS[21])) - 1.00000))*( STATES[6]*( STATES[7]*( STATES[8]*STATES[9])));
ALGEBRAIC[18] = 1.00000/(1.00000+exp(- ( STATES[0]*1000.00+9.10000)/7.00000));
ALGEBRAIC[19] = 0.250000+1.40000/(1.00000+exp(( - STATES[0]*1000.00 - 35.0000)/13.0000));
ALGEBRAIC[20] = 1.40000/(1.00000+exp(( STATES[0]*1000.00+5.00000)/5.00000));
ALGEBRAIC[21] = 1.00000/(1.00000+exp(( - STATES[0]*1000.00+50.0000)/20.0000));
ALGEBRAIC[22] = ( ( ALGEBRAIC[19]*ALGEBRAIC[20]+ALGEBRAIC[21])*1.00000)/1000.00;
ALGEBRAIC[23] = 1.00000/(1.00000+exp(( STATES[0]*1000.00+26.0000)/3.00000));
ALGEBRAIC[25] = 0.330000+0.670000/(1.00000+exp(( STATES[0]*1000.00+35.0000)/4.00000));
ALGEBRAIC[26] = ( ( 600.000*exp(- pow( STATES[0]*1000.00+25.0000, 2.00000)/170.000)+(31.0000/(1.00000+exp((25.0000 - STATES[0]*1000.00)/10.0000))+16.0000/(1.00000+exp((30.0000+ STATES[0]*1000.00)/10.0000))))*CONSTANTS[69])/1000.00;
ALGEBRAIC[27] = 1.00000/(1.00000+pow(STATES[2]/0.000600000, 8.00000));
ALGEBRAIC[28] = 0.100000/(1.00000+exp((STATES[2] - 0.000900000)/0.000100000));
ALGEBRAIC[29] = 0.300000/(1.00000+exp((STATES[2] - 0.000750000)/0.000800000));
ALGEBRAIC[30] = (ALGEBRAIC[27]+(ALGEBRAIC[28]+ALGEBRAIC[29]))/1.31560;
ALGEBRAIC[31] = (CONDVAR[9]>0.00000&&CONDVAR[10]>0.00000 ? 0.00000 : 1.00000);
ALGEBRAIC[32] = CONSTANTS[70]*( CONSTANTS[31]*( (STATES[0] - CONSTANTS[66])*( STATES[10]*( STATES[11]* pow((CONSTANTS[26]/5.40000), 1.0 / 2)))));
ALGEBRAIC[33] = 1.00000/(1.00000+exp((CONSTANTS[71] - STATES[0]*1000.00)/4.90000));
ALGEBRAIC[34] = 450.000/(1.00000+exp((- 45.0000 - STATES[0]*1000.00)/10.0000));
ALGEBRAIC[35] = 6.00000/(1.00000+exp((30.0000+ STATES[0]*1000.00)/11.5000));
ALGEBRAIC[36] = ( 1.00000*( ALGEBRAIC[34]*ALGEBRAIC[35]))/1000.00;
ALGEBRAIC[37] = 1.00000/(1.00000+exp(( STATES[0]*1000.00+88.0000)/50.0000));
ALGEBRAIC[38] = 3.00000/(1.00000+exp((- 60.0000 - STATES[0]*1000.00)/20.0000));
ALGEBRAIC[39] = 1.12000/(1.00000+exp((- 60.0000+ STATES[0]*1000.00)/20.0000));
ALGEBRAIC[40] = ( 1.00000*( ALGEBRAIC[38]*ALGEBRAIC[39]))/1000.00;
ALGEBRAIC[2] = (( CONSTANTS[20]*CONSTANTS[21])/CONSTANTS[22])*log((CONSTANTS[26]+ CONSTANTS[27]*CONSTANTS[23])/(CONSTANTS[25]+ CONSTANTS[27]*STATES[1]));
ALGEBRAIC[41] = CONSTANTS[72]*( CONSTANTS[34]*( (STATES[0] - ALGEBRAIC[2])*( pow(STATES[12], 2.00000)*(1.00000+0.600000/(1.00000+pow(( 3.80000*1.00000e-05)/STATES[2], 1.40000))))));
ALGEBRAIC[42] = 1.00000/(1.00000+exp(( - STATES[0]*1000.00 - 20.0000)/16.0000));
ALGEBRAIC[43] = 1100.00/ pow((1.00000+exp((- 10.0000 - STATES[0]*1000.00)/6.00000)), 1.0 / 2);
ALGEBRAIC[44] = 1.00000/(1.00000+exp((- 60.0000+ STATES[0]*1000.00)/20.0000));
ALGEBRAIC[45] = ( 1.00000*( ALGEBRAIC[43]*ALGEBRAIC[44]))/1000.00;
ALGEBRAIC[46] = 3.91000/(1.00000+exp( 0.594200*(( STATES[0]*1000.00 - CONSTANTS[66]*1000.00) - 200.000)));
ALGEBRAIC[47] = ( - 1.50900*exp( 0.000200000*(( STATES[0]*1000.00 - CONSTANTS[66]*1000.00)+100.000))+exp( 0.588600*(( STATES[0]*1000.00 - CONSTANTS[66]*1000.00) - 10.0000)))/(1.00000+exp( 0.454700*( STATES[0]*1000.00 - CONSTANTS[66]*1000.00)));
ALGEBRAIC[48] = ALGEBRAIC[46]/(ALGEBRAIC[46]+ALGEBRAIC[47]);
ALGEBRAIC[49] = CONSTANTS[35]*( ALGEBRAIC[48]*( (STATES[0] - CONSTANTS[66])* pow((CONSTANTS[26]/5.40000), 1.0 / 2)));
ALGEBRAIC[50] = CONSTANTS[36]*( STATES[13]*(STATES[0] - CONSTANTS[37]));
ALGEBRAIC[51] = 1.00000/(1.00000+exp(( STATES[0]*1000.00+77.8500)/5.00000));
ALGEBRAIC[52] = (1900.00/(1.00000+exp(( STATES[0]*1000.00+15.0000)/10.0000)))/1000.00;
ALGEBRAIC[53] = CONSTANTS[38]*(STATES[0] - ALGEBRAIC[1]);
ALGEBRAIC[3] = (( 0.500000*( CONSTANTS[20]*CONSTANTS[21]))/CONSTANTS[22])*log(CONSTANTS[24]/STATES[2]);
ALGEBRAIC[54] = CONSTANTS[39]*(STATES[0] - ALGEBRAIC[3]);
ALGEBRAIC[55] = (( (( CONSTANTS[42]*CONSTANTS[26])/(CONSTANTS[26]+CONSTANTS[40]))*STATES[1])/(STATES[1]+CONSTANTS[41]))/(1.00000+( 0.124500*exp(( - 0.100000*( STATES[0]*CONSTANTS[22]))/( CONSTANTS[20]*CONSTANTS[21]))+ 0.0353000*exp(( - STATES[0]*CONSTANTS[22])/( CONSTANTS[20]*CONSTANTS[21]))));
ALGEBRAIC[56] = ( CONSTANTS[43]*( exp(( CONSTANTS[45]*( STATES[0]*CONSTANTS[22]))/( CONSTANTS[20]*CONSTANTS[21]))*( pow(STATES[1], 3.00000)*CONSTANTS[24]) - exp(( (CONSTANTS[45] - 1.00000)*( STATES[0]*CONSTANTS[22]))/( CONSTANTS[20]*CONSTANTS[21]))*( pow(CONSTANTS[23], 3.00000)*( STATES[2]*CONSTANTS[44]))))/( (pow(CONSTANTS[48], 3.00000)+pow(CONSTANTS[23], 3.00000))*( (CONSTANTS[47]+CONSTANTS[24])*(1.00000+ CONSTANTS[46]*exp(( (CONSTANTS[45] - 1.00000)*( STATES[0]*CONSTANTS[22]))/( CONSTANTS[20]*CONSTANTS[21])))));
ALGEBRAIC[57] = ( CONSTANTS[49]*STATES[2])/(STATES[2]+CONSTANTS[50]);
ALGEBRAIC[58] = CONSTANTS[51]*( (STATES[0] - CONSTANTS[66])*( STATES[14]*STATES[15]));
ALGEBRAIC[59] = 1.00000/(1.00000+exp(( STATES[0]*1000.00+53.0000)/13.0000));
ALGEBRAIC[60] = (6.06000+39.1020/( 0.570000*exp( - 0.0800000*( STATES[0]*1000.00+44.0000))+ 0.0650000*exp( 0.100000*( STATES[0]*1000.00+45.9300))))/1000.00;
ALGEBRAIC[61] = 1.00000/(1.00000+exp(- ( STATES[0]*1000.00 - 22.3000)/18.7500));
ALGEBRAIC[62] = (2.75352+14.4052/( 1.03700*exp( 0.0900000*( STATES[0]*1000.00+30.6100))+ 0.369000*exp( - 0.120000*( STATES[0]*1000.00+23.8400))))/1000.00;
ALGEBRAIC[63] = (CONSTANTS[56]+( CONSTANTS[54]*pow(STATES[16], 2.00000))/(pow(CONSTANTS[55], 2.00000)+pow(STATES[16], 2.00000)))*( STATES[6]*( STATES[17]*0.0411000));
ALGEBRAIC[64] = CONSTANTS[63]/(1.00000+pow(CONSTANTS[58], 2.00000)/pow(STATES[2], 2.00000));
ALGEBRAIC[65] = (STATES[16] - STATES[2])*CONSTANTS[64];
ALGEBRAIC[66] = (CONDVAR[11]<=0.00000 ? 1.00000/(1.00000+pow(STATES[2]/0.000350000, 6.00000)) : 1.00000/(1.00000+pow(STATES[2]/0.000350000, 16.0000)));
ALGEBRAIC[67] = (CONDVAR[12]>0.00000&&CONDVAR[13]>0.00000 ? 0.00000 : 1.00000);
ALGEBRAIC[68] = 1.00000/(1.00000+( CONSTANTS[59]*CONSTANTS[61])/pow(STATES[2]+CONSTANTS[61], 2.00000));
ALGEBRAIC[69] = 1.00000/(1.00000+( CONSTANTS[60]*CONSTANTS[62])/pow(STATES[16]+CONSTANTS[62], 2.00000));
ALGEBRAIC[24] = (CONDVAR[8]>0.00000 ? 1.00000+ 1433.00*(STATES[2] - 50.0000*1.00000e-06) : 1.00000);
ALGEBRAIC[70] = ( (20.0000+( 1102.50*exp(- pow(pow( STATES[0]*1000.00+27.0000, 2.00000)/15.0000, 2.00000))+(200.000/(1.00000+exp((13.0000 - STATES[0]*1000.00)/10.0000))+180.000/(1.00000+exp((30.0000+ STATES[0]*1000.00)/10.0000)))))*ALGEBRAIC[24])/1000.00;
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
SI[8] = 1.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[13] = 1.0;
SI[14] = 1.0;
SI[15] = 1.0;
SI[16] = 1.0;
SI[17] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - CONSTANTS[2];
CONDVAR[1] = VOI - CONSTANTS[3];
CONDVAR[2] = ((VOI - CONSTANTS[2]) - floor((VOI - CONSTANTS[2])/CONSTANTS[65])*CONSTANTS[65]) - CONSTANTS[6];
CONDVAR[3] = STATES[0] - - 0.0400000;
CONDVAR[4] = STATES[0] - - 0.0400000;
CONDVAR[5] = STATES[0] - - 0.0400000;
CONDVAR[6] = STATES[0] - - 0.0400000;
CONDVAR[7] = STATES[0] - - 0.0400000;
CONDVAR[8] = (ALGEBRAIC[23] - STATES[7]) - 0.00000;
CONDVAR[9] = STATES[0] - - 0.0600000;
CONDVAR[10] = ALGEBRAIC[30] - STATES[9];
CONDVAR[11] = STATES[2] - 0.000350000;
CONDVAR[12] = ALGEBRAIC[66] - STATES[17];
CONDVAR[13] = STATES[0] - - 0.0600000;
}