/*
   There are a total of 16 entries in the algebraic variable array.
   There are a total of 6 entries in each of the rate and state variable arrays.
   There are a total of 19 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * CONSTANTS[0] is q_PKACI_init in component environment (fmol).
 * CONSTANTS[1] is q_Ip_init in component environment (fmol).
 * CONSTANTS[2] is q_Ip_PKACI_init in component environment (fmol).
 * CONSTANTS[3] is q_PP2a_init in component environment (fmol).
 * CONSTANTS[4] is q_Ip_PP2a_init in component environment (fmol).
 * CONSTANTS[5] is q_I_init in component environment (fmol).
 * ALGEBRAIC[0] is q_PKACI in component environment (fmol).
 * ALGEBRAIC[1] is q_Ip in component environment (fmol).
 * ALGEBRAIC[2] is q_Ip_PKACI in component environment (fmol).
 * ALGEBRAIC[3] is q_PP2a in component environment (fmol).
 * ALGEBRAIC[4] is q_Ip_PP2a in component environment (fmol).
 * ALGEBRAIC[5] is q_I in component environment (fmol).
 * STATES[0] is q_PKACI in component Inhib1 (fmol).
 * STATES[1] is q_Ip in component Inhib1 (fmol).
 * STATES[2] is q_Ip_PKACI in component Inhib1 (fmol).
 * STATES[3] is q_PP2a in component Inhib1 (fmol).
 * STATES[4] is q_Ip_PP2a in component Inhib1 (fmol).
 * STATES[5] is q_I in component Inhib1 (fmol).
 * CONSTANTS[6] is kappa_PKAi1 in component Inhib1_parameters (fmol_per_sec).
 * CONSTANTS[7] is kappa_PKAi2 in component Inhib1_parameters (fmol_per_sec).
 * CONSTANTS[8] is kappa_PP2ai1 in component Inhib1_parameters (fmol_per_sec).
 * CONSTANTS[9] is kappa_PP2ai2 in component Inhib1_parameters (fmol_per_sec).
 * CONSTANTS[10] is K_PKACI in component Inhib1_parameters (per_fmol).
 * CONSTANTS[11] is K_Ip in component Inhib1_parameters (per_fmol).
 * CONSTANTS[12] is K_Ip_PKACI in component Inhib1_parameters (per_fmol).
 * CONSTANTS[13] is K_PP2a in component Inhib1_parameters (per_fmol).
 * CONSTANTS[14] is K_Ip_PP2a in component Inhib1_parameters (per_fmol).
 * CONSTANTS[15] is K_I in component Inhib1_parameters (per_fmol).
 * CONSTANTS[16] is R in component constants (J_per_K_per_mol).
 * CONSTANTS[17] is T in component constants (kelvin).
 * ALGEBRAIC[12] is vPKAi1 in component Inhib1 (fmol_per_sec).
 * ALGEBRAIC[13] is vPKAi2 in component Inhib1 (fmol_per_sec).
 * ALGEBRAIC[14] is vPP2ai1 in component Inhib1 (fmol_per_sec).
 * ALGEBRAIC[15] is vPP2ai2 in component Inhib1 (fmol_per_sec).
 * ALGEBRAIC[6] is mu_PKACI in component Inhib1 (J_per_mol).
 * ALGEBRAIC[7] is mu_Ip in component Inhib1 (J_per_mol).
 * ALGEBRAIC[8] is mu_Ip_PKACI in component Inhib1 (J_per_mol).
 * ALGEBRAIC[9] is mu_PP2a in component Inhib1 (J_per_mol).
 * ALGEBRAIC[10] is mu_Ip_PP2a in component Inhib1 (J_per_mol).
 * ALGEBRAIC[11] is mu_I in component Inhib1 (J_per_mol).
 * CONSTANTS[18] is F in component constants (C_per_mol).
 * RATES[0] is d/dt q_PKACI in component Inhib1 (fmol).
 * RATES[1] is d/dt q_Ip in component Inhib1 (fmol).
 * RATES[2] is d/dt q_Ip_PKACI in component Inhib1 (fmol).
 * RATES[3] is d/dt q_PP2a in component Inhib1 (fmol).
 * RATES[4] is d/dt q_Ip_PP2a in component Inhib1 (fmol).
 * RATES[5] is d/dt q_I in component Inhib1 (fmol).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1.117E-03;
CONSTANTS[1] = 9.994E-04;
CONSTANTS[2] = 1e-13;
CONSTANTS[3] = 8.512E-03;
CONSTANTS[4] = 1e-13;
CONSTANTS[5] = 1.140E-02;
STATES[0] = 1e-16;
STATES[1] = 1e-16;
STATES[2] = 1e-16;
STATES[3] = 1e-16;
STATES[4] = 1e-16;
STATES[5] = 1e-16;
CONSTANTS[6] = 20.9606;
CONSTANTS[7] = 22.2706;
CONSTANTS[8] = 5.39341;
CONSTANTS[9] = 5.30013;
CONSTANTS[10] = 0.284433;
CONSTANTS[11] = 0.133405;
CONSTANTS[12] = 1.38688;
CONSTANTS[13] = 0.0704691;
CONSTANTS[14] = 0.342795;
CONSTANTS[15] = 0.150247;
CONSTANTS[16] = 8.31;
CONSTANTS[17] = 310;
CONSTANTS[18] = 96485;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = STATES[0]+CONSTANTS[0];
ALGEBRAIC[6] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[10]*ALGEBRAIC[0]);
ALGEBRAIC[2] = STATES[2]+CONSTANTS[2];
ALGEBRAIC[8] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[12]*ALGEBRAIC[2]);
ALGEBRAIC[5] = STATES[5]+CONSTANTS[5];
ALGEBRAIC[11] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[15]*ALGEBRAIC[5]);
ALGEBRAIC[12] =  CONSTANTS[6]*(exp((ALGEBRAIC[6]+ALGEBRAIC[11])/( CONSTANTS[16]*CONSTANTS[17])) - exp(ALGEBRAIC[8]/( CONSTANTS[16]*CONSTANTS[17])));
ALGEBRAIC[1] = STATES[1]+CONSTANTS[1];
ALGEBRAIC[7] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[11]*ALGEBRAIC[1]);
ALGEBRAIC[13] =  CONSTANTS[7]*(exp(ALGEBRAIC[8]/( CONSTANTS[16]*CONSTANTS[17])) - exp((ALGEBRAIC[6]+ALGEBRAIC[7])/( CONSTANTS[16]*CONSTANTS[17])));
RATES[0] = ALGEBRAIC[13] - ALGEBRAIC[12];
RATES[2] = ALGEBRAIC[12] - ALGEBRAIC[13];
ALGEBRAIC[3] = STATES[3]+CONSTANTS[3];
ALGEBRAIC[9] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[13]*ALGEBRAIC[3]);
ALGEBRAIC[4] = STATES[4]+CONSTANTS[4];
ALGEBRAIC[10] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[14]*ALGEBRAIC[4]);
ALGEBRAIC[14] =  CONSTANTS[8]*(exp((ALGEBRAIC[9]+ALGEBRAIC[7])/( CONSTANTS[16]*CONSTANTS[17])) - exp(ALGEBRAIC[10]/( CONSTANTS[16]*CONSTANTS[17])));
RATES[1] = ALGEBRAIC[13] - ALGEBRAIC[14];
ALGEBRAIC[15] =  CONSTANTS[9]*(exp(ALGEBRAIC[10]/( CONSTANTS[16]*CONSTANTS[17])) - exp((ALGEBRAIC[9]+ALGEBRAIC[11])/( CONSTANTS[16]*CONSTANTS[17])));
RATES[3] = ALGEBRAIC[15] - ALGEBRAIC[14];
RATES[4] = ALGEBRAIC[14] - ALGEBRAIC[15];
RATES[5] = ALGEBRAIC[15] - ALGEBRAIC[12];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = STATES[0]+CONSTANTS[0];
ALGEBRAIC[6] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[10]*ALGEBRAIC[0]);
ALGEBRAIC[2] = STATES[2]+CONSTANTS[2];
ALGEBRAIC[8] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[12]*ALGEBRAIC[2]);
ALGEBRAIC[5] = STATES[5]+CONSTANTS[5];
ALGEBRAIC[11] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[15]*ALGEBRAIC[5]);
ALGEBRAIC[12] =  CONSTANTS[6]*(exp((ALGEBRAIC[6]+ALGEBRAIC[11])/( CONSTANTS[16]*CONSTANTS[17])) - exp(ALGEBRAIC[8]/( CONSTANTS[16]*CONSTANTS[17])));
ALGEBRAIC[1] = STATES[1]+CONSTANTS[1];
ALGEBRAIC[7] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[11]*ALGEBRAIC[1]);
ALGEBRAIC[13] =  CONSTANTS[7]*(exp(ALGEBRAIC[8]/( CONSTANTS[16]*CONSTANTS[17])) - exp((ALGEBRAIC[6]+ALGEBRAIC[7])/( CONSTANTS[16]*CONSTANTS[17])));
ALGEBRAIC[3] = STATES[3]+CONSTANTS[3];
ALGEBRAIC[9] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[13]*ALGEBRAIC[3]);
ALGEBRAIC[4] = STATES[4]+CONSTANTS[4];
ALGEBRAIC[10] =  CONSTANTS[16]*CONSTANTS[17]*log( CONSTANTS[14]*ALGEBRAIC[4]);
ALGEBRAIC[14] =  CONSTANTS[8]*(exp((ALGEBRAIC[9]+ALGEBRAIC[7])/( CONSTANTS[16]*CONSTANTS[17])) - exp(ALGEBRAIC[10]/( CONSTANTS[16]*CONSTANTS[17])));
ALGEBRAIC[15] =  CONSTANTS[9]*(exp(ALGEBRAIC[10]/( CONSTANTS[16]*CONSTANTS[17])) - exp((ALGEBRAIC[9]+ALGEBRAIC[11])/( CONSTANTS[16]*CONSTANTS[17])));
}