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 50 entries in the algebraic variable array.
There are a total of 11 entries in each of the rate and state variable arrays.
There are a total of 71 entries in the constant variable array.
*/
/*
* VOI is time in component environment (millisecond).
* ALGEBRAIC[3] is SOVFThick in component sarcomere_geometry (dimensionless).
* ALGEBRAIC[4] is SOVFThin in component sarcomere_geometry (dimensionless).
* ALGEBRAIC[0] is sovr_ze in component sarcomere_geometry (micrometre).
* ALGEBRAIC[1] is sovr_cle in component sarcomere_geometry (micrometre).
* ALGEBRAIC[2] is len_sovr in component sarcomere_geometry (micrometre).
* CONSTANTS[0] is len_thin in component model_parameters (micrometre).
* CONSTANTS[1] is len_thick in component model_parameters (micrometre).
* CONSTANTS[2] is len_hbare in component model_parameters (micrometre).
* STATES[0] is SL in component normalised_active_and_passive_force (micrometre).
* STATES[1] is TRPNCaL in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* STATES[2] is TRPNCaH in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* ALGEBRAIC[5] is dTRPNCaL in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* ALGEBRAIC[6] is dTRPNCaH in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* ALGEBRAIC[42] is kn_pT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* ALGEBRAIC[43] is kp_nT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* CONSTANTS[60] is H in component Ca_binding_to_troponin_to_thin_filament_regulation (micromolar).
* CONSTANTS[62] is H_cons in component Ca_binding_to_troponin_to_thin_filament_regulation (micromolar).
* CONSTANTS[63] is konT in component Ca_binding_to_troponin_to_thin_filament_regulation (second_order_rate_constant).
* CONSTANTS[58] is koffLT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* CONSTANTS[59] is koffHT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* CONSTANTS[3] is Qkon in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[4] is Qkoff in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[5] is Qkn_p in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[6] is Qkp_n in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[7] is kon in component Ca_binding_to_troponin_to_thin_filament_regulation (second_order_rate_constant).
* CONSTANTS[8] is koffL in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* CONSTANTS[9] is koffH in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* CONSTANTS[10] is perm50 in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[11] is nperm in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[12] is kn_p in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* CONSTANTS[13] is kp_n in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
* CONSTANTS[14] is koffmod in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* ALGEBRAIC[7] is Tropreg in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* ALGEBRAIC[36] is permtot in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* ALGEBRAIC[37] is inprmt in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[15] is pH in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[16] is m in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* CONSTANTS[17] is kdHCa in component Ca_binding_to_troponin_to_thin_filament_regulation (micromolar).
* CONSTANTS[18] is TmpC in component model_parameters (celsius).
* CONSTANTS[19] is Cai in component model_parameters (micromolar).
* CONSTANTS[64] is fappT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* ALGEBRAIC[9] is gappT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* ALGEBRAIC[12] is hfT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* ALGEBRAIC[13] is hbT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* ALGEBRAIC[15] is gxbT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* CONSTANTS[20] is fapp in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* CONSTANTS[21] is gapp in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* CONSTANTS[22] is hf in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* CONSTANTS[23] is hb in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* CONSTANTS[24] is gxb in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
* CONSTANTS[25] is gslmod in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* ALGEBRAIC[10] is hfmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* ALGEBRAIC[11] is hbmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[26] is hfmdc in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[27] is hbmdc in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[28] is sigmap in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[29] is sigman in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[30] is xbmodsp in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[31] is Qfapp in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[32] is Qgapp in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[33] is Qhf in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[34] is Qhb in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[35] is Qgxb in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* ALGEBRAIC[14] is gxbmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* ALGEBRAIC[8] is gapslmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
* CONSTANTS[36] is x_0 in component model_parameters (micrometre).
* STATES[3] is xXBpostr in component mean_strain_of_strongly_bound_states (micrometre).
* STATES[4] is xXBprer in component mean_strain_of_strongly_bound_states (micrometre).
* STATES[5] is XBpostr in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* STATES[6] is XBprer in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* ALGEBRAIC[44] is dXBpostr in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* ALGEBRAIC[38] is dXBprer in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* CONSTANTS[65] is alpha1_plus in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* ALGEBRAIC[17] is alpha2_plus in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* ALGEBRAIC[18] is alpha3_plus in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* ALGEBRAIC[19] is alpha1_minus in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* ALGEBRAIC[20] is alpha2_minus in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* ALGEBRAIC[39] is alpha3_minus in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* CONSTANTS[37] is kMgATP in component regulation_and_crossbridge_cycling_state_equations (micromolar2).
* CONSTANTS[38] is kdADP in component regulation_and_crossbridge_cycling_state_equations (micromolar).
* CONSTANTS[39] is xPi_cons in component regulation_and_crossbridge_cycling_state_equations (micromolar).
* CONSTANTS[40] is MgATP_cons in component regulation_and_crossbridge_cycling_state_equations (micromolar).
* ALGEBRAIC[21] is fxbT in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
* STATES[7] is N_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* STATES[8] is P_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* STATES[9] is P in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* ALGEBRAIC[16] is N in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* CONSTANTS[41] is MgADP_cons in component regulation_and_crossbridge_cycling_state_equations (micromolar).
* CONSTANTS[42] is xPi in component regulation_and_crossbridge_cycling_state_equations (micromolar).
* CONSTANTS[43] is MgATP in component regulation_and_crossbridge_cycling_state_equations (micromolar).
* CONSTANTS[44] is MgADP in component regulation_and_crossbridge_cycling_state_equations (micromolar).
* ALGEBRAIC[46] is dxXBpostr in component mean_strain_of_strongly_bound_states (micrometre_per_millisecond).
* ALGEBRAIC[45] is dxXBprer in component mean_strain_of_strongly_bound_states (micrometre_per_millisecond).
* CONSTANTS[45] is xPsi in component mean_strain_of_strongly_bound_states (dimensionless).
* ALGEBRAIC[40] is dutyprer in component mean_strain_of_strongly_bound_states (dimensionless).
* ALGEBRAIC[41] is dutypostr in component mean_strain_of_strongly_bound_states (dimensionless).
* CONSTANTS[61] is dSL in component normalised_active_and_passive_force (micrometre_per_millisecond).
* CONSTANTS[67] is SSXBpostr in component normalised_active_and_passive_force (dimensionless).
* ALGEBRAIC[22] is SSXBprer in component normalised_active_and_passive_force (dimensionless).
* CONSTANTS[46] is kxb in component normalised_active_and_passive_force (millinewton_per_millimetre2).
* CONSTANTS[68] is Fnordv in component normalised_active_and_passive_force (millinewton_micrometre_per_millimetre2).
* ALGEBRAIC[23] is force in component normalised_active_and_passive_force (millinewton_micrometre_per_millimetre2).
* ALGEBRAIC[24] is active in component normalised_active_and_passive_force (unit_normalised_force).
* ALGEBRAIC[27] is ppforce in component normalised_active_and_passive_force (unit_normalised_force).
* ALGEBRAIC[25] is ppforce_t in component normalised_active_and_passive_force (unit_normalised_force).
* ALGEBRAIC[26] is ppforce_c in component normalised_active_and_passive_force (unit_normalised_force).
* CONSTANTS[69] is preload in component normalised_active_and_passive_force (unit_normalised_force).
* ALGEBRAIC[28] is afterload in component normalised_active_and_passive_force (unit_normalised_force).
* STATES[10] is intf in component normalised_active_and_passive_force (unit_normalised_force_millisecond).
* CONSTANTS[47] is SL_c in component normalised_active_and_passive_force (micrometre).
* CONSTANTS[48] is SLrest in component normalised_active_and_passive_force (micrometre).
* CONSTANTS[49] is SLset in component normalised_active_and_passive_force (micrometre).
* CONSTANTS[50] is PCon_t in component normalised_active_and_passive_force (unit_normalised_force).
* CONSTANTS[51] is PExp_t in component normalised_active_and_passive_force (per_micrometre).
* CONSTANTS[52] is PCon_c in component normalised_active_and_passive_force (unit_normalised_force).
* CONSTANTS[53] is PExp_c in component normalised_active_and_passive_force (per_micrometre).
* CONSTANTS[54] is KSE in component normalised_active_and_passive_force (unit_normalised_force_per_micrometre).
* CONSTANTS[66] is fxb in component normalised_active_and_passive_force (first_order_rate_constant).
* CONSTANTS[55] is SEon in component normalised_active_and_passive_force (dimensionless).
* ALGEBRAIC[29] is FrSBXB in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (dimensionless).
* ALGEBRAIC[47] is dFrSBXB in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (first_order_rate_constant).
* ALGEBRAIC[30] is dsovr_ze in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micrometre_per_millisecond).
* ALGEBRAIC[31] is dsovr_cle in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micrometre_per_millisecond).
* ALGEBRAIC[32] is dlen_sovr in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micrometre_per_millisecond).
* ALGEBRAIC[34] is dSOVFThick in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (first_order_rate_constant).
* ALGEBRAIC[33] is dSOVFThin in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (first_order_rate_constant).
* CONSTANTS[56] is kxb in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (millinewton_per_millimetre2).
* ALGEBRAIC[49] is dforce in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (millinewton_micrometre_per_millimetre2_per_millisecond).
* CONSTANTS[57] is Trop_conc in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micromolar).
* ALGEBRAIC[35] is TropTot in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micromolar).
* ALGEBRAIC[48] is dTropTot in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micromolar_per_millisecond).
* RATES[1] is d/dt TRPNCaL in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* RATES[2] is d/dt TRPNCaH in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
* RATES[7] is d/dt N_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* RATES[8] is d/dt P_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* RATES[9] is d/dt P in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* RATES[6] is d/dt XBprer in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* RATES[5] is d/dt XBpostr in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
* RATES[4] is d/dt xXBprer in component mean_strain_of_strongly_bound_states (micrometre).
* RATES[3] is d/dt xXBpostr in component mean_strain_of_strongly_bound_states (micrometre).
* RATES[0] is d/dt SL in component normalised_active_and_passive_force (micrometre).
* RATES[10] is d/dt intf in component normalised_active_and_passive_force (unit_normalised_force_millisecond).
* There are a total of 7 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1.2;
CONSTANTS[1] = 1.65;
CONSTANTS[2] = 0.1;
STATES[0] = 2.2;
STATES[1] = 0.0147730085063734;
STATES[2] = 0.13066096561522;
CONSTANTS[3] = 1.5;
CONSTANTS[4] = 1.3;
CONSTANTS[5] = 1.6;
CONSTANTS[6] = 1.6;
CONSTANTS[7] = 0.05;
CONSTANTS[8] = 0.25;
CONSTANTS[9] = 0.025;
CONSTANTS[10] = 0.5;
CONSTANTS[11] = 15;
CONSTANTS[12] = 0.5;
CONSTANTS[13] = 0.05;
CONSTANTS[14] = 1;
CONSTANTS[15] = 7.15;
CONSTANTS[16] = 1;
CONSTANTS[17] = 2e-2;
CONSTANTS[18] = 22;
CONSTANTS[19] = 200.0;
CONSTANTS[20] = 0.5;
CONSTANTS[21] = 0.07;
CONSTANTS[22] = 2;
CONSTANTS[23] = 0.4;
CONSTANTS[24] = 0.07;
CONSTANTS[25] = 6;
CONSTANTS[26] = 5;
CONSTANTS[27] = 0;
CONSTANTS[28] = 8;
CONSTANTS[29] = 1;
CONSTANTS[30] = 0.2;
CONSTANTS[31] = 6.25;
CONSTANTS[32] = 2.5;
CONSTANTS[33] = 6.25;
CONSTANTS[34] = 6.25;
CONSTANTS[35] = 6.25;
CONSTANTS[36] = 0.007;
STATES[3] = 0.00700005394873882;
STATES[4] = 3.41212828972468e-8;
STATES[5] = 1.81017564383744e-6;
STATES[6] = 3.0494964880038e-7;
CONSTANTS[37] = 15400e6;
CONSTANTS[38] = 4;
CONSTANTS[39] = 2e3;
CONSTANTS[40] = 5e3;
STATES[7] = 0.999999959256274;
STATES[8] = 4.07437173988636e-8;
STATES[9] = 0.999997834540066;
CONSTANTS[41] = 36;
CONSTANTS[42] = 2e3;
CONSTANTS[43] = 5e3;
CONSTANTS[44] = 36.3;
CONSTANTS[45] = 2;
CONSTANTS[46] = 120;
STATES[10] = -4.5113452510363e-6;
CONSTANTS[47] = 2.25;
CONSTANTS[48] = 1.85;
CONSTANTS[49] = 1.9;
CONSTANTS[50] = 0.002;
CONSTANTS[51] = 10;
CONSTANTS[52] = 0.02;
CONSTANTS[53] = 70;
CONSTANTS[54] = 1;
CONSTANTS[55] = 1;
CONSTANTS[56] = 120;
CONSTANTS[57] = 70;
CONSTANTS[58] = CONSTANTS[8]*CONSTANTS[14]*pow(CONSTANTS[4], (CONSTANTS[18] - 37.0000)/10.0000);
CONSTANTS[59] = CONSTANTS[9]*CONSTANTS[14]*pow(CONSTANTS[4], (CONSTANTS[18] - 37.0000)/10.0000);
CONSTANTS[60] = 1.00000e+07*pow(10.0000, - CONSTANTS[15]);
CONSTANTS[61] = 0.00000;
CONSTANTS[62] = 1.00000e+07*pow(10.0000, -7.15000);
CONSTANTS[70] = CONSTANTS[61];
CONSTANTS[63] = ((pow(CONSTANTS[17], CONSTANTS[16])+pow(CONSTANTS[62], CONSTANTS[16]))/(pow(CONSTANTS[17], CONSTANTS[16])+pow(CONSTANTS[60], CONSTANTS[16])))*( CONSTANTS[7]*pow(CONSTANTS[3], (CONSTANTS[18] - 37.0000)/10.0000));
CONSTANTS[64] = CONSTANTS[20]*CONSTANTS[30]*pow(CONSTANTS[31], (CONSTANTS[18] - 37.0000)/10.0000);
CONSTANTS[65] = CONSTANTS[64];
CONSTANTS[66] = ( CONSTANTS[38]*CONSTANTS[20]*CONSTANTS[22]*(CONSTANTS[24]/CONSTANTS[40]))/( (CONSTANTS[21]/CONSTANTS[39])*(CONSTANTS[23]/CONSTANTS[62])*CONSTANTS[37]);
CONSTANTS[67] = ( CONSTANTS[20]*CONSTANTS[22]+ CONSTANTS[66]*CONSTANTS[21]+ CONSTANTS[66]*CONSTANTS[23])/( CONSTANTS[22]*CONSTANTS[24]+ CONSTANTS[23]*CONSTANTS[21]+ CONSTANTS[24]*CONSTANTS[21]+ CONSTANTS[66]*CONSTANTS[23]+ CONSTANTS[24]*CONSTANTS[20]+ CONSTANTS[23]*CONSTANTS[20]+ CONSTANTS[66]*CONSTANTS[21]+ CONSTANTS[20]*CONSTANTS[22]+ CONSTANTS[66]*CONSTANTS[23]);
CONSTANTS[68] = CONSTANTS[46]*CONSTANTS[36]*CONSTANTS[67];
CONSTANTS[69] = (fabs(CONSTANTS[49] - CONSTANTS[48])/(CONSTANTS[49] - CONSTANTS[48]))*CONSTANTS[50]*(exp( CONSTANTS[51]*fabs(CONSTANTS[49] - CONSTANTS[48])) - 1.00000);
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[6] = 0.1001;
RATES[5] = 0.1001;
RATES[4] = 0.1001;
RATES[3] = 0.1001;
RATES[10] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[1] - ALGEBRAIC[5];
resid[1] = RATES[2] - ALGEBRAIC[6];
resid[2] = RATES[7] - ALGEBRAIC[43]*STATES[8] - ALGEBRAIC[42]*STATES[7];
resid[3] = RATES[8] - ALGEBRAIC[42]*STATES[7] - ALGEBRAIC[43]*STATES[8];
resid[4] = RATES[9] - ALGEBRAIC[42]*ALGEBRAIC[16] - (ALGEBRAIC[43]+CONSTANTS[65])*STATES[9];
resid[5] = RATES[6] - ALGEBRAIC[38];
resid[6] = RATES[5] - ALGEBRAIC[44];
resid[7] = RATES[4] - ALGEBRAIC[45];
resid[8] = RATES[3] - ALGEBRAIC[46];
resid[9] = RATES[10] - (CONSTANTS[69]+ALGEBRAIC[28]) - (ALGEBRAIC[27]+ALGEBRAIC[24]);
resid[10] = RATES[0] - CONSTANTS[70];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[22] = ( CONSTANTS[66]*CONSTANTS[23]+ CONSTANTS[24]*CONSTANTS[20]+ ALGEBRAIC[13]*CONSTANTS[20])/( CONSTANTS[22]*CONSTANTS[24]+ CONSTANTS[23]*CONSTANTS[21]+ CONSTANTS[24]*CONSTANTS[21]+ CONSTANTS[66]*CONSTANTS[23]+ CONSTANTS[24]*CONSTANTS[20]+ CONSTANTS[23]*CONSTANTS[20]+ CONSTANTS[66]*CONSTANTS[21]+ CONSTANTS[20]*CONSTANTS[22]+ CONSTANTS[66]*CONSTANTS[23]);
ALGEBRAIC[29] = (STATES[5]+STATES[6])/(CONSTANTS[67]+ALGEBRAIC[22]);
ALGEBRAIC[30] = (CONDVAR[5]<0.00000 ? - 0.500000*CONSTANTS[61] : 0.00000);
ALGEBRAIC[31] = (CONDVAR[6]>0.00000 ? - 0.500000*CONSTANTS[61] : 0.00000);
ALGEBRAIC[32] = ALGEBRAIC[30] - ALGEBRAIC[31];
ALGEBRAIC[33] = ALGEBRAIC[32]/CONSTANTS[0];
ALGEBRAIC[34] = ( 2.00000*ALGEBRAIC[32])/(CONSTANTS[1] - CONSTANTS[2]);
ALGEBRAIC[35] = CONSTANTS[57]*( (1.00000 - ALGEBRAIC[4])*STATES[1]+ ALGEBRAIC[4]*( ALGEBRAIC[29]*STATES[2]+ (1.00000 - ALGEBRAIC[29])*STATES[1]));
ALGEBRAIC[47] = (ALGEBRAIC[44]+ALGEBRAIC[38])/(CONSTANTS[67]+ALGEBRAIC[22]);
ALGEBRAIC[48] = CONSTANTS[57]*( - ALGEBRAIC[33]*STATES[1]+ (1.00000 - ALGEBRAIC[4])*ALGEBRAIC[5]+ ALGEBRAIC[33]*( ALGEBRAIC[29]*STATES[2]+ (1.00000 - ALGEBRAIC[29])*STATES[1])+ ALGEBRAIC[4]*(( ALGEBRAIC[47]*STATES[2]+ ALGEBRAIC[29]*ALGEBRAIC[6]+ (1.00000 - ALGEBRAIC[29])*ALGEBRAIC[5]) - ALGEBRAIC[47]*STATES[1]));
ALGEBRAIC[49] = CONSTANTS[56]*ALGEBRAIC[34]*( STATES[3]*STATES[5]+ STATES[4]*STATES[6])+ CONSTANTS[56]*ALGEBRAIC[3]*( ALGEBRAIC[46]*STATES[5]+ STATES[3]*ALGEBRAIC[44]+ ALGEBRAIC[45]*STATES[6]+ STATES[4]*ALGEBRAIC[38]);
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[5] = CONSTANTS[63]*CONSTANTS[19]*(1.00000 - STATES[1]) - CONSTANTS[58]*STATES[1];
ALGEBRAIC[6] = CONSTANTS[63]*CONSTANTS[19]*(1.00000 - STATES[2]) - CONSTANTS[59]*STATES[2];
ALGEBRAIC[16] = 1.00000 - (STATES[9]+STATES[6]+STATES[5]);
ALGEBRAIC[0] = (CONDVAR[0]<0.00000 ? CONSTANTS[1]/2.00000 : STATES[0]/2.00000);
ALGEBRAIC[1] = (CONDVAR[1]>0.00000 ? STATES[0]/2.00000 - (STATES[0] - CONSTANTS[0]) : CONSTANTS[2]/2.00000);
ALGEBRAIC[2] = ALGEBRAIC[0] - ALGEBRAIC[1];
ALGEBRAIC[3] = ( ALGEBRAIC[2]*2.00000)/(CONSTANTS[1] - CONSTANTS[2]);
ALGEBRAIC[23] = CONSTANTS[46]*ALGEBRAIC[3]*( STATES[3]*STATES[5]+ STATES[4]*STATES[6]);
ALGEBRAIC[24] = ( 1.00000*ALGEBRAIC[23])/CONSTANTS[68];
ALGEBRAIC[25] = ((STATES[0] - CONSTANTS[48])/fabs(STATES[0] - CONSTANTS[48]))*CONSTANTS[50]*(exp( CONSTANTS[51]*fabs(STATES[0] - CONSTANTS[48])) - 1.00000);
ALGEBRAIC[26] = (CONDVAR[4]>0.00000 ? CONSTANTS[52]*(exp( CONSTANTS[53]*fabs(STATES[0] - CONSTANTS[47])) - 1.00000) : 0.00000);
ALGEBRAIC[27] = ALGEBRAIC[25]+ALGEBRAIC[26];
ALGEBRAIC[28] = (CONSTANTS[55]==1.00000 ? CONSTANTS[54]*(CONSTANTS[49] - STATES[0]) : 0.00000);
ALGEBRAIC[10] = exp( (- STATES[4]/fabs(STATES[4]))*CONSTANTS[26]*pow(STATES[4]/CONSTANTS[36], 2.00000));
ALGEBRAIC[12] = CONSTANTS[22]*ALGEBRAIC[10]*CONSTANTS[30]*pow(CONSTANTS[33], (CONSTANTS[18] - 37.0000)/10.0000);
ALGEBRAIC[17] = ALGEBRAIC[12];
ALGEBRAIC[8] = 1.00000+ (1.00000 - ALGEBRAIC[3])*CONSTANTS[25];
ALGEBRAIC[9] = CONSTANTS[21]*ALGEBRAIC[8]*CONSTANTS[30]*pow(CONSTANTS[32], (CONSTANTS[18] - 37.0000)/10.0000);
ALGEBRAIC[19] = CONSTANTS[42]*(ALGEBRAIC[9]/CONSTANTS[39]);
ALGEBRAIC[11] = exp( ((STATES[3] - CONSTANTS[36])/fabs(STATES[3] - CONSTANTS[36]))*CONSTANTS[27]*pow((STATES[3] - CONSTANTS[36])/CONSTANTS[36], 2.00000));
ALGEBRAIC[13] = CONSTANTS[23]*ALGEBRAIC[11]*CONSTANTS[30]*pow(CONSTANTS[34], (CONSTANTS[18] - 37.0000)/10.0000);
ALGEBRAIC[20] = CONSTANTS[60]*(ALGEBRAIC[13]/CONSTANTS[62])*((CONSTANTS[38]+CONSTANTS[41])/CONSTANTS[41])*(CONSTANTS[44]/(CONSTANTS[38]+CONSTANTS[44]));
ALGEBRAIC[38] = ( CONSTANTS[65]*STATES[9]+ ALGEBRAIC[20]*STATES[5]) - (ALGEBRAIC[19]+ALGEBRAIC[17])*STATES[6];
ALGEBRAIC[4] = ALGEBRAIC[2]/CONSTANTS[0];
ALGEBRAIC[7] = (1.00000 - ALGEBRAIC[4])*STATES[1]+ ALGEBRAIC[4]*STATES[2];
ALGEBRAIC[36] = pow(fabs(1.00000/(1.00000+pow(CONSTANTS[10]/ALGEBRAIC[7], CONSTANTS[11]))), 1.0 / 2);
ALGEBRAIC[42] = CONSTANTS[12]*ALGEBRAIC[36]*pow(CONSTANTS[5], (CONSTANTS[18] - 37.0000)/10.0000);
ALGEBRAIC[37] = (CONDVAR[2]<0.00000 ? 1.00000/ALGEBRAIC[36] : 100.000);
ALGEBRAIC[43] = CONSTANTS[13]*ALGEBRAIC[37]*pow(CONSTANTS[6], (CONSTANTS[18] - 37.0000)/10.0000);
ALGEBRAIC[14] = (CONDVAR[3]<0.00000 ? exp( CONSTANTS[28]*pow((CONSTANTS[36] - STATES[3])/CONSTANTS[36], 2.00000)) : exp( CONSTANTS[29]*pow((STATES[3] - CONSTANTS[36])/CONSTANTS[36], 2.00000)));
ALGEBRAIC[15] = CONSTANTS[24]*ALGEBRAIC[14]*CONSTANTS[30]*pow(CONSTANTS[35], (CONSTANTS[18] - 37.0000)/10.0000);
ALGEBRAIC[18] = CONSTANTS[43]*(ALGEBRAIC[15]/CONSTANTS[40])*((CONSTANTS[38]+CONSTANTS[41])/(CONSTANTS[38]+CONSTANTS[44]));
ALGEBRAIC[21] = ( CONSTANTS[38]*CONSTANTS[64]*ALGEBRAIC[12]*(ALGEBRAIC[15]/CONSTANTS[40]))/( (ALGEBRAIC[9]/CONSTANTS[39])*(ALGEBRAIC[13]/CONSTANTS[62])*CONSTANTS[37]);
ALGEBRAIC[39] = ALGEBRAIC[21];
ALGEBRAIC[44] = ( ALGEBRAIC[39]*STATES[9]+ ALGEBRAIC[17]*STATES[6]) - (ALGEBRAIC[20]+ALGEBRAIC[18])*STATES[5];
ALGEBRAIC[40] = ( ALGEBRAIC[39]*ALGEBRAIC[20]+ ALGEBRAIC[18]*CONSTANTS[65]+ ALGEBRAIC[20]*CONSTANTS[65])/( CONSTANTS[65]*ALGEBRAIC[17]+ ALGEBRAIC[39]*ALGEBRAIC[19]+ ALGEBRAIC[39]*ALGEBRAIC[17]+ ALGEBRAIC[39]*ALGEBRAIC[20]+ ALGEBRAIC[18]*CONSTANTS[65]+ ALGEBRAIC[20]*CONSTANTS[65]+ ALGEBRAIC[17]*ALGEBRAIC[18]+ ALGEBRAIC[39]*ALGEBRAIC[19]+ ALGEBRAIC[18]*ALGEBRAIC[19]);
ALGEBRAIC[45] = CONSTANTS[61]/2.00000+ (CONSTANTS[45]/ALGEBRAIC[40])*(- ( CONSTANTS[65]*STATES[4])+ ALGEBRAIC[20]*(STATES[3] - (CONSTANTS[36]+STATES[4])));
ALGEBRAIC[41] = ( CONSTANTS[65]*ALGEBRAIC[17]+ ALGEBRAIC[39]*ALGEBRAIC[19]+ ALGEBRAIC[39]*ALGEBRAIC[17])/( CONSTANTS[65]*ALGEBRAIC[17]+ ALGEBRAIC[39]*ALGEBRAIC[19]+ ALGEBRAIC[39]*ALGEBRAIC[17]+ ALGEBRAIC[39]*ALGEBRAIC[20]+ ALGEBRAIC[18]*CONSTANTS[65]+ ALGEBRAIC[20]*CONSTANTS[65]+ ALGEBRAIC[17]*ALGEBRAIC[18]+ ALGEBRAIC[39]*ALGEBRAIC[19]+ ALGEBRAIC[18]*ALGEBRAIC[19]);
ALGEBRAIC[46] = CONSTANTS[61]/2.00000+ (CONSTANTS[45]/ALGEBRAIC[41])*( ALGEBRAIC[17]*((STATES[4]+CONSTANTS[36]) - STATES[3]));
}
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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = CONSTANTS[1]/2.00000 - STATES[0]/2.00000;
CONDVAR[1] = (STATES[0]/2.00000 - (STATES[0] - CONSTANTS[0])) - CONSTANTS[2]/2.00000;
CONDVAR[2] = 1.00000/ALGEBRAIC[36] - 100.000;
CONDVAR[3] = STATES[3] - CONSTANTS[36];
CONDVAR[4] = STATES[0] - CONSTANTS[47];
CONDVAR[5] = STATES[0] - CONSTANTS[1];
CONDVAR[6] = ( 2.00000*CONSTANTS[0] - STATES[0]) - CONSTANTS[2];
}