This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Jacobian matrix for ODEs
- From: Slaven Peles <peles at cns dot physics dot gatech dot edu>
- To: gsl-discuss at sources dot redhat dot com
- Date: Thu, 27 Jun 2002 10:15:10 -0400
- Subject: Jacobian matrix for ODEs
Hi,
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,
Slaven