www.neuron.yale.edu

Dear developers,

how expensive is one Jacobian evaluation in the adaptive solver in NEURON implementation? Can it be compared to one ode call or is it rather multiplied by the dimension of the system?

thank you in advance,
Brocke

It is the same order as evaluation of an y' = f(y) call. That is because solving the linear system using the (approximate)
jacobian has exactly the
same number of arithmetic operations as solving a tridiagonal system with the same rank.