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
