Jacobian matrix for ODEs

Slaven Peles peles@cns.physics.gatech.edu
Thu Jun 27 10:55:00 GMT 2002


I am looking at the example for adaptive step ODE integrators in GSL 
Reference Manual (http://sources.redhat.com/gsl/ref/gsl-ref_25.html#SEC377) 
and I am trying to figure out how to read Jacobian matrix elements during the 
integration. For the fixed step integration array dfdy[4] with elements of 
the Jacobi matrix is explicitely defined in main function, so accessing its 
elements is straightforward. I still don't know how to do it when using 
adaptive stepsize integration. Only solution that I can think of right now is 
to invoke function jac whenever Jacobian values are needed, but it seems to 
me like a duplicating the work. Any suggestions would be greatly appreciated,

More information about the Gsl-discuss mailing list