www.neuron.yale.edu

Thank you for your response and coding example... This really helps! I did have one follow-up question concerning the way the node is set up. From the code you have: node = na[cyt].nodes(soma(0.5))[0] na_cur = rxd_flux_mA(na_current, node) k_cur = rxd_flux_mA(k_current, node) leak_cur = rxd_flux_mA(...

Hello, I was wondering if there was a way to record the contribution of individual mechanisms that all act on the same species. For example, intracellular calcium in the cytosol (which would be defined as an rxd.Species()) has multiple mechanisms that affect its concentration. The mechanisms can be ...

Hello, I wanted to ask if it is possible for rxd to keep track of currents that arise from the movement of ions from one compartment to another. For example, if there was an rxd.Rate() that defined a mechanism for calcium influx (extracellular to intracellular), would the movement of these calcium i...

Could you elaborate on what the indexing on the nodes represents? What would the difference between using Glu.nodes[0] versus Glu.nodes[1] be, for example?

Hello, I'd like to ask if two separate functionalities are supported in RxD: (1) Assignment of a variable by an explicit algebraic equation (rather than a derivative) at every timestep (2) Assignment of a variable (or rate) by a time-dependent function To illustrate, I have a snippet of code that ac...

No worries! Thank you for the help! I will update if I find a way to accomplish this.

For this particular problem, what do you want your array of precalculated values to drive? Probability of transmitter release as a function of time? Transmitter concentration in the synaptic cleft? The detailed time course of synaptic conductance? If none of these, then what? I am trying to get the...

Thanks for the reply! Presumably you are doing that with Vector play(). I am actually using pointer manipulation such that an NMODL mechanism is able to access the contents of a numpy array (see https://www.neuron.yale.edu/phpBB/viewtopic.php?t=4167 ). The aim is to have the NMODL mechanism "pl...

Hello, I am wondering if there is a way to implement superposition for waveforms that are "saved" in the form of a constant, predetermined array. For example, the basic exp2syn.mod mechanism defines conductance waveforms that are governed by a double exponential function with respective DE...

Hello, I am working on creating my own mechanism .mod file and running into large value errors such as: exp(1.86734e+184) out of range, returning exp(700) No more errno warnings during this execution Without getting into too much detail, the mechanism works completely fine without any "out of r...

Thank you Ted and Robert for your help! The numpy_element_ref and h.setpointer functions were exactly what I was looking for!

I greatly appreciate the efforts. For my purposes, they are simply a set of indexed values. However, I am also curious to see the method that would solve the problem in the breakpoint case as well for future applications. Thank you!

Thanks for the reply, Ted. If I understand correctly, the Vector.play() changes a variable's value based on an array of values with a corresponding time vector which designates WHEN the variable's value is changed. In other words, a SINGLE value is being changed over time. What I believe you were tr...

Hello, I am still pretty new to NEURON. I've been trying to look for/come up with a way to set the values of an NMODL array that has the same exact values across many instances of the same NMODL mechanism. The reason is because there is a large array that I would like to use in NMODL but setting the...

Thank you for the reply Ted, appreciate it! Will keep all the tips you suggested in mind moving forward and will have a look at the NEURON guide as well. The use of the BREAKPOINT block to affect charge balance is exactly what I needed.