/*
   There are a total of 15 entries in the algebraic variable array.
   There are a total of 10 entries in each of the rate and state variable arrays.
   There are a total of 38 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * STATES[0] is ATP_cyt in component ATP_cyt (micromolar).
 * ALGEBRAIC[5] is J_hyd in component J_hyd (flux).
 * ALGEBRAIC[12] is J_CKMM in component J_CKMM (flux).
 * ALGEBRAIC[7] is J_diff_ATP in component J_diff_ATP (flux).
 * CONSTANTS[0] is V_cyt in component fractional_volumes (dimensionless).
 * STATES[1] is ADP_cyt in component ADP_cyt (micromolar).
 * ALGEBRAIC[8] is J_diff_ADP in component J_diff_ADP (flux).
 * STATES[2] is PCr_cyt in component PCr_cyt (micromolar).
 * ALGEBRAIC[9] is J_diff_PCr in component J_diff_PCr (flux).
 * STATES[3] is Cr_cyt in component Cr_cyt (micromolar).
 * ALGEBRAIC[10] is J_diff_Cr in component J_diff_Cr (flux).
 * STATES[4] is Pi_cyt in component Pi_cyt (micromolar).
 * ALGEBRAIC[11] is J_diff_Pi in component J_diff_Pi (flux).
 * STATES[5] is ATP_ims in component ATP_ims (micromolar).
 * ALGEBRAIC[13] is J_CKMi in component J_CKMi (flux).
 * ALGEBRAIC[14] is J_syn in component J_syn (flux).
 * CONSTANTS[1] is V_ims in component fractional_volumes (dimensionless).
 * STATES[6] is ADP_ims in component ADP_ims (micromolar).
 * STATES[7] is PCr_ims in component PCr_ims (micromolar).
 * STATES[8] is Cr_ims in component Cr_ims (micromolar).
 * STATES[9] is Pi_ims in component Pi_ims (micromolar).
 * ALGEBRAIC[0] is Den_MMCK in component J_CKMM (dimensionless).
 * CONSTANTS[2] is Kia in component J_CKMM (micromolar).
 * CONSTANTS[3] is Kb in component J_CKMM (micromolar).
 * CONSTANTS[4] is Kib in component J_CKMM (micromolar).
 * CONSTANTS[31] is KIb in component J_CKMM (micromolar).
 * CONSTANTS[32] is Kc in component J_CKMM (micromolar).
 * CONSTANTS[5] is Kic in component J_CKMM (micromolar).
 * CONSTANTS[6] is Kd in component J_CKMM (micromolar).
 * CONSTANTS[7] is Kid in component J_CKMM (micromolar).
 * CONSTANTS[8] is Vmax_MM_f in component J_CKMM (flux).
 * CONSTANTS[9] is Vmax_MM_b in component J_CKMM (flux).
 * ALGEBRAIC[1] is Den_MiCK in component J_CKMi (dimensionless).
 * CONSTANTS[10] is Kia in component J_CKMi (micromolar).
 * CONSTANTS[11] is Kb in component J_CKMi (micromolar).
 * CONSTANTS[12] is Kib in component J_CKMi (micromolar).
 * CONSTANTS[33] is KIb in component J_CKMi (micromolar).
 * CONSTANTS[34] is Kc in component J_CKMi (micromolar).
 * CONSTANTS[13] is Kic in component J_CKMi (micromolar).
 * CONSTANTS[14] is Kd in component J_CKMi (micromolar).
 * CONSTANTS[15] is Kid in component J_CKMi (micromolar).
 * CONSTANTS[16] is Vmax_Mi_f in component J_CKMi (flux).
 * CONSTANTS[17] is Vmax_Mi_b in component J_CKMi (flux).
 * ALGEBRAIC[4] is H_ATPmax in component J_hyd (flux).
 * CONSTANTS[18] is J_hyd_basis_1 in component J_hyd (flux).
 * CONSTANTS[19] is J_hyd_basis_2 in component J_hyd (flux).
 * CONSTANTS[20] is freq_1 in component J_hyd (dimensionless).
 * CONSTANTS[21] is freq_2 in component J_hyd (dimensionless).
 * CONSTANTS[35] is t_cycle_1 in component J_hyd (second).
 * CONSTANTS[36] is t_cycle_2 in component J_hyd (second).
 * ALGEBRAIC[3] is t_cycle in component J_hyd (second).
 * CONSTANTS[22] is nb_of_cycles_1 in component J_hyd (dimensionless).
 * CONSTANTS[37] is duration_1 in component J_hyd (second).
 * ALGEBRAIC[2] is ltime in component J_hyd (second).
 * ALGEBRAIC[6] is Den_syn in component J_syn (dimensionless).
 * CONSTANTS[23] is KPi in component J_syn (micromolar).
 * CONSTANTS[24] is KADP in component J_syn (micromolar).
 * CONSTANTS[25] is V_max_syn in component J_syn (flux).
 * CONSTANTS[26] is PS_tot_ATP in component J_diff_ATP (first_order_rate_constant).
 * CONSTANTS[27] is PS_tot_ADP in component J_diff_ADP (first_order_rate_constant).
 * CONSTANTS[28] is PS_tot_PCr in component J_diff_PCr (first_order_rate_constant).
 * CONSTANTS[29] is PS_tot_Cr in component J_diff_Cr (first_order_rate_constant).
 * CONSTANTS[30] is PS_tot_Pi in component J_diff_Pi (first_order_rate_constant).
 * RATES[0] is d/dt ATP_cyt in component ATP_cyt (micromolar).
 * RATES[1] is d/dt ADP_cyt in component ADP_cyt (micromolar).
 * RATES[2] is d/dt PCr_cyt in component PCr_cyt (micromolar).
 * RATES[3] is d/dt Cr_cyt in component Cr_cyt (micromolar).
 * RATES[4] is d/dt Pi_cyt in component Pi_cyt (micromolar).
 * RATES[5] is d/dt ATP_ims in component ATP_ims (micromolar).
 * RATES[6] is d/dt ADP_ims in component ADP_ims (micromolar).
 * RATES[7] is d/dt PCr_ims in component PCr_ims (micromolar).
 * RATES[8] is d/dt Cr_ims in component Cr_ims (micromolar).
 * RATES[9] is d/dt Pi_ims in component Pi_ims (micromolar).
 * There are a total of 9 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 5912.77;
CONSTANTS[0] = 0.75;
STATES[1] = 64;
STATES[2] = 5000;
STATES[3] = 10500;
STATES[4] = 913;
STATES[5] = 5912.77;
CONSTANTS[1] = 0.0625;
STATES[6] = 39;
STATES[7] = 5000;
STATES[8] = 10500;
STATES[9] = 910;
CONSTANTS[2] = 9.0E2;
CONSTANTS[3] = 1.55E4;
CONSTANTS[4] = 3.49E4;
CONSTANTS[5] = 2.224E2;
CONSTANTS[6] = 1.67E3;
CONSTANTS[7] = 4.73E3;
CONSTANTS[8] = 1.144E4;
CONSTANTS[9] = 4.804E4;
CONSTANTS[10] = 7.5E2;
CONSTANTS[11] = 5.2E3;
CONSTANTS[12] = 2.88E4;
CONSTANTS[13] = 2.048E2;
CONSTANTS[14] = 5.0E2;
CONSTANTS[15] = 1.6E3;
CONSTANTS[16] = 8.82E2;
CONSTANTS[17] = 3.704E3;
CONSTANTS[18] = 4.865e2;
CONSTANTS[19] = 6.276e2;
CONSTANTS[20] = 135;
CONSTANTS[21] = 220;
CONSTANTS[22] = 5;
CONSTANTS[23] = 8E2;
CONSTANTS[24] = 25;
CONSTANTS[25] = 1.504E4;
CONSTANTS[26] = 13.3;
CONSTANTS[27] = 13.3;
CONSTANTS[28] = 155.0;
CONSTANTS[29] = 155.0;
CONSTANTS[30] = 194.0;
CONSTANTS[31] = CONSTANTS[4];
CONSTANTS[32] = ( CONSTANTS[5]*CONSTANTS[6])/CONSTANTS[7];
CONSTANTS[33] = CONSTANTS[12];
CONSTANTS[34] = ( CONSTANTS[13]*CONSTANTS[14])/CONSTANTS[15];
CONSTANTS[35] = 60.0000/CONSTANTS[20];
CONSTANTS[36] = 60.0000/CONSTANTS[21];
CONSTANTS[37] =  CONSTANTS[22]*CONSTANTS[35];
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 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;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - (ALGEBRAIC[7] - (ALGEBRAIC[5]+ALGEBRAIC[12]))/CONSTANTS[0];
resid[1] = RATES[1] - (ALGEBRAIC[8]+ALGEBRAIC[5]+ALGEBRAIC[12])/CONSTANTS[0];
resid[2] = RATES[2] - (ALGEBRAIC[9]+ALGEBRAIC[12])/CONSTANTS[0];
resid[3] = RATES[3] - (ALGEBRAIC[10] - ALGEBRAIC[12])/CONSTANTS[0];
resid[4] = RATES[4] - (ALGEBRAIC[11]+ALGEBRAIC[5])/CONSTANTS[0];
resid[5] = RATES[5] - (ALGEBRAIC[14] - (ALGEBRAIC[7]+ALGEBRAIC[13]))/CONSTANTS[1];
resid[6] = RATES[6] - (ALGEBRAIC[13] - (ALGEBRAIC[14]+ALGEBRAIC[8]))/CONSTANTS[1];
resid[7] = RATES[7] - (ALGEBRAIC[13] - ALGEBRAIC[9])/CONSTANTS[1];
resid[8] = RATES[8] - - (ALGEBRAIC[13]+ALGEBRAIC[10])/CONSTANTS[1];
resid[9] = RATES[9] - - (ALGEBRAIC[14]+ALGEBRAIC[11])/CONSTANTS[1];
}
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[4] = (CONDVAR[2]<=0.00000 ?  6.00000*CONSTANTS[18] :  6.00000*CONSTANTS[19]);
ALGEBRAIC[3] = (CONDVAR[1]<=0.00000 ? CONSTANTS[35] : CONSTANTS[36]);
ALGEBRAIC[2] = (CONDVAR[0]<=0.00000 ? VOI -  CONSTANTS[35]*floor(VOI/CONSTANTS[35]) : (VOI - CONSTANTS[37]) -  CONSTANTS[36]*floor((VOI - CONSTANTS[37])/CONSTANTS[36]));
ALGEBRAIC[5] = (CONDVAR[3]>=0.00000&&CONDVAR[4]<0.00000 ?  (( ALGEBRAIC[4]*ALGEBRAIC[2])/ALGEBRAIC[3])*6.00000 : CONDVAR[5]>=0.00000&&CONDVAR[6]<0.00000 ?  ALGEBRAIC[4]*(1.00000 -  6.00000*(ALGEBRAIC[2]/ALGEBRAIC[3] - 1.00000/6.00000)) : CONDVAR[7]>=0.00000&&CONDVAR[8]<0.00000 ? 0.00000 : 0.0/0.0);
ALGEBRAIC[7] =  CONSTANTS[26]*(STATES[5] - STATES[0]);
ALGEBRAIC[8] =  CONSTANTS[27]*(STATES[6] - STATES[1]);
ALGEBRAIC[9] =  CONSTANTS[28]*(STATES[7] - STATES[2]);
ALGEBRAIC[10] =  CONSTANTS[29]*(STATES[8] - STATES[3]);
ALGEBRAIC[11] =  CONSTANTS[30]*(STATES[9] - STATES[4]);
ALGEBRAIC[0] = 1.00000+STATES[3]/CONSTANTS[4]+STATES[2]/CONSTANTS[7]+ STATES[0]*(1.00000/CONSTANTS[2]+STATES[3]/( CONSTANTS[2]*CONSTANTS[3]))+ STATES[1]*(1.00000/CONSTANTS[5]+STATES[2]/( CONSTANTS[7]*CONSTANTS[32])+STATES[3]/( CONSTANTS[5]*CONSTANTS[31]));
ALGEBRAIC[12] = ( CONSTANTS[8]*(( STATES[0]*STATES[3])/( CONSTANTS[2]*CONSTANTS[3])) -  CONSTANTS[9]*(( STATES[1]*STATES[2])/( CONSTANTS[5]*CONSTANTS[6])))/ALGEBRAIC[0];
ALGEBRAIC[1] = 1.00000+STATES[8]/CONSTANTS[12]+STATES[7]/CONSTANTS[15]+ STATES[5]*(1.00000/CONSTANTS[10]+STATES[8]/( CONSTANTS[10]*CONSTANTS[11]))+ STATES[6]*(1.00000/CONSTANTS[13]+STATES[7]/( CONSTANTS[15]*CONSTANTS[34])+STATES[8]/( CONSTANTS[13]*CONSTANTS[33]));
ALGEBRAIC[13] = ( CONSTANTS[16]*(( STATES[5]*STATES[8])/( CONSTANTS[10]*CONSTANTS[11])) -  CONSTANTS[17]*(( STATES[6]*STATES[7])/( CONSTANTS[13]*CONSTANTS[14])))/ALGEBRAIC[1];
ALGEBRAIC[6] = 1.00000+STATES[6]/CONSTANTS[24]+STATES[9]/CONSTANTS[23]+( STATES[6]*STATES[9])/( CONSTANTS[24]*CONSTANTS[23]);
ALGEBRAIC[14] =  CONSTANTS[25]*(( STATES[6]*STATES[9])/( CONSTANTS[23]*CONSTANTS[24]*ALGEBRAIC[6]));
}
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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - CONSTANTS[37];
CONDVAR[1] = VOI - CONSTANTS[37];
CONDVAR[2] = VOI - CONSTANTS[37];
CONDVAR[3] = ALGEBRAIC[2] - 0.00000;
CONDVAR[4] = ALGEBRAIC[2] -  (1.00000/6.00000)*ALGEBRAIC[3];
CONDVAR[5] = ALGEBRAIC[2] -  (1.00000/6.00000)*ALGEBRAIC[3];
CONDVAR[6] = ALGEBRAIC[2] -  (1.00000/3.00000)*ALGEBRAIC[3];
CONDVAR[7] = ALGEBRAIC[2] -  (1.00000/3.00000)*ALGEBRAIC[3];
CONDVAR[8] = ALGEBRAIC[2] - ALGEBRAIC[3];
}