initial=kwargs[key],
atolscale=1 * nM)
def recorder(ptr):
"""return a vector that records the pointer"""
vec = h.Vector()
vec.record(ptr)
return vec
rxd.set_solve_type(dimension=3)
cell = h.Section(name="cell")
cell.diam = cell.L = 5
r = rxd.Region(h.allsec())
rxd.nthread(4)
declare_parameters(
vsP=1.1 * nM / hour,
vmP=1.0 * nM / hour,
KmP=0.2 * nM,
KIP=1.0 * nM,
ksP=0.9 / hour,
vdP=2.2 * nM / hour,
KdP=0.2 * nM,
vsT=1.0 * nM / hour,
vmT=0.7 * nM / hour,
KmT=0.2 * nM,
KIT=1.0 * nM,
ksT=0.9 / hour,
vdT=3.0 * nM / hour,
from neuron import h, rxd
from neuron.rxd import v
from neuron.rxd.rxdmath import exp, log
import numpy
h.load_file("stdrun.hoc")
# set the number of rxd threads
rxd.nthread(2)
# set the seed
numpy.random.seed(63245)
# enable extracellular rxd (only necessary for NEURON prior to 7.7)
rxd.options.enable.extracellular = True
# Model parameters
e = 1.60217662e-19
scale = 1e-14 / e
h.celsius = 37
# Persistent Na current parameters
gnap = 0.4e-3 * scale
thmp = -40
thhp = -48
sigmp = 6
sighp = -6
nap_minf = 1.0 / (1.0 + exp(-(v - thmp) / sigmp))
nap_hinf = 1.0 / (1.0 + exp(-(v - thhp) / sighp))
# Fast Na current parameters