Generated Code
The following is c code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
There are a total of 42 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 47 entries in the constant variable array.
*/
/*
* VOI is time in component environment (millisecond).
* STATES[0] is Gi in component Gi (millimolar).
* ALGEBRAIC[4] is Jglut in component Jglut (millimolar_per_millisecond).
* ALGEBRAIC[6] is Jgk in component Jgk (millimolar_per_millisecond).
* STATES[1] is Ge in component Ge (millimolar).
* ALGEBRAIC[1] is G_infinity in component Ge (millimolar).
* CONSTANTS[0] is tau_G in component Ge (millisecond).
* CONSTANTS[1] is Gmax in component Ge (millimolar).
* CONSTANTS[2] is Gmin in component Ge (millimolar).
* CONSTANTS[3] is SG in component Ge (dimensionless).
* CONSTANTS[4] is I_ in component Ge (dimensionless).
* STATES[2] is I in component I (dimensionless).
* STATES[3] is G6P in component G6P (millimolar).
* CONSTANTS[5] is kappa in component model_parameters (dimensionless).
* ALGEBRAIC[41] is JPFK in component JPFK (micromolar_per_millisecond).
* STATES[4] is FBP in component FBP (micromolar).
* ALGEBRAIC[7] is JGPDH in component JGPDH (micromolar_per_millisecond).
* ALGEBRAIC[0] is F6P in component F6P (millimolar).
* CONSTANTS[6] is Kglut in component Jglut (millimolar).
* CONSTANTS[7] is Vglut in component Jglut (millimolar_per_millisecond).
* CONSTANTS[8] is Kgk in component Jgk (millimolar).
* CONSTANTS[9] is Vgk in component Jgk (millimolar_per_millisecond).
* CONSTANTS[10] is ngk in component Jgk (dimensionless).
* CONSTANTS[11] is lambda in component JPFK (dimensionless).
* CONSTANTS[12] is Vmax in component JPFK (micromolar_per_millisecond).
* CONSTANTS[46] is w0000 in component w (dimensionless).
* ALGEBRAIC[25] is w1000 in component w (dimensionless).
* ALGEBRAIC[8] is w0100 in component w (dimensionless).
* ALGEBRAIC[9] is w0010 in component w (dimensionless).
* ALGEBRAIC[27] is w0001 in component w (dimensionless).
* ALGEBRAIC[29] is w1100 in component w (dimensionless).
* ALGEBRAIC[31] is w1010 in component w (dimensionless).
* ALGEBRAIC[33] is w1001 in component w (dimensionless).
* ALGEBRAIC[34] is w0101 in component w (dimensionless).
* ALGEBRAIC[35] is w0011 in component w (dimensionless).
* ALGEBRAIC[10] is w0110 in component w (dimensionless).
* ALGEBRAIC[36] is w1110 in component w (dimensionless).
* ALGEBRAIC[37] is w0111 in component w (dimensionless).
* ALGEBRAIC[38] is w1011 in component w (dimensionless).
* ALGEBRAIC[39] is w1101 in component w (dimensionless).
* ALGEBRAIC[40] is w1111 in component w (dimensionless).
* CONSTANTS[13] is f13 in component w (dimensionless).
* CONSTANTS[14] is f23 in component w (dimensionless).
* CONSTANTS[15] is f41 in component w (dimensionless).
* CONSTANTS[16] is f42 in component w (dimensionless).
* CONSTANTS[17] is f43 in component w (dimensionless).
* CONSTANTS[18] is K1 in component w (micromolar).
* CONSTANTS[19] is K2 in component w (micromolar).
* CONSTANTS[20] is K3 in component w (micromolar).
* CONSTANTS[21] is K4 in component w (micromolar).
* ALGEBRAIC[23] is AMP in component AMP (micromolar).
* ALGEBRAIC[24] is ATP in component ATP (micromolar).
* CONSTANTS[22] is Atot in component ATP (micromolar).
* STATES[5] is ADP in component ADP (micromolar).
* CONSTANTS[23] is tau_a in component ADP (millisecond).
* CONSTANTS[24] is r in component ADP (dimensionless).
* CONSTANTS[25] is r1 in component ADP (micromolar).
* ALGEBRAIC[11] is gamma in component gamma (dimensionless).
* STATES[6] is Ca in component Ca (micromolar).
* CONSTANTS[26] is v_gamma in component gamma (dimensionless).
* CONSTANTS[27] is k_gamma in component gamma (micromolar_per_millisecond).
* STATES[7] is v in component membrane (millivolt).
* CONSTANTS[28] is cm in component membrane (femtofarad).
* ALGEBRAIC[14] is I_Ca in component I_Ca (picoampere).
* ALGEBRAIC[12] is I_K in component I_K (picoampere).
* ALGEBRAIC[17] is I_K_Ca in component I_K_Ca (picoampere).
* ALGEBRAIC[32] is I_K_ATP in component I_K_ATP (picoampere).
* CONSTANTS[29] is gK_ in component I_K (picosiemens).
* CONSTANTS[30] is vK in component model_parameters (millivolt).
* STATES[8] is n in component n (dimensionless).
* ALGEBRAIC[5] is n_infinity in component n (dimensionless).
* ALGEBRAIC[2] is tau_n in component n (millisecond).
* CONSTANTS[31] is gCa_ in component I_Ca (picosiemens).
* CONSTANTS[32] is vCa in component model_parameters (millivolt).
* ALGEBRAIC[13] is m_infinity in component m (dimensionless).
* CONSTANTS[33] is gkCa_ in component I_K_Ca (picosiemens).
* ALGEBRAIC[15] is gkCa in component I_K_Ca (picosiemens).
* CONSTANTS[34] is KD in component I_K_Ca (micromolar).
* CONSTANTS[35] is gkATP_ in component I_K_ATP (picosiemens).
* ALGEBRAIC[30] is gkATP in component I_K_ATP (picosiemens).
* ALGEBRAIC[28] is o_infinity in component I_K_ATP (dimensionless).
* ALGEBRAIC[19] is MgADP in component I_K_ATP (micromolar).
* ALGEBRAIC[21] is ADP3_ in component I_K_ATP (micromolar).
* ALGEBRAIC[26] is ATP4_ in component I_K_ATP (micromolar).
* CONSTANTS[36] is fcyt in component Ca (dimensionless).
* ALGEBRAIC[16] is Jmem in component Jmem (micromolar_per_millisecond).
* ALGEBRAIC[22] is Jer in component Jer (micromolar_per_millisecond).
* STATES[9] is Caer in component Caer (micromolar).
* CONSTANTS[37] is fer in component Caer (dimensionless).
* CONSTANTS[38] is Vcyt_Ver in component Caer (dimensionless).
* CONSTANTS[39] is kPMCA in component Jmem (first_order_rate_constant).
* CONSTANTS[40] is alpha in component Jmem (micromolar_per_millisecond).
* ALGEBRAIC[20] is Jleak in component Jleak (micromolar_per_millisecond).
* ALGEBRAIC[18] is JSERCA in component JSERCA (micromolar_per_millisecond).
* CONSTANTS[41] is kSERCA in component JSERCA (first_order_rate_constant).
* CONSTANTS[42] is pleak in component Jleak (first_order_rate_constant).
* ALGEBRAIC[3] is I_infinity in component I (dimensionless).
* CONSTANTS[43] is tau_I in component I (millisecond).
* CONSTANTS[44] is I_slope in component I (per_micromolar).
* CONSTANTS[45] is Canull in component I (micromolar).
* RATES[0] is d/dt Gi in component Gi (millimolar).
* RATES[1] is d/dt Ge in component Ge (millimolar).
* RATES[3] is d/dt G6P in component G6P (millimolar).
* RATES[4] is d/dt FBP in component FBP (micromolar).
* RATES[5] is d/dt ADP in component ADP (micromolar).
* RATES[7] is d/dt v in component membrane (millivolt).
* RATES[8] is d/dt n in component n (dimensionless).
* RATES[6] is d/dt Ca in component Ca (micromolar).
* RATES[9] is d/dt Caer in component Caer (micromolar).
* RATES[2] is d/dt I in component I (dimensionless).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 0;
STATES[1] = 7;
CONSTANTS[0] = 450000;
CONSTANTS[1] = 15;
CONSTANTS[2] = 1;
CONSTANTS[3] = 1;
CONSTANTS[4] = 5;
STATES[2] = 0;
STATES[3] = 200;
CONSTANTS[5] = 0.005;
STATES[4] = 40;
CONSTANTS[6] = 7;
CONSTANTS[7] = 8;
CONSTANTS[8] = 7;
CONSTANTS[9] = 0.8;
CONSTANTS[10] = 4;
CONSTANTS[11] = 0.06;
CONSTANTS[12] = 2;
CONSTANTS[13] = 0.02;
CONSTANTS[14] = 0.2;
CONSTANTS[15] = 20;
CONSTANTS[16] = 20;
CONSTANTS[17] = 20;
CONSTANTS[18] = 30;
CONSTANTS[19] = 1;
CONSTANTS[20] = 50000;
CONSTANTS[21] = 1000;
CONSTANTS[22] = 3000;
STATES[5] = 780;
CONSTANTS[23] = 300000;
CONSTANTS[24] = 0.5;
CONSTANTS[25] = 0.35;
STATES[6] = 0.1;
CONSTANTS[26] = 2.2;
CONSTANTS[27] = 0.1;
STATES[7] = -60;
CONSTANTS[28] = 5300;
CONSTANTS[29] = 2700;
CONSTANTS[30] = -75;
STATES[8] = 0;
CONSTANTS[31] = 1000;
CONSTANTS[32] = 25;
CONSTANTS[33] = 400;
CONSTANTS[34] = 0.5;
CONSTANTS[35] = 40000;
CONSTANTS[36] = 0.01;
STATES[9] = 185;
CONSTANTS[37] = 0.01;
CONSTANTS[38] = 31;
CONSTANTS[39] = 0.18;
CONSTANTS[40] = 4.5e-6;
CONSTANTS[41] = 0.4;
CONSTANTS[42] = 0.0002;
CONSTANTS[43] = 10000;
CONSTANTS[44] = 210;
CONSTANTS[45] = 0.055;
CONSTANTS[46] = 1.00000;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = CONSTANTS[2]+(CONSTANTS[1] - CONSTANTS[2])/(1.00000+exp((STATES[2] - CONSTANTS[4])/CONSTANTS[3]));
RATES[1] = (ALGEBRAIC[1] - STATES[1])/CONSTANTS[0];
ALGEBRAIC[3] = (STATES[6]>=CONSTANTS[45] ? CONSTANTS[44]*(STATES[6] - CONSTANTS[45]) : STATES[6]<CONSTANTS[45] ? 0.00000 : 0.0/0.0);
RATES[2] = (ALGEBRAIC[3] - STATES[2])/CONSTANTS[43];
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[7]+16.0000)/5.60000));
ALGEBRAIC[2] = 1.00000/( 0.0350000*cosh((STATES[7]+16.0000)/22.4000));
RATES[8] = (ALGEBRAIC[5] - STATES[8])/ALGEBRAIC[2];
ALGEBRAIC[4] = ( CONSTANTS[7]*(STATES[1] - STATES[0])*CONSTANTS[6])/( (CONSTANTS[6]+STATES[1])*(CONSTANTS[6]+STATES[0]));
ALGEBRAIC[6] = ( CONSTANTS[9]*pow(STATES[0], CONSTANTS[10]))/(pow(CONSTANTS[8], CONSTANTS[10])+pow(STATES[0], CONSTANTS[10]));
RATES[0] = ALGEBRAIC[4] - ALGEBRAIC[6];
ALGEBRAIC[13] = 1.00000/(1.00000+exp(- (STATES[7]+20.0000)/12.0000));
ALGEBRAIC[14] = CONSTANTS[31]*ALGEBRAIC[13]*(STATES[7] - CONSTANTS[32]);
ALGEBRAIC[16] = - ( CONSTANTS[40]*ALGEBRAIC[14]+ CONSTANTS[39]*STATES[6]);
ALGEBRAIC[20] = CONSTANTS[42]*(STATES[9] - STATES[6]);
ALGEBRAIC[18] = CONSTANTS[41]*STATES[6];
ALGEBRAIC[22] = ALGEBRAIC[20] - ALGEBRAIC[18];
RATES[6] = CONSTANTS[36]*(ALGEBRAIC[16]+ALGEBRAIC[22]);
RATES[9] = - CONSTANTS[37]*CONSTANTS[38]*ALGEBRAIC[22];
rootfind_0(VOI, CONSTANTS, RATES, STATES, ALGEBRAIC, pret);
ALGEBRAIC[7] = 0.200000* pow((STATES[4]/1.00000), 1.0 / 2);
ALGEBRAIC[11] = ( CONSTANTS[26]*ALGEBRAIC[7])/(CONSTANTS[27]+ALGEBRAIC[7]);
RATES[5] = (1.00000/CONSTANTS[23])*(ALGEBRAIC[24] - STATES[5]*exp( (CONSTANTS[24]+ALGEBRAIC[11])*(1.00000 - STATES[6]/CONSTANTS[25])));
ALGEBRAIC[12] = CONSTANTS[29]*STATES[8]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[15] = ( CONSTANTS[33]*pow(STATES[6], 2.00000))/(pow(CONSTANTS[34], 2.00000)+pow(STATES[6], 2.00000));
ALGEBRAIC[17] = ALGEBRAIC[15]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[19] = 0.165000*STATES[5];
ALGEBRAIC[21] = 0.135000*STATES[5];
ALGEBRAIC[26] = 0.00500000*ALGEBRAIC[24];
ALGEBRAIC[28] = ( 0.0800000*(1.00000+( 2.00000*ALGEBRAIC[19])/17.0000)+ 0.890000*pow(ALGEBRAIC[19]/17.0000, 2.00000))/( pow(1.00000+ALGEBRAIC[19]/17.0000, 2.00000)*(1.00000+ALGEBRAIC[21]/26.0000+ALGEBRAIC[26]/1.00000));
ALGEBRAIC[30] = CONSTANTS[35]*ALGEBRAIC[28];
ALGEBRAIC[32] = ALGEBRAIC[30]*(STATES[7] - CONSTANTS[30]);
RATES[7] = - (ALGEBRAIC[12]+ALGEBRAIC[14]+ALGEBRAIC[17]+ALGEBRAIC[32])/CONSTANTS[28];
ALGEBRAIC[25] = pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000);
ALGEBRAIC[8] = pow(STATES[4]/CONSTANTS[19], 1.00000);
ALGEBRAIC[0] = 0.300000*STATES[3];
ALGEBRAIC[9] = pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[27] = pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[29] = pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000);
ALGEBRAIC[31] = (1.00000/pow(CONSTANTS[13], 1.00000*1.00000))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[33] = (1.00000/pow(CONSTANTS[15], 1.00000*1.00000))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[34] = (1.00000/pow(CONSTANTS[16], 1.00000*1.00000))*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[35] = (1.00000/pow(CONSTANTS[17], 1.00000*1.00000))*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[10] = (1.00000/pow(CONSTANTS[14], 1.00000*1.00000))*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[36] = (1.00000/( pow(CONSTANTS[13], 1.00000*1.00000)*pow(CONSTANTS[14], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[37] = (1.00000/( pow(CONSTANTS[14], 1.00000*1.00000)*pow(CONSTANTS[16], 1.00000*1.00000)*pow(CONSTANTS[17], 1.00000*1.00000)))*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[38] = (1.00000/( pow(CONSTANTS[13], 1.00000*1.00000)*pow(CONSTANTS[15], 1.00000*1.00000)*pow(CONSTANTS[17], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[39] = (1.00000/( pow(CONSTANTS[15], 1.00000*1.00000)*pow(CONSTANTS[16], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[40] = (1.00000/( pow(CONSTANTS[13], 1.00000*1.00000)*pow(CONSTANTS[14], 1.00000*1.00000)*pow(CONSTANTS[15], 1.00000*1.00000)*pow(CONSTANTS[16], 1.00000*1.00000)*pow(CONSTANTS[17], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[41] = ( CONSTANTS[12]*( (1.00000 - CONSTANTS[11])*ALGEBRAIC[36]+ CONSTANTS[11]*(ALGEBRAIC[9]+ALGEBRAIC[31]+ALGEBRAIC[35]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38]+ALGEBRAIC[10]+ALGEBRAIC[40])))/(CONSTANTS[46]+ALGEBRAIC[25]+ALGEBRAIC[8]+ALGEBRAIC[9]+ALGEBRAIC[27]+ALGEBRAIC[29]+ALGEBRAIC[31]+ALGEBRAIC[33]+ALGEBRAIC[34]+ALGEBRAIC[35]+ALGEBRAIC[10]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38]+ALGEBRAIC[39]+ALGEBRAIC[40]);
RATES[3] = CONSTANTS[5]*(ALGEBRAIC[6] - ALGEBRAIC[41]);
RATES[4] = CONSTANTS[5]*(ALGEBRAIC[41] - 0.500000*ALGEBRAIC[7]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = CONSTANTS[2]+(CONSTANTS[1] - CONSTANTS[2])/(1.00000+exp((STATES[2] - CONSTANTS[4])/CONSTANTS[3]));
ALGEBRAIC[3] = (STATES[6]>=CONSTANTS[45] ? CONSTANTS[44]*(STATES[6] - CONSTANTS[45]) : STATES[6]<CONSTANTS[45] ? 0.00000 : 0.0/0.0);
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[7]+16.0000)/5.60000));
ALGEBRAIC[2] = 1.00000/( 0.0350000*cosh((STATES[7]+16.0000)/22.4000));
ALGEBRAIC[4] = ( CONSTANTS[7]*(STATES[1] - STATES[0])*CONSTANTS[6])/( (CONSTANTS[6]+STATES[1])*(CONSTANTS[6]+STATES[0]));
ALGEBRAIC[6] = ( CONSTANTS[9]*pow(STATES[0], CONSTANTS[10]))/(pow(CONSTANTS[8], CONSTANTS[10])+pow(STATES[0], CONSTANTS[10]));
ALGEBRAIC[13] = 1.00000/(1.00000+exp(- (STATES[7]+20.0000)/12.0000));
ALGEBRAIC[14] = CONSTANTS[31]*ALGEBRAIC[13]*(STATES[7] - CONSTANTS[32]);
ALGEBRAIC[16] = - ( CONSTANTS[40]*ALGEBRAIC[14]+ CONSTANTS[39]*STATES[6]);
ALGEBRAIC[20] = CONSTANTS[42]*(STATES[9] - STATES[6]);
ALGEBRAIC[18] = CONSTANTS[41]*STATES[6];
ALGEBRAIC[22] = ALGEBRAIC[20] - ALGEBRAIC[18];
ALGEBRAIC[7] = 0.200000* pow((STATES[4]/1.00000), 1.0 / 2);
ALGEBRAIC[11] = ( CONSTANTS[26]*ALGEBRAIC[7])/(CONSTANTS[27]+ALGEBRAIC[7]);
ALGEBRAIC[12] = CONSTANTS[29]*STATES[8]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[15] = ( CONSTANTS[33]*pow(STATES[6], 2.00000))/(pow(CONSTANTS[34], 2.00000)+pow(STATES[6], 2.00000));
ALGEBRAIC[17] = ALGEBRAIC[15]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[19] = 0.165000*STATES[5];
ALGEBRAIC[21] = 0.135000*STATES[5];
ALGEBRAIC[26] = 0.00500000*ALGEBRAIC[24];
ALGEBRAIC[28] = ( 0.0800000*(1.00000+( 2.00000*ALGEBRAIC[19])/17.0000)+ 0.890000*pow(ALGEBRAIC[19]/17.0000, 2.00000))/( pow(1.00000+ALGEBRAIC[19]/17.0000, 2.00000)*(1.00000+ALGEBRAIC[21]/26.0000+ALGEBRAIC[26]/1.00000));
ALGEBRAIC[30] = CONSTANTS[35]*ALGEBRAIC[28];
ALGEBRAIC[32] = ALGEBRAIC[30]*(STATES[7] - CONSTANTS[30]);
ALGEBRAIC[25] = pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000);
ALGEBRAIC[8] = pow(STATES[4]/CONSTANTS[19], 1.00000);
ALGEBRAIC[0] = 0.300000*STATES[3];
ALGEBRAIC[9] = pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[27] = pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[29] = pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000);
ALGEBRAIC[31] = (1.00000/pow(CONSTANTS[13], 1.00000*1.00000))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[33] = (1.00000/pow(CONSTANTS[15], 1.00000*1.00000))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[34] = (1.00000/pow(CONSTANTS[16], 1.00000*1.00000))*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[35] = (1.00000/pow(CONSTANTS[17], 1.00000*1.00000))*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[10] = (1.00000/pow(CONSTANTS[14], 1.00000*1.00000))*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[36] = (1.00000/( pow(CONSTANTS[13], 1.00000*1.00000)*pow(CONSTANTS[14], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000);
ALGEBRAIC[37] = (1.00000/( pow(CONSTANTS[14], 1.00000*1.00000)*pow(CONSTANTS[16], 1.00000*1.00000)*pow(CONSTANTS[17], 1.00000*1.00000)))*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[38] = (1.00000/( pow(CONSTANTS[13], 1.00000*1.00000)*pow(CONSTANTS[15], 1.00000*1.00000)*pow(CONSTANTS[17], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[39] = (1.00000/( pow(CONSTANTS[15], 1.00000*1.00000)*pow(CONSTANTS[16], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[40] = (1.00000/( pow(CONSTANTS[13], 1.00000*1.00000)*pow(CONSTANTS[14], 1.00000*1.00000)*pow(CONSTANTS[15], 1.00000*1.00000)*pow(CONSTANTS[16], 1.00000*1.00000)*pow(CONSTANTS[17], 1.00000*1.00000)))*pow(ALGEBRAIC[23]/CONSTANTS[18], 1.00000)*pow(STATES[4]/CONSTANTS[19], 1.00000)*pow(pow(ALGEBRAIC[0], 2.00000)/CONSTANTS[20], 1.00000)*pow(pow(ALGEBRAIC[24], 2.00000)/CONSTANTS[21], 1.00000);
ALGEBRAIC[41] = ( CONSTANTS[12]*( (1.00000 - CONSTANTS[11])*ALGEBRAIC[36]+ CONSTANTS[11]*(ALGEBRAIC[9]+ALGEBRAIC[31]+ALGEBRAIC[35]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38]+ALGEBRAIC[10]+ALGEBRAIC[40])))/(CONSTANTS[46]+ALGEBRAIC[25]+ALGEBRAIC[8]+ALGEBRAIC[9]+ALGEBRAIC[27]+ALGEBRAIC[29]+ALGEBRAIC[31]+ALGEBRAIC[33]+ALGEBRAIC[34]+ALGEBRAIC[35]+ALGEBRAIC[10]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38]+ALGEBRAIC[39]+ALGEBRAIC[40]);
}
void objfunc_0(double *p, double *hx, int m, int n, void *adata)
{
struct rootfind_info* rfi = (struct rootfind_info*)adata;
#define VOI rfi->aVOI
#define CONSTANTS rfi->aCONSTANTS
#define RATES rfi->aRATES
#define STATES rfi->aSTATES
#define ALGEBRAIC rfi->aALGEBRAIC
#define pret rfi->aPRET
ALGEBRAIC[23] = p[0];
ALGEBRAIC[24] = p[1];
hx[0] = ALGEBRAIC[24] - (CONSTANTS[22] - (STATES[5]+ALGEBRAIC[23]));
hx[1] = ALGEBRAIC[23] - pow(STATES[5], 2.00000)/ALGEBRAIC[24];
#undef VOI
#undef CONSTANTS
#undef RATES
#undef STATES
#undef ALGEBRAIC
#undef pret
}
void rootfind_0(double VOI, double* CONSTANTS, double* RATES,
double* STATES, double* ALGEBRAIC, int* pret)
{
static double p[2] = {0.1,0.1};
double bp[2], work[LM_DIF_WORKSZ(2, 2)];
struct rootfind_info rfi;
rfi.aVOI = VOI;
rfi.aCONSTANTS = CONSTANTS;
rfi.aRATES = RATES;
rfi.aSTATES = STATES;
rfi.aALGEBRAIC = ALGEBRAIC;
rfi.aPRET = pret;
do_levmar(objfunc_0, p, bp, work, pret, 2, &rfi);
ALGEBRAIC[23] = p[0];
ALGEBRAIC[24] = p[1];
}