This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

Re: Providing the Jacobian in odeiv


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


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]