/* 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; }