From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 14387 invoked by alias); 1 Mar 2009 14:37:50 -0000 Received: (qmail 14377 invoked by uid 22791); 1 Mar 2009 14:37:49 -0000 X-SWARE-Spam-Status: No, hits=-1.9 required=5.0 tests=AWL,BAYES_00 X-Spam-Check-By: sourceware.org Received: from smtp6.welho.com (HELO smtp6.welho.com) (213.243.153.40) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Sun, 01 Mar 2009 14:37:43 +0000 Received: from [10.0.0.133] (cs27048245.pp.htv.fi [89.27.48.245]) by smtp6.welho.com (Postfix) with ESMTP id D17D85BC00F for ; Sun, 1 Mar 2009 16:37:40 +0200 (EET) Message-ID: <49AA9DB5.6030908@iki.fi> Date: Sun, 01 Mar 2009 14:37:00 -0000 From: Tuomo Keskitalo User-Agent: Icedove 1.5.0.14eol (X11/20090105) MIME-Version: 1.0 To: 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> <498727E5.6080407@iki.fi> In-Reply-To: <498727E5.6080407@iki.fi> 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/msg00015.txt.bz2 Hello, On 02/02/2009 07:05 PM, Tuomo Keskitalo wrote: > 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. Ok, I've made a new beta release which includes new msbdf stepper and modifications to modnewton1 and the related steppers, among other changes. Please look at ode-initval/ChangeLog-Keskitalo for more information. Comments are welcome. git source can be cloned to current directory with command git clone http://iki.fi/tuomo.keskitalo/gsl/ode-initval-additions/.git > 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! This still applies. msadams and msbdf are somewhat complex steppers (multistep & multiorder, Nordsieck form) and need to be tested. I think that if I later on have time to implement Radau II A, it will not introduce any further changes to the framework. Stability detection should still be added to msbdf, otherwise it is feature-complete. Thereby, I hope to hear your opinion about this development branch. Can it be incorporated into GSL, or should it be e.g. an extension library? -- Tuomo.Keskitalo@iki.fi http://iki.fi/tuomo.keskitalo