HHNeuron.cpp

Go to the documentation of this file.

#include "HHNeuron.h"

// HHNeuron

void HHNeuron::derivatives(double t, const double y[], double f[]) 
{
    // y[0] ... Vm, y[1] ... m, y[2] ... h, y[3] ... n, y[4] ... nM
    double V = y[0];
/*
    if (nStepsInRefr>0) {
       f[0] = 0.0;
    }
    else
*/
    f[0] = 1.0/Cm*(- gbar_Na*pow(y[1], 3)*y[2]*(V-E_Na) 
                   - gbar_K*pow(y[3], 4)*(V-E_K) 
                   - gbar_M*y[4]*(V-E_M) 
                   + I0 
                   - G0*V);

    double v1 = V*1000.0;
    double v2 = v1 + 63.0;

    f[1] = alpha_m(v2) * (1.0 - y[1]) - beta_m(v2) * y[1];
    f[2] = alpha_h(v2) * (1.0 - y[2]) - beta_h(v2) * y[2];
    f[3] = alpha_n(v2) * (1.0 - y[3]) - beta_n(v2) * y[3];
    f[4] = alpha_nM(v1) * (1.0 - y[4]) - beta_nM(v1) * y[4];
}


int HHNeuron::reset(double dt) 
{
    double y0[5];
    double v1 = Vinit * 1000.0;
    double v2 = v1 + 63.0;

    y0[0] = Vinit;
    y0[1] = infty_m(v2); 
    y0[2] = infty_h(v2); 
    y0[3] = infty_n(v2);
    y0[4] = infty_nM(v1);

    double gNa = gbar_Na * pow(infty_m(v2), 3) * infty_h(v2);
    double gK  = gbar_K  * pow(infty_n(v2), 4);
    double gM  = gbar_M  * infty_nM(v1);

    double Gtot_init = 1.0/Rm + gNa + gK + gM;
    double Ich_init = gNa*E_Na + gK*E_K + gM*E_M;

    Em = Rm * (Vresting * Gtot_init - Ich_init);

    ODESystemBasedSpikingNeuron::reset(dt, y0);

    return 0;
}


// HHINeuron

void HHINeuron::derivatives(double t, const double y[], double f[]) 
{
    // y[0] ... Vm, y[1] ... m, y[2] ... h, y[3] ... n
    double V = y[0];

/*
    if (nStepsInRefr>0) {
       f[0] = 0.0;
    }
    else
*/
    f[0] = 1.0/Cm*(- gbar_Na*pow(y[1], 3)*y[2]*(V-E_Na) 
                   - gbar_K*pow(y[3], 4)*(V-E_K) 
                   + I0 
                   - G0*V);

    double v2 = V*1000.0 + 63.0;

    f[1] = alpha_m(v2) * (1 - y[1]) - beta_m(v2) * y[1];
    f[2] = alpha_h(v2) * (1 - y[2]) - beta_h(v2) * y[2];
    f[3] = alpha_n(v2) * (1 - y[3]) - beta_n(v2) * y[3];
}


int HHINeuron::reset(double dt) 
{
    double y0[4];
    double v2 = Vinit * 1000.0 + 63.0;

    y0[0] = Vinit;
    y0[1] = infty_m(v2);
    y0[2] = infty_h(v2);
    y0[3] = infty_n(v2);

    double gNa = gbar_Na * pow(infty_m(v2), 3) * infty_h(v2);
    double gK  = gbar_K  * pow(infty_n(v2), 4);

    double Gtot_init = 1.0/Rm + gNa + gK;
    double Ich_init = gNa*E_Na + gK*E_K;

    Em = Rm * (Vresting * Gtot_init - Ich_init);

    ODESystemBasedSpikingNeuron::reset(dt, y0);

    return 0;
}

---