From: Tuomo Keskitalo <Tuomo.Keskitalo@iki.fi>
To: Brian Gough <bjg@gnu.org>
Cc: GSL Discuss Mailing List <gsl-discuss@sourceware.org>
Subject: Re: ode-initval implicit solvers and development
Date: Mon, 02 Feb 2009 17:05:00 -0000 [thread overview]
Message-ID: <498727E5.6080407@iki.fi> (raw)
In-Reply-To: <m38woqp3je.wl%bjg@network-theory.co.uk>
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
next prev parent reply other threads:[~2009-02-02 17:05 UTC|newest]
Thread overview: 64+ messages / expand[flat|nested] mbox.gz Atom feed top
2008-09-30 17:07 Tuomo Keskitalo
2008-10-01 18:29 ` Brian Gough
2008-10-09 13:22 ` Brian Gough
2008-11-02 17:35 ` Tuomo Keskitalo
2008-11-03 18:09 ` Brian Gough
2009-01-24 11:52 ` Tuomo Keskitalo
2009-02-01 17:01 ` Brian Gough
2009-02-02 17:05 ` Tuomo Keskitalo [this message]
2009-03-01 14:37 ` Tuomo Keskitalo
2009-03-03 16:34 ` Brian Gough
2009-03-05 19:47 ` Tuomo Keskitalo
2009-03-05 19:54 ` Heikki Orsila
2009-03-06 20:03 ` Brian Gough
2009-04-05 12:28 ` Tuomo Keskitalo
2009-05-01 14:05 ` Tuomo Keskitalo
2009-05-04 11:23 ` Brian Gough
2009-05-08 10:51 ` Brian Gough
2009-08-06 13:51 ` GSL 2.0 roadmap Tuomo Keskitalo
2009-08-21 20:42 ` Brian Gough
2009-08-27 11:42 ` Tuomo Keskitalo
2009-08-27 12:51 ` Robert G. Brown
2009-08-28 13:57 ` Jordi Burguet Castell
2009-08-27 17:13 ` Robert G. Brown
2009-08-28 13:58 ` Brian Gough
2009-08-27 23:10 ` Gerard Jungman
2009-08-27 23:13 ` GSL 2.0 roadmap (one man's view) Gerard Jungman
2009-08-28 13:58 ` Brian Gough
2009-09-16 0:43 ` Gerard Jungman
2009-09-03 19:37 ` Brian Gough
2009-09-16 0:44 ` Gerard Jungman
2009-09-07 15:10 ` Brian Gough
2009-09-16 0:46 ` Gerard Jungman
2009-09-16 2:48 ` Robert G. Brown
2009-09-17 19:14 ` Brian Gough
2009-09-07 15:10 ` Brian Gough
2009-09-16 0:47 ` Gerard Jungman
2009-09-27 8:03 ` new double precision data structure? Tuomo Keskitalo
2009-09-28 8:44 ` James Bergstra
2009-09-28 15:48 ` Tuomo Keskitalo
2009-10-16 13:59 ` Brian Gough
2009-09-29 18:38 ` Gerard Jungman
2009-09-07 15:10 ` GSL 2.0 roadmap (one man's view) Brian Gough
2009-09-16 0:46 ` Gerard Jungman
2009-09-17 20:12 ` Brian Gough
2009-09-07 15:10 ` Brian Gough
2009-09-16 0:44 ` Gerard Jungman
2009-09-17 20:12 ` Brian Gough
[not found] ` <645d17210909090818u474f32f0q19a6334578b9f02c@mail.gmail.com>
2009-09-17 19:14 ` Brian Gough
2009-09-07 15:10 ` Brian Gough
2009-09-07 15:34 ` Rhys Ulerich
2009-09-07 18:21 ` Robert G. Brown
2009-09-16 0:47 ` Gerard Jungman
2009-09-18 3:51 ` column-major Z F
2009-09-21 12:08 ` column-major Brian Gough
2009-09-07 15:10 ` GSL 2.0 roadmap (one man's view) Brian Gough
2009-09-16 0:45 ` Gerard Jungman
2009-09-20 9:36 ` Tuomo Keskitalo
2009-09-20 13:23 ` Robert G. Brown
2009-09-20 15:31 ` Rhys Ulerich
2009-09-20 16:19 ` Robert G. Brown
2009-09-21 15:13 ` Brian Gough
2009-09-20 15:08 ` Rhys Ulerich
2009-09-21 12:08 ` Brian Gough
2009-08-28 13:58 ` GSL 2.0 roadmap Brian Gough
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=498727E5.6080407@iki.fi \
--to=tuomo.keskitalo@iki.fi \
--cc=bjg@gnu.org \
--cc=gsl-discuss@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).