From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 32608 invoked by alias); 2 Feb 2009 17:05:51 -0000 Received: (qmail 32598 invoked by uid 22791); 2 Feb 2009 17:05:50 -0000 X-SWARE-Spam-Status: No, hits=-1.7 required=5.0 tests=AWL,BAYES_00 X-Spam-Check-By: sourceware.org Received: from smtp5.welho.com (HELO smtp5.welho.com) (213.243.153.39) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Mon, 02 Feb 2009 17:05:44 +0000 Received: from [10.0.0.133] (cs27048245.pp.htv.fi [89.27.48.245]) by smtp5.welho.com (Postfix) with ESMTP id 5308F5BC021; Mon, 2 Feb 2009 19:05:41 +0200 (EET) Message-ID: <498727E5.6080407@iki.fi> Date: Mon, 02 Feb 2009 17:05:00 -0000 From: Tuomo Keskitalo User-Agent: Icedove 1.5.0.14eol (X11/20090105) MIME-Version: 1.0 To: Brian Gough CC: GSL Discuss Mailing List Subject: Re: ode-initval implicit solvers and development References: <48E25CA9.6080306@iki.fi> <490DE4BD.7070907@iki.fi> <497B00F6.2080400@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: 2009-q1/txt/msg00004.txt.bz2 Hello, On 02/01/2009 05:42 PM, Brian Gough wrote: > - Currently the control object is only used by the multistep adams > stepper. What other types of methods might need access to the error > control? All steppers that do some kind of iterative solution internally. For example the methods that use Newton iteration. > - In modnewton1.c there is the comment > > /* Stopping tolerance for Newton iteration */ > /* FIXME: tol should be user defined */ > > How would you handle the newton tolerance--do we need a way to pass > that in as well? Yes, I think that is needed. This would be done similarly to msadams, e.g. modnewton1 asks the tolerance via gsl_odeiv_control_errlev. I plan to implement the Newton iteration for BDF methods first, and then modify modnewton1. I'd also like to implement the implicit 5th order Radau method, which also needs Newton iteration. This will take some time. It might be best to wait until then to see what other changes I'll have to make to the ode-initval interfaces. Hopefully none. > Valgrind detects an uninitialised variable in msadams.c, which is > flagged with -Wall as well. I'll add initialization to those, to get rid of the warnings, in the next release. Thanks for your comments! I'd appreciate if anyone has the chance to test and report how msadams compares to other steppers in practice and if there are any problems with it! -- Tuomo.Keskitalo@iki.fi http://iki.fi/tuomo.keskitalo