/*
   There are a total of 9 entries in the algebraic variable array.
   There are a total of 18 entries in each of the rate and state variable arrays.
   There are a total of 40 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * CONSTANTS[0] is R_mt in component heart_parameters (kPa_second_per_liter).
 * CONSTANTS[1] is R_av in component heart_parameters (kPa_second_per_liter).
 * CONSTANTS[2] is R_tc in component heart_parameters (kPa_second_per_liter).
 * CONSTANTS[3] is R_pv in component heart_parameters (kPa_second_per_liter).
 * CONSTANTS[4] is R_pul in component heart_parameters (kPa_second_per_liter).
 * CONSTANTS[5] is R_sys in component heart_parameters (kPa_second_per_liter).
 * CONSTANTS[6] is HR in component heart_parameters (dimensionless).
 * CONSTANTS[7] is V_tot in component heart_parameters (liter).
 * CONSTANTS[8] is P_pl in component heart_parameters (kPa).
 * ALGEBRAIC[1] is e_t in component driver_function (dimensionless).
 * CONSTANTS[9] is A in component driver_function (dimensionless).
 * CONSTANTS[10] is B in component driver_function (dimensionless).
 * CONSTANTS[11] is C in component driver_function (dimensionless).
 * ALGEBRAIC[0] is tau in component driver_function (second).
 * CONSTANTS[12] is period in component driver_function (dimensionless).
 * ALGEBRAIC[2] is V_pcd in component pericardium (liter).
 * ALGEBRAIC[3] is P_pcd in component pericardium (kPa).
 * ALGEBRAIC[4] is P_peri in component pericardium (kPa).
 * STATES[0] is V_lv in component left_ventricle (liter).
 * STATES[1] is V_rv in component right_ventricle (liter).
 * CONSTANTS[13] is P_0_pcd in component pericardium (kPa).
 * CONSTANTS[14] is V_0_pcd in component pericardium (liter).
 * CONSTANTS[15] is lambda_pcd in component pericardium (per_liter).
 * STATES[6] is V_lvf in component left_ventricle (liter).
 * STATES[7] is P_lvf in component left_ventricle (kPa).
 * STATES[8] is P_lv in component left_ventricle (kPa).
 * STATES[9] is V_spt in component septum (liter).
 * STATES[10] is P_es_lvf in component lvf_calculator (kPa).
 * STATES[11] is P_ed_lvf in component lvf_calculator (kPa).
 * ALGEBRAIC[6] is P_pu in component pulmonary_vein (kPa).
 * ALGEBRAIC[7] is P_ao in component aorta (kPa).
 * CONSTANTS[16] is E_es_lvf in component lvf_calculator (kPa_per_liter).
 * CONSTANTS[17] is V_d_lvf in component lvf_calculator (liter).
 * CONSTANTS[18] is P_0_lvf in component lvf_calculator (kPa).
 * CONSTANTS[19] is lambda_lvf in component lvf_calculator (per_liter).
 * CONSTANTS[20] is V_0_lvf in component lvf_calculator (liter).
 * STATES[12] is V_rvf in component right_ventricle (liter).
 * STATES[13] is P_rvf in component right_ventricle (kPa).
 * STATES[14] is P_rv in component right_ventricle (kPa).
 * STATES[15] is P_es_rvf in component rvf_calculator (kPa).
 * STATES[16] is P_ed_rvf in component rvf_calculator (kPa).
 * ALGEBRAIC[5] is P_pa in component pulmonary_artery (kPa).
 * ALGEBRAIC[8] is P_vc in component vena_cava (kPa).
 * CONSTANTS[21] is E_es_rvf in component rvf_calculator (kPa_per_liter).
 * CONSTANTS[22] is V_d_rvf in component rvf_calculator (liter).
 * CONSTANTS[23] is P_0_rvf in component rvf_calculator (kPa).
 * CONSTANTS[24] is lambda_rvf in component rvf_calculator (per_liter).
 * CONSTANTS[25] is V_0_rvf in component rvf_calculator (liter).
 * STATES[17] is P_sept in component septum (kPa).
 * CONSTANTS[26] is E_es_spt in component septum (kPa_per_liter).
 * CONSTANTS[27] is V_d_spt in component septum (liter).
 * CONSTANTS[28] is P_0_spt in component septum (kPa).
 * CONSTANTS[29] is lambda_spt in component septum (per_liter).
 * CONSTANTS[30] is V_0_spt in component septum (liter).
 * CONSTANTS[31] is one in component septum (dimensionless).
 * CONSTANTS[32] is E_es_pa in component pulmonary_artery (kPa_per_liter).
 * STATES[2] is V_pa in component pulmonary_artery (liter).
 * CONSTANTS[33] is V_d_pa in component pulmonary_artery (liter).
 * CONSTANTS[34] is E_es_pu in component pulmonary_vein (kPa_per_liter).
 * STATES[3] is V_pu in component pulmonary_vein (liter).
 * CONSTANTS[35] is V_d_pu in component pulmonary_vein (liter).
 * CONSTANTS[36] is E_es_ao in component aorta (kPa_per_liter).
 * STATES[4] is V_ao in component aorta (liter).
 * CONSTANTS[37] is V_d_ao in component aorta (liter).
 * CONSTANTS[38] is E_es_vc in component vena_cava (kPa_per_liter).
 * STATES[5] is V_vc in component vena_cava (liter).
 * CONSTANTS[39] is V_d_vc in component vena_cava (liter).
 * RATES[0] is d/dt V_lv in component left_ventricle (liter).
 * RATES[1] is d/dt V_rv in component right_ventricle (liter).
 * RATES[2] is d/dt V_pa in component pulmonary_artery (liter).
 * RATES[3] is d/dt V_pu in component pulmonary_vein (liter).
 * RATES[4] is d/dt V_ao in component aorta (liter).
 * RATES[5] is d/dt V_vc in component vena_cava (liter).
 * There are a total of 24 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 0.06;
CONSTANTS[1] = 1.4;
CONSTANTS[2] = 0.18;
CONSTANTS[3] = 0.48;
CONSTANTS[4] = 19;
CONSTANTS[5] = 140;
CONSTANTS[6] = 80;
CONSTANTS[7] = 5.5;
CONSTANTS[8] = -0.533289474;
CONSTANTS[9] = 1;
CONSTANTS[10] = 80;
CONSTANTS[11] = 0.27;
CONSTANTS[12] = 0.405;
STATES[0] = 0.005;
STATES[1] = 0.005;
CONSTANTS[13] = 0.067;
CONSTANTS[14] = 0.2;
CONSTANTS[15] = 30;
CONSTANTS[16] = 454;
CONSTANTS[17] = 0.005;
CONSTANTS[18] = 0.17;
CONSTANTS[19] = 15;
CONSTANTS[20] = 0.005;
CONSTANTS[21] = 87;
CONSTANTS[22] = 0.005;
CONSTANTS[23] = 0.16;
CONSTANTS[24] = 15;
CONSTANTS[25] = 0.005;
CONSTANTS[26] = 6500;
CONSTANTS[27] = 0.002;
CONSTANTS[28] = 0.148;
CONSTANTS[29] = 435;
CONSTANTS[30] = 0.002;
CONSTANTS[31] = 1;
CONSTANTS[32] = 45;
STATES[2] = 0.16;
CONSTANTS[33] = 0.16;
CONSTANTS[34] = 0.8;
STATES[3] = 0.2;
CONSTANTS[35] = 0.2;
CONSTANTS[36] = 94;
STATES[4] = 0.8;
CONSTANTS[37] = 0.8;
CONSTANTS[38] = 1.5;
STATES[5] = 2.83;
CONSTANTS[39] = 2.83;
STATES[6] = 0.1001;
STATES[7] = 0.1001;
STATES[8] = 0.1001;
STATES[9] = 0.1001;
STATES[10] = 0.1001;
STATES[11] = 0.1001;
STATES[12] = 0.1001;
STATES[13] = 0.1001;
STATES[14] = 0.1001;
STATES[15] = 0.1001;
STATES[16] = 0.1001;
STATES[17] = 0.1001;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = STATES[6] - STATES[0] - STATES[9];
resid[1] = STATES[7] -  ALGEBRAIC[1]*STATES[10]+ (1.00000 - ALGEBRAIC[1])*STATES[11];
resid[2] = STATES[8] - STATES[7]+ALGEBRAIC[4];
resid[3] = RATES[0] - (CONDVAR[0]<0.00000&&CONDVAR[1]<0.00000 ? 0.00000 : CONDVAR[2]<0.00000 ? - (STATES[8] - ALGEBRAIC[7])/CONSTANTS[1] : CONDVAR[3]<0.00000 ? (ALGEBRAIC[6] - STATES[8])/CONSTANTS[0] : (ALGEBRAIC[6] - STATES[8])/CONSTANTS[0] - (STATES[8] - ALGEBRAIC[7])/CONSTANTS[1]);
resid[4] = STATES[10] -  CONSTANTS[16]*(STATES[6] - CONSTANTS[17]);
resid[5] = STATES[11] -  CONSTANTS[18]*(exp( CONSTANTS[19]*(STATES[6] - CONSTANTS[20])) - 1.00000);
resid[6] = STATES[12] - STATES[1]+STATES[9];
resid[7] = STATES[13] -  ALGEBRAIC[1]*STATES[15]+ (1.00000 - ALGEBRAIC[1])*STATES[16];
resid[8] = STATES[14] - STATES[13]+ALGEBRAIC[4];
resid[9] = RATES[1] - (CONDVAR[4]<0.00000&&CONDVAR[5]<0.00000 ? 0.00000 : CONDVAR[6]<0.00000 ? - (STATES[14] - ALGEBRAIC[5])/CONSTANTS[3] : CONDVAR[7]<0.00000 ? (ALGEBRAIC[8] - STATES[14])/CONSTANTS[2] : (ALGEBRAIC[8] - STATES[14])/CONSTANTS[2] - (STATES[14] - ALGEBRAIC[5])/CONSTANTS[3]);
resid[10] = STATES[15] -  CONSTANTS[21]*(STATES[12] - CONSTANTS[22]);
resid[11] = STATES[16] -  CONSTANTS[23]*(exp( CONSTANTS[24]*(STATES[12] - CONSTANTS[25])) - 1.00000);
resid[12] = STATES[17] - STATES[8] - STATES[14];
resid[13] = STATES[7] - ( ALGEBRAIC[1]*CONSTANTS[26]*(STATES[9] - CONSTANTS[27])+ (CONSTANTS[31] - ALGEBRAIC[1])*CONSTANTS[28]*(exp( CONSTANTS[29]*(STATES[9] - CONSTANTS[30])) - CONSTANTS[31])) - STATES[13];
resid[14] = RATES[2] - (CONDVAR[8]<0.00000&&CONDVAR[9]<0.00000 ? 0.00000 : CONDVAR[10]<0.00000 ? - (ALGEBRAIC[5] - ALGEBRAIC[6])/CONSTANTS[4] : CONDVAR[11]<0.00000 ? (STATES[14] - ALGEBRAIC[5])/CONSTANTS[3] : (STATES[14] - ALGEBRAIC[5])/CONSTANTS[3] - (ALGEBRAIC[5] - ALGEBRAIC[6])/CONSTANTS[4]);
resid[15] = RATES[3] - (CONDVAR[12]<0.00000&&CONDVAR[13]<0.00000 ? 0.00000 : CONDVAR[14]<0.00000 ? - (ALGEBRAIC[6] - STATES[8])/CONSTANTS[0] : CONDVAR[15]<0.00000 ? (ALGEBRAIC[5] - ALGEBRAIC[6])/CONSTANTS[4] : (ALGEBRAIC[5] - ALGEBRAIC[6])/CONSTANTS[4] - (ALGEBRAIC[6] - STATES[8])/CONSTANTS[0]);
resid[16] = RATES[4] - (CONDVAR[16]<0.00000&&CONDVAR[17]<0.00000 ? 0.00000 : CONDVAR[18]<0.00000 ? - (ALGEBRAIC[7] - ALGEBRAIC[8])/CONSTANTS[5] : CONDVAR[19]<0.00000 ? (STATES[8] - ALGEBRAIC[7])/CONSTANTS[1] : (STATES[8] - ALGEBRAIC[7])/CONSTANTS[1] - (ALGEBRAIC[7] - ALGEBRAIC[8])/CONSTANTS[5]);
resid[17] = RATES[5] - (CONDVAR[20]<0.00000&&CONDVAR[21]<0.00000 ? 0.00000 : CONDVAR[22]<0.00000 ? - (ALGEBRAIC[8] - STATES[14])/CONSTANTS[2] : CONDVAR[23]<0.00000 ? (ALGEBRAIC[7] - ALGEBRAIC[8])/CONSTANTS[5] : (ALGEBRAIC[7] - ALGEBRAIC[8])/CONSTANTS[5] - (ALGEBRAIC[8] - STATES[14])/CONSTANTS[2]);
}
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] =  (int)(VOI) % (int)(CONSTANTS[12]);
ALGEBRAIC[1] =  CONSTANTS[9]*exp( - CONSTANTS[10]*pow(( ALGEBRAIC[0]*CONSTANTS[6])/60.0000 - CONSTANTS[11], 2.00000));
ALGEBRAIC[2] = STATES[0]+STATES[1];
ALGEBRAIC[3] =  CONSTANTS[13]*(exp( CONSTANTS[15]*(ALGEBRAIC[2] - CONSTANTS[14])) - 1.00000);
ALGEBRAIC[4] = ALGEBRAIC[3]+CONSTANTS[8];
ALGEBRAIC[5] =  CONSTANTS[32]*(STATES[2] - CONSTANTS[33]);
ALGEBRAIC[6] =  CONSTANTS[34]*(STATES[3] - CONSTANTS[35]);
ALGEBRAIC[7] =  CONSTANTS[36]*(STATES[4] - CONSTANTS[37]);
ALGEBRAIC[8] =  CONSTANTS[38]*(STATES[5] - CONSTANTS[39]);
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[6] = 0.0;
SI[7] = 0.0;
SI[8] = 0.0;
SI[9] = 0.0;
SI[10] = 0.0;
SI[11] = 0.0;
SI[12] = 0.0;
SI[13] = 0.0;
SI[14] = 0.0;
SI[15] = 0.0;
SI[16] = 0.0;
SI[17] = 0.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = (ALGEBRAIC[6] - STATES[8]) - 0.00000;
CONDVAR[1] = (STATES[8] - ALGEBRAIC[7]) - 0.00000;
CONDVAR[2] = (ALGEBRAIC[6] - STATES[8]) - 0.00000;
CONDVAR[3] = (STATES[8] - ALGEBRAIC[7]) - 0.00000;
CONDVAR[4] = (ALGEBRAIC[8] - STATES[14]) - 0.00000;
CONDVAR[5] = (STATES[14] - ALGEBRAIC[5]) - 0.00000;
CONDVAR[6] = (ALGEBRAIC[8] - STATES[14]) - 0.00000;
CONDVAR[7] = (STATES[14] - ALGEBRAIC[5]) - 0.00000;
CONDVAR[8] = (STATES[14] - ALGEBRAIC[5]) - 0.00000;
CONDVAR[9] = (ALGEBRAIC[5] - ALGEBRAIC[6]) - 0.00000;
CONDVAR[10] = (STATES[14] - ALGEBRAIC[5]) - 0.00000;
CONDVAR[11] = (ALGEBRAIC[5] - ALGEBRAIC[6]) - 0.00000;
CONDVAR[12] = (ALGEBRAIC[5] - ALGEBRAIC[6]) - 0.00000;
CONDVAR[13] = (ALGEBRAIC[6] - STATES[8]) - 0.00000;
CONDVAR[14] = (ALGEBRAIC[5] - ALGEBRAIC[6]) - 0.00000;
CONDVAR[15] = (ALGEBRAIC[6] - STATES[8]) - 0.00000;
CONDVAR[16] = (STATES[8] - ALGEBRAIC[7]) - 0.00000;
CONDVAR[17] = (ALGEBRAIC[7] - ALGEBRAIC[8]) - 0.00000;
CONDVAR[18] = (STATES[8] - ALGEBRAIC[7]) - 0.00000;
CONDVAR[19] = (ALGEBRAIC[7] - ALGEBRAIC[8]) - 0.00000;
CONDVAR[20] = (ALGEBRAIC[7] - ALGEBRAIC[8]) - 0.00000;
CONDVAR[21] = (ALGEBRAIC[8] - STATES[14]) - 0.00000;
CONDVAR[22] = (ALGEBRAIC[7] - ALGEBRAIC[8]) - 0.00000;
CONDVAR[23] = (ALGEBRAIC[8] - STATES[14]) - 0.00000;
}