Generated Code

/*
   There are a total of 21 entries in the algebraic variable array.
   There are a total of 11 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 (millisecond).
 * STATES[0] is V in component membrane (millivolt).
 * CONSTANTS[0] is Cm in component membrane (picoF).
 * CONSTANTS[1] is i_app in component membrane (picoA).
 * ALGEBRAIC[8] is i_Na in component sodium_current (picoA).
 * ALGEBRAIC[14] is i_Ca_L in component L_type_calcium_current (picoA).
 * ALGEBRAIC[15] is i_Ca_T in component T_type_calcium_current (picoA).
 * ALGEBRAIC[16] is i_K_dr in component delayed_rectifier_K_channel_current (picoA).
 * ALGEBRAIC[17] is i_M in component M_like_K_current (picoA).
 * ALGEBRAIC[19] is i_ir in component inward_rectifier_K_current (picoA).
 * ALGEBRAIC[20] is i_d in component inward_leak_current (picoA).
 * CONSTANTS[2] is g_Na in component sodium_current (nanoS).
 * CONSTANTS[3] is V_Na in component sodium_current (millivolt).
 * ALGEBRAIC[0] is O in component sodium_current (dimensionless).
 * STATES[1] is A in component A (dimensionless).
 * CONSTANTS[4] is k1 in component reaction_constants (first_order_rate_constant).
 * CONSTANTS[5] is k1_ in component reaction_constants (first_order_rate_constant).
 * ALGEBRAIC[1] is alpha in component reaction_constants (first_order_rate_constant).
 * ALGEBRAIC[9] is beta in component reaction_constants (first_order_rate_constant).
 * STATES[2] is D in component D (dimensionless).
 * STATES[3] is A_ in component A_ (dimensionless).
 * CONSTANTS[19] is a in component reaction_constants (dimensionless).
 * STATES[4] is D_ in component D_ (dimensionless).
 * CONSTANTS[6] is k2 in component reaction_constants (first_order_rate_constant).
 * CONSTANTS[7] is k2_ in component reaction_constants (first_order_rate_constant).
 * CONSTANTS[8] is V_Ca in component L_type_calcium_current (millivolt).
 * CONSTANTS[9] is g_Ca_L in component L_type_calcium_current (nanoS).
 * STATES[5] is m in component L_type_calcium_current_m_gate (dimensionless).
 * ALGEBRAIC[2] is m_infinity in component L_type_calcium_current_m_gate (dimensionless).
 * ALGEBRAIC[10] is tau_m in component L_type_calcium_current_m_gate (millisecond).
 * CONSTANTS[10] is Vh in component L_type_calcium_current_m_gate (millivolt).
 * CONSTANTS[11] is g_Ca_T in component T_type_calcium_current (nanoS).
 * STATES[6] is m in component T_type_calcium_current_m_gate (dimensionless).
 * STATES[7] is h in component T_type_calcium_current_h_gate (dimensionless).
 * ALGEBRAIC[3] is m_infinity in component T_type_calcium_current_m_gate (dimensionless).
 * ALGEBRAIC[11] is tau_m in component T_type_calcium_current_m_gate (millisecond).
 * ALGEBRAIC[4] is h_infinity in component T_type_calcium_current_h_gate (dimensionless).
 * CONSTANTS[12] is tau_h in component T_type_calcium_current_h_gate (millisecond).
 * CONSTANTS[13] is V_K in component delayed_rectifier_K_channel_current (millivolt).
 * CONSTANTS[14] is g_K_dr in component delayed_rectifier_K_channel_current (nanoS).
 * STATES[8] is n in component delayed_rectifier_K_channel_current_n_gate (dimensionless).
 * STATES[9] is h in component delayed_rectifier_K_channel_current_h_gate (dimensionless).
 * ALGEBRAIC[5] is n_infinity in component delayed_rectifier_K_channel_current_n_gate (dimensionless).
 * ALGEBRAIC[12] is tau_n in component delayed_rectifier_K_channel_current_n_gate (millisecond).
 * ALGEBRAIC[6] is h_infinity in component delayed_rectifier_K_channel_current_h_gate (dimensionless).
 * CONSTANTS[15] is tau_h in component delayed_rectifier_K_channel_current_h_gate (millisecond).
 * CONSTANTS[16] is g_M in component M_like_K_current (nanoS).
 * STATES[10] is n in component M_like_K_current_n_gate (dimensionless).
 * ALGEBRAIC[7] is n_infinity in component M_like_K_current_n_gate (dimensionless).
 * ALGEBRAIC[13] is tau_n in component M_like_K_current_n_gate (millisecond).
 * CONSTANTS[17] is g_ir in component inward_rectifier_K_current (nanoS).
 * ALGEBRAIC[18] is n_infinity in component inward_rectifier_K_current_n_gate (dimensionless).
 * CONSTANTS[18] is g_d in component inward_leak_current (nanoS).
 * RATES[0] is d/dt V in component membrane (millivolt).
 * RATES[1] is d/dt A in component A (dimensionless).
 * RATES[3] is d/dt A_ in component A_ (dimensionless).
 * RATES[2] is d/dt D in component D (dimensionless).
 * RATES[4] is d/dt D_ in component D_ (dimensionless).
 * RATES[5] is d/dt m in component L_type_calcium_current_m_gate (dimensionless).
 * RATES[6] is d/dt m in component T_type_calcium_current_m_gate (dimensionless).
 * RATES[7] is d/dt h in component T_type_calcium_current_h_gate (dimensionless).
 * RATES[8] is d/dt n in component delayed_rectifier_K_channel_current_n_gate (dimensionless).
 * RATES[9] is d/dt h in component delayed_rectifier_K_channel_current_h_gate (dimensionless).
 * RATES[10] is d/dt n in component M_like_K_current_n_gate (dimensionless).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -60;
CONSTANTS[0] = 14;
CONSTANTS[1] = 15;
CONSTANTS[2] = 60;
CONSTANTS[3] = 60;
STATES[1] = 1;
CONSTANTS[4] = 0.3;
CONSTANTS[5] = 0.03;
STATES[2] = 0;
STATES[3] = 0;
STATES[4] = 0;
CONSTANTS[6] = 0.001;
CONSTANTS[7] = 0.01;
CONSTANTS[8] = 100;
CONSTANTS[9] = 1.3;
STATES[5] = 0;
CONSTANTS[10] = 40;
CONSTANTS[11] = 0.94;
STATES[6] = 0;
STATES[7] = 0;
CONSTANTS[12] = 22;
CONSTANTS[13] = -80;
CONSTANTS[14] = 20;
STATES[8] = 0;
STATES[9] = 0;
CONSTANTS[15] = 1000;
CONSTANTS[16] = 4;
STATES[10] = 0;
CONSTANTS[17] = 1.71;
CONSTANTS[18] = 0.044;
CONSTANTS[19] =  pow((( CONSTANTS[4]*CONSTANTS[7])/( CONSTANTS[5]*CONSTANTS[6])), 1.0 / 2);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+86.4000)/4.70000));
RATES[7] = (ALGEBRAIC[4] - STATES[7])/CONSTANTS[12];
ALGEBRAIC[6] = 0.700000/(1.00000+exp(- (STATES[0]+35.0000)/10.0000))+0.300000;
RATES[9] = (ALGEBRAIC[6] - STATES[9])/CONSTANTS[15];
ALGEBRAIC[1] = 10.0000/(1.00000+exp(- (STATES[0]+6.00000)/10.0000));
ALGEBRAIC[9] = 10.0000/(1.00000+exp((STATES[0]+54.4000)/4.60000));
RATES[1] = ( ALGEBRAIC[1]*STATES[2]+ CONSTANTS[5]*STATES[3]) - ( ALGEBRAIC[9]*STATES[1]+ CONSTANTS[4]*STATES[1]);
RATES[3] = ( ALGEBRAIC[1]*CONSTANTS[19]*STATES[4]+ CONSTANTS[4]*STATES[1]) - ( (ALGEBRAIC[9]/CONSTANTS[19])*STATES[3]+ CONSTANTS[5]*STATES[3]);
RATES[2] = ( ALGEBRAIC[9]*STATES[1]+ CONSTANTS[7]*STATES[4]) - ( ALGEBRAIC[1]*STATES[2]+ CONSTANTS[6]*STATES[2]);
RATES[4] = ( (ALGEBRAIC[9]/CONSTANTS[19])*STATES[3]+ CONSTANTS[6]*STATES[2]) - ( ALGEBRAIC[1]*CONSTANTS[19]*STATES[4]+ CONSTANTS[7]*STATES[4]);
ALGEBRAIC[2] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[10])/12.0000));
ALGEBRAIC[10] = 5.00000/(exp((STATES[0]+15.0000)/25.0000)+exp(- (STATES[0]+15.0000)/25.0000));
RATES[5] = (ALGEBRAIC[2] - STATES[5])/ALGEBRAIC[10];
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (STATES[0] - 56.1000)/10.0000));
ALGEBRAIC[11] = 7.00000/(exp((STATES[0]+50.0000)/9.00000)+exp(- (STATES[0]+50.0000)/9.00000))+0.800000;
RATES[6] = (ALGEBRAIC[3] - STATES[6])/ALGEBRAIC[11];
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[0]+25.0000)/15.0000));
ALGEBRAIC[12] = 15.0000/(exp((STATES[0]+30.0000)/15.0000)+exp(- (STATES[0]+30.0000)/15.0000))+1.00000;
RATES[8] = (ALGEBRAIC[5] - STATES[8])/ALGEBRAIC[12];
ALGEBRAIC[7] = 1.00000/(1.00000+exp(- (STATES[0]+37.0000)/4.00000));
ALGEBRAIC[13] = 80.0000/(exp((STATES[0]+30.0000)/15.0000)+exp(- (STATES[0]+30.0000)/15.0000));
RATES[10] = (ALGEBRAIC[7] - STATES[10])/ALGEBRAIC[13];
ALGEBRAIC[0] = pow(STATES[1], 3.00000);
ALGEBRAIC[8] =  CONSTANTS[2]*ALGEBRAIC[0]*(STATES[0] - CONSTANTS[3]);
ALGEBRAIC[14] =  CONSTANTS[9]*pow(STATES[5], 2.00000)*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[15] =  CONSTANTS[11]*pow(STATES[6], 2.00000)*STATES[7]*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[16] =  CONSTANTS[14]*pow(STATES[8], 4.00000)*STATES[9]*(STATES[0] - CONSTANTS[13]);
ALGEBRAIC[17] =  CONSTANTS[16]*STATES[10]*(STATES[0] - CONSTANTS[13]);
ALGEBRAIC[18] = ( 0.800000*1.00000)/(1.00000+exp((STATES[0]+80.0000)/12.0000))+0.200000;
ALGEBRAIC[19] =  CONSTANTS[17]*ALGEBRAIC[18]*(STATES[0] - CONSTANTS[13]);
ALGEBRAIC[20] =  CONSTANTS[18]*(STATES[0] - CONSTANTS[8]);
RATES[0] = (CONSTANTS[1] - (ALGEBRAIC[8]+ALGEBRAIC[14]+ALGEBRAIC[15]+ALGEBRAIC[16]+ALGEBRAIC[17]+ALGEBRAIC[19]+ALGEBRAIC[20]))/CONSTANTS[0];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+86.4000)/4.70000));
ALGEBRAIC[6] = 0.700000/(1.00000+exp(- (STATES[0]+35.0000)/10.0000))+0.300000;
ALGEBRAIC[1] = 10.0000/(1.00000+exp(- (STATES[0]+6.00000)/10.0000));
ALGEBRAIC[9] = 10.0000/(1.00000+exp((STATES[0]+54.4000)/4.60000));
ALGEBRAIC[2] = 1.00000/(1.00000+exp(- (STATES[0] - CONSTANTS[10])/12.0000));
ALGEBRAIC[10] = 5.00000/(exp((STATES[0]+15.0000)/25.0000)+exp(- (STATES[0]+15.0000)/25.0000));
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (STATES[0] - 56.1000)/10.0000));
ALGEBRAIC[11] = 7.00000/(exp((STATES[0]+50.0000)/9.00000)+exp(- (STATES[0]+50.0000)/9.00000))+0.800000;
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[0]+25.0000)/15.0000));
ALGEBRAIC[12] = 15.0000/(exp((STATES[0]+30.0000)/15.0000)+exp(- (STATES[0]+30.0000)/15.0000))+1.00000;
ALGEBRAIC[7] = 1.00000/(1.00000+exp(- (STATES[0]+37.0000)/4.00000));
ALGEBRAIC[13] = 80.0000/(exp((STATES[0]+30.0000)/15.0000)+exp(- (STATES[0]+30.0000)/15.0000));
ALGEBRAIC[0] = pow(STATES[1], 3.00000);
ALGEBRAIC[8] =  CONSTANTS[2]*ALGEBRAIC[0]*(STATES[0] - CONSTANTS[3]);
ALGEBRAIC[14] =  CONSTANTS[9]*pow(STATES[5], 2.00000)*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[15] =  CONSTANTS[11]*pow(STATES[6], 2.00000)*STATES[7]*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[16] =  CONSTANTS[14]*pow(STATES[8], 4.00000)*STATES[9]*(STATES[0] - CONSTANTS[13]);
ALGEBRAIC[17] =  CONSTANTS[16]*STATES[10]*(STATES[0] - CONSTANTS[13]);
ALGEBRAIC[18] = ( 0.800000*1.00000)/(1.00000+exp((STATES[0]+80.0000)/12.0000))+0.200000;
ALGEBRAIC[19] =  CONSTANTS[17]*ALGEBRAIC[18]*(STATES[0] - CONSTANTS[13]);
ALGEBRAIC[20] =  CONSTANTS[18]*(STATES[0] - CONSTANTS[8]);
}