Generated Code

/*
   There are a total of 4 entries in the algebraic variable array.
   There are a total of 4 entries in each of the rate and state variable arrays.
   There are a total of 14 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is r in component r (nanomolar).
 * ALGEBRAIC[0] is Ir in component r (flux).
 * CONSTANTS[0] is k6 in component r (first_order_rate_constant).
 * STATES[1] is s in component s (nanomolar).
 * ALGEBRAIC[1] is Is in component model_parameters (flux).
 * CONSTANTS[1] is k7 in component model_parameters (first_order_rate_constant).
 * STATES[2] is f in component f (dimensionless).
 * CONSTANTS[2] is k1 in component f (second_order_rate_constant).
 * CONSTANTS[3] is k2 in component f (first_order_rate_constant).
 * CONSTANTS[4] is k3 in component f (first_order_rate_constant).
 * ALGEBRAIC[2] is phi_b_s in component f (dimensionless).
 * CONSTANTS[5] is sb in component f (dimensionless).
 * CONSTANTS[6] is delta_b in component f (dimensionless).
 * CONSTANTS[7] is c in component model_parameters (nanomolar).
 * STATES[3] is h in component h (nanomolar).
 * CONSTANTS[8] is k4 in component h (first_order_rate_constant).
 * CONSTANTS[9] is k5 in component h (first_order_rate_constant).
 * ALGEBRAIC[3] is phi_r_s in component h (dimensionless).
 * CONSTANTS[10] is sr in component h (dimensionless).
 * CONSTANTS[11] is delta_r in component h (dimensionless).
 * CONSTANTS[12] is k8 in component model_parameters (first_order_rate_constant).
 * CONSTANTS[13] is j1 in component model_parameters (dimensionless).
 * RATES[0] is d/dt r in component r (nanomolar).
 * RATES[1] is d/dt s in component s (nanomolar).
 * RATES[2] is d/dt f in component f (dimensionless).
 * RATES[3] is d/dt h in component h (nanomolar).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 0.0;
CONSTANTS[0] = 5.0;
STATES[1] = 0.0;
CONSTANTS[1] = 5.0;
STATES[2] = 0.3;
CONSTANTS[2] = 0.1;
CONSTANTS[3] = 0.002;
CONSTANTS[4] = 0.018;
CONSTANTS[5] = 0.029;
CONSTANTS[6] = 0.3;
CONSTANTS[7] = 0.01;
STATES[3] = 0.0;
CONSTANTS[8] = 9.0;
CONSTANTS[9] = 71.0;
CONSTANTS[10] = -0.56;
CONSTANTS[11] = 0.3;
CONSTANTS[12] = 0.07;
CONSTANTS[13] = 10;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (VOI>=0.00000&&VOI<=90.0000 ? 0.00000 : VOI>=91.0000&&VOI<=92.0000 ? 10.0000 : VOI>=93.0000&&VOI<=113.000 ? 0.00000 : VOI>=114.000&&VOI<=115.000 ? 10.0000 : VOI>=116.000&&VOI<=136.000 ? 0.00000 : VOI>=137.000&&VOI<=138.000 ? 10.0000 : VOI>=139.000&&VOI<=159.000 ? 0.00000 : VOI>=160.000&&VOI<=161.000 ? 10.0000 : VOI>=162.000&&VOI<=252.000 ? 0.00000 : VOI>=253.000&&VOI<=254.000 ? 10.0000 : VOI>=255.000&&VOI<=275.000 ? 0.00000 : VOI>=276.000&&VOI<=277.000 ? 10.0000 : VOI>=278.000&&VOI<=298.000 ? 0.00000 : VOI>=299.000&&VOI<=300.000 ? 10.0000 : VOI>=301.000&&VOI<=321.000 ? 0.00000 : VOI>=322.000&&VOI<=323.000 ? 10.0000 : 0.00000);
RATES[0] = ALGEBRAIC[0] -  CONSTANTS[0]*STATES[0];
ALGEBRAIC[1] = (VOI>0.00000&&VOI<=90.0000 ? 10.0000 : VOI>90.0000&&VOI<=180.000 ? 0.00000 : VOI>180.000&&VOI<=270.000 ? 10.0000 : VOI>270.000&&VOI<=360.000 ? 0.00000 : 0.00000);
RATES[1] = ALGEBRAIC[1] -  CONSTANTS[1]*STATES[1];
ALGEBRAIC[2] = 1.00000/(1.00000+exp(- (arbitrary_log( 1.00000*STATES[1], 10) - CONSTANTS[5])/CONSTANTS[6]));
RATES[2] = - ( CONSTANTS[2]*(STATES[0]+CONSTANTS[7])*STATES[2])+ (CONSTANTS[3]+ CONSTANTS[4]*ALGEBRAIC[2])*(1.00000 - STATES[2]);
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (arbitrary_log( 1.00000*STATES[1], 10) - CONSTANTS[10])/CONSTANTS[11]));
RATES[3] =  CONSTANTS[13]*( (CONSTANTS[8]+ CONSTANTS[9]*(1.00000 - ALGEBRAIC[3]))*( (STATES[0]+CONSTANTS[7])*STATES[2]) -  CONSTANTS[12]*STATES[3]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (VOI>=0.00000&&VOI<=90.0000 ? 0.00000 : VOI>=91.0000&&VOI<=92.0000 ? 10.0000 : VOI>=93.0000&&VOI<=113.000 ? 0.00000 : VOI>=114.000&&VOI<=115.000 ? 10.0000 : VOI>=116.000&&VOI<=136.000 ? 0.00000 : VOI>=137.000&&VOI<=138.000 ? 10.0000 : VOI>=139.000&&VOI<=159.000 ? 0.00000 : VOI>=160.000&&VOI<=161.000 ? 10.0000 : VOI>=162.000&&VOI<=252.000 ? 0.00000 : VOI>=253.000&&VOI<=254.000 ? 10.0000 : VOI>=255.000&&VOI<=275.000 ? 0.00000 : VOI>=276.000&&VOI<=277.000 ? 10.0000 : VOI>=278.000&&VOI<=298.000 ? 0.00000 : VOI>=299.000&&VOI<=300.000 ? 10.0000 : VOI>=301.000&&VOI<=321.000 ? 0.00000 : VOI>=322.000&&VOI<=323.000 ? 10.0000 : 0.00000);
ALGEBRAIC[1] = (VOI>0.00000&&VOI<=90.0000 ? 10.0000 : VOI>90.0000&&VOI<=180.000 ? 0.00000 : VOI>180.000&&VOI<=270.000 ? 10.0000 : VOI>270.000&&VOI<=360.000 ? 0.00000 : 0.00000);
ALGEBRAIC[2] = 1.00000/(1.00000+exp(- (arbitrary_log( 1.00000*STATES[1], 10) - CONSTANTS[5])/CONSTANTS[6]));
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (arbitrary_log( 1.00000*STATES[1], 10) - CONSTANTS[10])/CONSTANTS[11]));
}