GNU Scientific Library - Bugs: bug #30540, please, add convergence checks in...
You are not allowed to post comments on this tracker with your current authentication level.
bug #30540: please, add convergence checks in rk4imp/rk2imp
Submitter: | Sergey B Kirpichev <skirpichev> | ||
Submitted: | Sat 24 Jul 2010 12:15:06 PM UTC | ||
Category: | Accuracy problem | Severity: | 3 - Normal |
Operating System: | Status: | Ready For Test | |
Assigned to: | None | Open/Closed: | Open |
Release: |
Thu 26 Aug 2010 12:05:28 PM UTC, comment #2: |
Sergey B Kirpichev <skirpichev> |
Fri 06 Aug 2010 01:09:11 PM UTC, comment #1: Thanks for the bug report.
|
-Deleted Account- <bjg> |
Sat 24 Jul 2010 12:15:06 PM UTC, original submission:
Quote from http://lists.gnu.org/archive/html/help-gsl/2009-12/msg00047.html
|
Sergey B Kirpichev <skirpichev> |
Depends on the following items: None found
Items that depend on this one: None found
Follow 6 latest changes.
Date | Changed by | Updated Field | Previous Value | => | Replaced by |
---|---|---|---|---|---|
2010-08-26 | skirpichev | Attached File | - | Added bug-ode2.c, #21314 | |
2010-08-06 | bjg | Status | None | Ready For Test | |
2010-07-24 | skirpichev | Attached File | - | Added bug.c, #21050 | |
Attached File | - | Added rk2imp.png, #21051 | |||
Attached File | - | Added rk4imp.png, #21052 | |||
Attached File | - | Added 0001-Add-convergence-checks-for-rk-imp.patch, #21053 |
Powered by Savane 3.13-f8d8.
Corresponding source code
I can't figure out how to use this code for fixed step size. imprk2/imprk4 seems to be not working at all for this case:
$ ./bug
Segmentation fault
Other steppers works (e.g. rk2) just fine. See bug-ode2.c.
Second, Newton-type iterations could be inefficient in general. See: Ernst Hairer, Christian Lubich, Gerhard Wanner, Geometric Numerical Integration. Structure-Preserving Algorithms for Ordinary Differential Equations. Springer Series in Comput. Mathematics, Vol. 31, Springer-Verlag 2002. pp. 330-332.
(file #21314)