Page 1 of 1

Jacobian calculation

Posted: Mon Apr 13, 2015 6:48 am
by brocke
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

Re: Jacobian calculation

Posted: Mon Apr 13, 2015 11:52 am
by hines
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.