This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Providing the Jacobian in odeiv
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Fleur Kelpin <fleur at bio dot vu dot nl>
- Cc: <gsl-discuss at sources dot redhat dot com>
- Date: Thu, 6 Dec 2001 20:44:36 +0000 (GMT)
- Subject: Re: Providing the Jacobian in odeiv
- References: <Pine.GSO.4.33.0112061728540.16991-100000@monodon>
Fleur Kelpin writes:
> I was wondering why a gsl_odeiv_system expressly needs a jacobian matrix
> to be defined, while only the bsimp solver appears to be using it. The
> current formulation makes it sound as if the gsl solvers cannot be used
> if no jacobian is provided.
Hi,
The interface is designed for the general case, so it includes the
jacobian. The algorithms are then all interchangeable.
Of course, you are free to pass a null pointer for the Jacobian
function if you know it won't be needed, e.g. if you use only RK4.
regards
Brian