/*
   There are a total of 5 entries in the algebraic variable array.
   There are a total of 9 entries in each of the rate and state variable arrays.
   There are a total of 31 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * CONSTANTS[0] is k1AA in component rate_constants (first_order_rate_constant).
 * CONSTANTS[1] is k3 in component rate_constants (first_order_rate_constant).
 * CONSTANTS[2] is kinh in component rate_constants (first_order_rate_constant).
 * CONSTANTS[3] is kcak in component rate_constants (first_order_rate_constant).
 * CONSTANTS[4] is kbPPase in component rate_constants (first_order_rate_constant).
 * CONSTANTS[5] is kfPPase in component rate_constants (first_order_rate_constant).
 * CONSTANTS[6] is khPPase in component rate_constants (first_order_rate_constant).
 * CONSTANTS[7] is kd_antiIE in component rate_constants (first_order_rate_constant).
 * CONSTANTS[8] is ka in component rate_constants (first_order_rate_constant).
 * CONSTANTS[9] is kc in component rate_constants (first_order_rate_constant).
 * CONSTANTS[10] is ke in component rate_constants (first_order_rate_constant).
 * CONSTANTS[11] is kg in component rate_constants (first_order_rate_constant).
 * CONSTANTS[12] is V2_ in component rate_constants (first_order_rate_constant).
 * CONSTANTS[13] is V2__ in component rate_constants (first_order_rate_constant).
 * CONSTANTS[14] is V25_ in component rate_constants (first_order_rate_constant).
 * CONSTANTS[15] is V25__ in component rate_constants (first_order_rate_constant).
 * CONSTANTS[16] is Vwee_ in component rate_constants (first_order_rate_constant).
 * CONSTANTS[17] is Vwee__ in component rate_constants (first_order_rate_constant).
 * CONSTANTS[18] is K_a in component rate_constants (dimensionless).
 * CONSTANTS[19] is K_b in component rate_constants (dimensionless).
 * CONSTANTS[20] is K_c in component rate_constants (dimensionless).
 * CONSTANTS[21] is K_d in component rate_constants (dimensionless).
 * CONSTANTS[22] is K_e in component rate_constants (dimensionless).
 * CONSTANTS[23] is K_f in component rate_constants (dimensionless).
 * CONSTANTS[24] is K_g in component rate_constants (dimensionless).
 * CONSTANTS[25] is K_h in component rate_constants (dimensionless).
 * CONSTANTS[26] is Cdc25_T in component concentration_variables (dimensionless).
 * STATES[0] is Cdc25_P in component Cdc25_P (dimensionless).
 * STATES[1] is Wee1_P in component Wee1_P (dimensionless).
 * CONSTANTS[27] is Wee1_T in component concentration_variables (dimensionless).
 * CONSTANTS[28] is UbE_T in component concentration_variables (dimensionless).
 * STATES[2] is UbE in component UbE (dimensionless).
 * ALGEBRAIC[0] is k25 in component rate_constants (first_order_rate_constant).
 * ALGEBRAIC[1] is kwee in component rate_constants (first_order_rate_constant).
 * ALGEBRAIC[2] is k2 in component rate_constants (first_order_rate_constant).
 * CONSTANTS[29] is Cdc2_T in component concentration_variables (dimensionless).
 * CONSTANTS[30] is IE_T in component concentration_variables (dimensionless).
 * STATES[3] is cyclin in component cyclin (dimensionless).
 * ALGEBRAIC[3] is Cdc2 in component Cdc2 (dimensionless).
 * STATES[4] is cyclin_Cdc2 in component cyclin_Cdc2 (dimensionless).
 * STATES[5] is MPF_active in component MPF_active (dimensionless).
 * STATES[6] is Tyr15P_dimer in component Tyr15P_dimer (dimensionless).
 * STATES[7] is MPF_pre in component MPF_pre (dimensionless).
 * ALGEBRAIC[4] is cyclin_T in component cyclin_T (dimensionless).
 * STATES[8] is IE_P in component IE_P (dimensionless).
 * RATES[3] is d/dt cyclin in component cyclin (dimensionless).
 * RATES[4] is d/dt cyclin_Cdc2 in component cyclin_Cdc2 (dimensionless).
 * RATES[6] is d/dt Tyr15P_dimer in component Tyr15P_dimer (dimensionless).
 * RATES[7] is d/dt MPF_pre in component MPF_pre (dimensionless).
 * RATES[5] is d/dt MPF_active in component MPF_active (dimensionless).
 * RATES[0] is d/dt Cdc25_P in component Cdc25_P (dimensionless).
 * RATES[1] is d/dt Wee1_P in component Wee1_P (dimensionless).
 * RATES[8] is d/dt IE_P in component IE_P (dimensionless).
 * RATES[2] is d/dt UbE in component UbE (dimensionless).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1;
CONSTANTS[1] = 0.01;
CONSTANTS[2] = 0.025;
CONSTANTS[3] = 0.25;
CONSTANTS[4] = 0.125;
CONSTANTS[5] = 0.1;
CONSTANTS[6] = 0.087;
CONSTANTS[7] = 0.095;
CONSTANTS[8] = 0.01;
CONSTANTS[9] = 0.1;
CONSTANTS[10] = 0.0133;
CONSTANTS[11] = 0.0065;
CONSTANTS[12] = 0.015;
CONSTANTS[13] = 1;
CONSTANTS[14] = 0.1;
CONSTANTS[15] = 2;
CONSTANTS[16] = 0.1;
CONSTANTS[17] = 1;
CONSTANTS[18] = 0.1;
CONSTANTS[19] = 0.1;
CONSTANTS[20] = 0.01;
CONSTANTS[21] = 0.01;
CONSTANTS[22] = 0.3;
CONSTANTS[23] = 0.3;
CONSTANTS[24] = 0.01;
CONSTANTS[25] = 0.01;
CONSTANTS[26] = 1;
STATES[0] = 0;
STATES[1] = 0;
CONSTANTS[27] = 1;
CONSTANTS[28] = 1;
STATES[2] = 0;
CONSTANTS[29] = 100;
CONSTANTS[30] = 1;
STATES[3] = 1;
STATES[4] = 0;
STATES[5] = 0;
STATES[6] = 0;
STATES[7] = 0;
STATES[8] = 0;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[5] = 0.1001;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[8] = 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[3] - (CONSTANTS[0] -  ALGEBRAIC[2]*STATES[3]) -  CONSTANTS[1]*STATES[3]*ALGEBRAIC[3];
resid[1] = RATES[4] - ( CONSTANTS[2]*STATES[5] -  (ALGEBRAIC[1]+CONSTANTS[3]+ALGEBRAIC[2])*STATES[4])+ ALGEBRAIC[0]*STATES[6]+ CONSTANTS[1]*STATES[3]*ALGEBRAIC[3];
resid[2] = RATES[6] - ( ALGEBRAIC[1]*STATES[4] -  (ALGEBRAIC[0]+CONSTANTS[3]+ALGEBRAIC[2])*STATES[6])+ CONSTANTS[2]*STATES[7];
resid[3] = RATES[7] - ( ALGEBRAIC[1]*STATES[5] -  (CONSTANTS[2]+ALGEBRAIC[0]+ALGEBRAIC[2])*STATES[7])+ CONSTANTS[3]*STATES[6];
resid[4] = RATES[5] - ( CONSTANTS[3]*STATES[4] -  (CONSTANTS[2]+ALGEBRAIC[1]+ALGEBRAIC[2])*STATES[5])+ ALGEBRAIC[0]*STATES[7];
resid[5] = RATES[0] - ( CONSTANTS[8]*STATES[5]*(CONSTANTS[26] - STATES[0]))/((CONSTANTS[18]+CONSTANTS[26]) - STATES[0]) - ( CONSTANTS[4]*STATES[0])/(CONSTANTS[19]+STATES[0]);
resid[6] = RATES[1] - ( CONSTANTS[10]*STATES[5]*(CONSTANTS[27] - STATES[1]))/((CONSTANTS[22]+CONSTANTS[27]) - STATES[1]) - ( CONSTANTS[5]*STATES[1])/(CONSTANTS[23]+STATES[1]);
resid[7] = RATES[8] - ( CONSTANTS[11]*STATES[5]*(CONSTANTS[30] - STATES[8]))/((CONSTANTS[24]+CONSTANTS[30]) - STATES[8]) - ( CONSTANTS[6]*STATES[8])/(CONSTANTS[25]+STATES[8]);
resid[8] = RATES[2] - ( CONSTANTS[9]*STATES[8]*(CONSTANTS[28] - STATES[2]))/((CONSTANTS[20]+CONSTANTS[28]) - STATES[2]) - ( CONSTANTS[7]*STATES[2])/(CONSTANTS[21]+STATES[2]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = STATES[3]+STATES[4]+STATES[5]+STATES[7]+STATES[6];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] =  CONSTANTS[14]*(CONSTANTS[26] - STATES[0])+ CONSTANTS[15]*STATES[0];
ALGEBRAIC[1] =  CONSTANTS[16]*STATES[1]+ CONSTANTS[17]*(CONSTANTS[27] - STATES[1]);
ALGEBRAIC[2] =  CONSTANTS[12]*(CONSTANTS[28] - STATES[2])+ CONSTANTS[13]*STATES[2];
ALGEBRAIC[3] = CONSTANTS[29] - (STATES[4]+STATES[5]+STATES[7]+STATES[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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}