This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: ODE NAN INF derivs
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Kevin Hobbs <kevin dot hobbs dot 1 at ohiou dot edu>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Mon, 12 Apr 2004 19:44:50 +0100
- Subject: Re: ODE NAN INF derivs
- References: <40758412.7080101@ohiou.edu>
Kevin Hobbs writes:
> I'd like advice on adding a way for my derivative function to reject the
> step, say by returning GSL_ERANGE, without terminating the model.
Your function(s) should return an error code instead of GSL_SUCCESS if
you want to avoid proceeding with NaN values. The library should
propagate the error code back up to the calling level, so you can
retry with a different step-size. If it doesn't, that is a bug.
--
Brian Gough