/*
   There are a total of 11 entries in the algebraic variable array.
   There are a total of 7 entries in each of the rate and state variable arrays.
   There are a total of 14 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is S1 in component S1 (millimolar).
 * CONSTANTS[0] is Jo in component glucose_influx_rate (flux).
 * ALGEBRAIC[3] is v1 in component v1 (flux).
 * STATES[1] is S2 in component S2 (millimolar).
 * ALGEBRAIC[4] is v2 in component v2 (flux).
 * ALGEBRAIC[8] is v6 in component v6 (flux).
 * STATES[2] is S3 in component S3 (millimolar).
 * ALGEBRAIC[5] is v3 in component v3 (flux).
 * STATES[3] is S4 in component S4 (millimolar).
 * ALGEBRAIC[6] is v4 in component v4 (flux).
 * ALGEBRAIC[10] is J in component S4_flux_rate_across_the_plasma_membrane (flux).
 * STATES[4] is S4_ex in component S4_ex (millimolar).
 * CONSTANTS[1] is phi in component S4_ex (dimensionless).
 * ALGEBRAIC[9] is v7 in component v7 (flux).
 * STATES[5] is A3 in component A3 (millimolar).
 * ALGEBRAIC[7] is v5 in component v5 (flux).
 * CONSTANTS[2] is A in component A2 (millimolar).
 * ALGEBRAIC[0] is A2 in component A2 (millimolar).
 * STATES[6] is N2 in component N2 (millimolar).
 * CONSTANTS[3] is N in component N1 (millimolar).
 * ALGEBRAIC[1] is N1 in component N1 (millimolar).
 * CONSTANTS[4] is K_I in component v1 (millimolar).
 * CONSTANTS[5] is k_1 in component v1 (second_order_rate_constant).
 * CONSTANTS[6] is q in component v1 (dimensionless).
 * ALGEBRAIC[2] is f_A3 in component v1 (dimensionless).
 * CONSTANTS[7] is k_2 in component v2 (second_order_rate_constant).
 * CONSTANTS[8] is k_3 in component v3 (second_order_rate_constant).
 * CONSTANTS[9] is k_4 in component v4 (second_order_rate_constant).
 * CONSTANTS[10] is k_5 in component v5 (first_order_rate_constant).
 * CONSTANTS[11] is k_6 in component v6 (second_order_rate_constant).
 * CONSTANTS[12] is k in component v7 (first_order_rate_constant).
 * CONSTANTS[13] is kappa in component S4_flux_rate_across_the_plasma_membrane (first_order_rate_constant).
 * RATES[0] is d/dt S1 in component S1 (millimolar).
 * RATES[1] is d/dt S2 in component S2 (millimolar).
 * RATES[2] is d/dt S3 in component S3 (millimolar).
 * RATES[3] is d/dt S4 in component S4 (millimolar).
 * RATES[4] is d/dt S4_ex in component S4_ex (millimolar).
 * RATES[5] is d/dt A3 in component A3 (millimolar).
 * RATES[6] is d/dt N2 in component N2 (millimolar).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 5.8;
CONSTANTS[0] = 3.0;
STATES[1] = 0.9;
STATES[2] = 0.2;
STATES[3] = 0.2;
STATES[4] = 0.1;
CONSTANTS[1] = 0.1;
STATES[5] = 2.4;
CONSTANTS[2] = 4.0;
STATES[6] = 0.1;
CONSTANTS[3] = 1.0;
CONSTANTS[4] = 0.52;
CONSTANTS[5] = 100.0;
CONSTANTS[6] = 4.0;
CONSTANTS[7] = 6.0;
CONSTANTS[8] = 16.0;
CONSTANTS[9] = 100.0;
CONSTANTS[10] = 1.28;
CONSTANTS[11] = 12.0;
CONSTANTS[12] = 1.3;
CONSTANTS[13] = 13.0;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = pow(1.00000+pow(STATES[5]/CONSTANTS[4], CONSTANTS[6]), -1.00000);
ALGEBRAIC[3] =  CONSTANTS[5]*STATES[0]*STATES[5]*ALGEBRAIC[2];
RATES[0] = CONSTANTS[0] - ALGEBRAIC[3];
ALGEBRAIC[1] = CONSTANTS[3] - STATES[6];
ALGEBRAIC[4] =  CONSTANTS[7]*STATES[1]*ALGEBRAIC[1];
ALGEBRAIC[0] = CONSTANTS[2] - STATES[5];
ALGEBRAIC[5] =  CONSTANTS[8]*STATES[2]*ALGEBRAIC[0];
RATES[2] = ALGEBRAIC[4] - ALGEBRAIC[5];
ALGEBRAIC[7] =  CONSTANTS[10]*STATES[5];
RATES[5] =  2.00000*ALGEBRAIC[5] - ( 2.00000*ALGEBRAIC[3]+ALGEBRAIC[7]);
ALGEBRAIC[8] =  CONSTANTS[11]*STATES[1]*STATES[6];
RATES[1] =  2.00000*ALGEBRAIC[3] - (ALGEBRAIC[4]+ALGEBRAIC[8]);
ALGEBRAIC[6] =  CONSTANTS[9]*STATES[3]*STATES[6];
RATES[6] = ALGEBRAIC[4] - (ALGEBRAIC[6]+ALGEBRAIC[8]);
ALGEBRAIC[10] =  CONSTANTS[13]*(STATES[3] - STATES[4]);
RATES[3] = ALGEBRAIC[5] - (ALGEBRAIC[6]+ALGEBRAIC[10]);
ALGEBRAIC[9] =  CONSTANTS[12]*STATES[4];
RATES[4] =  CONSTANTS[1]*ALGEBRAIC[10] - ALGEBRAIC[9];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = pow(1.00000+pow(STATES[5]/CONSTANTS[4], CONSTANTS[6]), -1.00000);
ALGEBRAIC[3] =  CONSTANTS[5]*STATES[0]*STATES[5]*ALGEBRAIC[2];
ALGEBRAIC[1] = CONSTANTS[3] - STATES[6];
ALGEBRAIC[4] =  CONSTANTS[7]*STATES[1]*ALGEBRAIC[1];
ALGEBRAIC[0] = CONSTANTS[2] - STATES[5];
ALGEBRAIC[5] =  CONSTANTS[8]*STATES[2]*ALGEBRAIC[0];
ALGEBRAIC[7] =  CONSTANTS[10]*STATES[5];
ALGEBRAIC[8] =  CONSTANTS[11]*STATES[1]*STATES[6];
ALGEBRAIC[6] =  CONSTANTS[9]*STATES[3]*STATES[6];
ALGEBRAIC[10] =  CONSTANTS[13]*(STATES[3] - STATES[4]);
ALGEBRAIC[9] =  CONSTANTS[12]*STATES[4];
}