Follow-up Comment #2, bug #30540 (project gsl):
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)
_______________________________________________________
Additional Item Attachment:
File name: bug-ode2.c Size:2 KB
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?30540>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
