I have found a strange behavior in CVode that I figured I would ask about. It seems like it might be a bug, but it's also quite possible I'm doing something wrong.
I have a six-compartment model, specified in the file lp-40.hoc, which uses a number of mechanisms that I've written. I've tested the mechanisms and the model in the NEURON GUI, so I'm reasonably confident they work properly.
However, if I start up NEURON (using nrniv, not nrngui), and execute
Code: Select all
load_file("lp-40.hoc")
objref cvode
cvode = new CVode(1)
finitialize(-50)
cvode.solve(2000)
If, on the other hand, I do:
Code: Select all
load_file("lp-40.hoc")
objref cvode
cvode = new CVode()
cvode.active(1)
finitialize(-50)
cvode.solve(2000)
This seems to imply that doing
Code: Select all
cvode = new CVode(1)
Code: Select all
cvode = new CVode()
cvode.active(1)
My version of NEURON is
Code: Select all
NEURON -- Version 5.8 2005-10-14 12:36:20 Main (88)
Is this behavior odd, or am I doing something wrong?
Thank you,
Adam