/*
   There are a total of 18 entries in the algebraic variable array.
   There are a total of 4 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 (second).
 * STATES[0] is HCO3_int in component concentrations (mM).
 * STATES[1] is HCO3_ext in component concentrations (mM).
 * STATES[2] is Cl_int in component concentrations (mM).
 * STATES[3] is Cl_ext in component concentrations (mM).
 * ALGEBRAIC[12] is J_AE1_HCO3 in component AE1 (mM_per_s).
 * ALGEBRAIC[17] is J_AE1_Cl in component AE1 (mM_per_s).
 * ALGEBRAIC[13] is J_HCO3_influx in component AE1 (mM_per_s).
 * ALGEBRAIC[14] is J_Cl_influx in component AE1 (mM_per_s).
 * CONSTANTS[0] is K_HCO3_ext in component AE1 (mM).
 * CONSTANTS[1] is K_HCO3_int in component AE1 (mM).
 * CONSTANTS[2] is K_Cl_ext in component AE1 (mM).
 * CONSTANTS[3] is K_Cl_int in component AE1 (mM).
 * CONSTANTS[4] is P_HCO3_ext in component AE1 (per_s).
 * CONSTANTS[5] is P_HCO3_int in component AE1 (per_s).
 * CONSTANTS[6] is P_Cl_ext in component AE1 (per_s).
 * CONSTANTS[7] is P_Cl_int in component AE1 (per_s).
 * ALGEBRAIC[1] is beta_ext in component AE1 (dimensionless).
 * ALGEBRAIC[8] is beta_int in component AE1 (dimensionless).
 * ALGEBRAIC[9] is gamma_ext in component AE1 (dimensionless).
 * ALGEBRAIC[10] is gamma_int in component AE1 (dimensionless).
 * ALGEBRAIC[11] is sigma in component AE1 (per_s).
 * CONSTANTS[8] is x_Tmax in component AE1 (mM).
 * CONSTANTS[9] is K_I in component AE1 (mM).
 * ALGEBRAIC[0] is x_T in component AE1 (mM).
 * ALGEBRAIC[15] is x_ext in component AE1 (mM).
 * ALGEBRAIC[16] is x_int in component AE1 (mM).
 * ALGEBRAIC[2] is Jo_bm in component AE1 (mM_per_s).
 * ALGEBRAIC[3] is Ji_bm in component AE1 (mM_per_s).
 * ALGEBRAIC[4] is Js_bm in component AE1 (mM_per_s).
 * ALGEBRAIC[5] is Jo_cm in component AE1 (mM_per_s).
 * ALGEBRAIC[6] is Ji_cm in component AE1 (mM_per_s).
 * ALGEBRAIC[7] is Js_cm in component AE1 (mM_per_s).
 * RATES[0] is d/dt HCO3_int in component concentrations (mM).
 * RATES[1] is d/dt HCO3_ext in component concentrations (mM).
 * RATES[2] is d/dt Cl_int in component concentrations (mM).
 * RATES[3] is d/dt Cl_ext in component concentrations (mM).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 26.0;
STATES[1] = 26.0;
STATES[2] = 1.0;
STATES[3] = 114.0;
CONSTANTS[0] = 198;
CONSTANTS[1] = 198;
CONSTANTS[2] = 50;
CONSTANTS[3] = 50;
CONSTANTS[4] = 1247;
CONSTANTS[5] = 135;
CONSTANTS[6] = 562;
CONSTANTS[7] = 61;
CONSTANTS[8] = 1;
CONSTANTS[9] = 172;
CONSTANTS[10] = 0.00000;
CONSTANTS[11] = 0.00000;
CONSTANTS[12] = 1.00000;
CONSTANTS[13] = 0.00000;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[0] = CONSTANTS[10];
RATES[1] = CONSTANTS[11];
RATES[2] = CONSTANTS[12];
RATES[3] = CONSTANTS[13];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = CONSTANTS[8]/(1.00000+STATES[0]/CONSTANTS[9]);
ALGEBRAIC[1] = STATES[1]/CONSTANTS[0];
ALGEBRAIC[2] = pow( (1.00000/ALGEBRAIC[0])*(1.00000/CONSTANTS[4]+1.00000/CONSTANTS[5]+CONSTANTS[1]/( CONSTANTS[5]*STATES[0])), - 1.00000);
ALGEBRAIC[3] = pow( (1.00000/ALGEBRAIC[0])*(1.00000/CONSTANTS[4]+1.00000/CONSTANTS[5]+CONSTANTS[0]/( CONSTANTS[4]*STATES[1])), - 1.00000);
ALGEBRAIC[4] = pow( (1.00000/ALGEBRAIC[0])*(1.00000/CONSTANTS[4]+1.00000/CONSTANTS[5]), - 1.00000);
ALGEBRAIC[5] = pow( (1.00000/ALGEBRAIC[0])*(1.00000/CONSTANTS[6]+1.00000/CONSTANTS[7]+CONSTANTS[3]/( CONSTANTS[7]*STATES[2])), - 1.00000);
ALGEBRAIC[6] = pow( (1.00000/ALGEBRAIC[0])*(1.00000/CONSTANTS[6]+1.00000/CONSTANTS[7]+CONSTANTS[2]/( CONSTANTS[6]*STATES[3])), - 1.00000);
ALGEBRAIC[7] = pow( (1.00000/ALGEBRAIC[0])*(1.00000/CONSTANTS[6]+1.00000/CONSTANTS[7]), - 1.00000);
ALGEBRAIC[8] = STATES[0]/CONSTANTS[1];
ALGEBRAIC[9] = STATES[3]/CONSTANTS[2];
ALGEBRAIC[10] = STATES[2]/CONSTANTS[3];
ALGEBRAIC[11] =  (1.00000+ALGEBRAIC[1]+ALGEBRAIC[9])*( CONSTANTS[5]*ALGEBRAIC[8]+ CONSTANTS[7]*ALGEBRAIC[10])+ (1.00000+ALGEBRAIC[8]+ALGEBRAIC[10])*( CONSTANTS[4]*ALGEBRAIC[1]+ CONSTANTS[6]*ALGEBRAIC[9]);
ALGEBRAIC[12] =  (ALGEBRAIC[0]/ALGEBRAIC[11])*( CONSTANTS[5]*ALGEBRAIC[8]*CONSTANTS[6]*ALGEBRAIC[9] -  CONSTANTS[4]*ALGEBRAIC[1]*CONSTANTS[7]*ALGEBRAIC[10]);
ALGEBRAIC[13] =  (ALGEBRAIC[0]/ALGEBRAIC[11])*CONSTANTS[4]*ALGEBRAIC[1]*( CONSTANTS[5]*ALGEBRAIC[8]+ CONSTANTS[7]*ALGEBRAIC[10]);
ALGEBRAIC[14] =  (ALGEBRAIC[0]/ALGEBRAIC[11])*CONSTANTS[6]*ALGEBRAIC[9]*( CONSTANTS[5]*ALGEBRAIC[8]+ CONSTANTS[7]*ALGEBRAIC[10]);
ALGEBRAIC[15] = ( ALGEBRAIC[0]*( CONSTANTS[5]*ALGEBRAIC[8]+ CONSTANTS[7]*ALGEBRAIC[10]))/ALGEBRAIC[11];
ALGEBRAIC[16] = ( ALGEBRAIC[0]*( CONSTANTS[4]*ALGEBRAIC[1]+ CONSTANTS[6]*ALGEBRAIC[9]))/ALGEBRAIC[11];
ALGEBRAIC[17] = - ALGEBRAIC[12];
}