This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: ode-initval implicit solvers and development
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Tuomo Keskitalo <Tuomo dot Keskitalo at iki dot fi>
- Cc: GSL Discuss Mailing List <gsl-discuss at sourceware dot org>
- Date: Thu, 09 Oct 2008 14:21:29 +0100
- Subject: Re: ode-initval implicit solvers and development
- References: <48E25CA9.6080306@iki.fi>
At Tue, 30 Sep 2008 20:06:49 +0300,
Tuomo Keskitalo wrote:
> my work with GSL ODE solvers is progressing (slowly, I am doing this on
> my own free time). I've written a modified Newton iteration method to
> solve the non-linear equation system that come up with implicit
> Runge-Kutta methods. Three very basic implicit solvers included use the
> solver. Also, I've made some changes and additions to ode-initval tests.
This is good, I would like to replace the existing rkimp solvers with
these in a future release. The main outstanding change needed is to
avoid recomputing the jacobian J(t,x) in the second call to
modnewton1_init. The parameters t and x are the same so it computed
once in the caller and passed to modnewton1_init as a parameter.
--
Brian Gough