/*
   There are a total of 0 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 5 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * CONSTANTS[0] is VVE in component stress_relaxation (litre).
 * STATES[0] is VV7 in component short_term_stress_relaxation (litre).
 * CONSTANTS[1] is SR in component parameter_values (dimensionless).
 * CONSTANTS[2] is SRK in component parameter_values (minute).
 * STATES[1] is VV6 in component long_term_stress_relaxation (litre).
 * CONSTANTS[3] is SR2 in component parameter_values (dimensionless).
 * CONSTANTS[4] is SRK2 in component parameter_values (minute).
 * RATES[0] is d/dt VV7 in component short_term_stress_relaxation (litre).
 * RATES[1] is d/dt VV6 in component long_term_stress_relaxation (litre).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 0.743224;
STATES[0] = 0.00366525;
CONSTANTS[1] = 1;
CONSTANTS[2] = 5;
STATES[1] = 0.0101913;
CONSTANTS[3] = 1;
CONSTANTS[4] = 10000;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[0] = ( (CONSTANTS[0] - 0.740000)*CONSTANTS[1] - STATES[0])/CONSTANTS[2];
RATES[1] = ( (CONSTANTS[0] - 0.740000)*CONSTANTS[3] - STATES[1])/CONSTANTS[4];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}