/*
   There are a total of 0 entries in the algebraic variable array.
   There are a total of 0 entries in each of the rate and state variable arrays.
   There are a total of 30 entries in the constant variable array.
 */
/*
 * CONSTANTS[28] is SHbO2 in component SHbO2 (dimensionless).
 * CONSTANTS[26] is KHbO2 in component KHbO2 (per_molar).
 * CONSTANTS[20] is O2 in component O2 (molar).
 * CONSTANTS[29] is SHbCO2 in component SHbCO2 (dimensionless).
 * CONSTANTS[27] is KHbCO2 in component KHbCO2 (per_molar).
 * CONSTANTS[22] is CO2 in component CO2 (molar).
 * CONSTANTS[24] is Hrbc in component model_parameters (molar).
 * CONSTANTS[0] is K2 in component model_parameters (per_molar).
 * CONSTANTS[1] is K2_ in component model_parameters (molar).
 * CONSTANTS[2] is K3 in component model_parameters (per_molar).
 * CONSTANTS[3] is K3_ in component model_parameters (molar).
 * CONSTANTS[25] is K4 in component K4 (per_molar).
 * CONSTANTS[4] is K5_ in component model_parameters (molar).
 * CONSTANTS[5] is K6_ in component model_parameters (molar).
 * CONSTANTS[6] is O2_S in component model_parameters (micromolar).
 * CONSTANTS[7] is H_S in component model_parameters (nanomolar).
 * CONSTANTS[8] is n1 in component model_parameters (dimensionless).
 * CONSTANTS[9] is n2 in component model_parameters (dimensionless).
 * CONSTANTS[10] is CO2_S in component model_parameters (millimolar).
 * CONSTANTS[11] is K4_ in component model_parameters (per_molar).
 * CONSTANTS[12] is n0 in component model_parameters (dimensionless).
 * CONSTANTS[19] is alpha_O2 in component alpha_O2 (M_mmHg).
 * CONSTANTS[13] is PO2 in component model_parameters (mmHg).
 * CONSTANTS[21] is alpha_CO2 in component alpha_CO2 (M_mmHg).
 * CONSTANTS[14] is PCO2 in component model_parameters (mmHg).
 * CONSTANTS[15] is Wpl in component model_parameters (ml_ml).
 * CONSTANTS[16] is T in component model_parameters (celsius).
 * CONSTANTS[17] is Rrbc in component model_parameters (dimensionless).
 * CONSTANTS[23] is Hpl in component model_parameters (molar).
 * CONSTANTS[18] is pHpl in component model_parameters (pH).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 29.5;
CONSTANTS[1] = 1E-6;
CONSTANTS[2] = 25.1;
CONSTANTS[3] = 1E-6;
CONSTANTS[4] = 2.63E-8;
CONSTANTS[5] = 1.91E-8;
CONSTANTS[6] = 146.0;
CONSTANTS[7] = 57.5;
CONSTANTS[8] = 1.06;
CONSTANTS[9] = 0.12;
CONSTANTS[10] = 1.31;
CONSTANTS[11] = 202123.0;
CONSTANTS[12] = 1.7;
CONSTANTS[13] = 100.0;
CONSTANTS[14] = 40.0;
CONSTANTS[15] = 0.94;
CONSTANTS[16] = 37.0;
CONSTANTS[17] = 0.69;
CONSTANTS[18] = 7.24;
CONSTANTS[19] =  ((1.37000 -  0.0137000*(CONSTANTS[16] - 37.0000))+ 0.000580000*pow(CONSTANTS[16] - 37.0000, 2.00000))*(1.00000e-06/CONSTANTS[15]);
CONSTANTS[20] =  CONSTANTS[19]*CONSTANTS[13];
CONSTANTS[21] =  ((3.07000 -  0.0570000*(CONSTANTS[16] - 37.0000))+ 0.00200000*pow(CONSTANTS[16] - 37.0000, 2.00000))*(1.00000e-05/CONSTANTS[15]);
CONSTANTS[22] =  CONSTANTS[21]*CONSTANTS[14];
CONSTANTS[23] = pow(10.0000, - CONSTANTS[18]);
CONSTANTS[24] = CONSTANTS[23]/CONSTANTS[17];
CONSTANTS[25] =  CONSTANTS[11]*pow(CONSTANTS[20]/CONSTANTS[6], CONSTANTS[12])*pow(CONSTANTS[24]/CONSTANTS[7], - CONSTANTS[8])*pow(CONSTANTS[22]/CONSTANTS[10], - CONSTANTS[9]);
CONSTANTS[26] = ( CONSTANTS[25]*( CONSTANTS[2]*CONSTANTS[22]*(1.00000+CONSTANTS[3]/CONSTANTS[24])+(1.00000+CONSTANTS[24]/CONSTANTS[5])))/( CONSTANTS[0]*CONSTANTS[22]*(1.00000+CONSTANTS[1]/CONSTANTS[24])+(1.00000+CONSTANTS[24]/CONSTANTS[4]));
CONSTANTS[27] = ( CONSTANTS[0]*(1.00000+CONSTANTS[1]/CONSTANTS[24])+ CONSTANTS[2]*CONSTANTS[25]*(1.00000+CONSTANTS[3]/CONSTANTS[24])*CONSTANTS[20])/(1.00000+CONSTANTS[24]/CONSTANTS[4]+ CONSTANTS[25]*(1.00000+CONSTANTS[24]/CONSTANTS[5])*CONSTANTS[20]);
CONSTANTS[28] = ( CONSTANTS[26]*CONSTANTS[20])/(1.00000+ CONSTANTS[26]*CONSTANTS[20]);
CONSTANTS[29] = ( CONSTANTS[27]*CONSTANTS[22])/(1.00000+ CONSTANTS[27]*CONSTANTS[22]);
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
getStateInformation(double* SI)
{
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}