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