GSL ode-initval development
Tuomo Keskitalo
Tuomo.Keskitalo@iki.fi
Sun Aug 31 07:26:00 GMT 2008
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
More information about the Gsl-discuss
mailing list