www.neuron.yale.edu

Here is the mod file Ted recommended for implementing vector.play in a memory-efficient way: COMMENT VecSyn, a "synaptic mechanism" whose conductance is driven by Vector.play. Actual synaptic conductance is gs. gs is the product of a scale factor k and gnorm, where gnorm is in microsiemens...

The first option: a sampled time course of synaptic conductance.

I have a simulation where I need to play hundreds of vectors into synaptic conductances. These vectors are all just scaled versions of one another. Here is an outline of the code: signal=np.loadtxt(#some txt file that defines trace for synapses gmax) syn_gmax_vec=h.Vector(signal) syn_veclist = [] fo...

Also, here is a previous thread with code for a minimal working model to calculate the LFP:
viewtopic.php?f=8&t=3656&p=15713#p15713

Okay, now I'm attempting to incorporate my firing rate model within a toy parallel network model. I started with the parallel ring network outlined here: https://www.neuron.yale.edu/phpBB/viewtopic.php?f=12&t=4270 I then added a dummy section (for the firing rate model) to the 0th host. In gener...

Never mind, I found my stupid mistake: I forgot to connect the dendrite to the soma. I was tempted to take this post down to save the embarrassment, but I figured I'd leave it up for posterity (since I couldn't find the code for this Python parallelized ring network anywhere online). Here's the work...

I'm trying to implement the code for a minimal parallel ring network, from the 2019 NEURON Course (I could not find the code anywhere online, so I adapted it from the course book). When I run the program (in serial, for now), for some reason the code to stimulate the 0th cell has no effect. When I p...

Perfect, thanks Robert. I see that your solution is exactly what I used before in the LinearMechanism implementation, but I thought it worked then only because I was passing a vector as the argument (i.e., I didn't realize that was the general way to dereference a pointer). Thanks again, and here's ...

I decided to try to implement the program I last posted using mod files, rather than LinearMechanism. I almost have it, except for an issue with the callback function similar to what I encountered a few posts ago. When I define the runtime_callback (to be used by extra_scatter_gather), I want to sen...

Ah, got it. The solution is to pass the vector 'y' to gnabar_callback, then access the 0th element within the function itself. The way I was doing it before, I was essentially just passing a constant value to gnabar_callback. Here's the updated code: from neuron import h from neuron.units import ms ...

Okay, thanks for the replies. I understand better what's going on with LinearMechanism now. I took the toy model from a few replies ago, and I'm trying to tweak it by adding an HH soma and having a dialogue between it and the firing rate model. I added a callback to LinearMechanism that gives the FR...

The LinearMechanism class looks very appealing, but I want to make sure I thoroughly understand it. A few questions: 1. I noticed in the documentation for LinearMechanism, there is the following warning: "Does not work with the CVODE integrator but does work with the differential-algebraic solv...

Thank you both very much! That will certainly get me started.

I am developing a thalamocortical network model of sleep. Currently, various cortical ionic conductances (e.g. potassium currents) are altered by a "hand of God" that sets them to different values to induce different stages of sleep. I would like to develop the model so that the ionic cond...

Never mind, I figured out the problem...I only printed out the LFP trace to three decimal places, and apparently the signal was so small that it required greater precision to see the fluctuations. Everything is working as it should.