/*
   There are a total of 4 entries in the algebraic variable array.
   There are a total of 3 entries in each of the rate and state variable arrays.
   There are a total of 18 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * CONSTANTS[0] is L_1 in component contraction (mm).
 * CONSTANTS[1] is L_2 in component contraction (mm).
 * CONSTANTS[2] is L_3 in component contraction (mm).
 * CONSTANTS[3] is L_4 in component contraction (mm).
 * CONSTANTS[4] is f_c in component contraction (newton).
 * CONSTANTS[5] is v_max in component contraction (mm_per_second).
 * CONSTANTS[6] is curv in component contraction (dimensionless).
 * CONSTANTS[7] is k_1 in component contraction (newton).
 * CONSTANTS[8] is k_2 in component contraction (per_mm).
 * CONSTANTS[9] is F_1 in component contraction (newton).
 * CONSTANTS[10] is d_LSEC1 in component contraction (mm).
 * CONSTANTS[11] is k_sh in component contraction (dimensionless).
 * CONSTANTS[12] is L_m in component contraction (mm).
 * CONSTANTS[13] is F_im in component contraction (newton).
 * CONSTANTS[14] is tau in component contraction (second).
 * STATES[1] is v_cc in component contraction (mm_per_second).
 * STATES[2] is f_v in component contraction (dimensionless).
 * ALGEBRAIC[0] is f_L in component contraction (newton).
 * ALGEBRAIC[3] is f_sec in component contraction (newton).
 * ALGEBRAIC[2] is delta_L_sec in component contraction (mm).
 * CONSTANTS[15] is delta_L_sec1 in component contraction (mm).
 * CONSTANTS[16] is k in component contraction (newton_per_mm).
 * STATES[0] is L_cc in component contraction (mm).
 * ALGEBRAIC[1] is L_mtc in component contraction (mm).
 * CONSTANTS[17] is A in component contraction (dimensionless).
 * RATES[0] is d/dt L_cc in component contraction (mm).
 * There are a total of 14 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = -24;
CONSTANTS[1] = -14;
CONSTANTS[2] = 4;
CONSTANTS[3] = 14;
CONSTANTS[4] = 0.51;
CONSTANTS[5] = -152;
CONSTANTS[6] = 9.5;
CONSTANTS[7] = 0.0064;
CONSTANTS[8] = 0.414;
CONSTANTS[9] = 8.9;
CONSTANTS[10] = 2.9;
CONSTANTS[11] = 2.5;
CONSTANTS[12] = 2.6;
CONSTANTS[13] = 20.6;
CONSTANTS[14] = 0.065;
CONSTANTS[15] = 4.1;
CONSTANTS[16] = 8.2;
STATES[0] = 0.2;
CONSTANTS[17] = 1;
STATES[1] = 0.1001;
STATES[2] = 0.1001;
RATES[0] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = STATES[1] -  (1.00000/STATES[2])*(ALGEBRAIC[3]/( CONSTANTS[17]*ALGEBRAIC[0]*CONSTANTS[13]));
resid[1] = RATES[0] - STATES[1];
resid[2] = STATES[2] - (CONSTANTS[5] - STATES[1])/(CONSTANTS[5]+ STATES[1]*CONSTANTS[6]);
}
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[3]>=0.00000&&CONDVAR[4]<=0.00000 ?  (CONSTANTS[4]/(CONSTANTS[1] - CONSTANTS[0]))*(STATES[0] - CONSTANTS[0]) : CONDVAR[5]>0.00000&&CONDVAR[6]<=0.00000 ?  ((1.00000 - CONSTANTS[4])/- CONSTANTS[1])*(STATES[0] - CONSTANTS[1]) : CONDVAR[7]>0.00000&&CONDVAR[8]<=0.00000 ? 1.00000 : CONDVAR[9]>0.00000&&CONDVAR[10]<=0.00000 ?  (- 1.00000/(CONSTANTS[3] - CONSTANTS[2]))*(STATES[0] - CONSTANTS[2]) : 0.0/0.0);
ALGEBRAIC[1] = (CONDVAR[11]<=0.00000 ? 0.290000 : CONDVAR[12]>0.00000&&CONDVAR[13]<0.00000 ? 0.220000 : 0.190000);
ALGEBRAIC[2] = (ALGEBRAIC[1] - STATES[0]) - CONSTANTS[12];
ALGEBRAIC[3] = (CONDVAR[0]>0.00000&&CONDVAR[1]<0.00000 ?  (CONSTANTS[9]/(exp(CONSTANTS[11]) - 1.00000))*(exp(( CONSTANTS[11]*ALGEBRAIC[2])/CONSTANTS[15]) - 1.00000) : CONDVAR[2]<=0.00000 ? CONSTANTS[9]+ CONSTANTS[16]*(ALGEBRAIC[2] - CONSTANTS[15]) : 0.0/0.0);
}
void
getStateInformation(double* SI)
{
SI[1] = 0.0;
SI[2] = 0.0;
SI[0] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = ALGEBRAIC[2] - 0.00000;
CONDVAR[1] = ALGEBRAIC[2] - CONSTANTS[15];
CONDVAR[2] = ALGEBRAIC[2] - CONSTANTS[15];
CONDVAR[3] = STATES[0] - CONSTANTS[0];
CONDVAR[4] = STATES[0] - CONSTANTS[1];
CONDVAR[5] = STATES[0] - CONSTANTS[1];
CONDVAR[6] = STATES[0] - 0.00000;
CONDVAR[7] = STATES[0] - 0.00000;
CONDVAR[8] = STATES[0] - CONSTANTS[2];
CONDVAR[9] = STATES[0] - CONSTANTS[2];
CONDVAR[10] = STATES[0] - CONSTANTS[3];
CONDVAR[11] = VOI - 1.00000;
CONDVAR[12] = VOI - 1.00000;
CONDVAR[13] = VOI - 5.00000;
}