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 42 entries in the algebraic variable array.
There are a total of 27 entries in each of the rate and state variable arrays.
There are a total of 64 entries in the constant variable array.
*/
/*
* VOI is time in component Environment (second).
* CONSTANTS[0] is O2_x in component Environment (molar).
* CONSTANTS[1] is X_AtC in component Environment (mole_per_second_per_l_cell).
* CONSTANTS[2] is RT in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[3] is F in component fixed_parameters (kilojoule_per_mole_per_millivolt).
* CONSTANTS[4] is dG_C1o in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[5] is dG_C3o in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[6] is dG_C4o in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[35] is dG_C1op in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[36] is dG_C3op in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[37] is dG_C4op in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[7] is dG_F1o in component fixed_parameters (kilojoule_per_mole).
* ALGEBRAIC[34] is dG_H in component fixed_parameters (kilojoule_per_mole).
* CONSTANTS[50] is k_dHPI in component fixed_parameters (molar).
* CONSTANTS[51] is k_dHatp in component fixed_parameters (molar).
* CONSTANTS[52] is k_dHadp in component fixed_parameters (molar).
* CONSTANTS[8] is K_DT in component fixed_parameters (molar).
* CONSTANTS[9] is K_DD in component fixed_parameters (molar).
* CONSTANTS[10] is K_AK in component fixed_parameters (dimensionless).
* CONSTANTS[11] is K_CK in component fixed_parameters (per_molar).
* CONSTANTS[56] is W_m in component fixed_parameters (l_water_per_l_mito).
* CONSTANTS[59] is W_c in component fixed_parameters (l_water_per_l_cyto).
* CONSTANTS[57] is W_x in component fixed_parameters (l_water_per_l_mito).
* CONSTANTS[58] is W_i in component fixed_parameters (l_water_per_l_mito).
* CONSTANTS[12] is V_cyto in component fixed_parameters (cyto_per_cell).
* CONSTANTS[13] is V_mito in component fixed_parameters (mito_per_cell).
* CONSTANTS[61] is Rm_cyto in component fixed_parameters (mito_per_cyto).
* CONSTANTS[54] is Rm_cell in component fixed_parameters (mito_per_cell).
* CONSTANTS[62] is Rc_cell in component fixed_parameters (cyto_per_cell).
* CONSTANTS[14] is n_A in component fixed_parameters (dimensionless).
* CONSTANTS[15] is C_tot in component fixed_parameters (molar).
* CONSTANTS[16] is CR_tot in component fixed_parameters (molar).
* CONSTANTS[17] is Q_tot in component fixed_parameters (molar).
* CONSTANTS[18] is NAD_tot in component fixed_parameters (molar).
* CONSTANTS[19] is pH_C in component fixed_parameters (dimensionless).
* CONSTANTS[53] is H_c in component fixed_parameters (molar).
* CONSTANTS[20] is K_c in component fixed_parameters (molar).
* CONSTANTS[60] is H_i in component fixed_parameters (molar).
* CONSTANTS[55] is K_i in component fixed_parameters (molar).
* ALGEBRAIC[5] is H_x in component dH_x_dt (molar).
* STATES[0] is deltaPsi in component ddeltaPsi_dt (millivolt).
* CONSTANTS[21] is k_PI1 in component adjustable_parameters (molar).
* CONSTANTS[22] is k_PI2 in component adjustable_parameters (molar).
* CONSTANTS[23] is k_PI3 in component adjustable_parameters (molar).
* CONSTANTS[24] is k_PI4 in component adjustable_parameters (molar).
* CONSTANTS[25] is k_PIH in component adjustable_parameters (molar).
* CONSTANTS[26] is r in component adjustable_parameters (dimensionless).
* CONSTANTS[27] is X_DH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
* CONSTANTS[28] is X_C1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[29] is X_C3 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_half_molar).
* CONSTANTS[30] is X_C4 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
* CONSTANTS[31] is X_F1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[32] is X_ANT in component adjustable_parameters (mole_per_second_per_l_mito).
* CONSTANTS[33] is X_PI1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
* CONSTANTS[34] is X_KH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[35] is X_Hle in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
* CONSTANTS[36] is X_K in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
* CONSTANTS[37] is k_mADP in component adjustable_parameters (molar).
* CONSTANTS[38] is X_AK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[39] is X_CK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[40] is X_MgA in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
* CONSTANTS[41] is k_O2 in component adjustable_parameters (molar).
* CONSTANTS[42] is x_buff in component adjustable_parameters (per_molar).
* CONSTANTS[43] is C_IM in component adjustable_parameters (mole_per_l_mito_per_millivolt).
* CONSTANTS[44] is X_A in component adjustable_parameters (micron_per_second).
* CONSTANTS[45] is X_PI2 in component adjustable_parameters (micron_per_second).
* CONSTANTS[46] is gamma in component adjustable_parameters (per_micron).
* ALGEBRAIC[0] is NAD_x in component dNAD_x_dt (molar).
* STATES[1] is NADH_x in component dNADH_x_dt (molar).
* ALGEBRAIC[1] is QH2_x in component dQH2_x_dt (molar).
* STATES[2] is Q_x in component dQ_x_dt (molar).
* ALGEBRAIC[2] is Cred_x in component dCred_x_dt (molar).
* STATES[3] is Cox_x in component dCox_x_dt (molar).
* ALGEBRAIC[3] is ATP_fx in component dATP_fx_dt (molar).
* STATES[4] is ATP_x in component dATP_x_dt (molar).
* STATES[5] is ATP_mx in component dATP_mx_dt (molar).
* ALGEBRAIC[4] is ADP_fx in component dADP_fx_dt (molar).
* STATES[6] is ADP_x in component dADP_x_dt (molar).
* STATES[7] is ADP_mx in component dADP_mx_dt (molar).
* STATES[8] is pH_x in component dH_x_dt (dimensionless).
* ALGEBRAIC[11] is J_DH in component J_DH (mole_per_second_per_l_mito).
* ALGEBRAIC[38] is J_C1 in component J_C1 (mole_per_second_per_l_mito).
* ALGEBRAIC[39] is J_C3 in component J_C3 (mole_per_second_per_l_mito).
* ALGEBRAIC[40] is J_C4 in component J_C4 (mole_per_second_per_l_mito).
* ALGEBRAIC[41] is J_F1 in component J_F1 (mole_per_second_per_l_mito).
* ALGEBRAIC[17] is J_PI1 in component J_PI1 (mole_per_second_per_l_mito).
* ALGEBRAIC[19] is J_KH in component J_KH (mole_per_second_per_l_mito).
* ALGEBRAIC[18] is J_Hle in component J_Hle (mole_per_second_per_l_mito).
* STATES[9] is K_x in component dK_x_dt (molar).
* ALGEBRAIC[20] is J_K in component J_K (mole_per_second_per_l_mito).
* STATES[10] is Mg_x in component dMg_x_dt (molar).
* ALGEBRAIC[25] is J_MgATP_x in component J_MgATP_x (mole_per_second_per_l_mito).
* ALGEBRAIC[26] is J_MgADP_x in component J_MgADP_x (mole_per_second_per_l_mito).
* ALGEBRAIC[14] is J_ANT in component J_ANT (mole_per_second_per_l_mito).
* STATES[11] is PI_x in component dPI_x_dt (molar).
* ALGEBRAIC[6] is ATP_fi in component dATP_fi_dt (molar).
* STATES[12] is ATP_i in component dATP_i_dt (molar).
* STATES[13] is ATP_mi in component dATP_mi_dt (molar).
* ALGEBRAIC[7] is ADP_fi in component dADP_fi_dt (molar).
* STATES[14] is ADP_i in component dADP_i_dt (molar).
* STATES[15] is ADP_mi in component dADP_mi_dt (molar).
* ALGEBRAIC[23] is J_ATP in component J_ATP (mole_per_second_per_l_mito).
* ALGEBRAIC[31] is J_AKi in component J_AKi (mole_per_second_per_l_mito).
* ALGEBRAIC[22] is J_ADP in component J_ADP (mole_per_second_per_l_mito).
* STATES[16] is AMP_i in component dAMP_i_dt (molar).
* ALGEBRAIC[21] is J_AMP in component J_AMP (mole_per_second_per_l_mito).
* ALGEBRAIC[27] is J_MgATP_i in component J_MgATP_i (mole_per_second_per_l_mito).
* ALGEBRAIC[28] is J_MgADP_i in component J_MgADP_i (mole_per_second_per_l_mito).
* STATES[17] is PI_i in component dPI_i_dt (molar).
* ALGEBRAIC[24] is J_PI2 in component J_PI2 (mole_per_second_per_l_mito).
* STATES[18] is Mg_i in component dMg_i_dt (molar).
* STATES[19] is ATP_c in component dATP_c_dt (molar).
* ALGEBRAIC[32] is J_AKc in component J_AKc (mole_per_second_per_l_cyto).
* CONSTANTS[63] is J_AtC in component J_AtC (mole_per_second_per_l_cell).
* ALGEBRAIC[33] is J_CKc in component J_CKc (mole_per_second_per_l_cyto).
* STATES[20] is ADP_c in component dADP_c_dt (molar).
* STATES[21] is AMP_c in component dAMP_c_dt (molar).
* STATES[22] is ATP_mc in component dATP_mc_dt (molar).
* ALGEBRAIC[29] is J_MgATP_c in component J_MgATP_c (mole_per_second_per_l_cyto).
* STATES[23] is ADP_mc in component dADP_mc_dt (molar).
* ALGEBRAIC[30] is J_MgADP_c in component J_MgADP_c (mole_per_second_per_l_cyto).
* STATES[24] is PI_c in component dPI_c_dt (molar).
* STATES[25] is Mg_c in component dMg_c_dt (molar).
* STATES[26] is PCr_c in component dPCr_c_dt (molar).
* ALGEBRAIC[8] is ATP_fc in component dATP_fc_dt (molar).
* ALGEBRAIC[9] is ADP_fc in component dADP_fc_dt (molar).
* ALGEBRAIC[10] is Cr_c in component dCr_c_dt (molar).
* CONSTANTS[47] is mincon in component J_ANT (molar).
* ALGEBRAIC[12] is Psi_x in component J_ANT (millivolt).
* ALGEBRAIC[13] is Psi_i in component J_ANT (millivolt).
* ALGEBRAIC[15] is H2PIi in component J_PI1 (molar).
* ALGEBRAIC[16] is H2PIx in component J_PI1 (molar).
* CONSTANTS[48] is mincond in component J_Hle (millivolt).
* CONSTANTS[49] is mincond in component J_K (millivolt).
* RATES[8] is d/dt pH_x in component dH_x_dt (dimensionless).
* RATES[9] is d/dt K_x in component dK_x_dt (molar).
* RATES[10] is d/dt Mg_x in component dMg_x_dt (molar).
* RATES[1] is d/dt NADH_x in component dNADH_x_dt (molar).
* RATES[2] is d/dt Q_x in component dQ_x_dt (molar).
* RATES[3] is d/dt Cox_x in component dCox_x_dt (molar).
* RATES[4] is d/dt ATP_x in component dATP_x_dt (molar).
* RATES[6] is d/dt ADP_x in component dADP_x_dt (molar).
* RATES[5] is d/dt ATP_mx in component dATP_mx_dt (molar).
* RATES[7] is d/dt ADP_mx in component dADP_mx_dt (molar).
* RATES[11] is d/dt PI_x in component dPI_x_dt (molar).
* RATES[12] is d/dt ATP_i in component dATP_i_dt (molar).
* RATES[14] is d/dt ADP_i in component dADP_i_dt (molar).
* RATES[16] is d/dt AMP_i in component dAMP_i_dt (molar).
* RATES[13] is d/dt ATP_mi in component dATP_mi_dt (molar).
* RATES[15] is d/dt ADP_mi in component dADP_mi_dt (molar).
* RATES[17] is d/dt PI_i in component dPI_i_dt (molar).
* RATES[18] is d/dt Mg_i in component dMg_i_dt (molar).
* RATES[19] is d/dt ATP_c in component dATP_c_dt (molar).
* RATES[20] is d/dt ADP_c in component dADP_c_dt (molar).
* RATES[21] is d/dt AMP_c in component dAMP_c_dt (molar).
* RATES[22] is d/dt ATP_mc in component dATP_mc_dt (molar).
* RATES[23] is d/dt ADP_mc in component dADP_mc_dt (molar).
* RATES[24] is d/dt PI_c in component dPI_c_dt (molar).
* RATES[25] is d/dt Mg_c in component dMg_c_dt (molar).
* RATES[26] is d/dt PCr_c in component dPCr_c_dt (molar).
* RATES[0] is d/dt deltaPsi in component ddeltaPsi_dt (millivolt).
* There are a total of 5 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 3.48e-5;
CONSTANTS[1] = 0e-4;
CONSTANTS[2] = 2.5775;
CONSTANTS[3] = 0.096484;
CONSTANTS[4] = -69.37;
CONSTANTS[5] = -32.53;
CONSTANTS[6] = -122.94;
CONSTANTS[7] = 36.03;
CONSTANTS[8] = 24e-6;
CONSTANTS[9] = 347e-6;
CONSTANTS[10] = 0.4331;
CONSTANTS[11] = 1.66e9;
CONSTANTS[12] = 0.894;
CONSTANTS[13] = 0.056;
CONSTANTS[14] = 3;
CONSTANTS[15] = 2.7e-3;
CONSTANTS[16] = 42.7e-3;
CONSTANTS[17] = 1.35e-3;
CONSTANTS[18] = 2.97e-3;
CONSTANTS[19] = 7;
CONSTANTS[20] = 0.15;
STATES[0] = 166.67;
CONSTANTS[21] = 0.1553e-3;
CONSTANTS[22] = 0.8222e-3;
CONSTANTS[23] = 0.3601e-3;
CONSTANTS[24] = 5.651e-3;
CONSTANTS[25] = 2.542e-4;
CONSTANTS[26] = 4.559;
CONSTANTS[27] = 0.0866;
CONSTANTS[28] = 4.405e3;
CONSTANTS[29] = 4.887;
CONSTANTS[30] = 6.766e-5;
CONSTANTS[31] = 1e3;
CONSTANTS[32] = 8.123e-3;
CONSTANTS[33] = 3.85e5;
CONSTANTS[34] = 5.651e7;
CONSTANTS[35] = 200;
CONSTANTS[36] = 0;
CONSTANTS[37] = 3.5e-6;
CONSTANTS[38] = 1e7;
CONSTANTS[39] = 1e7;
CONSTANTS[40] = 1e7;
CONSTANTS[41] = 1.2e-4;
CONSTANTS[42] = 100;
CONSTANTS[43] = 6.75e-6;
CONSTANTS[44] = 85;
CONSTANTS[45] = 327;
CONSTANTS[46] = 5.99;
STATES[1] = 0.0015723;
STATES[2] = 6.75e-4;
STATES[3] = 1.35e-3;
STATES[4] = 0.0026657;
STATES[5] = 0.0026046;
STATES[6] = 0.0073343;
STATES[7] = 0.0054765;
STATES[8] = 7;
STATES[9] = 0.14661;
STATES[10] = 1.0229e-3;
STATES[11] = 1.72e-4;
STATES[12] = 0.0065339;
STATES[13] = 0.0063812;
STATES[14] = 6.5773e-5;
STATES[15] = 4.8866e-5;
STATES[16] = 2.8837e-7;
STATES[17] = 1.72e-4;
STATES[18] = 0.0010029;
STATES[19] = 0.0081312;
STATES[20] = 6.85e-5;
STATES[21] = 3.0911e-7;
STATES[22] = 0.0079786;
STATES[23] = 5.1958e-5;
STATES[24] = 8.7702e-3;
STATES[25] = 0.001003;
STATES[26] = 23.589e-3;
CONSTANTS[47] = 1e-9;
CONSTANTS[48] = 1e-6;
CONSTANTS[49] = 1e-6;
CONSTANTS[50] = 1.00000*pow(10.0000, - 6.75000);
CONSTANTS[51] = 1.00000*pow(10.0000, - 6.48000);
CONSTANTS[52] = 1.00000*pow(10.0000, - 6.29000);
CONSTANTS[53] = 1.00000*pow(10.0000, - CONSTANTS[19]);
CONSTANTS[54] = CONSTANTS[13];
CONSTANTS[55] = CONSTANTS[20];
CONSTANTS[56] = 0.664000*1.09000;
CONSTANTS[57] = 0.900000*CONSTANTS[56];
CONSTANTS[58] = 0.100000*CONSTANTS[56];
CONSTANTS[59] = 0.807000*1.04400;
CONSTANTS[60] = CONSTANTS[53];
CONSTANTS[61] = CONSTANTS[13]/CONSTANTS[12];
CONSTANTS[62] = CONSTANTS[12];
CONSTANTS[63] = CONSTANTS[1];
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[6] = 0.1001;
RATES[5] = 0.1001;
RATES[7] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[14] = 0.1001;
RATES[16] = 0.1001;
RATES[13] = 0.1001;
RATES[15] = 0.1001;
RATES[17] = 0.1001;
RATES[18] = 0.1001;
RATES[19] = 0.1001;
RATES[20] = 0.1001;
RATES[21] = 0.1001;
RATES[22] = 0.1001;
RATES[23] = 0.1001;
RATES[24] = 0.1001;
RATES[25] = 0.1001;
RATES[26] = 0.1001;
RATES[0] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[8] - ( (- 1.00000/2.30300)*CONSTANTS[42]*((((((ALGEBRAIC[11] - 5.00000*ALGEBRAIC[38]) - 2.00000*ALGEBRAIC[39]) - 4.00000*ALGEBRAIC[40])+ (CONSTANTS[14] - 1.00000)*ALGEBRAIC[41]+ 2.00000*ALGEBRAIC[17]) - ALGEBRAIC[19])+ALGEBRAIC[18]))/CONSTANTS[57];
resid[1] = RATES[9] - (ALGEBRAIC[19]+ALGEBRAIC[20])/CONSTANTS[57];
resid[2] = RATES[10] - (- ALGEBRAIC[26] - ALGEBRAIC[25])/CONSTANTS[57];
resid[3] = RATES[1] - (ALGEBRAIC[11] - ALGEBRAIC[38])/CONSTANTS[57];
resid[4] = RATES[2] - - (ALGEBRAIC[38] - ALGEBRAIC[39])/CONSTANTS[57];
resid[5] = RATES[3] - - ( 2.00000*ALGEBRAIC[39] - 2.00000*ALGEBRAIC[40])/CONSTANTS[57];
resid[6] = RATES[4] - (ALGEBRAIC[41] - ALGEBRAIC[14])/CONSTANTS[57];
resid[7] = RATES[6] - (- ALGEBRAIC[41]+ALGEBRAIC[14])/CONSTANTS[57];
resid[8] = RATES[5] - ALGEBRAIC[25]/CONSTANTS[57];
resid[9] = RATES[7] - ALGEBRAIC[26]/CONSTANTS[57];
resid[10] = RATES[11] - (- ALGEBRAIC[41]+ALGEBRAIC[17])/CONSTANTS[57];
resid[11] = RATES[12] - (ALGEBRAIC[23]+ALGEBRAIC[14]+ALGEBRAIC[31])/CONSTANTS[58];
resid[12] = RATES[14] - ((ALGEBRAIC[22] - ALGEBRAIC[14]) - 2.00000*ALGEBRAIC[31])/CONSTANTS[58];
resid[13] = RATES[16] - (ALGEBRAIC[21]+ALGEBRAIC[31])/CONSTANTS[58];
resid[14] = RATES[13] - ALGEBRAIC[27]/CONSTANTS[58];
resid[15] = RATES[15] - ALGEBRAIC[28]/CONSTANTS[58];
resid[16] = RATES[17] - (- ALGEBRAIC[17]+ALGEBRAIC[24])/CONSTANTS[58];
resid[17] = RATES[18] - (- ALGEBRAIC[28] - ALGEBRAIC[27])/CONSTANTS[58];
resid[18] = RATES[19] - ((( - CONSTANTS[61]*ALGEBRAIC[23]+ALGEBRAIC[32]) - CONSTANTS[63]/CONSTANTS[62])+ALGEBRAIC[33])/CONSTANTS[59];
resid[19] = RATES[20] - ((( - CONSTANTS[61]*ALGEBRAIC[22] - 2.00000*ALGEBRAIC[32])+CONSTANTS[63]/CONSTANTS[62]) - ALGEBRAIC[33])/CONSTANTS[59];
resid[20] = RATES[21] - ( - CONSTANTS[61]*ALGEBRAIC[21]+ALGEBRAIC[32])/CONSTANTS[59];
resid[21] = RATES[22] - ALGEBRAIC[29]/CONSTANTS[59];
resid[22] = RATES[23] - ALGEBRAIC[30]/CONSTANTS[59];
resid[23] = RATES[24] - ( - ALGEBRAIC[24]*CONSTANTS[61]+CONSTANTS[63]/CONSTANTS[62])/CONSTANTS[59];
resid[24] = RATES[25] - (- ALGEBRAIC[29] - ALGEBRAIC[30])/CONSTANTS[59];
resid[25] = RATES[26] - - ALGEBRAIC[33]/CONSTANTS[59];
resid[26] = RATES[0] - ((((( 4.00000*ALGEBRAIC[38]+ 2.00000*ALGEBRAIC[39]+ 4.00000*ALGEBRAIC[40]) - CONSTANTS[14]*ALGEBRAIC[41]) - ALGEBRAIC[14]) - ALGEBRAIC[18]) - ALGEBRAIC[20])/CONSTANTS[43];
}
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] = CONSTANTS[18] - STATES[1];
ALGEBRAIC[11] = ( CONSTANTS[27]*( CONSTANTS[26]*ALGEBRAIC[0] - STATES[1])*(1.00000+STATES[11]/CONSTANTS[21]))/(1.00000+STATES[11]/CONSTANTS[22]);
ALGEBRAIC[3] = STATES[4] - STATES[5];
ALGEBRAIC[4] = STATES[6] - STATES[7];
ALGEBRAIC[6] = STATES[12] - STATES[13];
ALGEBRAIC[7] = STATES[14] - STATES[15];
ALGEBRAIC[12] = - 0.650000*STATES[0];
ALGEBRAIC[13] = 0.350000*STATES[0];
ALGEBRAIC[14] = (CONDVAR[1]>0.00000||CONDVAR[2]>0.00000 ? CONSTANTS[32]*(ALGEBRAIC[7]/(ALGEBRAIC[7]+ ALGEBRAIC[6]*exp(( - CONSTANTS[3]*ALGEBRAIC[13])/CONSTANTS[2])) - ALGEBRAIC[4]/(ALGEBRAIC[4]+ ALGEBRAIC[3]*exp(( - CONSTANTS[3]*ALGEBRAIC[12])/CONSTANTS[2])))*(ALGEBRAIC[7]/(ALGEBRAIC[7]+CONSTANTS[37])) : 0.00000);
ALGEBRAIC[5] = 1.00000*pow(10.0000, - STATES[8]);
ALGEBRAIC[15] = ( STATES[17]*CONSTANTS[60])/(CONSTANTS[60]+CONSTANTS[50]);
ALGEBRAIC[16] = ( STATES[11]*ALGEBRAIC[5])/(ALGEBRAIC[5]+CONSTANTS[50]);
ALGEBRAIC[17] = ( CONSTANTS[33]*( CONSTANTS[60]*ALGEBRAIC[15] - ALGEBRAIC[5]*ALGEBRAIC[16]))/(ALGEBRAIC[15]+CONSTANTS[25]);
ALGEBRAIC[18] = (CONDVAR[3]>0.00000 ? ( CONSTANTS[35]*STATES[0]*( CONSTANTS[60]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - ALGEBRAIC[5]))/(exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - 1.00000) : ( CONSTANTS[35]*CONSTANTS[2]*(CONSTANTS[60] - ALGEBRAIC[5]))/CONSTANTS[3]);
ALGEBRAIC[19] = CONSTANTS[34]*( CONSTANTS[55]*ALGEBRAIC[5] - STATES[9]*CONSTANTS[60]);
ALGEBRAIC[20] = (CONDVAR[4]>0.00000 ? ( CONSTANTS[36]*STATES[0]*( CONSTANTS[55]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - STATES[9]))/(exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - 1.00000) : ( CONSTANTS[36]*CONSTANTS[2]*(CONSTANTS[55] - STATES[9]))/CONSTANTS[3]);
ALGEBRAIC[21] = CONSTANTS[46]*CONSTANTS[44]*(STATES[21] - STATES[16]);
ALGEBRAIC[22] = CONSTANTS[46]*CONSTANTS[44]*(STATES[20] - STATES[14]);
ALGEBRAIC[23] = CONSTANTS[46]*CONSTANTS[44]*(STATES[19] - STATES[12]);
ALGEBRAIC[24] = CONSTANTS[46]*CONSTANTS[45]*(STATES[24] - STATES[17]);
ALGEBRAIC[25] = CONSTANTS[40]*( ALGEBRAIC[3]*STATES[10] - CONSTANTS[8]*STATES[5]);
ALGEBRAIC[26] = CONSTANTS[40]*( ALGEBRAIC[4]*STATES[10] - CONSTANTS[9]*STATES[7]);
ALGEBRAIC[27] = CONSTANTS[40]*( ALGEBRAIC[6]*STATES[18] - CONSTANTS[8]*STATES[13]);
ALGEBRAIC[28] = CONSTANTS[40]*( ALGEBRAIC[7]*STATES[18] - CONSTANTS[9]*STATES[15]);
ALGEBRAIC[8] = STATES[19] - STATES[22];
ALGEBRAIC[29] = CONSTANTS[40]*1.00000*( ALGEBRAIC[8]*STATES[25] - CONSTANTS[8]*STATES[22]);
ALGEBRAIC[9] = STATES[20] - STATES[23];
ALGEBRAIC[30] = CONSTANTS[40]*1.00000*( ALGEBRAIC[9]*STATES[25] - CONSTANTS[9]*STATES[23]);
ALGEBRAIC[31] = CONSTANTS[38]*( CONSTANTS[10]*STATES[14]*STATES[14] - STATES[16]*STATES[12]);
ALGEBRAIC[32] = CONSTANTS[38]*1.00000*( CONSTANTS[10]*STATES[20]*STATES[20] - STATES[21]*STATES[19]);
ALGEBRAIC[10] = CONSTANTS[16] - STATES[26];
ALGEBRAIC[33] = CONSTANTS[39]*1.00000*( CONSTANTS[11]*STATES[20]*STATES[26]*CONSTANTS[53] - STATES[19]*ALGEBRAIC[10]);
ALGEBRAIC[35] = CONSTANTS[4] - 1.00000*CONSTANTS[2]*log(ALGEBRAIC[5]/1.00000e-07);
ALGEBRAIC[34] = CONSTANTS[3]*STATES[0]+ CONSTANTS[2]*log(CONSTANTS[60]/ALGEBRAIC[5]);
ALGEBRAIC[1] = CONSTANTS[17] - STATES[2];
ALGEBRAIC[38] = CONSTANTS[28]*( exp(- (ALGEBRAIC[35]+ 4.00000*ALGEBRAIC[34])/CONSTANTS[2])*STATES[1]*STATES[2] - ALGEBRAIC[0]*ALGEBRAIC[1]);
ALGEBRAIC[36] = CONSTANTS[5]+ 2.00000*CONSTANTS[2]*log(ALGEBRAIC[5]/1.00000e-07);
ALGEBRAIC[2] = CONSTANTS[15] - STATES[3];
ALGEBRAIC[39] = (( CONSTANTS[29]*(1.00000+STATES[11]/CONSTANTS[23]))/(1.00000+STATES[11]/CONSTANTS[24]))*( exp(- ((ALGEBRAIC[36]+ 4.00000*ALGEBRAIC[34]) - 2.00000*CONSTANTS[3]*STATES[0])/( 2.00000*CONSTANTS[2]))*STATES[3]* pow(ALGEBRAIC[1], 1.0 / 2) - ALGEBRAIC[2]* pow(STATES[2], 1.0 / 2));
ALGEBRAIC[37] = CONSTANTS[6] - 2.00000*CONSTANTS[2]*log(ALGEBRAIC[5]/1.00000e-07);
ALGEBRAIC[40] = (CONDVAR[0]<0.00000 ? 0.00000 : (( CONSTANTS[30]*CONSTANTS[0])/(CONSTANTS[0]+CONSTANTS[41]))*(ALGEBRAIC[2]/CONSTANTS[15])*( exp(- (ALGEBRAIC[37]+ 2.00000*ALGEBRAIC[34])/( 2.00000*CONSTANTS[2]))*ALGEBRAIC[2]*pow(CONSTANTS[0]/1.00000, 0.250000) - STATES[3]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2])));
ALGEBRAIC[41] = CONSTANTS[31]*( (( exp(- (CONSTANTS[7] - CONSTANTS[14]*ALGEBRAIC[34])/CONSTANTS[2])*CONSTANTS[9])/CONSTANTS[8])*STATES[7]*STATES[11] - STATES[5]*1.00000);
}
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;
SI[18] = 1.0;
SI[19] = 1.0;
SI[20] = 1.0;
SI[21] = 1.0;
SI[22] = 1.0;
SI[23] = 1.0;
SI[24] = 1.0;
SI[25] = 1.0;
SI[26] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = CONSTANTS[0] - 1.00000e-12;
CONDVAR[1] = ALGEBRAIC[7] - CONSTANTS[47];
CONDVAR[2] = ALGEBRAIC[6] - CONSTANTS[47];
CONDVAR[3] = fabs(STATES[0]) - CONSTANTS[48];
CONDVAR[4] = fabs(STATES[0]) - CONSTANTS[49];
}