/*
   There are a total of 3 entries in the algebraic variable array.
   There are a total of 2 entries in each of the rate and state variable arrays.
   There are a total of 9 entries in the constant variable array.
 */
/*
 * VOI is t in component main (second).
 * STATES[0] is q_1 in component main (mole).
 * STATES[1] is q_2 in component main (mole).
 * CONSTANTS[0] is q in component main (mole).
 * CONSTANTS[1] is E_tot in component main (mole).
 * ALGEBRAIC[0] is v in component main (mol_per_s).
 * CONSTANTS[2] is kappa_3 in component main (mol_per_s).
 * CONSTANTS[3] is kappa_4 in component main (mol_per_s).
 * CONSTANTS[4] is K_s in component main (per_mol).
 * CONSTANTS[5] is K_m in component main (per_mol).
 * CONSTANTS[7] is k_m in component main (mole).
 * CONSTANTS[8] is v_max in component main (mol_per_s).
 * ALGEBRAIC[1] is v_MM in component main (mol_per_s).
 * CONSTANTS[6] is q20 in component main (mole).
 * ALGEBRAIC[2] is v_0 in component main (mol_per_s).
 * RATES[0] is d/dt q_1 in component main (mole).
 * RATES[1] is d/dt q_2 in component main (mole).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 10;
STATES[1] = 0;
CONSTANTS[0] = 0;
CONSTANTS[1] = 10;
CONSTANTS[2] = 1;
CONSTANTS[3] = 0.1;
CONSTANTS[4] = 0.1;
CONSTANTS[5] = 1;
CONSTANTS[6] = 0;
CONSTANTS[7] = 2.00000/( CONSTANTS[4]*(1.00000+CONSTANTS[2]/CONSTANTS[3]));
CONSTANTS[8] = ( CONSTANTS[1]*CONSTANTS[2]*CONSTANTS[5])/(1.00000+CONSTANTS[2]/CONSTANTS[3]);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = ( CONSTANTS[1]*CONSTANTS[2]*CONSTANTS[3]*CONSTANTS[5]*CONSTANTS[4]*(STATES[0] - STATES[1]))/( 2.00000*CONSTANTS[3]+ (CONSTANTS[2]+CONSTANTS[3])*CONSTANTS[4]*(STATES[0]+STATES[1])+ 2.00000*CONSTANTS[2]*CONSTANTS[4]*CONSTANTS[4]*STATES[0]*STATES[1]);
RATES[0] = - ALGEBRAIC[0];
RATES[1] = ALGEBRAIC[0] -  0.100000*STATES[1];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = ( CONSTANTS[1]*CONSTANTS[2]*CONSTANTS[3]*CONSTANTS[5]*CONSTANTS[4]*(STATES[0] - STATES[1]))/( 2.00000*CONSTANTS[3]+ (CONSTANTS[2]+CONSTANTS[3])*CONSTANTS[4]*(STATES[0]+STATES[1])+ 2.00000*CONSTANTS[2]*CONSTANTS[4]*CONSTANTS[4]*STATES[0]*STATES[1]);
ALGEBRAIC[1] = ( CONSTANTS[8]*STATES[0])/(CONSTANTS[7]+STATES[0]);
ALGEBRAIC[2] = ( CONSTANTS[1]*CONSTANTS[2]*CONSTANTS[3]*CONSTANTS[5]*CONSTANTS[4]*(STATES[0] - CONSTANTS[6]))/( 2.00000*CONSTANTS[3]+ (CONSTANTS[2]+CONSTANTS[3])*CONSTANTS[4]*(STATES[0]+CONSTANTS[6])+ 2.00000*CONSTANTS[2]*CONSTANTS[4]*CONSTANTS[4]*STATES[0]*CONSTANTS[6]);
}