Re: GSL ode-initval development

[Tuomo Keskitalo <Tuomo dot Keskitalo at iki dot fi>]

On 08/18/2008 09:09 PM, Brian Gough wrote:

> At Sun, 17 Aug 2008 11:38:15 +0300,
> Tuomo Keskitalo wrote:
> > Implicit solvers result in a group of algebraic non-linear
> > equations, for which there exists several solution strategies.
> If it is only a small number (and they don't need any additional
> parameters) then it's simplest to use the existing framework with
> explicit names for each combination, e.g. gsl_odeiv_step_foo_newt,
> gsl_odeiv_step_foo_imp, etc

I think there will not be too many (hopefully), and your suggestion is OK.

I have been writing a non-linear equation solver based on a modified 
Newton iteration method, and a question about failures has come up. What 
value should the non-linear equation solver return to the stepper (and 
from there to user) in the case that the method reaches a maximum number 
of Newton iterations and therefore fails? GSL_CONTINUE or GSL_FAILURE or 
something else?

--
Tuomo.Keskitalo@iki.fi
http://iki.fi/tuomo.keskitalo