Generated Code

/*
   There are a total of 22 entries in the algebraic variable array.
   There are a total of 10 entries in each of the rate and state variable arrays.
   There are a total of 43 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_mL).
 * CONSTANTS[1] is R_av in component heart_parameters (kPa_second_per_mL).
 * CONSTANTS[2] is R_tc in component heart_parameters (kPa_second_per_mL).
 * CONSTANTS[3] is R_pv in component heart_parameters (kPa_second_per_mL).
 * CONSTANTS[4] is R_pul in component heart_parameters (kPa_second_per_mL).
 * CONSTANTS[5] is R_sys in component heart_parameters (kPa_second_per_mL).
 * CONSTANTS[6] is L_tc in component heart_parameters (kPa_second2_per_mL).
 * CONSTANTS[7] is L_pv in component heart_parameters (kPa_second2_per_mL).
 * CONSTANTS[8] is L_mt in component heart_parameters (kPa_second2_per_mL).
 * CONSTANTS[9] is L_av in component heart_parameters (kPa_second2_per_mL).
 * CONSTANTS[10] is V_tot in component heart_parameters (mL).
 * CONSTANTS[11] is P_th in component heart_parameters (kPa).
 * ALGEBRAIC[1] is e_t in component driver_function (dimensionless).
 * CONSTANTS[12] is A in component driver_function (dimensionless).
 * CONSTANTS[13] is B in component driver_function (per_second2).
 * CONSTANTS[14] is C in component driver_function (second).
 * ALGEBRAIC[0] is tau in component driver_function (second).
 * CONSTANTS[15] is period in component driver_function (second).
 * ALGEBRAIC[2] is V_pcd in component pericardium (mL).
 * 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 (mL).
 * STATES[1] is V_rv in component right_ventricle (mL).
 * CONSTANTS[16] is P_0_pcd in component pericardium (kPa).
 * CONSTANTS[17] is V_0_pcd in component pericardium (mL).
 * CONSTANTS[18] is lambda_pcd in component pericardium (per_mL).
 * ALGEBRAIC[6] is V_lvf in component left_ventricle (mL).
 * ALGEBRAIC[9] is P_lvf in component left_ventricle (kPa).
 * ALGEBRAIC[10] is P_lv in component left_ventricle (kPa).
 * ALGEBRAIC[5] is V_spt in component septum (mL).
 * ALGEBRAIC[7] is P_es_lvf in component lvf_calculator (kPa).
 * ALGEBRAIC[8] is P_ed_lvf in component lvf_calculator (kPa).
 * ALGEBRAIC[18] is P_pu in component pulmonary_vein (kPa).
 * ALGEBRAIC[17] is P_ao in component aorta (kPa).
 * CONSTANTS[19] is E_es_lvf in component lvf_calculator (kPa_per_mL).
 * CONSTANTS[20] is lambda_lvf in component lvf_calculator (per_mL).
 * CONSTANTS[21] is P_0_lvf in component lvf_calculator (kPa).
 * STATES[2] is Q_mt in component flow (mL_per_second).
 * STATES[3] is Q_av in component flow (mL_per_second).
 * CONSTANTS[22] is V_d_lvf in component lvf_calculator (mL).
 * CONSTANTS[23] is V_0_lvf in component lvf_calculator (mL).
 * ALGEBRAIC[11] is V_rvf in component right_ventricle (mL).
 * ALGEBRAIC[14] is P_rvf in component right_ventricle (kPa).
 * ALGEBRAIC[15] is P_rv in component right_ventricle (kPa).
 * ALGEBRAIC[12] is P_es_rvf in component rvf_calculator (kPa).
 * ALGEBRAIC[13] is P_ed_rvf in component rvf_calculator (kPa).
 * ALGEBRAIC[16] is P_pa in component pulmonary_artery (kPa).
 * ALGEBRAIC[19] is P_vc in component vena_cava (kPa).
 * CONSTANTS[24] is E_es_rvf in component rvf_calculator (kPa_per_mL).
 * CONSTANTS[25] is lambda_rvf in component rvf_calculator (per_mL).
 * CONSTANTS[26] is P_0_rvf in component rvf_calculator (kPa).
 * STATES[4] is Q_tc in component flow (mL_per_second).
 * STATES[5] is Q_pv in component flow (mL_per_second).
 * CONSTANTS[27] is V_d_rvf in component rvf_calculator (mL).
 * CONSTANTS[28] is V_0_rvf in component rvf_calculator (mL).
 * CONSTANTS[29] is E_es_spt in component septum (kPa_per_mL).
 * CONSTANTS[30] is V_d_spt in component septum (mL).
 * CONSTANTS[31] is P_0_spt in component septum (kPa).
 * CONSTANTS[32] is lambda_spt in component septum (per_mL).
 * CONSTANTS[33] is V_0_spt in component septum (mL).
 * CONSTANTS[34] is one in component septum (dimensionless).
 * CONSTANTS[35] is E_es_pa in component pulmonary_artery (kPa_per_mL).
 * STATES[6] is V_pa in component pulmonary_artery (mL).
 * CONSTANTS[36] is V_d_pa in component pulmonary_artery (mL).
 * ALGEBRAIC[20] is Q_pul in component flow (mL_per_second).
 * CONSTANTS[37] is E_es_pu in component pulmonary_vein (kPa_per_mL).
 * STATES[7] is V_pu in component pulmonary_vein (mL).
 * CONSTANTS[38] is V_d_pu in component pulmonary_vein (mL).
 * CONSTANTS[39] is E_es_ao in component aorta (kPa_per_mL).
 * STATES[8] is V_ao in component aorta (mL).
 * CONSTANTS[40] is V_d_ao in component aorta (mL).
 * ALGEBRAIC[21] is Q_sys in component flow (mL_per_second).
 * CONSTANTS[41] is E_es_vc in component vena_cava (kPa_per_mL).
 * STATES[9] is V_vc in component vena_cava (mL).
 * CONSTANTS[42] is V_d_vc in component vena_cava (mL).
 * RATES[0] is d/dt V_lv in component left_ventricle (mL).
 * RATES[1] is d/dt V_rv in component right_ventricle (mL).
 * RATES[6] is d/dt V_pa in component pulmonary_artery (mL).
 * RATES[7] is d/dt V_pu in component pulmonary_vein (mL).
 * RATES[8] is d/dt V_ao in component aorta (mL).
 * RATES[9] is d/dt V_vc in component vena_cava (mL).
 * RATES[2] is d/dt Q_mt in component flow (mL_per_second).
 * RATES[3] is d/dt Q_av in component flow (mL_per_second).
 * RATES[4] is d/dt Q_tc in component flow (mL_per_second).
 * RATES[5] is d/dt Q_pv in component flow (mL_per_second).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 0.0158;
CONSTANTS[1] = 0.0180;
CONSTANTS[2] = 0.0237;
CONSTANTS[3] = 0.0055;
CONSTANTS[4] = 0.1552;
CONSTANTS[5] = 1.0889;
CONSTANTS[6] = 8.0093e-5;
CONSTANTS[7] = 1.4868e-4;
CONSTANTS[8] = 7.6968e-5;
CONSTANTS[9] = 1.2189e-4;
CONSTANTS[10] = 5.5;
CONSTANTS[11] = -4;
CONSTANTS[12] = 1;
CONSTANTS[13] = 80;
CONSTANTS[14] = 0.375;
CONSTANTS[15] = 0.75;
STATES[0] = 94.6812;
STATES[1] = 90.7302;
CONSTANTS[16] = 0.5003;
CONSTANTS[17] = 200;
CONSTANTS[18] = 0.03;
CONSTANTS[19] = 2.8798;
CONSTANTS[20] = 0.033;
CONSTANTS[21] = 0.1203;
STATES[2] = 245.5813;
STATES[3] = 0;
CONSTANTS[22] = 0;
CONSTANTS[23] = 0;
CONSTANTS[24] = 0.585;
CONSTANTS[25] = 0.023;
CONSTANTS[26] = 0.2157;
STATES[4] = 190.0661;
STATES[5] = 0;
CONSTANTS[27] = 0;
CONSTANTS[28] = 0;
CONSTANTS[29] = 48.754;
CONSTANTS[30] = 2;
CONSTANTS[31] = 1.1101;
CONSTANTS[32] = 0.435;
CONSTANTS[33] = 2;
CONSTANTS[34] = 1;
CONSTANTS[35] = 0.369;
STATES[6] = 43.0123;
CONSTANTS[36] = 0;
CONSTANTS[37] = 0.0073;
STATES[7] = 808.4579;
CONSTANTS[38] = 0;
CONSTANTS[39] = 0.6913;
STATES[8] = 133.3381;
CONSTANTS[40] = 0;
CONSTANTS[41] = 0.0059;
STATES[9] = 329.7803;
CONSTANTS[42] = 0;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[0] = (STATES[2]<0.00000&&STATES[3]<0.00000 ? 0.00000 : STATES[2]<0.00000 ? - STATES[3] : STATES[3]<0.00000 ? STATES[2] : STATES[2] - STATES[3]);
RATES[1] = (STATES[4]<0.00000&&STATES[5]<0.00000 ? 0.00000 : STATES[4]<0.00000 ? - STATES[5] : STATES[5]<0.00000 ? STATES[4] : STATES[4] - STATES[5]);
ALGEBRAIC[2] = STATES[0]+STATES[1];
ALGEBRAIC[3] =  CONSTANTS[16]*(exp( CONSTANTS[18]*(ALGEBRAIC[2] - CONSTANTS[17])) - 1.00000);
ALGEBRAIC[4] = ALGEBRAIC[3]+CONSTANTS[11];
ALGEBRAIC[0] = (VOI<=CONSTANTS[15] ? VOI : VOI<= CONSTANTS[15]*2.00000 ? VOI - CONSTANTS[15] : VOI<= CONSTANTS[15]*3.00000 ? VOI -  CONSTANTS[15]*2.00000 : VOI<= CONSTANTS[15]*4.00000 ? VOI -  CONSTANTS[15]*3.00000 : VOI<= CONSTANTS[15]*5.00000 ? VOI -  CONSTANTS[15]*4.00000 : VOI<= CONSTANTS[15]*6.00000 ? VOI -  CONSTANTS[15]*5.00000 : VOI<= CONSTANTS[15]*7.00000 ? VOI -  CONSTANTS[15]*6.00000 : VOI<= CONSTANTS[15]*8.00000 ? VOI -  CONSTANTS[15]*7.00000 : VOI<= CONSTANTS[15]*9.00000 ? VOI -  CONSTANTS[15]*8.00000 : VOI<= CONSTANTS[15]*10.0000 ? VOI -  CONSTANTS[15]*9.00000 : VOI<= CONSTANTS[15]*11.0000 ? VOI -  CONSTANTS[15]*10.0000 : VOI<= CONSTANTS[15]*12.0000 ? VOI -  CONSTANTS[15]*11.0000 : VOI<= CONSTANTS[15]*13.0000 ? VOI -  CONSTANTS[15]*12.0000 : 0.0/0.0);
ALGEBRAIC[1] =  CONSTANTS[12]*exp( - CONSTANTS[13]*pow(ALGEBRAIC[0] - CONSTANTS[14], 2.00000));
rootfind_0(VOI, CONSTANTS, RATES, STATES, ALGEBRAIC, pret);
ALGEBRAIC[6] = STATES[0] - ALGEBRAIC[5];
ALGEBRAIC[7] =  CONSTANTS[19]*(ALGEBRAIC[6] - CONSTANTS[22]);
ALGEBRAIC[8] =  CONSTANTS[21]*(exp( CONSTANTS[20]*(ALGEBRAIC[6] - CONSTANTS[23])) - 1.00000);
ALGEBRAIC[9] =  ALGEBRAIC[1]*ALGEBRAIC[7]+ (1.00000 - ALGEBRAIC[1])*ALGEBRAIC[8];
ALGEBRAIC[10] = ALGEBRAIC[9]+ALGEBRAIC[4];
ALGEBRAIC[17] =  CONSTANTS[39]*(STATES[8] - CONSTANTS[40]);
RATES[3] = (ALGEBRAIC[10] - ALGEBRAIC[17]<0.00000&&STATES[3]<0.00000 ? 0.00000 : ((ALGEBRAIC[10] - ALGEBRAIC[17]) -  STATES[3]*CONSTANTS[1])/CONSTANTS[9]);
ALGEBRAIC[11] = STATES[1]+ALGEBRAIC[5];
ALGEBRAIC[12] =  CONSTANTS[24]*(ALGEBRAIC[11] - CONSTANTS[27]);
ALGEBRAIC[13] =  CONSTANTS[26]*(exp( CONSTANTS[25]*(ALGEBRAIC[11] - CONSTANTS[28])) - 1.00000);
ALGEBRAIC[14] =  ALGEBRAIC[1]*ALGEBRAIC[12]+ (1.00000 - ALGEBRAIC[1])*ALGEBRAIC[13];
ALGEBRAIC[15] = ALGEBRAIC[14]+ALGEBRAIC[4];
ALGEBRAIC[16] =  CONSTANTS[35]*(STATES[6] - CONSTANTS[36])+CONSTANTS[11];
RATES[5] = (ALGEBRAIC[15] - ALGEBRAIC[16]<0.00000&&STATES[5]<0.00000 ? 0.00000 : ((ALGEBRAIC[15] - ALGEBRAIC[16]) -  STATES[5]*CONSTANTS[3])/CONSTANTS[7]);
ALGEBRAIC[18] =  CONSTANTS[37]*(STATES[7] - CONSTANTS[38])+CONSTANTS[11];
RATES[2] = (ALGEBRAIC[18] - ALGEBRAIC[10]<0.00000&&STATES[2]<0.00000 ? 0.00000 : ((ALGEBRAIC[18] - ALGEBRAIC[10]) -  STATES[2]*CONSTANTS[0])/CONSTANTS[8]);
ALGEBRAIC[19] =  CONSTANTS[41]*(STATES[9] - CONSTANTS[42]);
RATES[4] = (ALGEBRAIC[19] - ALGEBRAIC[15]<0.00000&&STATES[4]<0.00000 ? 0.00000 : ((ALGEBRAIC[19] - ALGEBRAIC[15]) -  STATES[4]*CONSTANTS[2])/CONSTANTS[6]);
ALGEBRAIC[20] = (ALGEBRAIC[16] - ALGEBRAIC[18])/CONSTANTS[4];
RATES[6] = (STATES[5]<0.00000 ? - ALGEBRAIC[20] : STATES[5] - ALGEBRAIC[20]);
RATES[7] = (STATES[2]<0.00000 ? ALGEBRAIC[20] : ALGEBRAIC[20] - STATES[2]);
ALGEBRAIC[21] = (ALGEBRAIC[17] - ALGEBRAIC[19])/CONSTANTS[5];
RATES[8] = (STATES[3]<0.00000 ? - ALGEBRAIC[21] : STATES[3] - ALGEBRAIC[21]);
RATES[9] = (STATES[4]<0.00000 ? ALGEBRAIC[21] : ALGEBRAIC[21] - STATES[4]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = STATES[0]+STATES[1];
ALGEBRAIC[3] =  CONSTANTS[16]*(exp( CONSTANTS[18]*(ALGEBRAIC[2] - CONSTANTS[17])) - 1.00000);
ALGEBRAIC[4] = ALGEBRAIC[3]+CONSTANTS[11];
ALGEBRAIC[0] = (VOI<=CONSTANTS[15] ? VOI : VOI<= CONSTANTS[15]*2.00000 ? VOI - CONSTANTS[15] : VOI<= CONSTANTS[15]*3.00000 ? VOI -  CONSTANTS[15]*2.00000 : VOI<= CONSTANTS[15]*4.00000 ? VOI -  CONSTANTS[15]*3.00000 : VOI<= CONSTANTS[15]*5.00000 ? VOI -  CONSTANTS[15]*4.00000 : VOI<= CONSTANTS[15]*6.00000 ? VOI -  CONSTANTS[15]*5.00000 : VOI<= CONSTANTS[15]*7.00000 ? VOI -  CONSTANTS[15]*6.00000 : VOI<= CONSTANTS[15]*8.00000 ? VOI -  CONSTANTS[15]*7.00000 : VOI<= CONSTANTS[15]*9.00000 ? VOI -  CONSTANTS[15]*8.00000 : VOI<= CONSTANTS[15]*10.0000 ? VOI -  CONSTANTS[15]*9.00000 : VOI<= CONSTANTS[15]*11.0000 ? VOI -  CONSTANTS[15]*10.0000 : VOI<= CONSTANTS[15]*12.0000 ? VOI -  CONSTANTS[15]*11.0000 : VOI<= CONSTANTS[15]*13.0000 ? VOI -  CONSTANTS[15]*12.0000 : 0.0/0.0);
ALGEBRAIC[1] =  CONSTANTS[12]*exp( - CONSTANTS[13]*pow(ALGEBRAIC[0] - CONSTANTS[14], 2.00000));
ALGEBRAIC[6] = STATES[0] - ALGEBRAIC[5];
ALGEBRAIC[7] =  CONSTANTS[19]*(ALGEBRAIC[6] - CONSTANTS[22]);
ALGEBRAIC[8] =  CONSTANTS[21]*(exp( CONSTANTS[20]*(ALGEBRAIC[6] - CONSTANTS[23])) - 1.00000);
ALGEBRAIC[9] =  ALGEBRAIC[1]*ALGEBRAIC[7]+ (1.00000 - ALGEBRAIC[1])*ALGEBRAIC[8];
ALGEBRAIC[10] = ALGEBRAIC[9]+ALGEBRAIC[4];
ALGEBRAIC[17] =  CONSTANTS[39]*(STATES[8] - CONSTANTS[40]);
ALGEBRAIC[11] = STATES[1]+ALGEBRAIC[5];
ALGEBRAIC[12] =  CONSTANTS[24]*(ALGEBRAIC[11] - CONSTANTS[27]);
ALGEBRAIC[13] =  CONSTANTS[26]*(exp( CONSTANTS[25]*(ALGEBRAIC[11] - CONSTANTS[28])) - 1.00000);
ALGEBRAIC[14] =  ALGEBRAIC[1]*ALGEBRAIC[12]+ (1.00000 - ALGEBRAIC[1])*ALGEBRAIC[13];
ALGEBRAIC[15] = ALGEBRAIC[14]+ALGEBRAIC[4];
ALGEBRAIC[16] =  CONSTANTS[35]*(STATES[6] - CONSTANTS[36])+CONSTANTS[11];
ALGEBRAIC[18] =  CONSTANTS[37]*(STATES[7] - CONSTANTS[38])+CONSTANTS[11];
ALGEBRAIC[19] =  CONSTANTS[41]*(STATES[9] - CONSTANTS[42]);
ALGEBRAIC[20] = (ALGEBRAIC[16] - ALGEBRAIC[18])/CONSTANTS[4];
ALGEBRAIC[21] = (ALGEBRAIC[17] - ALGEBRAIC[19])/CONSTANTS[5];
}

void objfunc_0(double *p, double *hx, int m, int n, void *adata)
{
  struct rootfind_info* rfi = (struct rootfind_info*)adata;
#define VOI rfi->aVOI
#define CONSTANTS rfi->aCONSTANTS
#define RATES rfi->aRATES
#define STATES rfi->aSTATES
#define ALGEBRAIC rfi->aALGEBRAIC
#define pret rfi->aPRET
  ALGEBRAIC[5] = *p;
  *hx = (0.00000) - (((( ALGEBRAIC[1]*CONSTANTS[29]*(ALGEBRAIC[5] - CONSTANTS[30])+ (CONSTANTS[34] - ALGEBRAIC[1])*CONSTANTS[31]*(exp( CONSTANTS[32]*(ALGEBRAIC[5] - CONSTANTS[33])) - CONSTANTS[34])) -  ALGEBRAIC[1]*CONSTANTS[19]*(STATES[0] - ALGEBRAIC[5])) -  (1.00000 - ALGEBRAIC[1])*CONSTANTS[21]*(exp( CONSTANTS[20]*(STATES[0] - ALGEBRAIC[5])) - 1.00000))+ ALGEBRAIC[1]*CONSTANTS[24]*(STATES[1]+ALGEBRAIC[5])+ (1.00000 - ALGEBRAIC[1])*CONSTANTS[26]*(exp( CONSTANTS[25]*(STATES[1]+ALGEBRAIC[5])) - 1.00000));
#undef VOI
#undef CONSTANTS
#undef RATES
#undef STATES
#undef ALGEBRAIC
#undef pret
}
void rootfind_0(double VOI, double* CONSTANTS, double* RATES, 
double* STATES, double* ALGEBRAIC, int* pret)
{
  static double val = 0.1;
  double bp, work[LM_DIF_WORKSZ(1, 1)];
  struct rootfind_info rfi;
  rfi.aVOI = VOI;
  rfi.aCONSTANTS = CONSTANTS;
  rfi.aRATES = RATES;
  rfi.aSTATES = STATES;
  rfi.aALGEBRAIC = ALGEBRAIC;
  rfi.aPRET = pret;
  do_levmar(objfunc_0, &val, &bp, work, pret, 1, &rfi);
  ALGEBRAIC[5] = val;
}