Generated Code

/*
   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[8] 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[6] 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[9] 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[7] is PPO in component lung_lymphatic_protein_flow (gram_per_minute).
 * ALGEBRAIC[5] 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).
 * There are a total of 2 condition variables.
 */
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];
RATES[0] = 0.1001;
RATES[1] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - ALGEBRAIC[10];
resid[1] = RATES[1] - ALGEBRAIC[12];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (CONDVAR[0]<0.00000 ? 0.00100000 : STATES[0]);
ALGEBRAIC[2] = (CONDVAR[1]<0.00000 ? 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[8] =  (((CONSTANTS[7] - ALGEBRAIC[1])+ALGEBRAIC[4]) - CONSTANTS[0])*CONSTANTS[6];
ALGEBRAIC[6] =  (ALGEBRAIC[1]+11.0000)*0.000300000;
ALGEBRAIC[9] = ALGEBRAIC[8] - ALGEBRAIC[6];
ALGEBRAIC[10] = ALGEBRAIC[9];
ALGEBRAIC[7] =  ALGEBRAIC[6]*ALGEBRAIC[3];
ALGEBRAIC[5] =  (CONSTANTS[3] - ALGEBRAIC[3])*0.000225000;
ALGEBRAIC[11] = ALGEBRAIC[5] - ALGEBRAIC[7];
ALGEBRAIC[12] = ALGEBRAIC[11];
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = STATES[0] - 0.00100000;
CONDVAR[1] = STATES[1] - 0.0250000;
}