/*
   There are a total of 15 entries in the algebraic variable array.
   There are a total of 7 entries in each of the rate and state variable arrays.
   There are a total of 40 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * CONSTANTS[0] is cap in component parameters (nanofarad).
 * CONSTANTS[1] is cc1lin in component parameters (per_second).
 * CONSTANTS[2] is cc_2 in component parameters (per_second).
 * CONSTANTS[3] is ck1lin in component parameters (per_second).
 * CONSTANTS[4] is ck_2 in component parameters (per_second).
 * CONSTANTS[5] is clmax in component parameters (nanosiemens).
 * CONSTANTS[6] is cnmax in component parameters (nanosiemens).
 * CONSTANTS[7] is cx1lin in component parameters (per_second).
 * CONSTANTS[8] is cx2 in component parameters (per_second).
 * CONSTANTS[9] is ef in component parameters (per_second).
 * CONSTANTS[10] is gl in component parameters (nanosiemens).
 * CONSTANTS[11] is hmc_1 in component parameters (uM).
 * CONSTANTS[12] is hmc_2 in component parameters (uM).
 * CONSTANTS[13] is inf in component parameters (uM_per_picocoulomb).
 * CONSTANTS[14] is inhmax in component parameters (dimensionless).
 * CONSTANTS[15] is k_1 in component parameters (per_uM_per_second).
 * CONSTANTS[16] is k_2 in component parameters (per_second).
 * CONSTANTS[17] is kI in component parameters (uM).
 * CONSTANTS[18] is kinh in component parameters (uM).
 * CONSTANTS[19] is kinhcng in component parameters (uM).
 * CONSTANTS[20] is n_1 in component parameters (dimensionless).
 * CONSTANTS[21] is n_2 in component parameters (dimensionless).
 * CONSTANTS[22] is nI in component parameters (dimensionless).
 * CONSTANTS[23] is ninh in component parameters (dimensionless).
 * CONSTANTS[24] is ninhcng in component parameters (dimensionless).
 * CONSTANTS[25] is pd in component parameters (per_second).
 * CONSTANTS[26] is r_1 in component parameters (per_second).
 * CONSTANTS[27] is r_2 in component parameters (per_second).
 * CONSTANTS[28] is smax in component parameters (uM_per_second).
 * CONSTANTS[29] is V_Cl in component parameters (millivolt).
 * CONSTANTS[30] is V_cng in component parameters (millivolt).
 * CONSTANTS[31] is V_l in component parameters (millivolt).
 * CONSTANTS[39] is F_vol in component parameters (picocoulomb_per_uM).
 * CONSTANTS[32] is F in component parameters (coulombs_per_mole).
 * CONSTANTS[33] is C_vol in component parameters (liter).
 * ALGEBRAIC[2] is O_stim in component O_stim (uM).
 * CONSTANTS[34] is od in component O_stim (uM).
 * CONSTANTS[35] is t_0 in component O_stim (second).
 * CONSTANTS[36] is t_1 in component O_stim (second).
 * ALGEBRAIC[0] is H_0 in component O_stim (dimensionless).
 * ALGEBRAIC[1] is H_1 in component O_stim (dimensionless).
 * STATES[0] is bLR in component bLR (dimensionless).
 * CONSTANTS[37] is R_tot in component bLR (dimensionless).
 * STATES[1] is aG in component aG (dimensionless).
 * CONSTANTS[38] is G_tot in component aG (dimensionless).
 * ALGEBRAIC[3] is k_G in component k_G (per_second).
 * ALGEBRAIC[4] is r_G in component r_G (per_second).
 * STATES[2] is cAMP in component cAMP (uM).
 * ALGEBRAIC[5] is synth in component synth (uM_per_second).
 * ALGEBRAIC[6] is degrad in component degrad (uM_per_second).
 * STATES[3] is aCaMK in component aCaMK (uM).
 * STATES[4] is Ca in component Ca (uM).
 * ALGEBRAIC[13] is I_CNG in component I_CNG (nanoampere).
 * ALGEBRAIC[11] is J_NCX in component J_NCX (uM_per_second).
 * ALGEBRAIC[7] is cc_1 in component cc_1 (uM_per_second).
 * STATES[5] is CaCaM in component CaCaM (uM).
 * ALGEBRAIC[8] is ck_1 in component ck_1 (uM_per_second).
 * STATES[6] is V in component V (millivolt).
 * ALGEBRAIC[10] is I_ClCa in component I_ClCa (nanoampere).
 * ALGEBRAIC[14] is I_NCX in component I_NCX (nanoampere).
 * ALGEBRAIC[12] is I_other in component I_other (nanoampere).
 * ALGEBRAIC[9] is inhcng in component inhcng (dimensionless).
 * RATES[0] is d/dt bLR in component bLR (dimensionless).
 * RATES[1] is d/dt aG in component aG (dimensionless).
 * RATES[2] is d/dt cAMP in component cAMP (uM).
 * RATES[4] is d/dt Ca in component Ca (uM).
 * RATES[5] is d/dt CaCaM in component CaCaM (uM).
 * RATES[3] is d/dt aCaMK in component aCaMK (uM).
 * RATES[6] is d/dt V in component V (millivolt).
 * There are a total of 2 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 0.004;
CONSTANTS[1] = 0.88;
CONSTANTS[2] = 26;
CONSTANTS[3] = 13;
CONSTANTS[4] = 0.9;
CONSTANTS[5] = 1;
CONSTANTS[6] = 1;
CONSTANTS[7] = 1;
CONSTANTS[8] = 13;
CONSTANTS[9] = 2;
CONSTANTS[10] = 6;
CONSTANTS[11] = 2;
CONSTANTS[12] = 3;
CONSTANTS[13] = 1.9;
CONSTANTS[14] = 5;
CONSTANTS[15] = 0.06;
CONSTANTS[16] = 20;
CONSTANTS[17] = 0.7;
CONSTANTS[18] = 2;
CONSTANTS[19] = 1;
CONSTANTS[20] = 2;
CONSTANTS[21] = 2;
CONSTANTS[22] = 2;
CONSTANTS[23] = 1.5;
CONSTANTS[24] = 1.3;
CONSTANTS[25] = 20;
CONSTANTS[26] = 10;
CONSTANTS[27] = 5;
CONSTANTS[28] = 71;
CONSTANTS[29] = -50;
CONSTANTS[30] = 0;
CONSTANTS[31] = -70;
CONSTANTS[32] = 9.649e4;
CONSTANTS[33] = 1e-13;
CONSTANTS[34] = 20;
CONSTANTS[35] = 0.5;
CONSTANTS[36] = 1.5;
STATES[0] = 0;
CONSTANTS[37] = 1;
STATES[1] = 0;
CONSTANTS[38] = 1;
STATES[2] = 1.35648992164649e-88;
STATES[3] = 6.60756525051462e-8;
STATES[4] = 5.09073088043779e-12;
STATES[5] = 1.86113118246926e-13;
STATES[6] = -70;
CONSTANTS[39] =  (1.00000e+12/1.00000)*(1.00000/1000.00)*CONSTANTS[32]*CONSTANTS[33];
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[3] = 0.1001;
RATES[6] = 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[15]*ALGEBRAIC[2]*(CONSTANTS[37] - STATES[0]) -  CONSTANTS[26]*STATES[0];
resid[1] = RATES[1] -  ALGEBRAIC[3]*(CONSTANTS[38] - STATES[1]) - ALGEBRAIC[4];
resid[2] = RATES[2] - ALGEBRAIC[5] - ALGEBRAIC[6];
resid[3] = RATES[4] - ( (1000.00/1.00000)*CONSTANTS[13]*ALGEBRAIC[13] - ALGEBRAIC[11]) - (ALGEBRAIC[7] -  CONSTANTS[2]*STATES[5]);
resid[4] = RATES[5] - ALGEBRAIC[7] -  CONSTANTS[2]*STATES[5];
resid[5] = RATES[3] - ALGEBRAIC[8] -  CONSTANTS[4]*STATES[3];
resid[6] = RATES[6] -  (1000.00/1.00000)*(1.00000/CONSTANTS[0])*(ALGEBRAIC[13]+ALGEBRAIC[10]+ALGEBRAIC[14]+ALGEBRAIC[12]);
}
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 ? 0.00000 : 1.00000);
ALGEBRAIC[1] = (CONDVAR[1]<0.00000 ? 0.00000 : 1.00000);
ALGEBRAIC[2] =  CONSTANTS[34]*(ALGEBRAIC[0] - ALGEBRAIC[1]);
ALGEBRAIC[3] =  CONSTANTS[16]*STATES[0];
ALGEBRAIC[4] =  CONSTANTS[27]*STATES[1];
ALGEBRAIC[5] = ( STATES[1]*CONSTANTS[28])/(1.00000+pow(STATES[3]/CONSTANTS[18], CONSTANTS[23]));
ALGEBRAIC[6] =  CONSTANTS[25]*STATES[2];
ALGEBRAIC[7] =  CONSTANTS[1]*STATES[4];
ALGEBRAIC[8] =  CONSTANTS[3]*STATES[5];
ALGEBRAIC[10] =  (( CONSTANTS[5]*pow(STATES[4], CONSTANTS[21]))/(pow(STATES[4], CONSTANTS[21])+pow(CONSTANTS[12], CONSTANTS[21])))*(1.00000/1000.00)*(CONSTANTS[29] - STATES[6]);
ALGEBRAIC[11] =  CONSTANTS[9]*STATES[4];
ALGEBRAIC[12] =  CONSTANTS[10]*(1.00000/1000.00)*(CONSTANTS[31] - STATES[6]);
ALGEBRAIC[9] = 1.00000+( (CONSTANTS[14] - 1.00000)*pow(STATES[5], CONSTANTS[24]))/(pow(STATES[5], CONSTANTS[24])+pow(CONSTANTS[19], CONSTANTS[24]));
ALGEBRAIC[13] =  (( CONSTANTS[6]*pow(STATES[2], CONSTANTS[20]))/(pow(STATES[2], CONSTANTS[20])+pow( ALGEBRAIC[9]*CONSTANTS[11], CONSTANTS[20])))*(1.00000/1000.00)*(CONSTANTS[30] - STATES[6]);
ALGEBRAIC[14] =  (1.00000/1000.00)*CONSTANTS[39]*ALGEBRAIC[11];
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - CONSTANTS[35];
CONDVAR[1] = VOI - CONSTANTS[36];
}