[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] Jacobian evaluation in ODE suite

From: Claude Lacoursiere
Subject: [Help-gsl] Jacobian evaluation in ODE suite
Date: Thu, 23 Sep 2010 17:05:59 +0200


I'm using the ODE module in a class I teach and I was surprised that the
Jacobian function is unused for all integration methods with the exception of

I tried all methods on the BZ reaction which is very stiff, and none of


used the provided Jacobian function a single time.

Using either the Matlab toolbox or the Octave odepkg, all implicit methods need
to evaluate the Jacobian to produce a decent solution for this
equation.  True, neither
of these have the Gauss'  colocation method available in GSL.

Does anyone know how to force the evaluation of Jacobians?  Is this a bug?

Attached is the code if anyone is curious.


Attachment: bz.c
Description: Text Data

reply via email to

[Prev in Thread] Current Thread [Next in Thread]