/* There are a total of 11 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 31 entries in the constant variable array. */ /* * VOI is time in component environment (second). * STATES[0] is ATP in component ATP (micromolar). * CONSTANTS[0] is V_hyd in component V_hyd (flux). * ALGEBRAIC[1] is V_MMCK in component V_MMCK (flux). * ALGEBRAIC[8] is J_ATP in component J_ATP (flux). * CONSTANTS[1] is V_cyt in component fractional_volumes (dimensionless). * STATES[1] is ADP in component ADP (micromolar). * ALGEBRAIC[9] is J_ADP in component J_ADP (flux). * STATES[2] is PCr in component PCr (micromolar). * ALGEBRAIC[5] is J_PCr in component J_PCr (flux). * STATES[3] is Cr in component Cr (micromolar). * ALGEBRAIC[6] is J_Cr in component J_Cr (flux). * STATES[4] is Pi in component Pi (micromolar). * ALGEBRAIC[10] is J_Pi in component J_Pi (flux). * STATES[5] is ATP_i in component ATP_i (micromolar). * ALGEBRAIC[3] is V_MiCK in component V_MiCK (flux). * ALGEBRAIC[7] is V_syn in component V_syn (flux). * CONSTANTS[2] is V_ims in component fractional_volumes (dimensionless). * STATES[6] is ADP_i in component ADP_i (micromolar). * STATES[7] is PCr_i in component PCr_i (micromolar). * STATES[8] is Cr_i in component Cr_i (micromolar). * STATES[9] is Pi_i in component Pi_i (micromolar). * ALGEBRAIC[0] is Den_MMCK in component V_MMCK (dimensionless). * CONSTANTS[3] is Kia in component V_MMCK (micromolar). * CONSTANTS[4] is Kb in component V_MMCK (micromolar). * CONSTANTS[5] is Kib in component V_MMCK (micromolar). * CONSTANTS[27] is KIb in component V_MMCK (micromolar). * CONSTANTS[28] is Kc in component V_MMCK (micromolar). * CONSTANTS[6] is Kic in component V_MMCK (micromolar). * CONSTANTS[7] is Kd in component V_MMCK (micromolar). * CONSTANTS[8] is Kid in component V_MMCK (micromolar). * CONSTANTS[9] is Vf in component V_MMCK (flux). * CONSTANTS[10] is Vb in component V_MMCK (flux). * ALGEBRAIC[2] is Den_MiCK in component V_MiCK (dimensionless). * CONSTANTS[11] is Kia in component V_MiCK (micromolar). * CONSTANTS[12] is Kb in component V_MiCK (micromolar). * CONSTANTS[13] is Kib in component V_MiCK (micromolar). * CONSTANTS[29] is KIb in component V_MiCK (micromolar). * CONSTANTS[30] is Kc in component V_MiCK (micromolar). * CONSTANTS[14] is Kic in component V_MiCK (micromolar). * CONSTANTS[15] is Kd in component V_MiCK (micromolar). * CONSTANTS[16] is Kid in component V_MiCK (micromolar). * CONSTANTS[17] is Vf in component V_MiCK (flux). * CONSTANTS[18] is Vb in component V_MiCK (flux). * ALGEBRAIC[4] is Den_syn in component V_syn (dimensionless). * CONSTANTS[19] is KPi in component V_syn (micromolar). * CONSTANTS[20] is KADP in component V_syn (micromolar). * CONSTANTS[21] is V_syn_max in component V_syn (flux). * CONSTANTS[22] is R_ATP in component J_ATP (first_order_rate_constant). * CONSTANTS[23] is R_ADP in component J_ADP (first_order_rate_constant). * CONSTANTS[24] is R_PCr in component J_PCr (first_order_rate_constant). * CONSTANTS[25] is R_Cr in component J_Cr (first_order_rate_constant). * CONSTANTS[26] is R_Pi in component J_Pi (first_order_rate_constant). * RATES[0] is d/dt ATP in component ATP (micromolar). * RATES[1] is d/dt ADP in component ADP (micromolar). * RATES[2] is d/dt PCr in component PCr (micromolar). * RATES[3] is d/dt Cr in component Cr (micromolar). * RATES[4] is d/dt Pi in component Pi (micromolar). * RATES[5] is d/dt ATP_i in component ATP_i (micromolar). * RATES[6] is d/dt ADP_i in component ADP_i (micromolar). * RATES[7] is d/dt PCr_i in component PCr_i (micromolar). * RATES[8] is d/dt Cr_i in component Cr_i (micromolar). * RATES[9] is d/dt Pi_i in component Pi_i (micromolar). */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = 9644.425; CONSTANTS[0] = 4.6E3; CONSTANTS[1] = 0.75; STATES[1] = 60.0; STATES[2] = 12500.0; STATES[3] = 13500.0; STATES[4] = 8000.0; STATES[5] = 9644.425; CONSTANTS[2] = 0.0625; STATES[6] = 2.5; STATES[7] = 12500.0; STATES[8] = 13500.0; STATES[9] = 8000.0; CONSTANTS[3] = 9.0E2; CONSTANTS[4] = 1.55E4; CONSTANTS[5] = 3.49E4; CONSTANTS[6] = 2.224E2; CONSTANTS[7] = 1.67E3; CONSTANTS[8] = 4.73E3; CONSTANTS[9] = 6.966E3; CONSTANTS[10] = 2.925E4; CONSTANTS[11] = 7.5E2; CONSTANTS[12] = 5.2E3; CONSTANTS[13] = 2.88E4; CONSTANTS[14] = 2.048E2; CONSTANTS[15] = 5.0E2; CONSTANTS[16] = 1.6E3; CONSTANTS[17] = 2.658E3; CONSTANTS[18] = 1.116E4; CONSTANTS[19] = 20.0; CONSTANTS[20] = 8.0E2; CONSTANTS[21] = 4.6E3; CONSTANTS[22] = 8.16; CONSTANTS[23] = 8.16; CONSTANTS[24] = 14.6; CONSTANTS[25] = 14.6; CONSTANTS[26] = 18.4; CONSTANTS[27] = CONSTANTS[5]; CONSTANTS[28] = ( CONSTANTS[6]*CONSTANTS[7])/CONSTANTS[8]; CONSTANTS[29] = CONSTANTS[13]; CONSTANTS[30] = ( CONSTANTS[14]*CONSTANTS[15])/CONSTANTS[16]; } void computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = 1.00000+STATES[3]/CONSTANTS[5]+STATES[2]/CONSTANTS[8]+ STATES[0]*(1.00000/CONSTANTS[3]+STATES[3]/( CONSTANTS[3]*CONSTANTS[4]))+ STATES[1]*(1.00000/CONSTANTS[6]+STATES[2]/( CONSTANTS[8]*CONSTANTS[28])+STATES[3]/( CONSTANTS[6]*CONSTANTS[27])); ALGEBRAIC[1] = ( CONSTANTS[9]*(( STATES[0]*STATES[3])/( CONSTANTS[3]*CONSTANTS[4])) - CONSTANTS[10]*(( STATES[1]*STATES[2])/( CONSTANTS[6]*CONSTANTS[7])))/ALGEBRAIC[0]; ALGEBRAIC[5] = CONSTANTS[24]*(STATES[7] - STATES[2]); RATES[2] = (ALGEBRAIC[5]+ALGEBRAIC[1])/CONSTANTS[1]; ALGEBRAIC[6] = CONSTANTS[25]*(STATES[8] - STATES[3]); RATES[3] = (ALGEBRAIC[6] - ALGEBRAIC[1])/CONSTANTS[1]; ALGEBRAIC[2] = 1.00000+STATES[8]/CONSTANTS[13]+STATES[7]/CONSTANTS[16]+ STATES[5]*(1.00000/CONSTANTS[11]+STATES[8]/( CONSTANTS[11]*CONSTANTS[12]))+ STATES[6]*(1.00000/CONSTANTS[14]+STATES[7]/( CONSTANTS[16]*CONSTANTS[30])+STATES[8]/( CONSTANTS[14]*CONSTANTS[29])); ALGEBRAIC[3] = ( CONSTANTS[17]*(( STATES[5]*STATES[8])/( CONSTANTS[11]*CONSTANTS[12])) - CONSTANTS[18]*(( STATES[6]*STATES[7])/( CONSTANTS[14]*CONSTANTS[15])))/ALGEBRAIC[2]; RATES[7] = (ALGEBRAIC[3] - ALGEBRAIC[5])/CONSTANTS[2]; RATES[8] = - (ALGEBRAIC[3]+ALGEBRAIC[6])/CONSTANTS[2]; ALGEBRAIC[8] = CONSTANTS[22]*(STATES[5] - STATES[0]); RATES[0] = (ALGEBRAIC[8] - (CONSTANTS[0]+ALGEBRAIC[1]))/CONSTANTS[1]; ALGEBRAIC[9] = CONSTANTS[23]*(STATES[6] - STATES[1]); RATES[1] = (ALGEBRAIC[9]+CONSTANTS[0]+ALGEBRAIC[1])/CONSTANTS[1]; ALGEBRAIC[10] = CONSTANTS[26]*(STATES[9] - STATES[4]); RATES[4] = (ALGEBRAIC[10]+CONSTANTS[0])/CONSTANTS[1]; ALGEBRAIC[4] = 1.00000+STATES[6]/CONSTANTS[20]+STATES[9]/CONSTANTS[19]+( STATES[6]*STATES[9])/( CONSTANTS[20]*CONSTANTS[19]); ALGEBRAIC[7] = CONSTANTS[21]*(( STATES[6]*STATES[9])/( CONSTANTS[19]*CONSTANTS[20]*ALGEBRAIC[4])); RATES[5] = - (ALGEBRAIC[8]+ALGEBRAIC[7]+ALGEBRAIC[3])/CONSTANTS[2]; RATES[6] = ((ALGEBRAIC[7]+ALGEBRAIC[3]) - ALGEBRAIC[9])/CONSTANTS[2]; RATES[9] = (ALGEBRAIC[7] - ALGEBRAIC[10])/CONSTANTS[2]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = 1.00000+STATES[3]/CONSTANTS[5]+STATES[2]/CONSTANTS[8]+ STATES[0]*(1.00000/CONSTANTS[3]+STATES[3]/( CONSTANTS[3]*CONSTANTS[4]))+ STATES[1]*(1.00000/CONSTANTS[6]+STATES[2]/( CONSTANTS[8]*CONSTANTS[28])+STATES[3]/( CONSTANTS[6]*CONSTANTS[27])); ALGEBRAIC[1] = ( CONSTANTS[9]*(( STATES[0]*STATES[3])/( CONSTANTS[3]*CONSTANTS[4])) - CONSTANTS[10]*(( STATES[1]*STATES[2])/( CONSTANTS[6]*CONSTANTS[7])))/ALGEBRAIC[0]; ALGEBRAIC[5] = CONSTANTS[24]*(STATES[7] - STATES[2]); ALGEBRAIC[6] = CONSTANTS[25]*(STATES[8] - STATES[3]); ALGEBRAIC[2] = 1.00000+STATES[8]/CONSTANTS[13]+STATES[7]/CONSTANTS[16]+ STATES[5]*(1.00000/CONSTANTS[11]+STATES[8]/( CONSTANTS[11]*CONSTANTS[12]))+ STATES[6]*(1.00000/CONSTANTS[14]+STATES[7]/( CONSTANTS[16]*CONSTANTS[30])+STATES[8]/( CONSTANTS[14]*CONSTANTS[29])); ALGEBRAIC[3] = ( CONSTANTS[17]*(( STATES[5]*STATES[8])/( CONSTANTS[11]*CONSTANTS[12])) - CONSTANTS[18]*(( STATES[6]*STATES[7])/( CONSTANTS[14]*CONSTANTS[15])))/ALGEBRAIC[2]; ALGEBRAIC[8] = CONSTANTS[22]*(STATES[5] - STATES[0]); ALGEBRAIC[9] = CONSTANTS[23]*(STATES[6] - STATES[1]); ALGEBRAIC[10] = CONSTANTS[26]*(STATES[9] - STATES[4]); ALGEBRAIC[4] = 1.00000+STATES[6]/CONSTANTS[20]+STATES[9]/CONSTANTS[19]+( STATES[6]*STATES[9])/( CONSTANTS[20]*CONSTANTS[19]); ALGEBRAIC[7] = CONSTANTS[21]*(( STATES[6]*STATES[9])/( CONSTANTS[19]*CONSTANTS[20]*ALGEBRAIC[4])); }