From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 32357 invoked by alias); 31 Aug 2008 07:26:27 -0000 Received: (qmail 32346 invoked by uid 22791); 31 Aug 2008 07:26:26 -0000 X-Spam-Check-By: sourceware.org Received: from smtp4.pp.htv.fi (HELO smtp4.pp.htv.fi) (213.243.153.38) by sourceware.org (qpsmtpd/0.31) with ESMTP; Sun, 31 Aug 2008 07:25:31 +0000 Received: from [10.0.0.133] (cs27054090.pp.htv.fi [89.27.54.90]) by smtp4.pp.htv.fi (Postfix) with ESMTP id 82EAB5BC01C; Sun, 31 Aug 2008 10:25:26 +0300 (EEST) Message-ID: <48BA4766.6070701@iki.fi> Date: Sun, 31 Aug 2008 07:26:00 -0000 From: Tuomo Keskitalo User-Agent: Icedove 1.5.0.14eol (X11/20080724) MIME-Version: 1.0 To: Brian Gough CC: gsl-discuss@sourceware.org Subject: Re: GSL ode-initval development References: <48A7E377.3020607@iki.fi> In-Reply-To: Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2008-q3/txt/msg00022.txt.bz2 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