Jacobian calculation

Anything that doesn't fit elsewhere.
Post Reply
brocke
Posts: 11
Joined: Fri Aug 16, 2013 9:21 am

Jacobian calculation

Post by brocke » Mon Apr 13, 2015 6:48 am

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

hines
Site Admin
Posts: 1580
Joined: Wed May 18, 2005 3:32 pm

Re: Jacobian calculation

Post by hines » Mon Apr 13, 2015 11:52 am

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.

Post Reply