initial value problems;
explicit Runge-Kutta methods;
numerical geometric integration;
preservation of invariants;
variable step-size codes;
D O I:
10.1137/04061979X
中图分类号:
O29 [应用数学];
学科分类号:
070104 ;
摘要:
A new strategy to preserve invariants in the numerical integration of initial value problems with explicit Runge-Kutta methods is presented. It is proved that this technique retains the order of the original method, has an easy and cheap implementation, and can be used in adaptive Runge-Kutta codes. Some numerical experiments with the classical code of Dormand and Prince, DoPri5(4), based on a pair of embedded methods with orders 5 and 4, are presented to show the behavior of the new method for several problems which possess invariants.
机构:
Imperial Coll, Dep of Mathematics,, London, Engl, Imperial Coll, Dep of Mathematics, London, EnglImperial Coll, Dep of Mathematics,, London, Engl, Imperial Coll, Dep of Mathematics, London, Engl