/*
   There are a total of 13 entries in the algebraic variable array.
   There are a total of 5 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 (second).
 * STATES[0] is V_s in component soma_compartment (mV).
 * CONSTANTS[0] is V_Na in component soma_compartment (mV).
 * CONSTANTS[1] is V_K in component soma_compartment (mV).
 * ALGEBRAIC[10] is I_Na in component soma_compartment (uA_per_cm2).
 * ALGEBRAIC[0] is I_K_DR in component soma_compartment (uA_per_cm2).
 * CONSTANTS[2] is g_K_DR in component soma_compartment (mS_per_cm2).
 * CONSTANTS[3] is g_Na in component soma_compartment (mS_per_cm2).
 * CONSTANTS[4] is g_c in component model_parameters (mS_per_cm2).
 * STATES[1] is V_D in component dendritic_compartment (mV).
 * CONSTANTS[5] is C_m in component model_parameters (uF_per_cm2).
 * CONSTANTS[6] is p in component model_parameters (dimensionless).
 * STATES[2] is n in component gating_variables (dimensionless).
 * STATES[3] is h in component gating_variables (dimensionless).
 * ALGEBRAIC[5] is m_infinity in component gating_variables (dimensionless).
 * CONSTANTS[7] is V_L in component dendritic_compartment (mV).
 * CONSTANTS[8] is V_NMDA in component dendritic_compartment (mV).
 * ALGEBRAIC[1] is I_L in component dendritic_compartment (uA_per_cm2).
 * ALGEBRAIC[12] is I_D in component dendritic_compartment (uA_per_cm2).
 * ALGEBRAIC[11] is I_pump in component dendritic_compartment (uA_per_cm2).
 * ALGEBRAIC[3] is I_NMDA in component dendritic_compartment (uA_per_cm2).
 * ALGEBRAIC[4] is I_Na_NMDA in component dendritic_compartment (uA_per_cm2).
 * CONSTANTS[9] is R_pump in component dendritic_compartment (uA_per_cm2).
 * CONSTANTS[10] is alpha in component dendritic_compartment (mMcm2_per_uAs).
 * CONSTANTS[11] is g_NMDA in component dendritic_compartment (mS_per_cm2).
 * CONSTANTS[12] is g_Na_NMDA in component dendritic_compartment (mS_per_cm2).
 * CONSTANTS[13] is g_L in component dendritic_compartment (mS_per_cm2).
 * STATES[4] is Na in component dendritic_compartment (mM).
 * CONSTANTS[14] is Na_eq in component dendritic_compartment (mM).
 * CONSTANTS[15] is K_p in component dendritic_compartment (mM).
 * CONSTANTS[16] is Mg_o in component dendritic_compartment (mM).
 * CONSTANTS[17] is K_Mg in component dendritic_compartment (mM).
 * CONSTANTS[18] is q in component dendritic_compartment (mV).
 * ALGEBRAIC[2] is phi_Na in component dendritic_compartment (dimensionless).
 * CONSTANTS[19] is phi_Na_eq in component dendritic_compartment (dimensionless).
 * ALGEBRAIC[6] is n_infinity in component gating_variables (dimensionless).
 * ALGEBRAIC[7] is h_infinity in component gating_variables (dimensionless).
 * ALGEBRAIC[8] is tau_h in component gating_variables (second).
 * ALGEBRAIC[9] is tau_n in component gating_variables (second).
 * RATES[0] is d/dt V_s in component soma_compartment (mV).
 * RATES[1] is d/dt V_D in component dendritic_compartment (mV).
 * RATES[4] is d/dt Na in component dendritic_compartment (mM).
 * RATES[3] is d/dt h in component gating_variables (dimensionless).
 * RATES[2] is d/dt n in component gating_variables (dimensionless).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -64;
CONSTANTS[0] = 55;
CONSTANTS[1] = -85;
CONSTANTS[2] = 3.2;
CONSTANTS[3] = 3.2;
CONSTANTS[4] = 0.0;
STATES[1] = -25.0;
CONSTANTS[5] = 1;
CONSTANTS[6] = 0.5;
STATES[2] = 0.002;
STATES[3] = 1;
CONSTANTS[7] = -50;
CONSTANTS[8] = 0;
CONSTANTS[9] = 18.0;
CONSTANTS[10] = 0.5;
CONSTANTS[11] = 1.25;
CONSTANTS[12] = 1.0;
CONSTANTS[13] = 0.18;
STATES[4] = 5.09;
CONSTANTS[14] = 8;
CONSTANTS[15] = 15;
CONSTANTS[16] = 1.4;
CONSTANTS[17] = 10.0;
CONSTANTS[18] = 12.5;
CONSTANTS[19] = pow(CONSTANTS[14], 3.00000)/(pow(CONSTANTS[14], 3.00000)+pow(CONSTANTS[15], 3.00000));
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[4] = 0.1001;
RATES[3] = 0.1001;
RATES[2] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - - ( 1000.00*(ALGEBRAIC[0]+ALGEBRAIC[10]+ (CONSTANTS[4]/CONSTANTS[6])*(STATES[0] - STATES[1])))/CONSTANTS[5];
resid[1] = RATES[1] - - ( 1000.00*(ALGEBRAIC[3]+ALGEBRAIC[11]+ALGEBRAIC[1]+ (CONSTANTS[4]/(1.00000 - CONSTANTS[6]))*(STATES[1] - STATES[0])))/CONSTANTS[5];
resid[2] = RATES[4] -  CONSTANTS[10]*(- ALGEBRAIC[4] -  3.00000*ALGEBRAIC[11]);
resid[3] = RATES[3] - (ALGEBRAIC[7] - STATES[3])/ALGEBRAIC[8];
resid[4] = RATES[2] - (ALGEBRAIC[6] - STATES[2])/ALGEBRAIC[9];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[12] = ALGEBRAIC[3]+ALGEBRAIC[11]+ALGEBRAIC[1];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] =  CONSTANTS[2]*(STATES[0] - CONSTANTS[1])*pow(STATES[2], 2.00000);
ALGEBRAIC[1] =  CONSTANTS[13]*(STATES[1] - CONSTANTS[7]);
ALGEBRAIC[3] =  (CONSTANTS[11]/(1.00000+ (CONSTANTS[16]/CONSTANTS[17])*exp(- (STATES[1]/CONSTANTS[18]))))*(STATES[1] - CONSTANTS[8]);
ALGEBRAIC[4] =  (CONSTANTS[12]/(1.00000+ (CONSTANTS[16]/CONSTANTS[17])*exp(- (STATES[1]/CONSTANTS[18]))))*(STATES[1] - CONSTANTS[0]);
ALGEBRAIC[6] = 1.00000/(1.00000+exp(- (STATES[0]+31.0000)/5.30000));
ALGEBRAIC[7] = 1.00000/(1.00000+exp((STATES[0]+30.0000)/8.30000));
ALGEBRAIC[8] =  0.400000*(1.00000+2.00000/(1.00000+exp((STATES[0]+25.0000)/5.00000)));
ALGEBRAIC[9] = (0.800000+1.60000/(1.00000+exp( 0.100000*(STATES[0]+25.0000))))/(1.00000+exp( - 0.100000*(STATES[0]+70.0000)));
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[0]+35.0000)/6.20000));
ALGEBRAIC[10] =  CONSTANTS[3]*STATES[3]*(STATES[0] - CONSTANTS[0])*pow(ALGEBRAIC[5], 3.00000);
ALGEBRAIC[2] = pow(STATES[4], 3.00000)/(pow(STATES[4], 3.00000)+pow(CONSTANTS[15], 3.00000));
ALGEBRAIC[11] =  CONSTANTS[9]*(ALGEBRAIC[2] - CONSTANTS[19]);
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}