www.neuron.yale.edu

Hi, I need a bit help with NEURON+Python. In hoc file, I can easily write new init procedure and it will be called from the run() function. But how can I substitute standard init() if I use Python? I've tried something like this: def myinit(): h.finitialize() h.fcurrent() h.frecord_init() h.init = m...

Thank you, Michael! It works. In other words h.stdinit() should be called just before pc.psolve(), otherwise pc.set_maxstep() won't be initialized and any spike won't be delivered. If this is correct, maybe it will be useful to add pc.pinit( maxstep) function in future version, for calling all requi...

In a parallelized implementation of a network model every spike source must be assigned a unique gid. This applies to instances of the VecStim class, which is a class of artificial spiking cell. This involves using ParallelContext's setgid2node method in order to associate the VecStim's gid with th...

Hi all, I'm writing a code for network which is activated by the input spike trains recorded experimentally or from another model. I load all presynaptic spike events into a vectors and then play these vectors by VecStim-s. Everything worked perfect until I ran simulation on parallel MPI cluster. Su...

Thank you, Ted. I just want set all segments without iterations like in hoc notation: soma g_pas = 0.05 However I've figured out how I can do this in python also: >>> soma = h.Section() >>> soma.insert('pas') <nrn.Section object at 0x96edc08> >>> soma.nseg=7 >>> print soma.g_pas 0.001 >>> for seg in...

Hi all, I have strange problem with using NEURON+Python. According Hines et.al. 2009 'NEURON and Python', parameters of inserted modules should be accessible through object's attribute with same name as SUFFIX. soma = h.Section() soma.insert ('pas') soma.pas.g = 0.0003 For example, I inserted a simp...

I don't try to connect two of them, but I've tried activate one neuron by Exp2Syn controlled by netstim. So it works. I wouldn't like publish this model on ModelDB until it is properly tested. So you can download it from my webpage http://nisms.krinc.ru/rth/izh/izhcur.zip and test it. I'll very appr...

correction: diam should be 1/PI

Bill, It seems it works. Here alpha version for testing. Please take a look. I didn't check all 21 modes, but what I've checked works perfect. I'll appreciate for any comments and tests. file izhcur.mod: TITLE IzhikevichCurrent COMMENT This dirty trick implements Izhikevich model as a non-specific c...

maybe we can insert two modules, one for current and the other one for threshold (POINT_PROCESS), but how the latter one will reach the u variable in current module, I don't know....

yes Bill, because NET_RECEIVE section can exist only in POINT_PROCESS module.

Well, here is a first shot. Module for current, izhcur.mod TITLE IzhikevichCurrent COMMENT It is a dirty trick for implementation of Izhikevich model as a non-specific current. This mod file treats Izhikevich model as v'=(Iizh(v)+i)/cm where Iizh=e*v^2+f*v+g-u; u is treated as gating variable Do not...

Thank you, Bill. Right now I have a question. You use vv, because v is embedded in standard segment. Why can we not use variable v and i in standard segment for Izh. model? We can rewrite his model as v'=(Izh(v)+I)/c and where Izh(v)=e*v*v+f*v+g+u - is the non-specific current, c is equal 1 and u is...

I'm trying to make a mixed model, where some multicompartment or single-compartment cells may be arbitrarily replaced by Izhikevich model(39948). Conditions of Izh. model should be the same as for conduction-based cells: the same Exp2Syn synapses, InNp noise current generator and so on. I'm not sure...

Hi Ted, Maybe a little bit off-top. I use NetStim->NetCon->Exp2Syn system to get a barrage of synaptic inputs on distal dendrites of 8 compartments pyramidal cell model. (The code is too long and too simple to post it here) It worked very nice and did what I want until I tried to run simulation with...