A question about parallel Neuron:
The parallel network model I simulate doesn't seem to like too big synaptic delays. I used synaptic delays up to 15 with set_maxstep(100).
To be precise, there seems to be a problem with a gaba mech that I use.
Below you find the error message and the respective gaba.mod.
[A test simulation running the exact same code only with less cells and shorter run time on a serial machine runs without any problems!]
Is there a problem with the mod file?
I'm using Neuron VERSION 6.0.867 (1738) 2007-05-10 on a beowulf cluster running SUSE Linux 10.
Many thanks for some help!
UB
The error message:
Code: Select all
Assertion failed: file gaba.c, line 242
23 /usr/local/Cluster-Apps/nrn/nrn-6.0-867/mx-noiv-mpi-gcc/x86_64/bin/nrniv: _tsav <= t
23 in p_DA3Run.hoc near line 218
23 prun()
^
23 ParallelContext[0].psolve(30000)
23 ParallelNetManager[0].psolve(30000)
23 prun()
[23] MPI Abort by user Aborting program !
MX:node012:send req(already completed):req status 8:Remote endpoint is closed
type: 1 (send_small)
state (0x10):
dead
mcp_handle : 59
seg:0x8edee0,8
dest:peer_index=5,eid=2,seqnum=25185
slength=8,xfer_length=8
caller: 0x518f80
This is the respective gaba.mod:
Code: Select all
TITLE gaba synapse
NEURON {
POINT_PROCESS gaba
NONSPECIFIC_CURRENT i
RANGE g,a,b,gGABAmax,tauD,tauF,util
}
UNITS {
(uS) = (microsiemens)
(nA) = (nanoamp)
(mV) = (millivolt)
}
PARAMETER {
tcon = .5 (ms)
tcoff = 5.0 (ms)
egaba = -75 (mV)
gGABAmax = 0 (uS)
tauD = 800 (ms)
tauF = 800 (ms)
util= .3
}
ASSIGNED {
v (mV)
i (nA)
g (uS)
factor
}
INITIAL {
a=0
b=0
factor=tcon*tcoff/(tcoff-tcon)
}
STATE {
a
b
}
BREAKPOINT {
SOLVE states METHOD cnexp
g = b-a
i = gGABAmax*g*(v-egaba)
}
DERIVATIVE states {
a' = -a/tcon
b' = -b/tcoff
}
NET_RECEIVE(wgt,R,u,tlast (ms),nspike) {
LOCAL x
if (nspike==0) { R=1 u=util }
else {
if (tauF>0) { u=util+(1-util)*u*exp(-(t-tlast)/tauF) }
if (tauD>0) { R=1+(R*(1-u)-1)*exp(-(t-tlast)/tauD) }
}
x=wgt*factor*R*u
state_discontinuity(a,a+x)
state_discontinuity(b,b+x)
tlast=t
nspike= nspike+1
}