Generated Code

/*
   There are a total of 6 entries in the algebraic variable array.
   There are a total of 2 entries in each of the rate and state variable arrays.
   There are a total of 15 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is G in component G (mg_per_litre).
 * ALGEBRAIC[0] is Gin in component G (mg_per_litre_minute).
 * ALGEBRAIC[2] is f2_G in component f2_G (mg_per_minute).
 * ALGEBRAIC[3] is f3_G in component f3_G (per_litre).
 * ALGEBRAIC[4] is f4_I in component f4_I (mg_per_minute).
 * ALGEBRAIC[5] is f5_I in component f5_I (mg_per_minute).
 * STATES[1] is I in component I (mU_per_litre).
 * ALGEBRAIC[1] is Iin in component I (mU_per_litre_minute).
 * CONSTANTS[0] is di in component model_parameters (first_order_rate_constant).
 * CONSTANTS[1] is C2 in component model_parameters (mg_per_litre).
 * CONSTANTS[2] is Ub in component model_parameters (mg_per_minute).
 * CONSTANTS[3] is Vg in component model_parameters (litre).
 * CONSTANTS[4] is C3 in component model_parameters (mg_per_litre).
 * CONSTANTS[5] is C4 in component model_parameters (mU_per_litre).
 * CONSTANTS[6] is Vi in component model_parameters (litre).
 * CONSTANTS[7] is U0 in component model_parameters (mg_per_minute).
 * CONSTANTS[8] is Um in component model_parameters (mg_per_minute).
 * CONSTANTS[9] is beta in component model_parameters (dimensionless).
 * CONSTANTS[10] is ti in component model_parameters (minute).
 * CONSTANTS[11] is Vp in component model_parameters (litre).
 * CONSTANTS[12] is C5 in component model_parameters (mU_per_litre).
 * CONSTANTS[13] is Rg in component model_parameters (mg_per_minute).
 * CONSTANTS[14] is alpha in component model_parameters (litre_per_mU).
 * RATES[0] is d/dt G in component G (mg_per_litre).
 * RATES[1] is d/dt I in component I (mU_per_litre).
 * There are a total of 40 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 150.0;
STATES[1] = 5.0;
CONSTANTS[0] = 0.0107;
CONSTANTS[1] = 144;
CONSTANTS[2] = 72.0;
CONSTANTS[3] = 10.0;
CONSTANTS[4] = 1000.0;
CONSTANTS[5] = 80.0;
CONSTANTS[6] = 11.0;
CONSTANTS[7] = 40.0;
CONSTANTS[8] = 940.0;
CONSTANTS[9] = 1.77;
CONSTANTS[10] = 100.0;
CONSTANTS[11] = 3.0;
CONSTANTS[12] = 26.0;
CONSTANTS[13] = 180.0;
CONSTANTS[14] = 0.29;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - (ALGEBRAIC[0]+ 1.00000*ALGEBRAIC[5]) - ( 1.00000*ALGEBRAIC[2]+ ALGEBRAIC[3]*ALGEBRAIC[4]);
resid[1] = RATES[1] - ALGEBRAIC[1] -  CONSTANTS[0]*STATES[1];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (CONDVAR[0]<0.00000&&CONDVAR[1]>=0.00000 ? 0.0500000+ (5.00000/15.0000)*VOI : CONDVAR[2]<0.00000&&CONDVAR[3]>=0.00000 ? 0.0500000+ 5.00000*((45.0000 - VOI)/(45.0000 - 15.0000)) : CONDVAR[4]<=0.00000&&CONDVAR[5]>=0.00000 ? 0.0500000 : CONDVAR[6]>0.00000&&CONDVAR[7]<0.00000 ? 0.0500000+ (5.00000/255.000)*VOI : CONDVAR[8]<0.00000&&CONDVAR[9]>=0.00000 ? 0.0500000+ 5.00000*((285.000 - VOI)/(285.000 - 255.000)) : CONDVAR[10]<=0.00000&&CONDVAR[11]>=0.00000 ? 0.0500000 : CONDVAR[12]<0.00000&&CONDVAR[13]>0.00000 ? 0.0500000+ (5.00000/495.000)*VOI : CONDVAR[14]<0.00000&&CONDVAR[15]>=0.00000 ? 0.0500000+ 5.00000*((525.000 - VOI)/(525.000 - 495.000)) : CONDVAR[16]<=0.00000&&CONDVAR[17]>=0.00000 ? 0.0500000 : CONDVAR[18]>0.00000&&CONDVAR[19]<0.00000 ? 0.0500000+ (5.00000/735.000)*VOI : CONDVAR[20]<0.00000&&CONDVAR[21]>=0.00000 ? 0.0500000+ 5.00000*((765.000 - VOI)/(765.000 - 735.000)) : CONDVAR[22]<=0.00000&&CONDVAR[23]>=0.00000 ? 0.0500000 : 0.0/0.0);
ALGEBRAIC[1] = (CONDVAR[24]<=0.00000&&CONDVAR[25]>=0.00000 ? 0.250000 : CONDVAR[26]<0.00000&&CONDVAR[27]>=0.00000 ? 0.250000+ 1.00000*(1.00000+(VOI - 120.000)/90.0000) : CONDVAR[28]<0.00000&&CONDVAR[29]>=0.00000 ? 0.250000+ 1.00000*(1.00000 -  0.500000*((VOI - 120.000)/120.000)) : CONDVAR[30]<=0.00000&&CONDVAR[31]>=0.00000 ? 0.250000+ 0.500000*(1.00000 - (VOI - 240.000)/240.000) : CONDVAR[32]<0.00000&&CONDVAR[33]>0.00000 ? 0.250000 : CONDVAR[34]<0.00000&&CONDVAR[35]>=0.00000 ? 0.250000+ 1.00000*(1.00000+(VOI - 620.000)/90.0000) : CONDVAR[36]<0.00000&&CONDVAR[37]>=0.00000 ? 0.250000+ 1.00000*(1.00000 -  0.500000*((VOI - 620.000)/120.000)) : CONDVAR[38]<=0.00000&&CONDVAR[39]>=0.00000 ? 0.250000+ 0.500000*(1.00000 - (VOI - 740.000)/240.000) : 0.0/0.0);
ALGEBRAIC[2] =  CONSTANTS[2]*(1.00000 - exp(- STATES[0]/( CONSTANTS[1]*CONSTANTS[3]*1.00000)));
ALGEBRAIC[3] = STATES[0]/( CONSTANTS[4]*CONSTANTS[3]);
ALGEBRAIC[4] = CONSTANTS[7]+(CONSTANTS[8] - CONSTANTS[7])/(1.00000+exp( - CONSTANTS[9]*log(STATES[1]/( CONSTANTS[5]*(1.00000/CONSTANTS[6]+1.00000/( 0.200000*CONSTANTS[10]))))));
ALGEBRAIC[5] = CONSTANTS[13]/(1.00000+exp( CONSTANTS[14]*(STATES[1]/( 1.00000*( 1.00000*CONSTANTS[11] - CONSTANTS[12])))));
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - 15.0000;
CONDVAR[1] = VOI - 0.00000;
CONDVAR[2] = VOI - 45.0000;
CONDVAR[3] = VOI - 15.0000;
CONDVAR[4] = VOI - 240.000;
CONDVAR[5] = VOI - 45.0000;
CONDVAR[6] = VOI - 240.000;
CONDVAR[7] = VOI - 255.000;
CONDVAR[8] = VOI - 285.000;
CONDVAR[9] = VOI - 255.000;
CONDVAR[10] = VOI - 480.000;
CONDVAR[11] = VOI - 285.000;
CONDVAR[12] = VOI - 495.000;
CONDVAR[13] = VOI - 480.000;
CONDVAR[14] = VOI - 525.000;
CONDVAR[15] = VOI - 495.000;
CONDVAR[16] = VOI - 720.000;
CONDVAR[17] = VOI - 525.000;
CONDVAR[18] = VOI - 720.000;
CONDVAR[19] = VOI - 735.000;
CONDVAR[20] = VOI - 765.000;
CONDVAR[21] = VOI - 735.000;
CONDVAR[22] = VOI - 960.000;
CONDVAR[23] = VOI - 765.000;
CONDVAR[24] = VOI - 30.0000;
CONDVAR[25] = VOI - 0.00000;
CONDVAR[26] = VOI - 120.000;
CONDVAR[27] = VOI - 30.0000;
CONDVAR[28] = VOI - 240.000;
CONDVAR[29] = VOI - 120.000;
CONDVAR[30] = VOI - 480.000;
CONDVAR[31] = VOI - 240.000;
CONDVAR[32] = VOI - 530.000;
CONDVAR[33] = VOI - 480.000;
CONDVAR[34] = VOI - 620.000;
CONDVAR[35] = VOI - 530.000;
CONDVAR[36] = VOI - 740.000;
CONDVAR[37] = VOI - 620.000;
CONDVAR[38] = VOI - 980.000;
CONDVAR[39] = VOI - 740.000;
}