/*
   There are a total of 24 entries in the algebraic variable array.
   There are a total of 22 entries in each of the rate and state variable arrays.
   There are a total of 60 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * CONSTANTS[0] is Glc_x_0 in component Glc_x_0 (millimolar).
 * STATES[0] is Glc_x in component Glc_x (millimolar).
 * ALGEBRAIC[1] is GlcTrans in component GlcTrans (flux).
 * ALGEBRAIC[0] is inGlc in component inGlc (flux).
 * STATES[1] is Glc in component Glc (millimolar).
 * ALGEBRAIC[2] is HK in component HK (flux).
 * STATES[2] is G6P in component G6P (millimolar).
 * ALGEBRAIC[3] is PGI in component PGI (flux).
 * ALGEBRAIC[21] is storage in component storage (flux).
 * STATES[3] is F6P in component F6P (millimolar).
 * ALGEBRAIC[4] is PFK in component PFK (flux).
 * STATES[4] is FBP in component FBP (millimolar).
 * ALGEBRAIC[5] is ALD in component ALD (flux).
 * STATES[5] is GAP in component GAP (millimolar).
 * ALGEBRAIC[7] is GAPDH in component GAPDH (flux).
 * ALGEBRAIC[6] is TIM in component TIM (flux).
 * STATES[6] is DHAP in component DHAP (millimolar).
 * ALGEBRAIC[14] is lpGlyc in component lpGlyc (flux).
 * STATES[7] is BPG in component BPG (millimolar).
 * ALGEBRAIC[8] is lpPEP in component lpPEP (flux).
 * STATES[8] is PEP in component PEP (millimolar).
 * ALGEBRAIC[9] is PK in component PK (flux).
 * STATES[9] is Pyr in component Pyr (millimolar).
 * ALGEBRAIC[10] is PDC in component PDC (flux).
 * STATES[10] is ATP in component ATP (millimolar).
 * ALGEBRAIC[22] is consum in component consum (flux).
 * ALGEBRAIC[23] is AK in component AK (flux).
 * STATES[11] is ADP in component ADP (millimolar).
 * STATES[12] is AMP in component AMP (millimolar).
 * CONSTANTS[1] is CN_x_0 in component CN_x_0 (millimolar).
 * STATES[13] is CN_x in component CN_x (millimolar).
 * ALGEBRAIC[19] is lacto in component lacto (flux).
 * ALGEBRAIC[20] is inCN in component inCN (flux).
 * STATES[14] is ACA in component ACA (millimolar).
 * ALGEBRAIC[18] is difACA in component difACA (flux).
 * ALGEBRAIC[11] is ADH in component ADH (flux).
 * STATES[15] is ACA_x in component ACA_x (millimolar).
 * ALGEBRAIC[17] is outACA in component outACA (flux).
 * STATES[16] is EtOH in component EtOH (millimolar).
 * ALGEBRAIC[13] is difEtOH in component difEtOH (flux).
 * STATES[17] is EtOH_x in component EtOH_x (millimolar).
 * ALGEBRAIC[12] is outEtOH in component outEtOH (flux).
 * STATES[18] is Glyc in component Glyc (millimolar).
 * ALGEBRAIC[16] is difGlyc in component difGlyc (flux).
 * STATES[19] is Glyc_x in component Glyc_x (millimolar).
 * ALGEBRAIC[15] is outGlyc in component outGlyc (flux).
 * STATES[20] is NADH in component NADH (millimolar).
 * STATES[21] is NAD in component NAD (millimolar).
 * CONSTANTS[2] is k0 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[3] is K2_Glc in component GlcTrans (millimolar).
 * CONSTANTS[4] is K2_IG6P in component GlcTrans (millimolar).
 * CONSTANTS[5] is K2_IIG6P in component GlcTrans (millimolar).
 * CONSTANTS[6] is V_2m in component GlcTrans (flux).
 * CONSTANTS[7] is P2 in component GlcTrans (dimensionless).
 * CONSTANTS[8] is y_vol in component model_parameters (dimensionless).
 * CONSTANTS[9] is K3_DGlc in component HK (millimolar).
 * CONSTANTS[10] is K3_Glc in component HK (millimolar).
 * CONSTANTS[11] is K3_ATP in component HK (millimolar).
 * CONSTANTS[12] is V_3m in component HK (flux).
 * CONSTANTS[13] is V_4m in component PGI (flux).
 * CONSTANTS[14] is K4_G6P in component PGI (millimolar).
 * CONSTANTS[15] is K4_F6P in component PGI (millimolar).
 * CONSTANTS[16] is K4_eq in component PGI (dimensionless).
 * CONSTANTS[17] is K5 in component PFK (millimolar2).
 * CONSTANTS[18] is kappa5 in component PFK (dimensionless).
 * CONSTANTS[19] is V_5m in component PFK (flux).
 * CONSTANTS[20] is K6_eq in component ALD (millimolar).
 * CONSTANTS[21] is K6_FBP in component ALD (millimolar).
 * CONSTANTS[22] is K6_DHAP in component ALD (millimolar).
 * CONSTANTS[23] is K6_GAP in component ALD (millimolar).
 * CONSTANTS[24] is K6_IGAP in component ALD (millimolar).
 * CONSTANTS[25] is V_6r in component ALD (flux).
 * CONSTANTS[26] is V_6f in component ALD (flux).
 * CONSTANTS[27] is K7_eq in component TIM (dimensionless).
 * CONSTANTS[28] is K7_DHAP in component TIM (millimolar).
 * CONSTANTS[29] is K7_GAP in component TIM (millimolar).
 * CONSTANTS[30] is V_7m in component TIM (flux).
 * CONSTANTS[31] is K8_NAD in component GAPDH (millimolar).
 * CONSTANTS[32] is K8_NADH in component GAPDH (millimolar).
 * CONSTANTS[33] is K8_GAP in component GAPDH (millimolar).
 * CONSTANTS[34] is K8_BPG in component GAPDH (millimolar).
 * CONSTANTS[35] is K8_eq in component GAPDH (dimensionless).
 * CONSTANTS[36] is V_8m in component GAPDH (flux).
 * CONSTANTS[37] is k9f in component lpPEP (second_order_rate_constant).
 * CONSTANTS[38] is k9r in component lpPEP (second_order_rate_constant).
 * CONSTANTS[39] is K10_PEP in component PK (millimolar).
 * CONSTANTS[40] is K10_ADP in component PK (millimolar).
 * CONSTANTS[41] is V_10m in component PK (flux).
 * CONSTANTS[42] is K11 in component PDC (millimolar).
 * CONSTANTS[43] is V_11m in component PDC (flux).
 * CONSTANTS[44] is K12_NADH in component ADH (millimolar).
 * CONSTANTS[45] is K12_ACA in component ADH (millimolar).
 * CONSTANTS[46] is V_12m in component ADH (flux).
 * CONSTANTS[47] is k13 in component difEtOH (first_order_rate_constant).
 * CONSTANTS[48] is K15_NADH in component lpGlyc (millimolar).
 * CONSTANTS[49] is K15_INADH in component lpGlyc (millimolar).
 * CONSTANTS[50] is K15_INAD in component lpGlyc (millimolar).
 * CONSTANTS[51] is K15_DHAP in component lpGlyc (millimolar).
 * CONSTANTS[52] is V_15m in component lpGlyc (flux).
 * CONSTANTS[53] is k16 in component difGlyc (first_order_rate_constant).
 * CONSTANTS[54] is k18 in component difACA (first_order_rate_constant).
 * CONSTANTS[55] is k20 in component lacto (second_order_rate_constant).
 * CONSTANTS[56] is k22 in component storage (second_order_rate_constant).
 * CONSTANTS[57] is k23 in component consum (first_order_rate_constant).
 * CONSTANTS[58] is k24f in component AK (second_order_rate_constant).
 * CONSTANTS[59] is k24r in component AK (second_order_rate_constant).
 * RATES[0] is d/dt Glc_x in component Glc_x (millimolar).
 * RATES[1] is d/dt Glc in component Glc (millimolar).
 * RATES[2] is d/dt G6P in component G6P (millimolar).
 * RATES[3] is d/dt F6P in component F6P (millimolar).
 * RATES[4] is d/dt FBP in component FBP (millimolar).
 * RATES[5] is d/dt GAP in component GAP (millimolar).
 * RATES[6] is d/dt DHAP in component DHAP (millimolar).
 * RATES[7] is d/dt BPG in component BPG (millimolar).
 * RATES[8] is d/dt PEP in component PEP (millimolar).
 * RATES[9] is d/dt Pyr in component Pyr (millimolar).
 * RATES[10] is d/dt ATP in component ATP (millimolar).
 * RATES[11] is d/dt ADP in component ADP (millimolar).
 * RATES[12] is d/dt AMP in component AMP (millimolar).
 * RATES[13] is d/dt CN_x in component CN_x (millimolar).
 * RATES[14] is d/dt ACA in component ACA (millimolar).
 * RATES[15] is d/dt ACA_x in component ACA_x (millimolar).
 * RATES[16] is d/dt EtOH in component EtOH (millimolar).
 * RATES[17] is d/dt EtOH_x in component EtOH_x (millimolar).
 * RATES[18] is d/dt Glyc in component Glyc (millimolar).
 * RATES[19] is d/dt Glyc_x in component Glyc_x (millimolar).
 * RATES[20] is d/dt NADH in component NADH (millimolar).
 * RATES[21] is d/dt NAD in component NAD (millimolar).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 24.0;
STATES[0] = 6.7;
STATES[1] = 0.573074;
STATES[2] = 4.2;
STATES[3] = 0.49;
STATES[4] = 4.64;
STATES[5] = 0.115;
STATES[6] = 2.95;
STATES[7] = 0.00027;
STATES[8] = 0.04;
STATES[9] = 8.7;
STATES[10] = 2.1;
STATES[11] = 1.5;
STATES[12] = 0.33;
CONSTANTS[1] = 5.60;
STATES[13] = 5.20358;
STATES[14] = 1.48153;
STATES[15] = 1.28836;
STATES[16] = 19.2379;
STATES[17] = 16.4514;
STATES[18] = 4.196;
STATES[19] = 1.68478;
STATES[20] = 0.33;
STATES[21] = 0.65;
CONSTANTS[2] = 0.048;
CONSTANTS[3] = 1.7;
CONSTANTS[4] = 1.2;
CONSTANTS[5] = 7.2;
CONSTANTS[6] = 1014.96;
CONSTANTS[7] = 1.0;
CONSTANTS[8] = 59.0;
CONSTANTS[9] = 0.37;
CONSTANTS[10] = 0.0;
CONSTANTS[11] = 0.1;
CONSTANTS[12] = 51.7547;
CONSTANTS[13] = 496.042;
CONSTANTS[14] = 0.8;
CONSTANTS[15] = 0.15;
CONSTANTS[16] = 0.13;
CONSTANTS[17] = 0.021;
CONSTANTS[18] = 0.15;
CONSTANTS[19] = 45.4327;
CONSTANTS[20] = 0.081;
CONSTANTS[21] = 0.3;
CONSTANTS[22] = 2.0;
CONSTANTS[23] = 4.0;
CONSTANTS[24] = 10.0;
CONSTANTS[25] = 1.10391E4;
CONSTANTS[26] = 2.20782E3;
CONSTANTS[27] = 0.055;
CONSTANTS[28] = 1.23;
CONSTANTS[29] = 1.27;
CONSTANTS[30] = 1.16365E2;
CONSTANTS[31] = 0.1;
CONSTANTS[32] = 0.06;
CONSTANTS[33] = 0.6;
CONSTANTS[34] = 0.01;
CONSTANTS[35] = 0.0055;
CONSTANTS[36] = 8.33858E2;
CONSTANTS[37] = 4.43866E5;
CONSTANTS[38] = 1.52862E3;
CONSTANTS[39] = 0.2;
CONSTANTS[40] = 0.17;
CONSTANTS[41] = 3.43096E2;
CONSTANTS[42] = 0.3;
CONSTANTS[43] = 5.31328E1;
CONSTANTS[44] = 0.1;
CONSTANTS[45] = 0.71;
CONSTANTS[46] = 8.98023E1;
CONSTANTS[47] = 16.72;
CONSTANTS[48] = 0.13;
CONSTANTS[49] = 0.034;
CONSTANTS[50] = 0.13;
CONSTANTS[51] = 25.0;
CONSTANTS[52] = 8.14797E1;
CONSTANTS[53] = 1.9;
CONSTANTS[54] = 24.7;
CONSTANTS[55] = 2.83828E-3;
CONSTANTS[56] = 2.25932;
CONSTANTS[57] = 3.20760;
CONSTANTS[58] = 4.32900E2;
CONSTANTS[59] = 1.33333E2;
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;
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[16] = 0.1001;
RATES[17] = 0.1001;
RATES[18] = 0.1001;
RATES[19] = 0.1001;
RATES[20] = 0.1001;
RATES[21] = 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[0] - ALGEBRAIC[1];
resid[1] = RATES[1] -  59.0000*ALGEBRAIC[1] - ALGEBRAIC[2];
resid[2] = RATES[2] - ALGEBRAIC[2] - (ALGEBRAIC[3]+ALGEBRAIC[21]);
resid[3] = RATES[3] - ALGEBRAIC[3] - ALGEBRAIC[4];
resid[4] = RATES[4] - ALGEBRAIC[4] - ALGEBRAIC[5];
resid[5] = RATES[5] - (ALGEBRAIC[5]+ALGEBRAIC[6]) - ALGEBRAIC[7];
resid[6] = RATES[6] - ALGEBRAIC[5] - (ALGEBRAIC[6]+ALGEBRAIC[14]);
resid[7] = RATES[7] - ALGEBRAIC[7] - ALGEBRAIC[8];
resid[8] = RATES[8] - ALGEBRAIC[8] - ALGEBRAIC[9];
resid[9] = RATES[9] - ALGEBRAIC[9] - ALGEBRAIC[10];
resid[10] = RATES[10] - (ALGEBRAIC[9]+ALGEBRAIC[8]) - (ALGEBRAIC[4]+ALGEBRAIC[21]+ALGEBRAIC[2]+ALGEBRAIC[22]+ALGEBRAIC[23]);
resid[11] = RATES[11] - (ALGEBRAIC[4]+ALGEBRAIC[21]+ALGEBRAIC[2]+ALGEBRAIC[22]+ 2.00000*ALGEBRAIC[23]) - (ALGEBRAIC[9]+ALGEBRAIC[8]);
resid[12] = RATES[12] - - ALGEBRAIC[23];
resid[13] = RATES[13] - ALGEBRAIC[20] - ALGEBRAIC[19];
resid[14] = RATES[14] - ALGEBRAIC[10] - ( 59.0000*ALGEBRAIC[18]+ALGEBRAIC[11]);
resid[15] = RATES[15] - ALGEBRAIC[18] - (ALGEBRAIC[19]+ALGEBRAIC[17]);
resid[16] = RATES[16] - ALGEBRAIC[11] -  59.0000*ALGEBRAIC[13];
resid[17] = RATES[17] - ALGEBRAIC[13] - ALGEBRAIC[12];
resid[18] = RATES[18] - ALGEBRAIC[14] -  59.0000*ALGEBRAIC[16];
resid[19] = RATES[19] - ALGEBRAIC[16] - ALGEBRAIC[15];
resid[20] = RATES[20] - ALGEBRAIC[7] - (ALGEBRAIC[11]+ALGEBRAIC[14]);
resid[21] = RATES[21] - (ALGEBRAIC[11]+ALGEBRAIC[14]) - ALGEBRAIC[7];
}
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[2]*(CONSTANTS[0] - STATES[0]);
ALGEBRAIC[1] =  (CONSTANTS[6]/CONSTANTS[8])*((STATES[0]/CONSTANTS[3])/(1.00000+STATES[0]/CONSTANTS[3]+ (( CONSTANTS[7]*(STATES[0]/CONSTANTS[3])+1.00000)/( CONSTANTS[7]*(STATES[1]/CONSTANTS[3])+1.00000))*(1.00000+STATES[1]/CONSTANTS[3]+STATES[2]/CONSTANTS[4]+( STATES[1]*STATES[2])/( CONSTANTS[3]*CONSTANTS[5])))) -  (CONSTANTS[6]/CONSTANTS[8])*((STATES[1]/CONSTANTS[3])/(1.00000+STATES[1]/CONSTANTS[3]+ (( CONSTANTS[7]*(STATES[1]/CONSTANTS[3])+1.00000)/( CONSTANTS[7]*(STATES[0]/CONSTANTS[3])+1.00000))*(1.00000+STATES[0]/CONSTANTS[3])+STATES[2]/CONSTANTS[4]+( STATES[1]*STATES[2])/( CONSTANTS[3]*CONSTANTS[5])));
ALGEBRAIC[2] = ( CONSTANTS[12]*STATES[10]*STATES[1])/( CONSTANTS[9]*CONSTANTS[11]+ CONSTANTS[10]*STATES[10]+ CONSTANTS[11]*STATES[1]+ STATES[1]*STATES[10]);
ALGEBRAIC[3] = ( CONSTANTS[13]*STATES[2])/(CONSTANTS[14]+STATES[2]+ (CONSTANTS[14]/CONSTANTS[15])*STATES[3]) - ( CONSTANTS[13]*(STATES[3]/CONSTANTS[16]))/(CONSTANTS[14]+STATES[2]+ (CONSTANTS[14]/CONSTANTS[15])*STATES[3]);
ALGEBRAIC[4] = ( CONSTANTS[19]*pow(STATES[3], 2.00000))/( CONSTANTS[17]*(1.00000+ CONSTANTS[18]*(STATES[10]/STATES[12])*(STATES[10]/STATES[12]))+pow(STATES[3], 2.00000));
ALGEBRAIC[5] = ( CONSTANTS[26]*STATES[4])/(CONSTANTS[21]+STATES[4]+( STATES[5]*CONSTANTS[22]*CONSTANTS[26])/( CONSTANTS[20]*CONSTANTS[25])+( STATES[6]*CONSTANTS[23]*CONSTANTS[26])/( CONSTANTS[20]*CONSTANTS[25])+( STATES[4]*STATES[5])/CONSTANTS[24]+( STATES[5]*STATES[6]*CONSTANTS[26])/( CONSTANTS[20]*CONSTANTS[25])) - ( CONSTANTS[26]*(( STATES[5]*STATES[6])/CONSTANTS[20]))/(CONSTANTS[21]+STATES[4]+( STATES[5]*CONSTANTS[22]*CONSTANTS[26])/( CONSTANTS[20]*CONSTANTS[25])+( STATES[6]*CONSTANTS[23]*CONSTANTS[26])/( CONSTANTS[20]*CONSTANTS[25])+( STATES[4]*STATES[5])/CONSTANTS[24]+( STATES[5]*STATES[6]*CONSTANTS[26])/( CONSTANTS[20]*CONSTANTS[25]));
ALGEBRAIC[6] = ( CONSTANTS[30]*STATES[6])/(CONSTANTS[28]+STATES[6]+ (CONSTANTS[28]/CONSTANTS[29])*STATES[5]) - ( CONSTANTS[30]*(STATES[5]/CONSTANTS[27]))/(CONSTANTS[28]+STATES[6]+ (CONSTANTS[28]/CONSTANTS[29])*STATES[5]);
ALGEBRAIC[7] = ( CONSTANTS[36]*STATES[5]*STATES[21])/( CONSTANTS[33]*CONSTANTS[31]*(1.00000+STATES[5]/CONSTANTS[33]+STATES[7]/CONSTANTS[34])*(1.00000+STATES[21]/CONSTANTS[31]+STATES[20]/CONSTANTS[32])) - ( CONSTANTS[36]*(( STATES[7]*STATES[20])/CONSTANTS[35]))/( CONSTANTS[33]*CONSTANTS[31]*(1.00000+STATES[5]/CONSTANTS[33]+STATES[7]/CONSTANTS[34])*(1.00000+STATES[21]/CONSTANTS[31]+STATES[20]/CONSTANTS[32]));
ALGEBRAIC[8] =  CONSTANTS[37]*STATES[7]*STATES[11] -  CONSTANTS[38]*STATES[8]*STATES[10];
ALGEBRAIC[9] = ( CONSTANTS[41]*STATES[11]*STATES[8])/( (CONSTANTS[39]+STATES[8])*(CONSTANTS[40]+STATES[11]));
ALGEBRAIC[10] = ( CONSTANTS[43]*STATES[9])/(CONSTANTS[42]+STATES[9]);
ALGEBRAIC[11] = ( CONSTANTS[46]*STATES[14]*STATES[20])/( (CONSTANTS[44]+STATES[20])*(CONSTANTS[45]+STATES[14]));
ALGEBRAIC[12] =  CONSTANTS[2]*STATES[17];
ALGEBRAIC[13] =  (CONSTANTS[47]/CONSTANTS[8])*(STATES[16] - STATES[17]);
ALGEBRAIC[14] = ( CONSTANTS[52]*STATES[6])/( CONSTANTS[51]*(1.00000+ (CONSTANTS[49]/STATES[20])*(1.00000+STATES[21]/CONSTANTS[50]))+ STATES[6]*(1.00000+ (CONSTANTS[48]/STATES[20])*(1.00000+STATES[21]/CONSTANTS[50])));
ALGEBRAIC[15] =  CONSTANTS[2]*STATES[19];
ALGEBRAIC[16] =  (CONSTANTS[53]/CONSTANTS[8])*(STATES[18] - STATES[19]);
ALGEBRAIC[17] =  CONSTANTS[2]*STATES[15];
ALGEBRAIC[18] =  (CONSTANTS[54]/CONSTANTS[8])*(STATES[14] - STATES[15]);
ALGEBRAIC[19] =  CONSTANTS[55]*STATES[15]*STATES[13];
ALGEBRAIC[20] =  CONSTANTS[2]*(CONSTANTS[1] - STATES[13]);
ALGEBRAIC[21] =  CONSTANTS[56]*STATES[10]*STATES[2];
ALGEBRAIC[22] =  CONSTANTS[57]*STATES[10];
ALGEBRAIC[23] =  CONSTANTS[58]*STATES[12]*STATES[10] -  CONSTANTS[59]*pow(STATES[11], 2.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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}