Generated Code
The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
There are a total of 93 entries in the algebraic variable array.
There are a total of 33 entries in each of the rate and state variable arrays.
There are a total of 139 entries in the constant variable array.
*/
/*
* VOI is t in component main (second).
* CONSTANTS[0] is rho in component main (Js2_per_m5).
* CONSTANTS[1] is g in component main (m_per_s2).
* CONSTANTS[2] is theta_deg in component main (dim).
* CONSTANTS[138] is theta_rad in component main (radian).
* ALGEBRAIC[0] is z in component main (dim).
* CONSTANTS[3] is T in component main (second).
* ALGEBRAIC[1] is mt in component main (dim).
* ALGEBRAIC[2] is mta in component main (dim).
* CONSTANTS[4] is delay in component main (dim).
* CONSTANTS[5] is To_pk in component main (dim).
* ALGEBRAIC[3] is E_B_LV in component main (kPa_per_L).
* CONSTANTS[6] is E_B_LVmax in component main (kPa_per_L).
* CONSTANTS[7] is E_B_LVmin in component main (kPa_per_L).
* STATES[0] is q_B_LV in component main (litre).
* CONSTANTS[8] is q_B_LV_us in component main (litre).
* ALGEBRAIC[4] is u_B_LV in component main (kPa).
* ALGEBRAIC[5] is E_B_RV in component main (kPa_per_L).
* CONSTANTS[9] is E_B_RVmax in component main (kPa_per_L).
* CONSTANTS[10] is E_B_RVmin in component main (kPa_per_L).
* STATES[1] is q_B_RV in component main (litre).
* CONSTANTS[11] is q_B_RV_us in component main (litre).
* ALGEBRAIC[6] is u_B_RV in component main (kPa).
* CONSTANTS[12] is E_B_LA in component main (kPa_per_L).
* STATES[2] is q_B_LA in component main (litre).
* CONSTANTS[13] is q_B_LA_us in component main (litre).
* ALGEBRAIC[7] is u_B_LA in component main (kPa).
* CONSTANTS[14] is E_B_RA in component main (kPa_per_L).
* STATES[3] is q_B_RA in component main (litre).
* CONSTANTS[15] is q_B_RA_us in component main (litre).
* ALGEBRAIC[8] is u_B_RA in component main (kPa).
* CONSTANTS[16] is R_B_AV in component main (kPa_s_per_L).
* ALGEBRAIC[82] is v_B_AV in component main (L_per_s).
* CONSTANTS[17] is R_B_PV in component main (kPa_s_per_L).
* ALGEBRAIC[83] is v_B_PV in component main (L_per_s).
* CONSTANTS[18] is R_B_TV in component main (kPa_s_per_L).
* ALGEBRAIC[9] is v_B_TV in component main (L_per_s).
* CONSTANTS[19] is R_B_MV in component main (kPa_s_per_L).
* ALGEBRAIC[10] is v_B_MV in component main (L_per_s).
* STATES[4] is q_B_PA in component main (litre).
* CONSTANTS[20] is q_B_PA_us in component main (litre).
* ALGEBRAIC[11] is u_B_PA in component main (kPa).
* CONSTANTS[21] is E_B_PA in component main (kPa_per_L).
* STATES[5] is q_B_lung in component main (litre).
* CONSTANTS[22] is q_B_lung_us in component main (litre).
* ALGEBRAIC[12] is u_B_lung in component main (kPa).
* CONSTANTS[23] is E_B_lung in component main (kPa_per_L).
* CONSTANTS[24] is z_lung in component main (meter).
* ALGEBRAIC[13] is v_B_lung1 in component main (L_per_s).
* ALGEBRAIC[84] is v_B_lung2 in component main (L_per_s).
* CONSTANTS[25] is R_B_lung1 in component main (kPa_s_per_L).
* CONSTANTS[26] is R_B_lung2 in component main (kPa_s_per_L).
* CONSTANTS[27] is E_B_pulmVein in component main (kPa_per_L).
* STATES[6] is q_B_pulmVein in component main (litre).
* CONSTANTS[28] is q_B_pulmVein_us in component main (litre).
* ALGEBRAIC[14] is u_B_pulmVein in component main (kPa).
* ALGEBRAIC[15] is v_B_pulmVein in component main (L_per_s).
* CONSTANTS[29] is R_B_pulmVein in component main (kPa_s_per_L).
* STATES[7] is q_B_brain in component main (litre).
* CONSTANTS[30] is q_B_brain_us in component main (litre).
* STATES[8] is q_B_brainVein in component main (litre).
* CONSTANTS[31] is q_B_brainVein_us in component main (litre).
* ALGEBRAIC[16] is u_B_brain in component main (kPa).
* ALGEBRAIC[17] is u_B_brainVein in component main (kPa).
* CONSTANTS[32] is E_B_brain in component main (kPa_per_L).
* CONSTANTS[33] is E_B_brainVein in component main (kPa_per_L).
* ALGEBRAIC[85] is v_B_brain1 in component main (L_per_s).
* ALGEBRAIC[18] is v_B_brain2 in component main (L_per_s).
* ALGEBRAIC[19] is v_B_brain3 in component main (L_per_s).
* CONSTANTS[34] is R_B_brain1 in component main (kPa_s_per_L).
* CONSTANTS[35] is R_B_brain2 in component main (kPa_s_per_L).
* CONSTANTS[36] is R_B_brain3 in component main (kPa_s_per_L).
* STATES[9] is q_B_AA in component main (litre).
* CONSTANTS[37] is q_B_AA_us in component main (litre).
* ALGEBRAIC[20] is u_B_AA in component main (kPa).
* CONSTANTS[38] is E_B_AA in component main (kPa_per_L).
* STATES[10] is q_B_celiac in component main (litre).
* CONSTANTS[39] is q_B_celiac_us in component main (litre).
* ALGEBRAIC[21] is u_B_celiac in component main (kPa).
* CONSTANTS[40] is E_B_celiac in component main (kPa_per_L).
* CONSTANTS[41] is z_celiac in component main (meter).
* ALGEBRAIC[22] is v_B_celiac in component main (L_per_s).
* CONSTANTS[42] is R_B_celiac in component main (kPa_s_per_L).
* STATES[11] is q_B_supMes in component main (litre).
* CONSTANTS[43] is q_B_supMes_us in component main (litre).
* ALGEBRAIC[23] is u_B_supMes in component main (kPa).
* CONSTANTS[44] is E_B_supMes in component main (kPa_per_L).
* CONSTANTS[45] is z_supMes in component main (meter).
* ALGEBRAIC[24] is v_B_supMes in component main (L_per_s).
* CONSTANTS[46] is R_B_supMes in component main (kPa_s_per_L).
* ALGEBRAIC[86] is v_B_infMes in component main (L_per_s).
* CONSTANTS[47] is R_B_infMes in component main (kPa_s_per_L).
* STATES[12] is q_B_stomach in component main (litre).
* CONSTANTS[48] is q_B_stomach_us in component main (litre).
* ALGEBRAIC[25] is u_B_stomach in component main (kPa).
* CONSTANTS[49] is E_B_stomach in component main (kPa_per_L).
* CONSTANTS[50] is z_stomach in component main (meter).
* ALGEBRAIC[26] is v_B_stomach1 in component main (L_per_s).
* ALGEBRAIC[87] is v_B_stomach2 in component main (L_per_s).
* CONSTANTS[51] is R_B_stomach1 in component main (kPa_s_per_L).
* CONSTANTS[52] is R_B_stomach2 in component main (kPa_s_per_L).
* STATES[13] is q_B_spleen in component main (litre).
* CONSTANTS[53] is q_B_spleen_us in component main (litre).
* ALGEBRAIC[27] is u_B_spleen in component main (kPa).
* CONSTANTS[54] is E_B_spleen in component main (kPa_per_L).
* CONSTANTS[55] is z_spleen in component main (meter).
* ALGEBRAIC[28] is v_B_spleen1 in component main (L_per_s).
* ALGEBRAIC[88] is v_B_spleen2 in component main (L_per_s).
* CONSTANTS[56] is R_B_spleen1 in component main (kPa_s_per_L).
* CONSTANTS[57] is R_B_spleen2 in component main (kPa_s_per_L).
* STATES[14] is q_B_pancreas in component main (litre).
* CONSTANTS[58] is q_B_pancreas_us in component main (litre).
* ALGEBRAIC[29] is u_B_pancreas in component main (kPa).
* CONSTANTS[59] is E_B_pancreas in component main (kPa_per_L).
* CONSTANTS[60] is z_pancreas in component main (meter).
* ALGEBRAIC[30] is v_B_pancreas1 in component main (L_per_s).
* ALGEBRAIC[31] is v_B_pancreas2 in component main (L_per_s).
* ALGEBRAIC[89] is v_B_pancreas3 in component main (L_per_s).
* CONSTANTS[61] is R_B_pancreas1 in component main (kPa_s_per_L).
* CONSTANTS[62] is R_B_pancreas2 in component main (kPa_s_per_L).
* CONSTANTS[63] is R_B_pancreas3 in component main (kPa_s_per_L).
* STATES[15] is q_B_intestine in component main (litre).
* CONSTANTS[64] is q_B_intestine_us in component main (litre).
* ALGEBRAIC[32] is u_B_intestine in component main (kPa).
* CONSTANTS[65] is E_B_intestine in component main (kPa_per_L).
* CONSTANTS[66] is z_ in component main (meter).
* ALGEBRAIC[33] is v_B_intestine1 in component main (L_per_s).
* ALGEBRAIC[90] is v_B_intestine2 in component main (L_per_s).
* CONSTANTS[67] is R_B_intestine1 in component main (kPa_s_per_L).
* CONSTANTS[68] is R_B_intestine2 in component main (kPa_s_per_L).
* STATES[16] is q_B_colon in component main (litre).
* CONSTANTS[69] is q_B_colon_us in component main (litre).
* ALGEBRAIC[34] is u_B_colon in component main (kPa).
* CONSTANTS[70] is E_B_colon in component main (kPa_per_L).
* CONSTANTS[71] is z_colon in component main (meter).
* ALGEBRAIC[35] is v_B_colon1 in component main (L_per_s).
* ALGEBRAIC[91] is v_B_colon2 in component main (L_per_s).
* CONSTANTS[72] is R_B_colon1 in component main (kPa_s_per_L).
* CONSTANTS[73] is R_B_colon2 in component main (kPa_s_per_L).
* STATES[17] is q_B_portalVein in component main (litre).
* CONSTANTS[74] is q_B_portalVein_us in component main (litre).
* ALGEBRAIC[36] is u_B_portalVein in component main (kPa).
* CONSTANTS[75] is E_B_portalVein in component main (kPa_per_L).
* CONSTANTS[76] is z_portalVein in component main (meter).
* ALGEBRAIC[92] is v_B_portalVein in component main (L_per_s).
* CONSTANTS[77] is R_B_portalVein in component main (kPa_s_per_L).
* STATES[18] is q_B_liver in component main (litre).
* CONSTANTS[78] is q_B_liver_us in component main (litre).
* ALGEBRAIC[37] is u_B_liver in component main (kPa).
* CONSTANTS[79] is E_B_liver in component main (kPa_per_L).
* CONSTANTS[80] is z_liver in component main (meter).
* ALGEBRAIC[38] is v_B_liver1 in component main (L_per_s).
* ALGEBRAIC[39] is v_B_liver2 in component main (L_per_s).
* CONSTANTS[81] is R_B_liver1 in component main (kPa_s_per_L).
* CONSTANTS[82] is R_B_liver2 in component main (kPa_s_per_L).
* STATES[19] is q_B_lGlomerulus in component main (litre).
* STATES[20] is q_B_rGlomerulus in component main (litre).
* STATES[21] is q_B_lPeritubular in component main (litre).
* STATES[22] is q_B_rPeritubular in component main (litre).
* CONSTANTS[83] is q_B_lGlomerulus_us in component main (litre).
* CONSTANTS[84] is q_B_rGlomerulus_us in component main (litre).
* CONSTANTS[85] is q_B_lPeritubular_us in component main (litre).
* CONSTANTS[86] is q_B_rPeritubular_us in component main (litre).
* ALGEBRAIC[40] is u_B_lGlomerulus in component main (kPa).
* ALGEBRAIC[41] is u_B_rGlomerulus in component main (kPa).
* ALGEBRAIC[42] is u_B_lPeritubular in component main (kPa).
* ALGEBRAIC[43] is u_B_rPeritubular in component main (kPa).
* CONSTANTS[87] is E_B_lGlomerulus in component main (kPa_per_L).
* CONSTANTS[88] is E_B_rGlomerulus in component main (kPa_per_L).
* CONSTANTS[89] is E_B_lPeritubular in component main (kPa_per_L).
* CONSTANTS[90] is E_B_rPeritubular in component main (kPa_per_L).
* CONSTANTS[91] is z_kidney in component main (meter).
* ALGEBRAIC[44] is v_B_lRenal1 in component main (L_per_s).
* ALGEBRAIC[45] is v_B_rRenal1 in component main (L_per_s).
* ALGEBRAIC[46] is v_B_lRenal2 in component main (L_per_s).
* ALGEBRAIC[47] is v_B_rRenal2 in component main (L_per_s).
* ALGEBRAIC[48] is v_B_lRenal3 in component main (L_per_s).
* ALGEBRAIC[49] is v_B_rRenal3 in component main (L_per_s).
* CONSTANTS[92] is R_B_lRenal1 in component main (kPa_s_per_L).
* CONSTANTS[93] is R_B_rRenal1 in component main (kPa_s_per_L).
* CONSTANTS[94] is R_B_lRenal2 in component main (kPa_s_per_L).
* CONSTANTS[95] is R_B_rRenal2 in component main (kPa_s_per_L).
* CONSTANTS[96] is R_B_lRenal3 in component main (kPa_s_per_L).
* CONSTANTS[97] is R_B_rRenal3 in component main (kPa_s_per_L).
* STATES[23] is q_B_musTorso in component main (litre).
* STATES[24] is q_B_musTorsoVein in component main (litre).
* CONSTANTS[98] is q_B_musTorso_us in component main (litre).
* CONSTANTS[99] is q_B_musTorsoVein_us in component main (litre).
* ALGEBRAIC[50] is u_B_musTorso in component main (kPa).
* ALGEBRAIC[51] is u_B_musTorsoVein in component main (kPa).
* CONSTANTS[100] is E_B_musTorso in component main (kPa_per_L).
* CONSTANTS[101] is E_B_musTorsoVein in component main (kPa_per_L).
* CONSTANTS[102] is z_musTorso in component main (meter).
* ALGEBRAIC[52] is v_B_musTorso1 in component main (L_per_s).
* ALGEBRAIC[53] is v_B_musTorso2 in component main (L_per_s).
* ALGEBRAIC[54] is v_B_musTorso3 in component main (L_per_s).
* CONSTANTS[103] is R_B_musTorso1 in component main (kPa_s_per_L).
* CONSTANTS[104] is R_B_musTorso2 in component main (kPa_s_per_L).
* CONSTANTS[105] is R_B_musTorso3 in component main (kPa_s_per_L).
* STATES[25] is q_B_lArm in component main (litre).
* STATES[26] is q_B_lArmVein in component main (litre).
* CONSTANTS[106] is q_B_lArm_us in component main (litre).
* CONSTANTS[107] is q_B_lArmVein_us in component main (litre).
* ALGEBRAIC[55] is u_B_lArm in component main (kPa).
* ALGEBRAIC[56] is u_B_lArmVein in component main (kPa).
* CONSTANTS[108] is E_B_lArm in component main (kPa_per_L).
* CONSTANTS[109] is E_B_lArmVein in component main (kPa_per_L).
* CONSTANTS[110] is z_lArm in component main (meter).
* ALGEBRAIC[57] is v_B_lArm1 in component main (L_per_s).
* ALGEBRAIC[58] is v_B_lArm2 in component main (L_per_s).
* ALGEBRAIC[59] is v_B_lArm3 in component main (L_per_s).
* CONSTANTS[111] is R_B_lArm1 in component main (kPa_s_per_L).
* CONSTANTS[112] is R_B_lArm2 in component main (kPa_s_per_L).
* CONSTANTS[113] is R_B_lArm3 in component main (kPa_s_per_L).
* STATES[27] is q_B_rArm in component main (litre).
* STATES[28] is q_B_rArmVein in component main (litre).
* CONSTANTS[114] is q_B_rArm_us in component main (litre).
* CONSTANTS[115] is q_B_rArmVein_us in component main (litre).
* ALGEBRAIC[60] is u_B_rArm in component main (kPa).
* ALGEBRAIC[61] is u_B_rArmVein in component main (kPa).
* CONSTANTS[116] is E_B_rArm in component main (kPa_per_L).
* CONSTANTS[117] is E_B_rArmVein in component main (kPa_per_L).
* CONSTANTS[118] is z_rArm in component main (meter).
* ALGEBRAIC[62] is v_B_rArm1 in component main (L_per_s).
* ALGEBRAIC[63] is v_B_rArm2 in component main (L_per_s).
* ALGEBRAIC[64] is v_B_rArm3 in component main (L_per_s).
* CONSTANTS[119] is R_B_rArm1 in component main (kPa_s_per_L).
* CONSTANTS[120] is R_B_rArm2 in component main (kPa_s_per_L).
* CONSTANTS[121] is R_B_rArm3 in component main (kPa_s_per_L).
* STATES[29] is q_B_lLeg in component main (litre).
* STATES[30] is q_B_lLegVein in component main (litre).
* CONSTANTS[122] is q_B_lLeg_us in component main (litre).
* CONSTANTS[123] is q_B_lLegVein_us in component main (litre).
* ALGEBRAIC[65] is u_B_lLeg in component main (kPa).
* ALGEBRAIC[66] is u_B_lLegVein in component main (kPa).
* CONSTANTS[124] is E_B_lLeg in component main (kPa_per_L).
* CONSTANTS[125] is E_B_lLegVein in component main (kPa_per_L).
* CONSTANTS[126] is z_lLeg in component main (meter).
* ALGEBRAIC[67] is v_B_lLeg1 in component main (L_per_s).
* ALGEBRAIC[68] is v_B_lLeg2 in component main (L_per_s).
* ALGEBRAIC[69] is v_B_lLeg3 in component main (L_per_s).
* CONSTANTS[127] is R_B_lLeg1 in component main (kPa_s_per_L).
* CONSTANTS[128] is R_B_lLeg2 in component main (kPa_s_per_L).
* CONSTANTS[129] is R_B_lLeg3 in component main (kPa_s_per_L).
* STATES[31] is q_B_rLeg in component main (litre).
* STATES[32] is q_B_rLegVein in component main (litre).
* CONSTANTS[130] is q_B_rLeg_us in component main (litre).
* CONSTANTS[131] is q_B_rLegVein_us in component main (litre).
* ALGEBRAIC[70] is u_B_rLeg in component main (kPa).
* ALGEBRAIC[71] is u_B_rLegVein in component main (kPa).
* CONSTANTS[132] is E_B_rLeg in component main (kPa_per_L).
* CONSTANTS[133] is E_B_rLegVein in component main (kPa_per_L).
* CONSTANTS[134] is z_rLeg in component main (meter).
* ALGEBRAIC[72] is v_B_rLeg1 in component main (L_per_s).
* ALGEBRAIC[73] is v_B_rLeg2 in component main (L_per_s).
* ALGEBRAIC[74] is v_B_rLeg3 in component main (L_per_s).
* CONSTANTS[135] is R_B_rLeg1 in component main (kPa_s_per_L).
* CONSTANTS[136] is R_B_rLeg2 in component main (kPa_s_per_L).
* CONSTANTS[137] is R_B_rLeg3 in component main (kPa_s_per_L).
* ALGEBRAIC[75] is q_B_heartTot in component main (litre).
* ALGEBRAIC[76] is q_B_lungTot in component main (litre).
* ALGEBRAIC[77] is q_B_renalTot in component main (litre).
* ALGEBRAIC[78] is q_B_gutTot in component main (litre).
* ALGEBRAIC[79] is q_B_brainTot in component main (litre).
* ALGEBRAIC[80] is q_B_limbTot in component main (litre).
* ALGEBRAIC[81] is q_Blood in component main (litre).
* RATES[0] is d/dt q_B_LV in component main (litre).
* RATES[1] is d/dt q_B_RV in component main (litre).
* RATES[2] is d/dt q_B_LA in component main (litre).
* RATES[3] is d/dt q_B_RA in component main (litre).
* RATES[4] is d/dt q_B_PA in component main (litre).
* RATES[5] is d/dt q_B_lung in component main (litre).
* RATES[6] is d/dt q_B_pulmVein in component main (litre).
* RATES[7] is d/dt q_B_brain in component main (litre).
* RATES[8] is d/dt q_B_brainVein in component main (litre).
* RATES[9] is d/dt q_B_AA in component main (litre).
* RATES[10] is d/dt q_B_celiac in component main (litre).
* RATES[11] is d/dt q_B_supMes in component main (litre).
* RATES[12] is d/dt q_B_stomach in component main (litre).
* RATES[13] is d/dt q_B_spleen in component main (litre).
* RATES[14] is d/dt q_B_pancreas in component main (litre).
* RATES[15] is d/dt q_B_intestine in component main (litre).
* RATES[16] is d/dt q_B_colon in component main (litre).
* RATES[17] is d/dt q_B_portalVein in component main (litre).
* RATES[18] is d/dt q_B_liver in component main (litre).
* RATES[19] is d/dt q_B_lGlomerulus in component main (litre).
* RATES[20] is d/dt q_B_rGlomerulus in component main (litre).
* RATES[21] is d/dt q_B_lPeritubular in component main (litre).
* RATES[22] is d/dt q_B_rPeritubular in component main (litre).
* RATES[23] is d/dt q_B_musTorso in component main (litre).
* RATES[24] is d/dt q_B_musTorsoVein in component main (litre).
* RATES[25] is d/dt q_B_lArm in component main (litre).
* RATES[26] is d/dt q_B_lArmVein in component main (litre).
* RATES[27] is d/dt q_B_rArm in component main (litre).
* RATES[28] is d/dt q_B_rArmVein in component main (litre).
* RATES[29] is d/dt q_B_lLeg in component main (litre).
* RATES[30] is d/dt q_B_lLegVein in component main (litre).
* RATES[31] is d/dt q_B_rLeg in component main (litre).
* RATES[32] is d/dt q_B_rLegVein in component main (litre).
* There are a total of 10 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1e3;
CONSTANTS[1] = 9.81;
CONSTANTS[2] = 90;
CONSTANTS[3] = 1;
CONSTANTS[4] = 0.1;
CONSTANTS[5] = 0.1;
CONSTANTS[6] = 320;
CONSTANTS[7] = 6.667;
STATES[0] = 1;
CONSTANTS[8] = 0.005;
CONSTANTS[9] = 70;
CONSTANTS[10] = 2;
STATES[1] = 0.2;
CONSTANTS[11] = 0.01;
CONSTANTS[12] = 10;
STATES[2] = 0.02;
CONSTANTS[13] = 0.004;
CONSTANTS[14] = 10;
STATES[3] = 0.02;
CONSTANTS[15] = 0.004;
CONSTANTS[16] = 1;
CONSTANTS[17] = 1;
CONSTANTS[18] = 1;
CONSTANTS[19] = 1;
STATES[4] = 0.3;
CONSTANTS[20] = 0.2;
CONSTANTS[21] = 200;
STATES[5] = 0.2;
CONSTANTS[22] = 0.1;
CONSTANTS[23] = 100;
CONSTANTS[24] = 0.20;
CONSTANTS[25] = 3;
CONSTANTS[26] = 10;
CONSTANTS[27] = 20;
STATES[6] = 0.02;
CONSTANTS[28] = 0.012;
CONSTANTS[29] = 1;
STATES[7] = 0.1;
CONSTANTS[30] = 0.06;
STATES[8] = 0.1;
CONSTANTS[31] = 0.06;
CONSTANTS[32] = 20;
CONSTANTS[33] = 20;
CONSTANTS[34] = 100;
CONSTANTS[35] = 600;
CONSTANTS[36] = 1;
STATES[9] = 0.2;
CONSTANTS[37] = 0.12;
CONSTANTS[38] = 80;
STATES[10] = 0.1;
CONSTANTS[39] = 0.06;
CONSTANTS[40] = 400;
CONSTANTS[41] = 0.20;
CONSTANTS[42] = 100;
STATES[11] = 0.1;
CONSTANTS[43] = 0.06;
CONSTANTS[44] = 400;
CONSTANTS[45] = 0.20;
CONSTANTS[46] = 100;
CONSTANTS[47] = 100;
STATES[12] = 0.1;
CONSTANTS[48] = 0.06;
CONSTANTS[49] = 100;
CONSTANTS[50] = 0.20;
CONSTANTS[51] = 100;
CONSTANTS[52] = 2400;
STATES[13] = 0.1;
CONSTANTS[53] = 0.06;
CONSTANTS[54] = 100;
CONSTANTS[55] = 0.20;
CONSTANTS[56] = 100;
CONSTANTS[57] = 3000;
STATES[14] = 0.1;
CONSTANTS[58] = 0.06;
CONSTANTS[59] = 100;
CONSTANTS[60] = 0.20;
CONSTANTS[61] = 100;
CONSTANTS[62] = 100;
CONSTANTS[63] = 4000;
STATES[15] = 0.1;
CONSTANTS[64] = 0.06;
CONSTANTS[65] = 100;
CONSTANTS[66] = 0.20;
CONSTANTS[67] = 100;
CONSTANTS[68] = 1200;
STATES[16] = 0.1;
CONSTANTS[69] = 0.06;
CONSTANTS[70] = 100;
CONSTANTS[71] = 0.20;
CONSTANTS[72] = 100;
CONSTANTS[73] = 2000;
STATES[17] = 0.1;
CONSTANTS[74] = 0.06;
CONSTANTS[75] = 50;
CONSTANTS[76] = 0.20;
CONSTANTS[77] = 100;
STATES[18] = 0.2;
CONSTANTS[78] = 0.06;
CONSTANTS[79] = 10;
CONSTANTS[80] = -0.50;
CONSTANTS[81] = 800;
CONSTANTS[82] = 1;
STATES[19] = 0.1;
STATES[20] = 0.1;
STATES[21] = 0.1;
STATES[22] = 0.1;
CONSTANTS[83] = 0.06;
CONSTANTS[84] = 0.06;
CONSTANTS[85] = 0.06;
CONSTANTS[86] = 0.06;
CONSTANTS[87] = 150;
CONSTANTS[88] = 150;
CONSTANTS[89] = 50;
CONSTANTS[90] = 50;
CONSTANTS[91] = -0.25;
CONSTANTS[92] = 100;
CONSTANTS[93] = 100;
CONSTANTS[94] = 2000;
CONSTANTS[95] = 2000;
CONSTANTS[96] = 1;
CONSTANTS[97] = 1;
STATES[23] = 0.2;
STATES[24] = 0.2;
CONSTANTS[98] = 0.12;
CONSTANTS[99] = 0.12;
CONSTANTS[100] = 60;
CONSTANTS[101] = 10;
CONSTANTS[102] = -0.50;
CONSTANTS[103] = 100;
CONSTANTS[104] = 1000;
CONSTANTS[105] = 1;
STATES[25] = 0.1;
STATES[26] = 0.1;
CONSTANTS[106] = 0.06;
CONSTANTS[107] = 0.06;
CONSTANTS[108] = 150;
CONSTANTS[109] = 60;
CONSTANTS[110] = -0.50;
CONSTANTS[111] = 100;
CONSTANTS[112] = 6000;
CONSTANTS[113] = 1;
STATES[27] = 0.1;
STATES[28] = 0.1;
CONSTANTS[114] = 0.06;
CONSTANTS[115] = 0.06;
CONSTANTS[116] = 150;
CONSTANTS[117] = 60;
CONSTANTS[118] = -0.50;
CONSTANTS[119] = 100;
CONSTANTS[120] = 6000;
CONSTANTS[121] = 1;
STATES[29] = 0.1;
STATES[30] = 0.1;
CONSTANTS[122] = 0.06;
CONSTANTS[123] = 0.06;
CONSTANTS[124] = 150;
CONSTANTS[125] = 60;
CONSTANTS[126] = -0.50;
CONSTANTS[127] = 100;
CONSTANTS[128] = 4000;
CONSTANTS[129] = 1;
STATES[31] = 0.1;
STATES[32] = 0.1;
CONSTANTS[130] = 0.06;
CONSTANTS[131] = 0.06;
CONSTANTS[132] = 150;
CONSTANTS[133] = 60;
CONSTANTS[134] = -0.50;
CONSTANTS[135] = 100;
CONSTANTS[136] = 4000;
CONSTANTS[137] = 1;
CONSTANTS[138] = ( CONSTANTS[2]* 3.14159265358979)/180.000;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[13] = 0.1001;
RATES[14] = 0.1001;
RATES[15] = 0.1001;
RATES[16] = 0.1001;
RATES[17] = 0.1001;
RATES[18] = 0.1001;
RATES[19] = 0.1001;
RATES[20] = 0.1001;
RATES[21] = 0.1001;
RATES[22] = 0.1001;
RATES[23] = 0.1001;
RATES[24] = 0.1001;
RATES[25] = 0.1001;
RATES[26] = 0.1001;
RATES[27] = 0.1001;
RATES[28] = 0.1001;
RATES[29] = 0.1001;
RATES[30] = 0.1001;
RATES[31] = 0.1001;
RATES[32] = 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] - ALGEBRAIC[82];
resid[1] = RATES[1] - ALGEBRAIC[9] - ALGEBRAIC[83];
resid[2] = RATES[2] - ALGEBRAIC[15] - ALGEBRAIC[10];
resid[3] = RATES[3] - (ALGEBRAIC[19]+ALGEBRAIC[39]+ALGEBRAIC[48]+ALGEBRAIC[49]+ALGEBRAIC[54]+ALGEBRAIC[59]+ALGEBRAIC[64]+ALGEBRAIC[69]+ALGEBRAIC[74]) - ALGEBRAIC[9];
resid[4] = RATES[4] - ALGEBRAIC[83] - ALGEBRAIC[13];
resid[5] = RATES[5] - ALGEBRAIC[13] - ALGEBRAIC[84];
resid[6] = RATES[6] - ALGEBRAIC[84] - ALGEBRAIC[15];
resid[7] = RATES[7] - ALGEBRAIC[85] - ALGEBRAIC[18];
resid[8] = RATES[8] - ALGEBRAIC[18] - ALGEBRAIC[19];
resid[9] = RATES[9] - ((((((((((ALGEBRAIC[82] - ALGEBRAIC[85]) - ALGEBRAIC[22]) - ALGEBRAIC[24]) - ALGEBRAIC[86]) - ALGEBRAIC[44]) - ALGEBRAIC[45]) - ALGEBRAIC[52]) - ALGEBRAIC[57]) - ALGEBRAIC[62]) - ALGEBRAIC[67]) - ALGEBRAIC[72];
resid[10] = RATES[10] - (((ALGEBRAIC[22] - ALGEBRAIC[38]) - ALGEBRAIC[26]) - ALGEBRAIC[28]) - ALGEBRAIC[30];
resid[11] = RATES[11] - ((ALGEBRAIC[24] - ALGEBRAIC[31]) - ALGEBRAIC[33]) - ALGEBRAIC[35];
resid[12] = RATES[12] - ALGEBRAIC[26] - ALGEBRAIC[87];
resid[13] = RATES[13] - ALGEBRAIC[28] - ALGEBRAIC[88];
resid[14] = RATES[14] - (ALGEBRAIC[30]+ALGEBRAIC[31]) - ALGEBRAIC[89];
resid[15] = RATES[15] - ALGEBRAIC[33] - ALGEBRAIC[90];
resid[16] = RATES[16] - (ALGEBRAIC[35]+ALGEBRAIC[86]) - ALGEBRAIC[91];
resid[17] = RATES[17] - (ALGEBRAIC[87]+ALGEBRAIC[88]+ALGEBRAIC[89]+ALGEBRAIC[90]+ALGEBRAIC[91]) - ALGEBRAIC[92];
resid[18] = RATES[18] - (ALGEBRAIC[38]+ALGEBRAIC[92]) - ALGEBRAIC[39];
resid[19] = RATES[19] - ALGEBRAIC[44] - ALGEBRAIC[46];
resid[20] = RATES[20] - ALGEBRAIC[44] - ALGEBRAIC[46];
resid[21] = RATES[21] - ALGEBRAIC[46] - ALGEBRAIC[48];
resid[22] = RATES[22] - ALGEBRAIC[46] - ALGEBRAIC[48];
resid[23] = RATES[23] - ALGEBRAIC[52] - ALGEBRAIC[53];
resid[24] = RATES[24] - ALGEBRAIC[53] - ALGEBRAIC[54];
resid[25] = RATES[25] - ALGEBRAIC[57] - ALGEBRAIC[58];
resid[26] = RATES[26] - ALGEBRAIC[58] - ALGEBRAIC[59];
resid[27] = RATES[27] - ALGEBRAIC[62] - ALGEBRAIC[63];
resid[28] = RATES[28] - ALGEBRAIC[63] - ALGEBRAIC[64];
resid[29] = RATES[29] - ALGEBRAIC[67] - ALGEBRAIC[68];
resid[30] = RATES[30] - ALGEBRAIC[68] - ALGEBRAIC[69];
resid[31] = RATES[31] - ALGEBRAIC[72] - ALGEBRAIC[73];
resid[32] = RATES[32] - ALGEBRAIC[73] - ALGEBRAIC[74];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (CONDVAR[0]>0.00000&&CONDVAR[1]<0.00000 ? 1.00000 : 0.00000);
ALGEBRAIC[2] = (VOI/CONSTANTS[3] - CONSTANTS[4]) - floor((VOI - CONSTANTS[4]*CONSTANTS[3])/CONSTANTS[3]);
ALGEBRAIC[47] = (ALGEBRAIC[41] - ALGEBRAIC[43])/CONSTANTS[95];
ALGEBRAIC[75] = STATES[0]+STATES[1]+STATES[2]+STATES[3]+STATES[9]+STATES[4];
ALGEBRAIC[76] = STATES[5]+STATES[6];
ALGEBRAIC[77] = STATES[19]+STATES[20]+STATES[21]+STATES[22];
ALGEBRAIC[78] = STATES[10]+STATES[11]+STATES[12]+STATES[13]+STATES[14]+STATES[15]+STATES[16]+STATES[17]+STATES[18]+STATES[23]+STATES[24];
ALGEBRAIC[79] = STATES[7]+STATES[8];
ALGEBRAIC[80] = STATES[25]+STATES[26]+STATES[27]+STATES[28]+STATES[29]+STATES[30]+STATES[31]+STATES[32];
ALGEBRAIC[81] = ALGEBRAIC[75]+ALGEBRAIC[76]+ALGEBRAIC[77]+ALGEBRAIC[78]+ALGEBRAIC[79]+ALGEBRAIC[80];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = VOI/CONSTANTS[3] - floor(VOI/CONSTANTS[3]);
ALGEBRAIC[5] = (CONDVAR[4]<0.00000 ? (( (CONSTANTS[9] - CONSTANTS[10])*ALGEBRAIC[1])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[1]/CONSTANTS[5], 2.00000)))+CONSTANTS[10] : CONDVAR[5]<0.00000 ? (( (CONSTANTS[9] - CONSTANTS[10])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[10] : CONSTANTS[10]);
ALGEBRAIC[6] = ALGEBRAIC[5]*(STATES[1] - CONSTANTS[11]);
ALGEBRAIC[8] = CONSTANTS[14]*(STATES[3] - CONSTANTS[15]);
ALGEBRAIC[9] = (CONDVAR[8]>0.00000 ? (ALGEBRAIC[8] - ALGEBRAIC[6])/CONSTANTS[18] : 0.00000);
ALGEBRAIC[3] = (CONDVAR[2]<0.00000 ? (( (CONSTANTS[6] - CONSTANTS[7])*ALGEBRAIC[1])/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(ALGEBRAIC[1]/CONSTANTS[5], 2.00000)))+CONSTANTS[7] : CONDVAR[3]<0.00000 ? (( (CONSTANTS[6] - CONSTANTS[7])*0.400000)/CONSTANTS[5])*exp( 0.500000*(1.00000 - pow(0.400000/CONSTANTS[5], 2.00000)))+CONSTANTS[7] : CONSTANTS[7]);
ALGEBRAIC[4] = ALGEBRAIC[3]*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[7] = CONSTANTS[12]*(STATES[2] - CONSTANTS[13]);
ALGEBRAIC[10] = (CONDVAR[9]>0.00000 ? (ALGEBRAIC[7] - ALGEBRAIC[4])/CONSTANTS[19] : 0.00000);
ALGEBRAIC[11] = CONSTANTS[21]*(STATES[4] - CONSTANTS[20]);
ALGEBRAIC[12] = CONSTANTS[23]*(STATES[5] - CONSTANTS[22]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[24]*cos(CONSTANTS[138]);
ALGEBRAIC[13] = (ALGEBRAIC[11] - ALGEBRAIC[12])/CONSTANTS[25];
ALGEBRAIC[14] = CONSTANTS[27]*(STATES[6] - CONSTANTS[28]);
ALGEBRAIC[15] = (ALGEBRAIC[14] - ALGEBRAIC[7])/CONSTANTS[29];
ALGEBRAIC[16] = CONSTANTS[32]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[17] = CONSTANTS[33]*(STATES[8] - CONSTANTS[31]);
ALGEBRAIC[18] = (ALGEBRAIC[16] - ALGEBRAIC[17])/CONSTANTS[35];
ALGEBRAIC[19] = (ALGEBRAIC[17] - ALGEBRAIC[8])/CONSTANTS[36];
ALGEBRAIC[20] = CONSTANTS[38]*(STATES[9] - CONSTANTS[37]);
ALGEBRAIC[21] = CONSTANTS[40]*(STATES[10] - CONSTANTS[39]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[41]*cos(CONSTANTS[138]);
ALGEBRAIC[22] = (ALGEBRAIC[20] - ALGEBRAIC[21])/CONSTANTS[42];
ALGEBRAIC[23] = CONSTANTS[44]*(STATES[11] - CONSTANTS[43]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[45]*cos(CONSTANTS[138]);
ALGEBRAIC[24] = (ALGEBRAIC[20] - ALGEBRAIC[23])/CONSTANTS[46];
ALGEBRAIC[25] = CONSTANTS[49]*(STATES[12] - CONSTANTS[48]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[50]*cos(CONSTANTS[138]);
ALGEBRAIC[26] = (ALGEBRAIC[21] - ALGEBRAIC[25])/CONSTANTS[51];
ALGEBRAIC[27] = CONSTANTS[70]*(STATES[13] - CONSTANTS[53]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[55]*cos(CONSTANTS[138]);
ALGEBRAIC[28] = (ALGEBRAIC[21] - ALGEBRAIC[27])/CONSTANTS[56];
ALGEBRAIC[29] = CONSTANTS[59]*(STATES[14] - CONSTANTS[58]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[60]*cos(CONSTANTS[138]);
ALGEBRAIC[30] = (ALGEBRAIC[21] - ALGEBRAIC[29])/CONSTANTS[61];
ALGEBRAIC[31] = (ALGEBRAIC[23] - ALGEBRAIC[29])/CONSTANTS[62];
ALGEBRAIC[32] = CONSTANTS[65]*(STATES[15] - CONSTANTS[64]);
ALGEBRAIC[33] = (ALGEBRAIC[23] - ALGEBRAIC[32])/CONSTANTS[67];
ALGEBRAIC[34] = CONSTANTS[70]*(STATES[16] - CONSTANTS[69]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[71]*cos(CONSTANTS[138]);
ALGEBRAIC[35] = (ALGEBRAIC[23] - ALGEBRAIC[34])/CONSTANTS[72];
ALGEBRAIC[37] = CONSTANTS[79]*(STATES[18] - CONSTANTS[78]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[80]*cos(CONSTANTS[138]);
ALGEBRAIC[38] = (ALGEBRAIC[21] - ALGEBRAIC[37])/CONSTANTS[81];
ALGEBRAIC[39] = (ALGEBRAIC[37] - ALGEBRAIC[8])/CONSTANTS[82];
ALGEBRAIC[40] = CONSTANTS[87]*(STATES[19] - CONSTANTS[83]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[44] = (ALGEBRAIC[20] - ALGEBRAIC[40])/CONSTANTS[92];
ALGEBRAIC[41] = CONSTANTS[88]*(STATES[20] - CONSTANTS[84]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[45] = (ALGEBRAIC[20] - ALGEBRAIC[41])/CONSTANTS[93];
ALGEBRAIC[42] = CONSTANTS[89]*(STATES[21] - CONSTANTS[85]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[46] = (ALGEBRAIC[40] - ALGEBRAIC[42])/CONSTANTS[94];
ALGEBRAIC[48] = (ALGEBRAIC[42] - ALGEBRAIC[8])/CONSTANTS[96];
ALGEBRAIC[43] = CONSTANTS[90]*(STATES[22] - CONSTANTS[86]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[91]*cos(CONSTANTS[138]);
ALGEBRAIC[49] = (ALGEBRAIC[43] - ALGEBRAIC[8])/CONSTANTS[97];
ALGEBRAIC[50] = CONSTANTS[100]*(STATES[23] - CONSTANTS[98]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[102]*cos(CONSTANTS[138]);
ALGEBRAIC[52] = (ALGEBRAIC[20] - ALGEBRAIC[50])/CONSTANTS[103];
ALGEBRAIC[51] = CONSTANTS[101]*(STATES[24] - CONSTANTS[99]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[102]*cos(CONSTANTS[138]);
ALGEBRAIC[53] = (ALGEBRAIC[50] - ALGEBRAIC[51])/CONSTANTS[104];
ALGEBRAIC[54] = (ALGEBRAIC[51] - ALGEBRAIC[8])/CONSTANTS[105];
ALGEBRAIC[55] = CONSTANTS[108]*(STATES[25] - CONSTANTS[106]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[110]*cos(CONSTANTS[138]);
ALGEBRAIC[57] = (ALGEBRAIC[20] - ALGEBRAIC[55])/CONSTANTS[111];
ALGEBRAIC[56] = CONSTANTS[109]*(STATES[26] - CONSTANTS[107]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[110]*cos(CONSTANTS[138]);
ALGEBRAIC[58] = (ALGEBRAIC[55] - ALGEBRAIC[56])/CONSTANTS[112];
ALGEBRAIC[59] = (ALGEBRAIC[56] - ALGEBRAIC[8])/CONSTANTS[113];
ALGEBRAIC[60] = CONSTANTS[116]*(STATES[27] - CONSTANTS[114]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[118]*cos(CONSTANTS[138]);
ALGEBRAIC[62] = (ALGEBRAIC[20] - ALGEBRAIC[60])/CONSTANTS[119];
ALGEBRAIC[61] = CONSTANTS[117]*(STATES[28] - CONSTANTS[114]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[118]*cos(CONSTANTS[138]);
ALGEBRAIC[63] = (ALGEBRAIC[60] - ALGEBRAIC[61])/CONSTANTS[120];
ALGEBRAIC[64] = (ALGEBRAIC[61] - ALGEBRAIC[8])/CONSTANTS[121];
ALGEBRAIC[65] = CONSTANTS[124]*(STATES[29] - CONSTANTS[122]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[126]*cos(CONSTANTS[138]);
ALGEBRAIC[67] = (ALGEBRAIC[20] - ALGEBRAIC[65])/CONSTANTS[127];
ALGEBRAIC[66] = CONSTANTS[125]*(STATES[30] - CONSTANTS[122]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[126]*cos(CONSTANTS[138]);
ALGEBRAIC[68] = (ALGEBRAIC[65] - ALGEBRAIC[66])/CONSTANTS[128];
ALGEBRAIC[69] = (ALGEBRAIC[66] - ALGEBRAIC[8])/CONSTANTS[129];
ALGEBRAIC[70] = CONSTANTS[132]*(STATES[31] - CONSTANTS[130]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[134]*cos(CONSTANTS[138]);
ALGEBRAIC[72] = (ALGEBRAIC[20] - ALGEBRAIC[70])/CONSTANTS[135];
ALGEBRAIC[71] = CONSTANTS[133]*(STATES[32] - CONSTANTS[130]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[134]*cos(CONSTANTS[138]);
ALGEBRAIC[73] = (ALGEBRAIC[70] - ALGEBRAIC[71])/CONSTANTS[136];
ALGEBRAIC[74] = (ALGEBRAIC[71] - ALGEBRAIC[8])/CONSTANTS[137];
ALGEBRAIC[82] = (CONDVAR[6]>0.00000 ? (ALGEBRAIC[4] - ALGEBRAIC[20])/CONSTANTS[16] : 0.00000);
ALGEBRAIC[83] = (CONDVAR[7]>0.00000 ? (ALGEBRAIC[6] - ALGEBRAIC[11])/CONSTANTS[17] : 0.00000);
ALGEBRAIC[84] = (ALGEBRAIC[12] - ALGEBRAIC[14])/CONSTANTS[26];
ALGEBRAIC[85] = (ALGEBRAIC[20] - ALGEBRAIC[16])/CONSTANTS[34];
ALGEBRAIC[86] = (ALGEBRAIC[20] - ALGEBRAIC[34])/CONSTANTS[47];
ALGEBRAIC[36] = CONSTANTS[75]*(STATES[17] - CONSTANTS[74]) - CONSTANTS[0]*CONSTANTS[1]*CONSTANTS[76]*cos(CONSTANTS[138]);
ALGEBRAIC[87] = (ALGEBRAIC[25] - ALGEBRAIC[36])/CONSTANTS[52];
ALGEBRAIC[88] = (ALGEBRAIC[27] - ALGEBRAIC[36])/CONSTANTS[57];
ALGEBRAIC[89] = (ALGEBRAIC[29] - ALGEBRAIC[36])/CONSTANTS[63];
ALGEBRAIC[90] = (ALGEBRAIC[32] - ALGEBRAIC[36])/CONSTANTS[68];
ALGEBRAIC[91] = (ALGEBRAIC[34] - ALGEBRAIC[36])/CONSTANTS[73];
ALGEBRAIC[92] = (ALGEBRAIC[36] - ALGEBRAIC[37])/CONSTANTS[77];
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
SI[8] = 1.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[13] = 1.0;
SI[14] = 1.0;
SI[15] = 1.0;
SI[16] = 1.0;
SI[17] = 1.0;
SI[18] = 1.0;
SI[19] = 1.0;
SI[20] = 1.0;
SI[21] = 1.0;
SI[22] = 1.0;
SI[23] = 1.0;
SI[24] = 1.0;
SI[25] = 1.0;
SI[26] = 1.0;
SI[27] = 1.0;
SI[28] = 1.0;
SI[29] = 1.0;
SI[30] = 1.0;
SI[31] = 1.0;
SI[32] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - 100.000;
CONDVAR[1] = VOI - 101.000;
CONDVAR[2] = ALGEBRAIC[1] - 0.400000;
CONDVAR[3] = ALGEBRAIC[1] - 0.600000;
CONDVAR[4] = ALGEBRAIC[1] - 0.400000;
CONDVAR[5] = ALGEBRAIC[1] - 0.600000;
CONDVAR[6] = ALGEBRAIC[4] - ALGEBRAIC[20];
CONDVAR[7] = ALGEBRAIC[6] - ALGEBRAIC[11];
CONDVAR[8] = ALGEBRAIC[8] - ALGEBRAIC[6];
CONDVAR[9] = ALGEBRAIC[7] - ALGEBRAIC[4];
}