Out of range error (when applying high-voltage pulses)
Posted: Mon Jan 22, 2024 3:39 am
Hi,
I am running a model by Cameron McIntyre of an axon (3810 from Model DB) which is stimulated by extracellular voltage pulses(Python + Neuron). Recently, when using large amplitude pulses I get the following error:
exp(703.768) out of range, returning exp(700)
exp(707.553) out of range, returning exp(700)
exp(711.301) out of range, returning exp(700)
exp(715.013) out of range, returning exp(700)
No more errno warnings during this execution
I added code in Python to specify voltage stimulation. When the voltage is greater than 100V, I will get the above error. This is similar to the issue in the following post.
FUNCTION vtrap(x,y) {
if (fabs(x/y) > 1e-6) {
vtrap = x/(exp(x/y) - 1)
} else {
vtrap = y*(1 - x/y/2) : from Taylor's series
}
}
The above error still exists.
Can anyone help?
I am running a model by Cameron McIntyre of an axon (3810 from Model DB) which is stimulated by extracellular voltage pulses(Python + Neuron). Recently, when using large amplitude pulses I get the following error:
exp(703.768) out of range, returning exp(700)
exp(707.553) out of range, returning exp(700)
exp(711.301) out of range, returning exp(700)
exp(715.013) out of range, returning exp(700)
No more errno warnings during this execution
I added code in Python to specify voltage stimulation. When the voltage is greater than 100V, I will get the above error. This is similar to the issue in the following post.
Especially, when I use the method in the post to modify the vtrap function to the following form:ny wrote: ↑Fri Feb 25, 2011 7:24 am Hi,
I am running a model by Cameron McIntyre of an axon which is stimulated by extracellular voltage pulses. Recently, when using large amplitude pulses I get the following error:
exp(700.68) out of range, returning exp(700)
exp(711.563) out of range, returning exp(700)
exp(721.082) out of range, returning exp(700)
exp(729.442) out of range, returning exp(700)
I have read other posts on the forum and I have checked the Axnode.mod mechanism and it uses "vtrap" functions to check for zero denominators.
Can you help?
Nada
FUNCTION vtrap(x,y) {
if (fabs(x/y) > 1e-6) {
vtrap = x/(exp(x/y) - 1)
} else {
vtrap = y*(1 - x/y/2) : from Taylor's series
}
}
The above error still exists.
Can anyone help?