www.neuron.yale.edu

I will assume that you are using the code in https://www.neuron.yale.edu/ftp/ted/neuron/extracellular_stim_and_rec.zip . The amplitude and time course of the extracellular stimulus is calculated by the procedure stim_waveform() in the file stim.hoc. The functionality of that procedure is very limite...

Hello NEURON community: This post talks about the underlying equation of AlphaSynapse. https://www.sas.upenn.edu/LabManuals/BBB251/NIA/NEUROLAB/HELP/alphasyn.htm I am wondering if the "v", postsynaptic voltage, is constant (and set as resting voltage) or should change during current injec...

I seem to have misplaced my scrying stone, and without that it is difficult to diagnose the unseen. Perhaps you could share what you mean by "add the 3D location" e.g. by providing code examples of what you did that worked, and what you did that didn't work.

Good Python suggestion, urid. Here are some NEURON-specific comments: First, just to make sure that everyone who reads this thread knows what we're writing about-- 1. NMODL's global variables. An NMODL-specified mechanism's PARAMETERs are global by default (have the same value in all segments in whi...

You're quite welcome, Shaina. Glad to help!

I've never seen this msg signify a discontinuity. It's usually caused by NMODL code with rate constants that involve exponential functions. Typically the exp is in a denominator where limiting the value returned by exp(x) to exp(700) will do no harm. Nothing will blow up numerically, and your simula...

Hi, Shaina. Sorry about the delay, but I've been traveling a bit after the course we did at ASU week before last. Open a terminal and cd into the the folder that holds your code. If there is a file called nrnmech.dll, delete it. Or if you see a folder called x86_64, delete that and all its contents....

The ODE for each segment is of the form transmembrane capacitive current + transmembrane ionic current = sum of all currents arriving in the segment from adjacent segments + sum of all currents that enter the segment from all point processes that generate currents Notice that current in these equati...

The primary utility of multisplit is to achieve load balance in network simulations, and that is useful only when the network contains a very few model cells that are very large. Otherwise the performance boost will not be worth the programmer's effort required to restructure your model. Few single ...

Runtime depends on model complexity. A crude but effective measure of complexity is the total number of ODEs. For every segment there is one ODE for the cable equation, plus one ODE for each state variable of every mechanism associated with that segment that requires numerical integration. This incl...

About this Forum: If you need to paste code, be sure to surround it with the proper markup (that's what the < / > formatting button is for). Otherwise things like leftsquarebracket i rightsquarebracket will be treated as formatting instructions (interpreted as "italic"). Here's how to fix ...

I asked because DBS protocols use biphasic, charge-balanced waveforms. Vector play works for arbitrary waveforms. If the conductive medium is nondispersive ("purely resistive"), then the time course of extracelluar potential vx(x,y,z) will be of the form A(x,y,z)*f(t) where f is a function...

Haven't been much help yet--you've done all the work so far. Implementing stimulation will be easy if you only need a single rectangular pulse--the IClamp class will do for that. But what if your stimulus current should have a more complex time course, e.g. two or more monophasic rectangular pulses,...

To provide specific advice, I will have to know what you're trying to do. All that is clear from the example code is that 1. it assumes that soma2 and dend2 have the same nseg value. This is unlikely to be true. 2. it assumes that this nseg value changes on every pass through the "for seg in se...

To give specific advice, I'll have to reproduce the error for myself. Can you email the mod file in question to ted dot carnevale at yale dot edu?