Generated Code

/*
   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[8] 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[4] 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[9] is I_pump in component dendritic_compartment (uA_per_cm2).
 * ALGEBRAIC[10] is I_NMDA in component dendritic_compartment (uA_per_cm2).
 * ALGEBRAIC[11] 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[5] is phi_Na in component dendritic_compartment (dimensionless).
 * CONSTANTS[19] is phi_Na_eq in component dendritic_compartment (dimensionless).
 * ALGEBRAIC[2] is n_infinity in component gating_variables (dimensionless).
 * ALGEBRAIC[3] is h_infinity in component gating_variables (dimensionless).
 * ALGEBRAIC[6] is tau_h in component gating_variables (second).
 * ALGEBRAIC[7] 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).
 */
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));
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[3] = 1.00000/(1.00000+exp((STATES[0]+30.0000)/8.30000));
ALGEBRAIC[6] =  0.400000*(1.00000+2.00000/(1.00000+exp((STATES[0]+25.0000)/5.00000)));
RATES[3] = (ALGEBRAIC[3] - STATES[3])/ALGEBRAIC[6];
ALGEBRAIC[2] = 1.00000/(1.00000+exp(- (STATES[0]+31.0000)/5.30000));
ALGEBRAIC[7] = (0.800000+1.60000/(1.00000+exp( 0.100000*(STATES[0]+25.0000))))/(1.00000+exp( - 0.100000*(STATES[0]+70.0000)));
RATES[2] = (ALGEBRAIC[2] - STATES[2])/ALGEBRAIC[7];
ALGEBRAIC[4] = 1.00000/(1.00000+exp(- (STATES[0]+35.0000)/6.20000));
ALGEBRAIC[8] =  CONSTANTS[3]*STATES[3]*(STATES[0] - CONSTANTS[0])*pow(ALGEBRAIC[4], 3.00000);
ALGEBRAIC[0] =  CONSTANTS[2]*(STATES[0] - CONSTANTS[1])*pow(STATES[2], 2.00000);
RATES[0] = - ( 1000.00*(ALGEBRAIC[0]+ALGEBRAIC[8]+ (CONSTANTS[4]/CONSTANTS[6])*(STATES[0] - STATES[1])))/CONSTANTS[5];
ALGEBRAIC[1] =  CONSTANTS[13]*(STATES[1] - CONSTANTS[7]);
ALGEBRAIC[5] = pow(STATES[4], 3.00000)/(pow(STATES[4], 3.00000)+pow(CONSTANTS[15], 3.00000));
ALGEBRAIC[9] =  CONSTANTS[9]*(ALGEBRAIC[5] - CONSTANTS[19]);
ALGEBRAIC[10] =  (CONSTANTS[11]/(1.00000+ (CONSTANTS[16]/CONSTANTS[17])*exp(- (STATES[1]/CONSTANTS[18]))))*(STATES[1] - CONSTANTS[8]);
RATES[1] = - ( 1000.00*(ALGEBRAIC[10]+ALGEBRAIC[9]+ALGEBRAIC[1]+ (CONSTANTS[4]/(1.00000 - CONSTANTS[6]))*(STATES[1] - STATES[0])))/CONSTANTS[5];
ALGEBRAIC[11] =  (CONSTANTS[12]/(1.00000+ (CONSTANTS[16]/CONSTANTS[17])*exp(- (STATES[1]/CONSTANTS[18]))))*(STATES[1] - CONSTANTS[0]);
RATES[4] =  CONSTANTS[10]*(- ALGEBRAIC[11] -  3.00000*ALGEBRAIC[9]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[3] = 1.00000/(1.00000+exp((STATES[0]+30.0000)/8.30000));
ALGEBRAIC[6] =  0.400000*(1.00000+2.00000/(1.00000+exp((STATES[0]+25.0000)/5.00000)));
ALGEBRAIC[2] = 1.00000/(1.00000+exp(- (STATES[0]+31.0000)/5.30000));
ALGEBRAIC[7] = (0.800000+1.60000/(1.00000+exp( 0.100000*(STATES[0]+25.0000))))/(1.00000+exp( - 0.100000*(STATES[0]+70.0000)));
ALGEBRAIC[4] = 1.00000/(1.00000+exp(- (STATES[0]+35.0000)/6.20000));
ALGEBRAIC[8] =  CONSTANTS[3]*STATES[3]*(STATES[0] - CONSTANTS[0])*pow(ALGEBRAIC[4], 3.00000);
ALGEBRAIC[0] =  CONSTANTS[2]*(STATES[0] - CONSTANTS[1])*pow(STATES[2], 2.00000);
ALGEBRAIC[1] =  CONSTANTS[13]*(STATES[1] - CONSTANTS[7]);
ALGEBRAIC[5] = pow(STATES[4], 3.00000)/(pow(STATES[4], 3.00000)+pow(CONSTANTS[15], 3.00000));
ALGEBRAIC[9] =  CONSTANTS[9]*(ALGEBRAIC[5] - CONSTANTS[19]);
ALGEBRAIC[10] =  (CONSTANTS[11]/(1.00000+ (CONSTANTS[16]/CONSTANTS[17])*exp(- (STATES[1]/CONSTANTS[18]))))*(STATES[1] - CONSTANTS[8]);
ALGEBRAIC[11] =  (CONSTANTS[12]/(1.00000+ (CONSTANTS[16]/CONSTANTS[17])*exp(- (STATES[1]/CONSTANTS[18]))))*(STATES[1] - CONSTANTS[0]);
ALGEBRAIC[12] = ALGEBRAIC[10]+ALGEBRAIC[9]+ALGEBRAIC[1];
}