/*
   There are a total of 13 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 8 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * CONSTANTS[0] is PPC in component pulmonary_fluid_dynamics (mmHg).
 * CONSTANTS[1] is PPA in component pulmonary_fluid_dynamics (mmHg).
 * CONSTANTS[2] is PLA in component pulmonary_fluid_dynamics (mmHg).
 * CONSTANTS[3] is CPP in component pulmonary_fluid_dynamics (gram_per_L).
 * CONSTANTS[4] is RPV in component pulmonary_fluid_dynamics (mmHg_minute_per_L).
 * CONSTANTS[5] is RPA in component pulmonary_fluid_dynamics (mmHg_minute_per_L).
 * CONSTANTS[7] is PCP in component pulmonary_capillary_pressure (mmHg).
 * ALGEBRAIC[4] is POS in component colloid_osmotic_pressure_of_pulmonary_interstitium (mmHg).
 * ALGEBRAIC[1] is PPI in component pulmonary_interstitial_fluid_pressure (mmHg).
 * ALGEBRAIC[5] is PFI in component fluid_filtration_into_pulmonary_interstitium (L_per_minute).
 * CONSTANTS[6] is CPF in component parameter_values (L_per_minute_per_mmHg).
 * ALGEBRAIC[7] is PLF in component lung_lymphatic_protein_flow (L_per_minute).
 * ALGEBRAIC[10] is DFP in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
 * ALGEBRAIC[0] is VPF in component pulmonary_interstitial_free_fluid_volume (litre).
 * ALGEBRAIC[8] is DFZ in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
 * STATES[0] is VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
 * ALGEBRAIC[9] is PPO in component lung_lymphatic_protein_flow (gram_per_minute).
 * ALGEBRAIC[6] is PPN in component protein_leakage_into_pulmonary_interstitium (gram_per_minute).
 * ALGEBRAIC[12] is PPD in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
 * ALGEBRAIC[3] is CPN in component concentration_of_protein_in_pulmonary_interstitium (gram_per_L).
 * ALGEBRAIC[11] is PPZ in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
 * STATES[1] is PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
 * ALGEBRAIC[2] is PPR in component concentration_of_protein_in_pulmonary_interstitium (gram).
 * RATES[0] is d/dt VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
 * RATES[1] is d/dt PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 29.9941;
CONSTANTS[1] = 15.6376;
CONSTANTS[2] = 2;
CONSTANTS[3] = 71.9719;
CONSTANTS[4] = 1.55719;
CONSTANTS[5] = 1.5683;
CONSTANTS[6] = 0.0003;
STATES[0] = 0.0123238;
STATES[1] = 0.419998;
CONSTANTS[7] = ( (CONSTANTS[1] - CONSTANTS[2])*CONSTANTS[4])/(CONSTANTS[4]+CONSTANTS[5])+CONSTANTS[2];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (STATES[0]<0.00100000 ? 0.00100000 : STATES[0]);
ALGEBRAIC[2] = (STATES[1]<0.0250000 ? 0.0250000 : STATES[1]);
ALGEBRAIC[3] = ALGEBRAIC[2]/ALGEBRAIC[0];
ALGEBRAIC[4] =  ALGEBRAIC[3]*0.400000;
ALGEBRAIC[1] = 2.00000 - 0.150000/ALGEBRAIC[0];
ALGEBRAIC[5] =  (((CONSTANTS[7] - ALGEBRAIC[1])+ALGEBRAIC[4]) - CONSTANTS[0])*CONSTANTS[6];
ALGEBRAIC[7] =  (ALGEBRAIC[1]+11.0000)*0.000300000;
ALGEBRAIC[8] = ALGEBRAIC[5] - ALGEBRAIC[7];
ALGEBRAIC[10] = ALGEBRAIC[8];
RATES[0] = ALGEBRAIC[10];
ALGEBRAIC[9] =  ALGEBRAIC[7]*ALGEBRAIC[3];
ALGEBRAIC[6] =  (CONSTANTS[3] - ALGEBRAIC[3])*0.000225000;
ALGEBRAIC[11] = ALGEBRAIC[6] - ALGEBRAIC[9];
ALGEBRAIC[12] = ALGEBRAIC[11];
RATES[1] = ALGEBRAIC[12];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (STATES[0]<0.00100000 ? 0.00100000 : STATES[0]);
ALGEBRAIC[2] = (STATES[1]<0.0250000 ? 0.0250000 : STATES[1]);
ALGEBRAIC[3] = ALGEBRAIC[2]/ALGEBRAIC[0];
ALGEBRAIC[4] =  ALGEBRAIC[3]*0.400000;
ALGEBRAIC[1] = 2.00000 - 0.150000/ALGEBRAIC[0];
ALGEBRAIC[5] =  (((CONSTANTS[7] - ALGEBRAIC[1])+ALGEBRAIC[4]) - CONSTANTS[0])*CONSTANTS[6];
ALGEBRAIC[7] =  (ALGEBRAIC[1]+11.0000)*0.000300000;
ALGEBRAIC[8] = ALGEBRAIC[5] - ALGEBRAIC[7];
ALGEBRAIC[10] = ALGEBRAIC[8];
ALGEBRAIC[9] =  ALGEBRAIC[7]*ALGEBRAIC[3];
ALGEBRAIC[6] =  (CONSTANTS[3] - ALGEBRAIC[3])*0.000225000;
ALGEBRAIC[11] = ALGEBRAIC[6] - ALGEBRAIC[9];
ALGEBRAIC[12] = ALGEBRAIC[11];
}