/* There are a total of 15 entries in the algebraic variable array. There are a total of 7 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 cap in component parameters (nanofarad). * CONSTANTS[1] is cc1lin in component parameters (per_second). * CONSTANTS[2] is cc_2 in component parameters (per_second). * CONSTANTS[3] is ck1lin in component parameters (per_second). * CONSTANTS[4] is ck_2 in component parameters (per_second). * CONSTANTS[5] is clmax in component parameters (nanosiemens). * CONSTANTS[6] is cnmax in component parameters (nanosiemens). * CONSTANTS[7] is cx1lin in component parameters (per_second). * CONSTANTS[8] is cx2 in component parameters (per_second). * CONSTANTS[9] is ef in component parameters (per_second). * CONSTANTS[10] is gl in component parameters (nanosiemens). * CONSTANTS[11] is hmc_1 in component parameters (uM). * CONSTANTS[12] is hmc_2 in component parameters (uM). * CONSTANTS[13] is inf in component parameters (uM_per_picocoulomb). * CONSTANTS[14] is inhmax in component parameters (dimensionless). * CONSTANTS[15] is k_1 in component parameters (per_uM_per_second). * CONSTANTS[16] is k_2 in component parameters (per_second). * CONSTANTS[17] is kI in component parameters (uM). * CONSTANTS[18] is kinh in component parameters (uM). * CONSTANTS[19] is kinhcng in component parameters (uM). * CONSTANTS[20] is n_1 in component parameters (dimensionless). * CONSTANTS[21] is n_2 in component parameters (dimensionless). * CONSTANTS[22] is nI in component parameters (dimensionless). * CONSTANTS[23] is ninh in component parameters (dimensionless). * CONSTANTS[24] is ninhcng in component parameters (dimensionless). * CONSTANTS[25] is pd in component parameters (per_second). * CONSTANTS[26] is r_1 in component parameters (per_second). * CONSTANTS[27] is r_2 in component parameters (per_second). * CONSTANTS[28] is smax in component parameters (uM_per_second). * CONSTANTS[29] is V_Cl in component parameters (millivolt). * CONSTANTS[30] is V_cng in component parameters (millivolt). * CONSTANTS[31] is V_l in component parameters (millivolt). * CONSTANTS[39] is F_vol in component parameters (picocoulomb_per_uM). * CONSTANTS[32] is F in component parameters (coulombs_per_mole). * CONSTANTS[33] is C_vol in component parameters (liter). * ALGEBRAIC[9] is O_stim in component O_stim (uM). * CONSTANTS[34] is od in component O_stim (uM). * CONSTANTS[35] is t_0 in component O_stim (second). * CONSTANTS[36] is t_1 in component O_stim (second). * ALGEBRAIC[0] is H_0 in component O_stim (dimensionless). * ALGEBRAIC[5] is H_1 in component O_stim (dimensionless). * STATES[0] is bLR in component bLR (dimensionless). * CONSTANTS[37] is R_tot in component bLR (dimensionless). * STATES[1] is aG in component aG (dimensionless). * CONSTANTS[38] is G_tot in component aG (dimensionless). * ALGEBRAIC[1] is k_G in component k_G (per_second). * ALGEBRAIC[6] is r_G in component r_G (per_second). * STATES[2] is cAMP in component cAMP (uM). * ALGEBRAIC[2] is synth in component synth (uM_per_second). * ALGEBRAIC[7] is degrad in component degrad (uM_per_second). * STATES[3] is aCaMK in component aCaMK (uM). * STATES[4] is Ca in component Ca (uM). * ALGEBRAIC[10] is I_CNG in component I_CNG (nanoampere). * ALGEBRAIC[12] is J_NCX in component J_NCX (uM_per_second). * ALGEBRAIC[3] is cc_1 in component cc_1 (uM_per_second). * STATES[5] is CaCaM in component CaCaM (uM). * ALGEBRAIC[4] is ck_1 in component ck_1 (uM_per_second). * STATES[6] is V in component V (millivolt). * ALGEBRAIC[11] is I_ClCa in component I_ClCa (nanoampere). * ALGEBRAIC[13] is I_NCX in component I_NCX (nanoampere). * ALGEBRAIC[14] is I_other in component I_other (nanoampere). * ALGEBRAIC[8] is inhcng in component inhcng (dimensionless). * RATES[0] is d/dt bLR in component bLR (dimensionless). * RATES[1] is d/dt aG in component aG (dimensionless). * RATES[2] is d/dt cAMP in component cAMP (uM). * RATES[4] is d/dt Ca in component Ca (uM). * RATES[5] is d/dt CaCaM in component CaCaM (uM). * RATES[3] is d/dt aCaMK in component aCaMK (uM). * RATES[6] is d/dt V in component V (millivolt). */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 0.004; CONSTANTS[1] = 0.88; CONSTANTS[2] = 26; CONSTANTS[3] = 13; CONSTANTS[4] = 0.9; CONSTANTS[5] = 1; CONSTANTS[6] = 1; CONSTANTS[7] = 1; CONSTANTS[8] = 13; CONSTANTS[9] = 2; CONSTANTS[10] = 6; CONSTANTS[11] = 2; CONSTANTS[12] = 3; CONSTANTS[13] = 1.9; CONSTANTS[14] = 5; CONSTANTS[15] = 0.06; CONSTANTS[16] = 20; CONSTANTS[17] = 0.7; CONSTANTS[18] = 2; CONSTANTS[19] = 1; CONSTANTS[20] = 2; CONSTANTS[21] = 2; CONSTANTS[22] = 2; CONSTANTS[23] = 1.5; CONSTANTS[24] = 1.3; CONSTANTS[25] = 20; CONSTANTS[26] = 10; CONSTANTS[27] = 5; CONSTANTS[28] = 71; CONSTANTS[29] = -50; CONSTANTS[30] = 0; CONSTANTS[31] = -70; CONSTANTS[32] = 9.649e4; CONSTANTS[33] = 1e-13; CONSTANTS[34] = 20; CONSTANTS[35] = 0.5; CONSTANTS[36] = 1.5; STATES[0] = 0; CONSTANTS[37] = 1; STATES[1] = 0; CONSTANTS[38] = 1; STATES[2] = 1.35648992164649e-88; STATES[3] = 6.60756525051462e-8; STATES[4] = 5.09073088043779e-12; STATES[5] = 1.86113118246926e-13; STATES[6] = -70; CONSTANTS[39] = (1.00000e+12/1.00000)*(1.00000/1000.00)*CONSTANTS[32]*CONSTANTS[33]; } void computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[3] = CONSTANTS[1]*STATES[4]; RATES[5] = ALGEBRAIC[3] - CONSTANTS[2]*STATES[5]; ALGEBRAIC[4] = CONSTANTS[3]*STATES[5]; RATES[3] = ALGEBRAIC[4] - CONSTANTS[4]*STATES[3]; ALGEBRAIC[1] = CONSTANTS[16]*STATES[0]; ALGEBRAIC[6] = CONSTANTS[27]*STATES[1]; RATES[1] = ALGEBRAIC[1]*(CONSTANTS[38] - STATES[1]) - ALGEBRAIC[6]; ALGEBRAIC[2] = ( STATES[1]*CONSTANTS[28])/(1.00000+pow(STATES[3]/CONSTANTS[18], CONSTANTS[23])); ALGEBRAIC[7] = CONSTANTS[25]*STATES[2]; RATES[2] = ALGEBRAIC[2] - ALGEBRAIC[7]; ALGEBRAIC[0] = (VOI