[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] [bug #30540] please, add convergence checks in rk4imp/rk2imp

From: Sergey B Kirpichev
Subject: [Bug-gsl] [bug #30540] please, add convergence checks in rk4imp/rk2imp
Date: Thu, 26 Aug 2010 12:05:28 +0000
User-agent: Mozilla/5.0 (X11; U; Linux x86_64; en; rv: Gecko/20080528 Epiphany/2.22

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:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]