/* There are a total of 30 entries in the algebraic variable array. There are a total of 36 entries in each of the rate and state variable arrays. There are a total of 136 entries in the constant variable array. */ /* * VOI is time in component environment (minute). * STATES[0] is mass in component mass (dimensionless). * CONSTANTS[0] is kg in component model_parameters (first_order_rate_constant). * STATES[1] is Cln2 in component Cln2 (dimensionless). * CONSTANTS[1] is ks_n2_ in component model_parameters (first_order_rate_constant). * CONSTANTS[2] is ks_n2__ in component model_parameters (first_order_rate_constant). * CONSTANTS[3] is kd_n2 in component model_parameters (first_order_rate_constant). * ALGEBRAIC[29] is SBF in component SBF_MBF (dimensionless). * STATES[2] is Clb5 in component Clb5 (dimensionless). * CONSTANTS[4] is ks_b5_ in component model_parameters (first_order_rate_constant). * CONSTANTS[5] is ks_b5__ in component model_parameters (first_order_rate_constant). * CONSTANTS[6] is kdi_f5 in component model_parameters (first_order_rate_constant). * CONSTANTS[7] is kdi_b5 in component model_parameters (first_order_rate_constant). * CONSTANTS[8] is kas_b5 in component model_parameters (first_order_rate_constant). * CONSTANTS[9] is kas_f5 in component model_parameters (first_order_rate_constant). * CONSTANTS[10] is kd3_f6 in component model_parameters (first_order_rate_constant). * CONSTANTS[11] is kd3_c1 in component model_parameters (first_order_rate_constant). * ALGEBRAIC[12] is Vd_b5 in component Vd_b5 (first_order_rate_constant). * ALGEBRAIC[28] is MBF in component SBF_MBF (dimensionless). * STATES[3] is C5P in component C5P (dimensionless). * STATES[4] is C5 in component C5 (dimensionless). * STATES[5] is F5P in component F5P (dimensionless). * STATES[6] is F5 in component F5 (dimensionless). * STATES[7] is Sic1 in component Sic1 (dimensionless). * STATES[8] is Cdc6 in component Cdc6 (dimensionless). * STATES[9] is Clb2 in component Clb2 (dimensionless). * CONSTANTS[12] is ks_b2_ in component model_parameters (first_order_rate_constant). * CONSTANTS[13] is ks_b2__ in component model_parameters (first_order_rate_constant). * CONSTANTS[14] is kdi_b2 in component model_parameters (first_order_rate_constant). * CONSTANTS[15] is kdi_f2 in component model_parameters (first_order_rate_constant). * CONSTANTS[16] is kas_b2 in component model_parameters (first_order_rate_constant). * CONSTANTS[17] is kas_f2 in component model_parameters (first_order_rate_constant). * ALGEBRAIC[13] is Vd_b2 in component Vd_b2 (first_order_rate_constant). * ALGEBRAIC[1] is Mcm1 in component Mcm1 (dimensionless). * STATES[10] is C2P in component C2P (dimensionless). * STATES[11] is C2 in component C2 (dimensionless). * STATES[12] is F2P in component F2P (dimensionless). * STATES[13] is F2 in component F2 (dimensionless). * CONSTANTS[18] is ks_c1_ in component model_parameters (first_order_rate_constant). * CONSTANTS[19] is ks_c1__ in component model_parameters (first_order_rate_constant). * CONSTANTS[20] is kpp_c1 in component model_parameters (first_order_rate_constant). * ALGEBRAIC[16] is Vkp_c1 in component Vkp_c1 (first_order_rate_constant). * STATES[14] is Swi5 in component Swi5 (dimensionless). * STATES[15] is Cdc14 in component Cdc14 (dimensionless). * STATES[16] is Sic1P in component Sic1P (dimensionless). * CONSTANTS[21] is ks_f6_ in component model_parameters (first_order_rate_constant). * CONSTANTS[22] is ks_f6__ in component model_parameters (first_order_rate_constant). * CONSTANTS[23] is ks_f6___ in component model_parameters (first_order_rate_constant). * CONSTANTS[24] is kpp_f6 in component model_parameters (first_order_rate_constant). * ALGEBRAIC[17] is Vkp_f6 in component Vkp_f6 (first_order_rate_constant). * STATES[17] is Cdc6P in component Cdc6P (dimensionless). * STATES[18] is Pds1 in component Pds1 (dimensionless). * CONSTANTS[25] is ks_pds_ in component model_parameters (first_order_rate_constant). * CONSTANTS[26] is ks1_pds__ in component model_parameters (first_order_rate_constant). * CONSTANTS[27] is ks2_pds__ in component model_parameters (first_order_rate_constant). * CONSTANTS[28] is kdi_esp in component model_parameters (first_order_rate_constant). * ALGEBRAIC[23] is Vd_pds in component Vd_pds (first_order_rate_constant). * CONSTANTS[29] is kas_esp in component model_parameters (first_order_rate_constant). * ALGEBRAIC[11] is PE in component PE (dimensionless). * STATES[19] is Esp1 in component Esp1 (dimensionless). * STATES[20] is ORI in component ORI (dimensionless). * CONSTANTS[30] is ks_ori in component model_parameters (first_order_rate_constant). * CONSTANTS[31] is kd_ori in component model_parameters (first_order_rate_constant). * CONSTANTS[32] is epsilon_ori_b5 in component model_parameters (dimensionless). * CONSTANTS[33] is epsilon_ori_b2 in component model_parameters (dimensionless). * STATES[21] is BUD in component BUD (dimensionless). * CONSTANTS[34] is ks_bud in component model_parameters (first_order_rate_constant). * CONSTANTS[35] is kd_bud in component model_parameters (first_order_rate_constant). * CONSTANTS[36] is epsilon_bud_n2 in component model_parameters (dimensionless). * CONSTANTS[37] is epsilon_bud_n3 in component model_parameters (dimensionless). * CONSTANTS[38] is epsilon_bud_b5 in component model_parameters (dimensionless). * ALGEBRAIC[2] is Cln3 in component Cln3 (dimensionless). * STATES[22] is SPN in component SPN (dimensionless). * CONSTANTS[39] is ks_spn in component model_parameters (first_order_rate_constant). * CONSTANTS[40] is kd_spn in component model_parameters (first_order_rate_constant). * CONSTANTS[41] is Jspn in component model_parameters (dimensionless). * ALGEBRAIC[27] is G_sbf in component G_sbf (dimensionless). * CONSTANTS[42] is Ji_sbf in component model_parameters (dimensionless). * CONSTANTS[43] is Ja_sbf in component model_parameters (dimensionless). * ALGEBRAIC[15] is Vi_sbf in component Vi_sbf (first_order_rate_constant). * ALGEBRAIC[14] is Va_sbf in component Va_sbf (first_order_rate_constant). * ALGEBRAIC[0] is G_mcm in component G_mcm (dimensionless). * CONSTANTS[44] is Ji_mcm in component model_parameters (dimensionless). * CONSTANTS[45] is Ja_mcm in component model_parameters (dimensionless). * CONSTANTS[46] is ki_mcm in component model_parameters (first_order_rate_constant). * CONSTANTS[47] is ka_mcm in component model_parameters (first_order_rate_constant). * CONSTANTS[48] is C0 in component model_parameters (dimensionless). * CONSTANTS[49] is Dn3 in component model_parameters (dimensionless). * CONSTANTS[50] is Jn3 in component model_parameters (dimensionless). * ALGEBRAIC[3] is Bck2 in component Bck2 (dimensionless). * CONSTANTS[51] is B0 in component model_parameters (dimensionless). * ALGEBRAIC[4] is Clb5_T in component Clb5_T (dimensionless). * ALGEBRAIC[5] is Clb2_T in component Clb2_T (dimensionless). * ALGEBRAIC[6] is Sic1_T in component Sic1_T (dimensionless). * STATES[23] is Swi5_T in component Swi5_T (dimensionless). * CONSTANTS[52] is ks_swi_ in component model_parameters (first_order_rate_constant). * CONSTANTS[53] is ks_swi__ in component model_parameters (first_order_rate_constant). * CONSTANTS[54] is kd_swi in component model_parameters (first_order_rate_constant). * CONSTANTS[55] is ka_swi in component model_parameters (first_order_rate_constant). * CONSTANTS[56] is ki_swi in component model_parameters (first_order_rate_constant). * STATES[24] is APC_P in component APC_P (dimensionless). * CONSTANTS[57] is ka_apc in component model_parameters (first_order_rate_constant). * CONSTANTS[58] is ki_apc in component model_parameters (first_order_rate_constant). * CONSTANTS[59] is Ja_apc in component model_parameters (dimensionless). * CONSTANTS[60] is Ji_apc in component model_parameters (dimensionless). * STATES[25] is Cdc20_T in component Cdc20_T (dimensionless). * CONSTANTS[61] is ks_20_ in component model_parameters (first_order_rate_constant). * CONSTANTS[62] is ks_20__ in component model_parameters (first_order_rate_constant). * CONSTANTS[63] is kd_20 in component model_parameters (first_order_rate_constant). * STATES[26] is Cdc20_A in component Cdc20_A (dimensionless). * CONSTANTS[64] is ka_20_ in component model_parameters (first_order_rate_constant). * CONSTANTS[65] is ka_20__ in component model_parameters (first_order_rate_constant). * ALGEBRAIC[24] is kmad2 in component model_parameters (first_order_rate_constant). * STATES[27] is Cdh1_T in component Cdh1_T (dimensionless). * CONSTANTS[66] is ks_cdh in component model_parameters (first_order_rate_constant). * CONSTANTS[67] is kd_cdh in component model_parameters (first_order_rate_constant). * STATES[28] is Cdh1 in component Cdh1 (dimensionless). * CONSTANTS[68] is Ja_cdh in component model_parameters (dimensionless). * CONSTANTS[69] is Ji_cdh in component model_parameters (dimensionless). * ALGEBRAIC[18] is Va_cdh in component Va_cdh (first_order_rate_constant). * ALGEBRAIC[19] is Vi_cdh in component Vi_cdh (first_order_rate_constant). * STATES[29] is Tem1 in component Tem1 (dimensionless). * ALGEBRAIC[25] is kbub2 in component model_parameters (first_order_rate_constant). * ALGEBRAIC[26] is klte1 in component model_parameters (first_order_rate_constant). * CONSTANTS[70] is Ja_tem in component model_parameters (dimensionless). * CONSTANTS[71] is Ji_tem in component model_parameters (dimensionless). * CONSTANTS[72] is Tem1_T in component model_parameters (dimensionless). * STATES[30] is Cdc15 in component Cdc15 (dimensionless). * CONSTANTS[73] is ka_15_ in component model_parameters (first_order_rate_constant). * CONSTANTS[74] is ka_15__ in component model_parameters (first_order_rate_constant). * CONSTANTS[75] is ka_15___ in component model_parameters (first_order_rate_constant). * CONSTANTS[76] is ki_15 in component model_parameters (first_order_rate_constant). * CONSTANTS[77] is Cdc15_T in component model_parameters (dimensionless). * STATES[31] is Cdc14_T in component Cdc14_T (dimensionless). * CONSTANTS[78] is ks_14 in component model_parameters (first_order_rate_constant). * CONSTANTS[79] is kd_14 in component model_parameters (first_order_rate_constant). * CONSTANTS[80] is kd_net in component model_parameters (first_order_rate_constant). * CONSTANTS[81] is kdi_rent in component model_parameters (first_order_rate_constant). * CONSTANTS[82] is kdi_rentp in component model_parameters (first_order_rate_constant). * CONSTANTS[83] is kas_rent in component model_parameters (first_order_rate_constant). * CONSTANTS[84] is kas_rentp in component model_parameters (first_order_rate_constant). * STATES[32] is RENT in component RENT (dimensionless). * ALGEBRAIC[9] is RENTP in component RENTP (dimensionless). * STATES[33] is Net1 in component Net1 (dimensionless). * ALGEBRAIC[10] is Net1P in component Net1P (dimensionless). * STATES[34] is Net1_T in component Net1_T (dimensionless). * CONSTANTS[85] is ks_net in component model_parameters (first_order_rate_constant). * ALGEBRAIC[7] is Cdc6_T in component Cdc6_T (dimensionless). * ALGEBRAIC[8] is CKI_T in component CKI_T (dimensionless). * CONSTANTS[86] is Esp1_T in component model_parameters (dimensionless). * CONSTANTS[87] is kd_b5_ in component model_parameters (first_order_rate_constant). * CONSTANTS[88] is kd_b5__ in component model_parameters (first_order_rate_constant). * CONSTANTS[89] is kd_b2_ in component model_parameters (first_order_rate_constant). * CONSTANTS[90] is kd_b2__ in component model_parameters (first_order_rate_constant). * CONSTANTS[91] is kd_b2p in component model_parameters (first_order_rate_constant). * CONSTANTS[92] is ka_sbf in component model_parameters (first_order_rate_constant). * CONSTANTS[93] is epsilon_sbf_n2 in component model_parameters (dimensionless). * CONSTANTS[94] is epsilon_sbf_n3 in component model_parameters (dimensionless). * CONSTANTS[95] is epsilon_sbf_b5 in component model_parameters (dimensionless). * CONSTANTS[96] is ki_sbf_ in component model_parameters (first_order_rate_constant). * CONSTANTS[97] is ki_sbf__ in component model_parameters (first_order_rate_constant). * CONSTANTS[98] is kd1_c1 in component model_parameters (first_order_rate_constant). * CONSTANTS[99] is kd2_c1 in component model_parameters (first_order_rate_constant). * CONSTANTS[100] is Jd2_c1 in component model_parameters (dimensionless). * CONSTANTS[101] is epsilon_c1_n2 in component model_parameters (dimensionless). * CONSTANTS[102] is epsilon_c1_n3 in component model_parameters (dimensionless). * CONSTANTS[103] is epsilon_c1_k2 in component model_parameters (dimensionless). * CONSTANTS[104] is epsilon_c1_b5 in component model_parameters (dimensionless). * CONSTANTS[105] is epsilon_c1_b2 in component model_parameters (dimensionless). * CONSTANTS[106] is Jd2_f6 in component model_parameters (dimensionless). * CONSTANTS[107] is kd1_f6 in component model_parameters (first_order_rate_constant). * CONSTANTS[108] is kd2_f6 in component model_parameters (first_order_rate_constant). * CONSTANTS[109] is epsilon_f6_n2 in component model_parameters (dimensionless). * CONSTANTS[110] is epsilon_f6_n3 in component model_parameters (dimensionless). * CONSTANTS[111] is epsilon_f6_k2 in component model_parameters (dimensionless). * CONSTANTS[112] is epsilon_f6_b5 in component model_parameters (dimensionless). * CONSTANTS[113] is epsilon_f6_b2 in component model_parameters (dimensionless). * CONSTANTS[114] is ka_cdh_ in component model_parameters (first_order_rate_constant). * CONSTANTS[115] is ka_cdh__ in component model_parameters (first_order_rate_constant). * CONSTANTS[116] is ki_cdh_ in component model_parameters (first_order_rate_constant). * CONSTANTS[117] is ki_cdh__ in component model_parameters (first_order_rate_constant). * CONSTANTS[118] is epsilon_cdh_n2 in component model_parameters (dimensionless). * CONSTANTS[119] is epsilon_cdh_n3 in component model_parameters (dimensionless). * CONSTANTS[120] is epsilon_cdh_b5 in component model_parameters (dimensionless). * CONSTANTS[121] is epsilon_cdh_b2 in component model_parameters (dimensionless). * ALGEBRAIC[20] is Vpp_net in component Vpp_net (first_order_rate_constant). * CONSTANTS[122] is kpp_net_ in component model_parameters (first_order_rate_constant). * CONSTANTS[123] is kpp_net__ in component model_parameters (first_order_rate_constant). * STATES[35] is PPX in component PPX (dimensionless). * ALGEBRAIC[21] is Vkp_net in component Vkp_net (first_order_rate_constant). * CONSTANTS[124] is kkp_net_ in component model_parameters (first_order_rate_constant). * CONSTANTS[125] is kkp_net__ in component model_parameters (first_order_rate_constant). * CONSTANTS[126] is ks_ppx in component model_parameters (first_order_rate_constant). * ALGEBRAIC[22] is Vd_ppx in component Vd_ppx (first_order_rate_constant). * CONSTANTS[127] is kd_ppx_ in component model_parameters (first_order_rate_constant). * CONSTANTS[128] is kd_ppx__ in component model_parameters (first_order_rate_constant). * CONSTANTS[129] is Jpds in component model_parameters (dimensionless). * CONSTANTS[130] is J20_ppx in component model_parameters (dimensionless). * CONSTANTS[131] is kd1_pds_ in component model_parameters (first_order_rate_constant). * CONSTANTS[132] is kd2_pds__ in component model_parameters (first_order_rate_constant). * CONSTANTS[133] is kd3_pds__ in component model_parameters (first_order_rate_constant). * CONSTANTS[134] is Kez in component model_parameters (dimensionless). * CONSTANTS[135] is Kez2 in component model_parameters (dimensionless). * RATES[0] is d/dt mass in component mass (dimensionless). * RATES[1] is d/dt Cln2 in component Cln2 (dimensionless). * RATES[2] is d/dt Clb5 in component Clb5 (dimensionless). * RATES[9] is d/dt Clb2 in component Clb2 (dimensionless). * RATES[7] is d/dt Sic1 in component Sic1 (dimensionless). * RATES[16] is d/dt Sic1P in component Sic1P (dimensionless). * RATES[11] is d/dt C2 in component C2 (dimensionless). * RATES[4] is d/dt C5 in component C5 (dimensionless). * RATES[10] is d/dt C2P in component C2P (dimensionless). * RATES[3] is d/dt C5P in component C5P (dimensionless). * RATES[8] is d/dt Cdc6 in component Cdc6 (dimensionless). * RATES[17] is d/dt Cdc6P in component Cdc6P (dimensionless). * RATES[13] is d/dt F2 in component F2 (dimensionless). * RATES[18] is d/dt Pds1 in component Pds1 (dimensionless). * RATES[19] is d/dt Esp1 in component Esp1 (dimensionless). * RATES[20] is d/dt ORI in component ORI (dimensionless). * RATES[21] is d/dt BUD in component BUD (dimensionless). * RATES[22] is d/dt SPN in component SPN (dimensionless). * RATES[6] is d/dt F5 in component F5 (dimensionless). * RATES[12] is d/dt F2P in component F2P (dimensionless). * RATES[5] is d/dt F5P in component F5P (dimensionless). * RATES[23] is d/dt Swi5_T in component Swi5_T (dimensionless). * RATES[14] is d/dt Swi5 in component Swi5 (dimensionless). * RATES[24] is d/dt APC_P in component APC_P (dimensionless). * RATES[25] is d/dt Cdc20_T in component Cdc20_T (dimensionless). * RATES[26] is d/dt Cdc20_A in component Cdc20_A (dimensionless). * RATES[27] is d/dt Cdh1_T in component Cdh1_T (dimensionless). * RATES[28] is d/dt Cdh1 in component Cdh1 (dimensionless). * RATES[29] is d/dt Tem1 in component Tem1 (dimensionless). * RATES[30] is d/dt Cdc15 in component Cdc15 (dimensionless). * RATES[31] is d/dt Cdc14_T in component Cdc14_T (dimensionless). * RATES[15] is d/dt Cdc14 in component Cdc14 (dimensionless). * RATES[34] is d/dt Net1_T in component Net1_T (dimensionless). * RATES[33] is d/dt Net1 in component Net1 (dimensionless). * RATES[32] is d/dt RENT in component RENT (dimensionless). * RATES[35] is d/dt PPX in component PPX (dimensionless). * There are a total of 6 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = 1.206; CONSTANTS[0] = 0.007702; STATES[1] = 0.0652; CONSTANTS[1] = 0; CONSTANTS[2] = 0.15; CONSTANTS[3] = 0.12; STATES[2] = 0.0518; CONSTANTS[4] = 0.0008; CONSTANTS[5] = 0.005; CONSTANTS[6] = 0.01; CONSTANTS[7] = 0.06; CONSTANTS[8] = 50; CONSTANTS[9] = 0.01; CONSTANTS[10] = 1; CONSTANTS[11] = 1; STATES[3] = 0.0069; STATES[4] = 0.0701; STATES[5] = 7.9e-6; STATES[6] = 7.2e-5; STATES[7] = 0.0229; STATES[8] = 0.1076; STATES[9] = 0.1469; CONSTANTS[12] = 0.001; CONSTANTS[13] = 0.04; CONSTANTS[14] = 0.05; CONSTANTS[15] = 0.5; CONSTANTS[16] = 50; CONSTANTS[17] = 15; STATES[10] = 0.024; STATES[11] = 0.2384; STATES[12] = 0.0274; STATES[13] = 0.2361; CONSTANTS[18] = 0.012; CONSTANTS[19] = 0.12; CONSTANTS[20] = 4; STATES[14] = 0.9562; STATES[15] = 0.4683; STATES[16] = 0.0064; CONSTANTS[21] = 0.024; CONSTANTS[22] = 0.12; CONSTANTS[23] = 0.004; CONSTANTS[24] = 4; STATES[17] = 0.0155; STATES[18] = 0.0256; CONSTANTS[25] = 0; CONSTANTS[26] = 0.03; CONSTANTS[27] = 0.055; CONSTANTS[28] = 0.5; CONSTANTS[29] = 50; STATES[19] = 0.3013; STATES[20] = 0.0009; CONSTANTS[30] = 2; CONSTANTS[31] = 0.06; CONSTANTS[32] = 0.9; CONSTANTS[33] = 0.45; STATES[21] = 0.0085; CONSTANTS[34] = 0.2; CONSTANTS[35] = 0.06; CONSTANTS[36] = 0.25; CONSTANTS[37] = 0.05; CONSTANTS[38] = 1; STATES[22] = 0.0305; CONSTANTS[39] = 0.1; CONSTANTS[40] = 0.06; CONSTANTS[41] = 0.14; CONSTANTS[42] = 0.01; CONSTANTS[43] = 0.01; CONSTANTS[44] = 0.1; CONSTANTS[45] = 0.1; CONSTANTS[46] = 0.15; CONSTANTS[47] = 1; CONSTANTS[48] = 0.4; CONSTANTS[49] = 1; CONSTANTS[50] = 6; CONSTANTS[51] = 0.054; STATES[23] = 0.9765; CONSTANTS[52] = 0.005; CONSTANTS[53] = 0.08; CONSTANTS[54] = 0.08; CONSTANTS[55] = 2; CONSTANTS[56] = 0.05; STATES[24] = 0.1015; CONSTANTS[57] = 0.1; CONSTANTS[58] = 0.15; CONSTANTS[59] = 0.1; CONSTANTS[60] = 0.1; STATES[25] = 1.9163; CONSTANTS[61] = 0.006; CONSTANTS[62] = 0.6; CONSTANTS[63] = 0.3; STATES[26] = 0.4443; CONSTANTS[64] = 0.05; CONSTANTS[65] = 0.2; STATES[27] = 1; CONSTANTS[66] = 0.01; CONSTANTS[67] = 0.01; STATES[28] = 0.9305; CONSTANTS[68] = 0.03; CONSTANTS[69] = 0.03; STATES[29] = 0.9039; CONSTANTS[70] = 0.1; CONSTANTS[71] = 0.1; CONSTANTS[72] = 1; STATES[30] = 0.6565; CONSTANTS[73] = 0.002; CONSTANTS[74] = 1; CONSTANTS[75] = 0.001; CONSTANTS[76] = 0.5; CONSTANTS[77] = 1; STATES[31] = 2; CONSTANTS[78] = 0.2; CONSTANTS[79] = 0.1; CONSTANTS[80] = 0.03; CONSTANTS[81] = 1; CONSTANTS[82] = 2; CONSTANTS[83] = 200; CONSTANTS[84] = 1; STATES[32] = 1.0495; STATES[33] = 0.0186; STATES[34] = 2.8; CONSTANTS[85] = 0.084; CONSTANTS[86] = 1; CONSTANTS[87] = 0.01; CONSTANTS[88] = 0.16; CONSTANTS[89] = 0.003; CONSTANTS[90] = 0.4; CONSTANTS[91] = 0.15; CONSTANTS[92] = 0.38; CONSTANTS[93] = 2; CONSTANTS[94] = 10; CONSTANTS[95] = 2; CONSTANTS[96] = 0.6; CONSTANTS[97] = 8; CONSTANTS[98] = 0.01; CONSTANTS[99] = 1; CONSTANTS[100] = 0.05; CONSTANTS[101] = 0.06; CONSTANTS[102] = 0.3; CONSTANTS[103] = 0.03; CONSTANTS[104] = 0.1; CONSTANTS[105] = 0.45; CONSTANTS[106] = 0.05; CONSTANTS[107] = 0.01; CONSTANTS[108] = 1; CONSTANTS[109] = 0.06; CONSTANTS[110] = 0.3; CONSTANTS[111] = 0.03; CONSTANTS[112] = 0.1; CONSTANTS[113] = 0.55; CONSTANTS[114] = 0.01; CONSTANTS[115] = 0.8; CONSTANTS[116] = 0.001; CONSTANTS[117] = 0.08; CONSTANTS[118] = 0.4; CONSTANTS[119] = 0.25; CONSTANTS[120] = 8; CONSTANTS[121] = 1.2; CONSTANTS[122] = 0.05; CONSTANTS[123] = 3; STATES[35] = 0.1232; CONSTANTS[124] = 0.01; CONSTANTS[125] = 0.6; CONSTANTS[126] = 0.1; CONSTANTS[127] = 0.17; CONSTANTS[128] = 2; CONSTANTS[129] = 0.04; CONSTANTS[130] = 0.15; CONSTANTS[131] = 0.01; CONSTANTS[132] = 0.2; CONSTANTS[133] = 0.04; CONSTANTS[134] = 0.3; CONSTANTS[135] = 0.2; RATES[0] = 0.1001; RATES[1] = 0.1001; RATES[2] = 0.1001; RATES[9] = 0.1001; RATES[7] = 0.1001; RATES[16] = 0.1001; RATES[11] = 0.1001; RATES[4] = 0.1001; RATES[10] = 0.1001; RATES[3] = 0.1001; RATES[8] = 0.1001; RATES[17] = 0.1001; RATES[13] = 0.1001; RATES[18] = 0.1001; RATES[19] = 0.1001; RATES[20] = 0.1001; RATES[21] = 0.1001; RATES[22] = 0.1001; RATES[6] = 0.1001; RATES[12] = 0.1001; RATES[5] = 0.1001; RATES[23] = 0.1001; RATES[14] = 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[15] = 0.1001; RATES[34] = 0.1001; RATES[33] = 0.1001; RATES[32] = 0.1001; RATES[35] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[0] - STATES[0]*CONSTANTS[0]; resid[1] = RATES[1] - (CONSTANTS[1]+ CONSTANTS[2]*ALGEBRAIC[29])*STATES[0] - CONSTANTS[3]*STATES[1]; resid[2] = RATES[2] - ( (CONSTANTS[4]+ CONSTANTS[5]*ALGEBRAIC[28])*STATES[0]+ CONSTANTS[11]*STATES[3]+ CONSTANTS[7]*STATES[4]+ CONSTANTS[10]*STATES[5]+ CONSTANTS[6]*STATES[6]) - (ALGEBRAIC[12]+ CONSTANTS[8]*STATES[7]+ CONSTANTS[9]*STATES[8])*STATES[2]; resid[3] = RATES[9] - ( (CONSTANTS[12]+ CONSTANTS[13]*ALGEBRAIC[1])*STATES[0]+ CONSTANTS[11]*STATES[10]+ CONSTANTS[14]*STATES[11]+ CONSTANTS[10]*STATES[12]+ CONSTANTS[15]*STATES[13]) - (ALGEBRAIC[13]+ CONSTANTS[16]*STATES[7]+ CONSTANTS[17]*STATES[8])*STATES[9]; resid[4] = RATES[7] - (CONSTANTS[18]+ CONSTANTS[19]*STATES[14]+ (ALGEBRAIC[13]+CONSTANTS[14])*STATES[11]+ (ALGEBRAIC[12]+CONSTANTS[7])*STATES[4]+ CONSTANTS[20]*STATES[15]*STATES[16]) - ( CONSTANTS[16]*STATES[9]+ CONSTANTS[8]*STATES[2]+ALGEBRAIC[16])*STATES[7]; resid[5] = RATES[16] - ( ALGEBRAIC[16]*STATES[7] - ( CONSTANTS[20]*STATES[15]+CONSTANTS[11])*STATES[16])+ ALGEBRAIC[13]*STATES[10]+ ALGEBRAIC[12]*STATES[3]; resid[6] = RATES[11] - ( CONSTANTS[16]*STATES[9]*STATES[7]+ CONSTANTS[20]*STATES[15]*STATES[10]) - (CONSTANTS[14]+ALGEBRAIC[13]+ALGEBRAIC[16])*STATES[11]; resid[7] = RATES[4] - ( CONSTANTS[8]*STATES[2]*STATES[7]+ CONSTANTS[20]*STATES[15]*STATES[3]) - (CONSTANTS[7]+ALGEBRAIC[12]+ALGEBRAIC[16])*STATES[4]; resid[8] = RATES[10] - ALGEBRAIC[16]*STATES[11] - ( CONSTANTS[20]*STATES[15]+CONSTANTS[11]+ALGEBRAIC[13])*STATES[10]; resid[9] = RATES[3] - ALGEBRAIC[16]*STATES[4] - ( CONSTANTS[20]*STATES[15]+CONSTANTS[11]+ALGEBRAIC[12])*STATES[3]; resid[10] = RATES[8] - (CONSTANTS[21]+ CONSTANTS[22]*STATES[14]+ CONSTANTS[23]*ALGEBRAIC[29]+ (ALGEBRAIC[13]+CONSTANTS[15])*STATES[13]+ (ALGEBRAIC[12]+CONSTANTS[6])*STATES[6]+ CONSTANTS[24]*STATES[15]*STATES[17]) - ( CONSTANTS[17]*STATES[9]+ CONSTANTS[9]*STATES[2]+ALGEBRAIC[17])*STATES[8]; resid[11] = RATES[17] - ( ALGEBRAIC[17]*STATES[8] - ( CONSTANTS[24]*STATES[15]+CONSTANTS[10])*STATES[17])+ ALGEBRAIC[13]*STATES[12]+ ALGEBRAIC[12]*STATES[5]; resid[12] = RATES[13] - ( CONSTANTS[17]*STATES[9]*STATES[8]+ CONSTANTS[24]*STATES[15]*STATES[12]) - (CONSTANTS[15]+ALGEBRAIC[13]+ALGEBRAIC[17])*STATES[13]; resid[13] = RATES[18] - (CONSTANTS[25]+ CONSTANTS[26]*ALGEBRAIC[29]+ CONSTANTS[27]*ALGEBRAIC[1]+ CONSTANTS[28]*ALGEBRAIC[11]) - (ALGEBRAIC[23]+ CONSTANTS[29]*STATES[19])*STATES[18]; resid[14] = RATES[19] - - CONSTANTS[29]*STATES[18]*STATES[19]+ (CONSTANTS[28]+ALGEBRAIC[23])*ALGEBRAIC[11]; resid[15] = RATES[20] - CONSTANTS[30]*( CONSTANTS[32]*STATES[2]+ CONSTANTS[33]*STATES[9]) - CONSTANTS[31]*STATES[20]; resid[16] = RATES[21] - CONSTANTS[34]*( CONSTANTS[36]*STATES[1]+ CONSTANTS[37]*ALGEBRAIC[2]+ CONSTANTS[38]*STATES[2]) - CONSTANTS[35]*STATES[21]; resid[17] = RATES[22] - ( CONSTANTS[39]*STATES[9])/(CONSTANTS[41]+STATES[9]) - CONSTANTS[40]*STATES[22]; resid[18] = RATES[6] - ( CONSTANTS[9]*STATES[2]*STATES[8]+ CONSTANTS[24]*STATES[15]*STATES[5]) - (CONSTANTS[6]+ALGEBRAIC[12]+ALGEBRAIC[17])*STATES[6]; resid[19] = RATES[12] - ALGEBRAIC[17]*STATES[13] - ( CONSTANTS[24]*STATES[15]+CONSTANTS[10]+ALGEBRAIC[13])*STATES[12]; resid[20] = RATES[5] - ALGEBRAIC[17]*STATES[6] - ( CONSTANTS[24]*STATES[15]+CONSTANTS[10]+ALGEBRAIC[12])*STATES[5]; resid[21] = RATES[23] - (CONSTANTS[52]+ CONSTANTS[53]*ALGEBRAIC[1]) - CONSTANTS[54]*STATES[23]; resid[22] = RATES[14] - (CONSTANTS[52]+ CONSTANTS[53]*ALGEBRAIC[1]+ CONSTANTS[55]*STATES[15]*(STATES[23] - STATES[14])) - (CONSTANTS[54]+ CONSTANTS[56]*STATES[9])*STATES[14]; resid[23] = RATES[24] - ( CONSTANTS[57]*STATES[9]*(1.00000 - STATES[24]))/((CONSTANTS[59]+1.00000) - STATES[24]) - ( CONSTANTS[58]*STATES[24])/(CONSTANTS[60]+STATES[24]); resid[24] = RATES[25] - (CONSTANTS[61]+ CONSTANTS[62]*ALGEBRAIC[1]) - CONSTANTS[63]*STATES[25]; resid[25] = RATES[26] - (CONSTANTS[64]+ CONSTANTS[65]*STATES[24])*(STATES[25] - STATES[26]) - (ALGEBRAIC[24]+CONSTANTS[63])*STATES[26]; resid[26] = RATES[27] - CONSTANTS[66] - CONSTANTS[67]*STATES[27]; resid[27] = RATES[28] - (CONSTANTS[66]+( ALGEBRAIC[18]*(STATES[27] - STATES[28]))/((CONSTANTS[68]+STATES[27]) - STATES[28])) - ( CONSTANTS[67]*STATES[28]+( ALGEBRAIC[19]*STATES[28])/(CONSTANTS[69]+STATES[28])); resid[28] = RATES[29] - ( ALGEBRAIC[26]*(CONSTANTS[72] - STATES[29]))/((CONSTANTS[70]+CONSTANTS[72]) - STATES[29]) - ( ALGEBRAIC[25]*STATES[29])/(CONSTANTS[71]+STATES[29]); resid[29] = RATES[30] - ( CONSTANTS[73]*(CONSTANTS[72] - STATES[29])+ CONSTANTS[74]*STATES[29]+ CONSTANTS[75]*STATES[15])*(CONSTANTS[77] - STATES[30]) - CONSTANTS[76]*STATES[30]; resid[30] = RATES[31] - CONSTANTS[78] - CONSTANTS[79]*STATES[31]; resid[31] = RATES[15] - (CONSTANTS[78]+ CONSTANTS[80]*(STATES[32]+ALGEBRAIC[9])+ CONSTANTS[81]*STATES[32]+ CONSTANTS[82]*ALGEBRAIC[9]) - ( CONSTANTS[79]*STATES[15]+ ( CONSTANTS[83]*STATES[33]+ CONSTANTS[84]*ALGEBRAIC[10])*STATES[15]); resid[32] = RATES[34] - CONSTANTS[85] - CONSTANTS[80]*STATES[34]; resid[33] = RATES[33] - (CONSTANTS[85]+ CONSTANTS[79]*STATES[32]+ CONSTANTS[81]*STATES[32]+ ALGEBRAIC[20]*ALGEBRAIC[10]) - ( CONSTANTS[80]*STATES[33]+ CONSTANTS[83]*STATES[15]*STATES[33]+ ALGEBRAIC[21]*STATES[33]); resid[34] = RATES[32] - ( CONSTANTS[83]*STATES[15]*STATES[33]+ ALGEBRAIC[20]*ALGEBRAIC[9]) - ( (CONSTANTS[79]+CONSTANTS[80])*STATES[32]+ CONSTANTS[81]*STATES[32]+ ALGEBRAIC[21]*STATES[32]); resid[35] = RATES[35] - CONSTANTS[126] - ALGEBRAIC[22]*STATES[35]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[4] = STATES[2]+STATES[4]+STATES[3]+STATES[6]+STATES[5]; ALGEBRAIC[5] = STATES[9]+STATES[11]+STATES[10]+STATES[13]+STATES[12]; ALGEBRAIC[8] = ALGEBRAIC[6]+ALGEBRAIC[7]; } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = ( 2.00000*CONSTANTS[44]*CONSTANTS[47]*STATES[9])/((CONSTANTS[46]+ CONSTANTS[45]*CONSTANTS[46]+ CONSTANTS[44]*CONSTANTS[47]*STATES[9]+ pow((pow((CONSTANTS[46]+ CONSTANTS[45]*CONSTANTS[46]+ CONSTANTS[44]*CONSTANTS[47]*STATES[9]) - CONSTANTS[47]*STATES[9], 2.00000) - 4.00000*(CONSTANTS[46] - CONSTANTS[47]*STATES[9])*CONSTANTS[44]*CONSTANTS[47]*STATES[9]), 1.0 / 2)) - CONSTANTS[47]*STATES[9]); ALGEBRAIC[1] = ALGEBRAIC[0]; ALGEBRAIC[2] = ( CONSTANTS[48]*CONSTANTS[49]*STATES[0])/(CONSTANTS[50]+ CONSTANTS[49]*STATES[0]); ALGEBRAIC[9] = STATES[31] - (STATES[32]+STATES[15]); ALGEBRAIC[10] = (STATES[34]+STATES[15]) - (STATES[33]+STATES[31]); ALGEBRAIC[11] = CONSTANTS[86] - STATES[19]; ALGEBRAIC[12] = CONSTANTS[87]+ CONSTANTS[88]*STATES[26]; ALGEBRAIC[13] = CONSTANTS[89]+ CONSTANTS[90]*STATES[28]+ CONSTANTS[91]*STATES[26]; ALGEBRAIC[3] = CONSTANTS[51]*STATES[0]; ALGEBRAIC[6] = STATES[7]+STATES[16]+STATES[11]+STATES[10]+STATES[4]+STATES[3]; ALGEBRAIC[16] = CONSTANTS[98]+( CONSTANTS[99]*( CONSTANTS[102]*ALGEBRAIC[2]+ CONSTANTS[103]*ALGEBRAIC[3]+ CONSTANTS[101]*STATES[1]+ CONSTANTS[104]*STATES[2]+ CONSTANTS[105]*STATES[9]))/(CONSTANTS[100]+ALGEBRAIC[6]); ALGEBRAIC[7] = STATES[8]+STATES[17]+STATES[13]+STATES[12]+STATES[6]+STATES[5]; ALGEBRAIC[17] = CONSTANTS[107]+( CONSTANTS[108]*( CONSTANTS[110]*ALGEBRAIC[2]+ CONSTANTS[111]*ALGEBRAIC[3]+ CONSTANTS[109]*STATES[1]+ CONSTANTS[112]*STATES[2]+ CONSTANTS[113]*STATES[9]))/(CONSTANTS[106]+ALGEBRAIC[7]); ALGEBRAIC[18] = CONSTANTS[114]+ CONSTANTS[115]*STATES[15]; ALGEBRAIC[19] = CONSTANTS[116]+ CONSTANTS[117]*( CONSTANTS[119]*ALGEBRAIC[2]+ CONSTANTS[118]*STATES[1]+ CONSTANTS[120]*STATES[2]+ CONSTANTS[121]*STATES[9]); ALGEBRAIC[20] = CONSTANTS[122]+ CONSTANTS[123]*STATES[35]; ALGEBRAIC[21] = (CONSTANTS[124]+ CONSTANTS[125]*STATES[30])*STATES[0]; ALGEBRAIC[22] = CONSTANTS[127]+( CONSTANTS[128]*(CONSTANTS[130]+STATES[26])*CONSTANTS[129])/(CONSTANTS[129]+STATES[18]); ALGEBRAIC[23] = CONSTANTS[131]+ CONSTANTS[132]*STATES[26]+ CONSTANTS[133]*STATES[28]; ALGEBRAIC[24] = (CONDVAR[0]>0.00000&&CONDVAR[1]<0.00000 ? 8.00000 : 0.0100000); ALGEBRAIC[25] = (CONDVAR[2]>0.00000&&CONDVAR[3]<0.00000 ? 1.00000 : 0.200000); ALGEBRAIC[26] = (CONDVAR[4]>0.00000&&CONDVAR[5]>0.00000 ? 1.00000 : 0.100000); ALGEBRAIC[15] = CONSTANTS[96]+ CONSTANTS[97]*STATES[9]; ALGEBRAIC[14] = CONSTANTS[92]*( CONSTANTS[93]*STATES[1]+ CONSTANTS[94]*(ALGEBRAIC[2]+ALGEBRAIC[3])+ CONSTANTS[95]*STATES[2]); ALGEBRAIC[27] = ( 2.00000*CONSTANTS[42]*ALGEBRAIC[14])/((ALGEBRAIC[15]+ CONSTANTS[43]*ALGEBRAIC[15]+ CONSTANTS[42]*ALGEBRAIC[14]+ pow((pow((ALGEBRAIC[15]+ CONSTANTS[43]*ALGEBRAIC[15]+ CONSTANTS[42]*ALGEBRAIC[14]) - ALGEBRAIC[14], 2.00000) - 4.00000*(ALGEBRAIC[15] - ALGEBRAIC[14])*CONSTANTS[42]*ALGEBRAIC[14]), 1.0 / 2)) - ALGEBRAIC[14]); ALGEBRAIC[28] = ALGEBRAIC[27]; ALGEBRAIC[29] = ALGEBRAIC[28]; } 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; SI[33] = 1.0; SI[34] = 1.0; SI[35] = 1.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = STATES[20] - 1.00000; CONDVAR[1] = STATES[22] - 1.00000; CONDVAR[2] = STATES[20] - 1.00000; CONDVAR[3] = STATES[22] - 1.00000; CONDVAR[4] = STATES[22] - 1.00000; CONDVAR[5] = STATES[9] - CONSTANTS[134]; }