Generated Code

/*
   There are a total of 5 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 20 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is G in component G (mg_l).
 * CONSTANTS[0] is Vg in component G (litre).
 * CONSTANTS[1] is Gin in component G (mg_l_min).
 * ALGEBRAIC[0] is f2_G in component G (mg_min).
 * ALGEBRAIC[1] is f3_G in component G (per_l).
 * ALGEBRAIC[2] is f4_I in component G (mg_min).
 * ALGEBRAIC[3] is f5_I in component G (mg_min).
 * CONSTANTS[2] is C2 in component G (mg_l).
 * CONSTANTS[3] is C3 in component G (mg_l).
 * CONSTANTS[4] is C4 in component G (mU_l).
 * CONSTANTS[5] is C5 in component G (mU_l).
 * CONSTANTS[6] is U0 in component G (mg_min).
 * CONSTANTS[7] is Um in component G (mg_min).
 * CONSTANTS[8] is Ub in component G (mg_min).
 * CONSTANTS[9] is beta in component G (dimensionless).
 * CONSTANTS[10] is Rg in component G (mg_min).
 * CONSTANTS[11] is alpha in component G (l_mU).
 * CONSTANTS[12] is ti in component G (minute).
 * CONSTANTS[13] is E in component G (l_min).
 * CONSTANTS[14] is Vp in component G (litre).
 * CONSTANTS[15] is Vi in component G (litre).
 * STATES[1] is I in component I (mU_l).
 * ALGEBRAIC[4] is f1_G in component I (mU_min).
 * CONSTANTS[16] is Rm in component I (mU_min).
 * CONSTANTS[17] is C1 in component I (mg_l).
 * CONSTANTS[18] is a1 in component I (mg_l).
 * CONSTANTS[19] is di in component I (first_order_rate_constant).
 * RATES[0] is d/dt G in component G (mg_l).
 * RATES[1] is d/dt I in component I (mU_l).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 1000.0;
CONSTANTS[0] = 10.0;
CONSTANTS[1] = 0.54;
CONSTANTS[2] = 144.0;
CONSTANTS[3] = 1000.0;
CONSTANTS[4] = 80.0;
CONSTANTS[5] = 26.0;
CONSTANTS[6] = 40.0;
CONSTANTS[7] = 940.0;
CONSTANTS[8] = 72.0;
CONSTANTS[9] = 1.77;
CONSTANTS[10] = 180.0;
CONSTANTS[11] = 0.29;
CONSTANTS[12] = 100;
CONSTANTS[13] = 0.2;
CONSTANTS[14] = 3.0;
CONSTANTS[15] = 11.0;
STATES[1] = 9000.0;
CONSTANTS[16] = 210.0;
CONSTANTS[17] = 2000.0;
CONSTANTS[18] = 300.0;
CONSTANTS[19] = 0.06;
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] - CONSTANTS[1]+ 1.00000*ALGEBRAIC[3]+- ( 1.00000*ALGEBRAIC[0]+ ALGEBRAIC[1]*ALGEBRAIC[2]);
resid[1] = RATES[1] -  1.00000*ALGEBRAIC[4] -  CONSTANTS[19]*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] =  CONSTANTS[8]*(1.00000 - exp( 1.00000*(- STATES[0]/( CONSTANTS[2]*CONSTANTS[0]))));
ALGEBRAIC[1] = STATES[0]/( CONSTANTS[3]*CONSTANTS[0]);
ALGEBRAIC[2] = CONSTANTS[6]+(CONSTANTS[7] - CONSTANTS[6])/(1.00000+exp( - CONSTANTS[9]*log( 1.00000*(STATES[1]/( CONSTANTS[4]*(1.00000/CONSTANTS[15]+1.00000/( CONSTANTS[13]*CONSTANTS[12])))))));
ALGEBRAIC[3] = CONSTANTS[10]/(1.00000+exp( CONSTANTS[11]*(STATES[1]/( 1.00000*CONSTANTS[14] -  1.00000*CONSTANTS[5]))));
ALGEBRAIC[4] = CONSTANTS[16]/(1.00000+exp(((CONSTANTS[17] - STATES[0])/( 1.00000*CONSTANTS[0]))/CONSTANTS[18]));
}
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)
{
}