from templates.golgi.Golgi_template import Goc
from neuron import h
from neuron import gui
import numpy as np
h.nrncontrolmenu()
#h.load_file("/home/bremen/nrn7.2alpha/nrn/share/lib/hoc/parcom.hoc")
##pc = h.ParallelContext()
#h.ParallelComputeTool[0].nthread(2)
nmf = 0 #100
npf = 00 #(30000)
nmli = 00 #300
cellula = Goc(record_all=True)
h.celsius = 37
vm = h.Vector()
tvec = h.Vector()
vm.record(cellula.soma(.5)._ref_v, sec=cellula.soma)
tvec.record(h._ref_t)
stim = [
h.IClamp(cellula.soma(.5)),
h.IClamp(cellula.soma(.5)),
h.IClamp(cellula.soma(.5)),
h.IClamp(cellula.soma(.5))
]
stim[0].delay = 10000
from matplotlib import pyplot as plt from templates.golgi.Golgi_template import Goc from templates.gapj.gapj import Gapj import toolbox as tb from neuron import h from neuron import gui import random as rnd h.nrncontrolmenu() cvode = h.CVode() cvode.active(0) gocs = Goc(np.array([0., 0., 0.]), record_all=1, Dt=0.1) stimdata = dict() stimdata['stim0del'] = 2000 stimdata['stim0dur'] = 1000 stimdata['stim0amp'] = 0.9 # Nano #~ stimdata['stim1del'] = 3100 #~ stimdata['stim1dur'] = 1000 #~ stimdata['stim1amp'] = 0.5 stimdata['timeglobal'] = 4000 h.dt = 0.025 h.celsius = 37 h.tstop = stimdata['timeglobal'] h.v_init = -65
sys.path.append('../../')
import h5py
import numpy as np
from matplotlib import pyplot as plt
from templates.golgi.Golgi_template import Goc
from templates.gapj.gapj import Gapj
import toolbox as tb
from neuron import h
from neuron import gui
h.nrn_load_dll('$i386/')
gocs = Goc(np.array([0.,0.,0.]),record_all=1)
# stim = h.IClamp(.5, sec=gocs[0][1][1].soma)
# stim.amp = -0.015
# stim.dur = 200
gocs.createsyn(ngoc=1,record_all=1)
vc = h.VClamp(0.5)
vc.amp[0] = 0
vc.dur[0] = 200
stim = h.NetStim(0.5)
stim.number = 1
stim.interval = 10
stim.start = 20
from templates.golgi.Golgi_template import Goc from templates.gapj.gapj import Gapj import toolbox as tb from neuron import h from neuron import gui import random as rnd h.nrncontrolmenu() cvode = h.CVode() cvode.active(0) gocs = Goc(np.array([0.,0.,0.]),record_all=1, Dt = 0.1) # stim = h.IClamp(.5, sec=gocs[0][1][1].soma) # stim.amp = -0.015 # stim.dur = 200 npf = 1 #npf = 1 #npf = 4257 #npf = 3912 #npf = 1997 ngrc = 0 #ngrc = 383 #ngrc = 173 #ngrc = 1
sys.path.append('../../')
import h5py
import numpy as np
from matplotlib import pyplot as plt
from templates.golgi.Golgi_template import Goc
from templates.gapj.gapj import Gapj
import toolbox as tb
from neuron import h
from neuron import gui
h.nrn_load_dll('$i386/')
gocs = Goc(np.array([0., 0., 0.]), record_all=1, Dt=0.025)
# stim = h.IClamp(.5, sec=gocs[0][1][1].soma)
# stim.amp = -0.015
# stim.dur = 200
gocs.createsyn(nmli=1, record_all=1, Dt=0.025)
vc = h.VClamp(0.5, sec=gocs.soma)
vc.amp[0] = -70
vc.dur[0] = 1000
# cc = h.IClamp(0.5)
# cc.amp = -0.065
# cc.dur = 2000
# cc.delay = 0
from templates.golgi.Golgi_template import Goc
from templates.gapj.gapj import Gapj
import toolbox as tb
from mpi4py import MPI
comm = MPI.COMM_WORLD
from neuron import h
from neuron import gui
pc = h.ParallelContext()
rank = int(pc.id())
size = int(pc.nhost())
h.nrn_load_dll('$i386/')
gocs = [tb.makecell(Goc(np.array([0., 0., 0.])), rank, rank, rank, pc)]
print gocs
if (rank == 1):
stim = h.IClamp(.5, sec=gocs[0][1][1].soma)
stim.amp = -0.015
stim.dur = 200
gocs[0][1][1].creategap(1, gocs[0][1][0], 50, pc, record_all=1)
print "ngapj ", len(gocs[0][1][1].GAP_L)
gocs[0][1][1].GAP_L[0].gapj.g = 130 * 9 # pS
for g in gocs[0][1][1].GAP_L:
print "gapj g", g.gapj.g
goc2goc_gap_nc = pc.target_var(gocs[0][1][1].GAP_L[0].gapj._ref_vgap, not rank)