/*
   There are a total of 6 entries in the algebraic variable array.
   There are a total of 3 entries in each of the rate and state variable arrays.
   There are a total of 28 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * CONSTANTS[27] is v_GK in component v_GK (flux).
 * CONSTANTS[0] is V_GK_min in component v_GK (enzyme_activity).
 * CONSTANTS[24] is V_GK in component v_GK (enzyme_activity).
 * CONSTANTS[1] is Sgk in component v_GK (millimolar).
 * CONSTANTS[2] is h_GK in component v_GK (dimensionless).
 * CONSTANTS[26] is Glc in component Glc (millimolar).
 * CONSTANTS[3] is min_to_sec in component model_parameters (dimensionless).
 * CONSTANTS[4] is dw_per_ml in component model_parameters (dimensionless).
 * ALGEBRAIC[2] is v_PFK in component v_PFK (flux).
 * CONSTANTS[5] is V_PFK_min in component v_PFK (enzyme_activity).
 * CONSTANTS[22] is V_PFK in component v_PFK (enzyme_activity).
 * CONSTANTS[6] is Spfk in component v_PFK (millimolar).
 * CONSTANTS[7] is Sfba in component v_PFK (millimolar).
 * CONSTANTS[8] is Xpfk in component v_PFK (millimolar).
 * CONSTANTS[9] is hx in component v_PFK (dimensionless).
 * CONSTANTS[10] is alpha in component v_PFK (dimensionless).
 * CONSTANTS[11] is h_PFK in component v_PFK (dimensionless).
 * CONSTANTS[12] is h_act in component v_PFK (dimensionless).
 * STATES[0] is FBP in component FBP (millimolar).
 * ALGEBRAIC[0] is F6P in component F6P (millimolar).
 * ALGEBRAIC[5] is v_FBA in component v_FBA (flux).
 * CONSTANTS[13] is V_FBA_min in component v_FBA (enzyme_activity).
 * CONSTANTS[23] is V_FBA in component v_FBA (enzyme_activity).
 * CONSTANTS[14] is Qfba in component v_FBA (millimolar).
 * CONSTANTS[15] is Sfba in component v_FBA (millimolar).
 * CONSTANTS[16] is Pfba in component v_FBA (millimolar).
 * CONSTANTS[17] is Keq_FBA in component v_FBA (millimolar).
 * ALGEBRAIC[1] is G3P in component G3P (millimolar).
 * ALGEBRAIC[4] is DHAP in component DHAP (millimolar).
 * ALGEBRAIC[3] is v_GAPDH in component v_GAPDH (flux).
 * CONSTANTS[18] is V_GAPDH_min in component v_GAPDH (enzyme_activity).
 * CONSTANTS[25] is V_GAPDH in component v_GAPDH (enzyme_activity).
 * CONSTANTS[19] is Sgapdh in component v_GAPDH (millimolar).
 * STATES[1] is G6P_F6P in component G6P_F6P (millimolar).
 * CONSTANTS[20] is Keq_GPI in component F6P (dimensionless).
 * STATES[2] is DHAP_G3P in component DHAP_G3P (millimolar).
 * CONSTANTS[21] is Keq_TPI in component G3P (dimensionless).
 * RATES[1] is d/dt G6P_F6P in component G6P_F6P (millimolar).
 * RATES[0] is d/dt FBP in component FBP (millimolar).
 * RATES[2] is d/dt DHAP_G3P in component DHAP_G3P (millimolar).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 10.0;
CONSTANTS[1] = 8.0;
CONSTANTS[2] = 1.7;
CONSTANTS[3] = 60.0;
CONSTANTS[4] = 0.3333;
CONSTANTS[5] = 100.0;
CONSTANTS[6] = 4.0;
CONSTANTS[7] = 0.005;
CONSTANTS[8] = 0.01;
CONSTANTS[9] = 2.5;
CONSTANTS[10] = 5.0;
CONSTANTS[11] = 2.5;
CONSTANTS[12] = 1.0;
STATES[0] = 0.00063612;
CONSTANTS[13] = 25.0;
CONSTANTS[14] = 0.275;
CONSTANTS[15] = 0.005;
CONSTANTS[16] = 0.5;
CONSTANTS[17] = 0.1;
CONSTANTS[18] = 250.0;
CONSTANTS[19] = 0.005;
STATES[1] = 3.71728;
CONSTANTS[20] = 0.3;
STATES[2] = 0.00262966;
CONSTANTS[21] = 0.045455;
CONSTANTS[22] = ( CONSTANTS[5]*CONSTANTS[4])/CONSTANTS[3];
CONSTANTS[23] = ( CONSTANTS[13]*CONSTANTS[4])/CONSTANTS[3];
CONSTANTS[24] = ( CONSTANTS[0]*CONSTANTS[4])/CONSTANTS[3];
CONSTANTS[25] = ( CONSTANTS[18]*CONSTANTS[4])/CONSTANTS[3];
CONSTANTS[26] = 10.0000;
CONSTANTS[27] = ( CONSTANTS[24]*pow(CONSTANTS[26]/CONSTANTS[1], CONSTANTS[2]))/(1.00000+pow(CONSTANTS[26]/CONSTANTS[1], CONSTANTS[2]));
RATES[1] = 0.1001;
RATES[0] = 0.1001;
RATES[2] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[1] - CONSTANTS[27] - ALGEBRAIC[2];
resid[1] = RATES[0] - ALGEBRAIC[2] - ALGEBRAIC[5];
resid[2] = RATES[2] -  2.00000*ALGEBRAIC[5] - ALGEBRAIC[3];
}
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] = ( STATES[1]*CONSTANTS[20])/(1.00000+CONSTANTS[20]);
ALGEBRAIC[2] = ( CONSTANTS[22]*pow(ALGEBRAIC[0]/CONSTANTS[6], CONSTANTS[11] -  (CONSTANTS[11] - CONSTANTS[12])*((STATES[0]/CONSTANTS[7])/(1.00000+STATES[0]/CONSTANTS[7]))))/(pow(ALGEBRAIC[0]/CONSTANTS[6], CONSTANTS[11] -  (CONSTANTS[11] - CONSTANTS[12])*((STATES[0]/CONSTANTS[7])/(1.00000+STATES[0]/CONSTANTS[7])))+(1.00000+pow(STATES[0]/CONSTANTS[8], CONSTANTS[9]))/(1.00000+ pow(CONSTANTS[10], CONSTANTS[11] -  (CONSTANTS[11] - CONSTANTS[12])*((STATES[0]/CONSTANTS[7])/(1.00000+STATES[0]/CONSTANTS[7])))*pow(STATES[0]/CONSTANTS[8], CONSTANTS[9])));
ALGEBRAIC[1] = ( STATES[2]*CONSTANTS[21])/(1.00000+CONSTANTS[21]);
ALGEBRAIC[3] = ( CONSTANTS[25]*ALGEBRAIC[1])/(CONSTANTS[19]+ALGEBRAIC[1]);
ALGEBRAIC[4] = STATES[2] - ALGEBRAIC[1];
ALGEBRAIC[5] = ( CONSTANTS[23]*(STATES[0]/CONSTANTS[15] - ( ALGEBRAIC[1]*ALGEBRAIC[4])/( CONSTANTS[16]*CONSTANTS[14]*CONSTANTS[17])))/(1.00000+STATES[0]/CONSTANTS[15]+ALGEBRAIC[4]/CONSTANTS[14]+( ALGEBRAIC[1]*ALGEBRAIC[4])/( CONSTANTS[16]*CONSTANTS[14]));
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}