/*
   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 17 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (hour).
 * ALGEBRAIC[0] is Qv in component v (per_second).
 * STATES[0] is Vv in component v (mV).
 * CONSTANTS[0] is tau_v in component v (second).
 * CONSTANTS[1] is v_vm in component v (mV_second).
 * ALGEBRAIC[1] is Qm in component m (per_second).
 * CONSTANTS[2] is Qmax in component model_parameters (per_second).
 * ALGEBRAIC[3] is D in component D (mV).
 * CONSTANTS[3] is theta in component model_parameters (mV).
 * CONSTANTS[4] is sigma in component model_parameters (mV).
 * CONSTANTS[16] is Qa in component a (per_second).
 * CONSTANTS[14] is Va in component a (mV).
 * CONSTANTS[5] is Vao in component a (mV).
 * STATES[1] is Vm in component m (mV).
 * CONSTANTS[6] is tau_m in component m (second).
 * CONSTANTS[7] is v_mv in component m (mV_second).
 * CONSTANTS[8] is v_maQao in component m (mV).
 * STATES[2] is H in component H (nM).
 * CONSTANTS[9] is chi in component H (hour).
 * CONSTANTS[10] is mu in component H (nM_second).
 * ALGEBRAIC[2] is C in component D (dimensionless).
 * CONSTANTS[11] is c0 in component D (dimensionless).
 * CONSTANTS[15] is omega in component D (per_hour).
 * CONSTANTS[12] is v_vc in component D (mV).
 * CONSTANTS[13] is v_vh in component D (mV_per_nM).
 * RATES[0] is d/dt Vv in component v (mV).
 * RATES[1] is d/dt Vm in component m (mV).
 * RATES[2] is d/dt H in component H (nM).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 0.0;
CONSTANTS[0] = 10.0;
CONSTANTS[1] = -1.9;
CONSTANTS[2] = 100.0;
CONSTANTS[3] = 10.0;
CONSTANTS[4] = 3.0;
CONSTANTS[5] = 1.0;
STATES[1] = 0.0;
CONSTANTS[6] = 10.0;
CONSTANTS[7] = -1.9;
CONSTANTS[8] = 1.0;
STATES[2] = 15.0;
CONSTANTS[9] = 10.8;
CONSTANTS[10] = 3.6;
CONSTANTS[11] = 1.0;
CONSTANTS[12] = -6.3;
CONSTANTS[13] = 0.19;
CONSTANTS[14] = CONSTANTS[5];
CONSTANTS[15] = ( 2.00000* 3.14159265358979)/24.0000;
CONSTANTS[16] = CONSTANTS[2]/(1.00000+exp(- (CONSTANTS[14] - CONSTANTS[3])/CONSTANTS[4]));
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = CONSTANTS[2]/(1.00000+exp(- (STATES[0] - CONSTANTS[3])/CONSTANTS[4]));
RATES[1] = ((CONSTANTS[8]+ CONSTANTS[7]*ALGEBRAIC[0]) - STATES[1])/(CONSTANTS[6]/3600.00);
ALGEBRAIC[1] = CONSTANTS[2]/(1.00000+exp(- (STATES[1] - CONSTANTS[3])/CONSTANTS[4]));
RATES[2] = ( CONSTANTS[10]*ALGEBRAIC[1] - STATES[2])/CONSTANTS[9];
ALGEBRAIC[2] = CONSTANTS[11]+cos( CONSTANTS[15]*VOI);
ALGEBRAIC[3] =  CONSTANTS[12]*ALGEBRAIC[2]+ CONSTANTS[13]*STATES[2];
RATES[0] = (( CONSTANTS[1]*ALGEBRAIC[1]+ALGEBRAIC[3]) - STATES[0])/(CONSTANTS[0]/3600.00);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = CONSTANTS[2]/(1.00000+exp(- (STATES[0] - CONSTANTS[3])/CONSTANTS[4]));
ALGEBRAIC[1] = CONSTANTS[2]/(1.00000+exp(- (STATES[1] - CONSTANTS[3])/CONSTANTS[4]));
ALGEBRAIC[2] = CONSTANTS[11]+cos( CONSTANTS[15]*VOI);
ALGEBRAIC[3] =  CONSTANTS[12]*ALGEBRAIC[2]+ CONSTANTS[13]*STATES[2];
}