FS_neuron::FS_neuron(const vector<double>& param)
{
    // set neuron params from vector
    V = param[0];
    ENa = param[1]; // mv Na reversal potential
    EK  = param[2]; // mv K reversal potential
    El  = param[3]; // mv Leakage reversal potential
    gbarNa = param[4]; // mS/cm^2 Na conductance
    gbarK  = param[5];  // mS/cm^2 K conductance
    gl = param[6];
    fi = param[7];
    Iextmean = param[8];
    variance = param[9];

    unsigned int seed = chrono::system_clock::now().time_since_epoch().count();
    generator =  new default_random_engine (seed);
    normRand = new normal_distribution <double> (Iextmean, variance);
    Iext = (*normRand)(*generator); // Iext generate from normal distribution

    m = alpha_m() / (alpha_m() + beta_m());
    n = alpha_n() / (alpha_n() + beta_n());
    h = alpha_h() / (alpha_h() + beta_h());

    gNa = gbarNa*m*m*m*h;
    gK = gbarK*n*n*n*n;

    Isyn = 0;
    countSp = true;
    th = -20;
}
Example #2
0
double tau_h(const double V_tilde)	{ return 1./(alpha_h(V_tilde)+beta_h(V_tilde)); }
Example #3
0
double h_inf(const double V_tilde)	{ return alpha_h(V_tilde)/(alpha_h(V_tilde)+beta_h(V_tilde)); }
Example #4
0
static inline double tau_h(double V) {
	return 1.0 / (alpha_h(V) + beta_h(V));
}
Example #5
0
static inline double h_inf(double V) {
	return alpha_h(V) / (alpha_h(V) + beta_h(V));
}
double FS_neuron::h_integrate(double dt)
{
    double h_0 = alpha_h() / (alpha_h() + beta_h());
    double tau_h = 1 / (alpha_h() + beta_h());
    return h_0 -(h_0 - h)*exp(-dt/tau_h);
}