From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 9515 invoked by alias); 31 Mar 2012 07:42:00 -0000 Received: (qmail 9507 invoked by uid 22791); 31 Mar 2012 07:41:58 -0000 X-SWARE-Spam-Status: No, hits=-2.6 required=5.0 tests=AWL,BAYES_00,KHOP_THREADED,RCVD_IN_DNSWL_NONE X-Spam-Check-By: sourceware.org Received: from gw02.mail.saunalahti.fi (HELO gw02.mail.saunalahti.fi) (195.197.172.116) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Sat, 31 Mar 2012 07:41:39 +0000 Received: from [192.168.0.104] (a91-156-39-232.elisa-laajakaista.fi [91.156.39.232]) by gw02.mail.saunalahti.fi (Postfix) with ESMTP id E58D1139797; Sat, 31 Mar 2012 10:41:33 +0300 (EEST) Message-ID: <4F76B52D.6000209@iki.fi> Date: Sat, 31 Mar 2012 07:42:00 -0000 From: Tuomo Keskitalo User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:11.0) Gecko/20120313 Thunderbird/11.0 MIME-Version: 1.0 To: Varad Deshmukh CC: gsl-discuss@sourceware.org Subject: Re: Contributing to GSL References: <4F728F55.2010404@cs.ucsb.edu> In-Reply-To: <4F728F55.2010404@cs.ucsb.edu> 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: 2012-q1/txt/msg00022.txt.bz2 Hello, the TODO in ode-initval2 directory is slightly outdated. A BDF method has already been implemented in msbdf stepper in gsl-1.15. The selection of steppers in ode-initval2 is pretty nice, I think, as far as dense general-purpose ode-iv methods go. I sometimes thought about trying to implement Radau 5 method described in the second book (Hairer, E., Wanner, G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems) but never did. If you want to fight with ODEs, you could try that. However, I cannot promise I'd have much time to contribute. DAEs (differential-algebraic equations) are a different story, even though they're close to ODEs. There's currently no DAE solver framework in GSL. It is something that would be very nice to have, though. If anyone is willing to try to design and implement one, feel free to do it! BR, Tuomo On 03/28/2012 07:11 AM, Varad Deshmukh wrote: > Hi, > I am Varad Deshmukh, a Master's student at UCSB majoring in Computer > Science. My research interests are in Scientific Computing and HPC, and > I am interested in contributing to GSL. I recently took up a course in > numerical simulation and I figured a good way to start would be looking > into ODE solvers. > > The TODO list mentions work to be done for adding other stepper > methods - especially the BDF methods for DAE's implemented in DASSL and > DASPK. Could you please elaborate more on which new methods you want to > focus on? I have been reading the code for the Runge-Kutta methods - the > driver, evolve and control formats - and have a rough idea about the > design. > > It would be great if I could get some help regarding GSL development in > general and ode-initval2 in particular. If you feel that I should work > with smaller bugs and fix them before going into a new implementation, I > would like to hear suggestions. > > Thank you, > Varad. > -- Tuomo.Keskitalo@iki.fi http://iki.fi/tuomo.keskitalo