Another approach to Runge-Kutta methods

The condition equations are derived by the introduction of a system of equivalent differential equations, avoiding the usual formalism with trees and elementary differentials. Solutions to the condition equations are found by direct optimization, avoiding the necessity to introduce simplifying assumptions upon the Runge-Kutta coefficients. More favourable coefficients, in view of rounding errors, are found

BOKASUN: a fast and precise numerical program to calculate the Master Integrals of the two-loop sunrise diagrams

We present the program BOKASUN for fast and precise evaluation of the Master Integrals of the two-loop self-mass sunrise diagram for arbitrary values of the internal masses and the external four-momentum. We use a combination of two methods: a Bernoulli accelerated series expansion and a Runge-Kutta numerical solution of a system of linear differential equations

Spatial Manifestations of Order Reduction in Runge-Kutta Methods for Initial Boundary Value Problems

This paper studies the spatial manifestations of order reduction that occur when time-stepping initial-boundary-value problems (IBVPs) with high-order Runge-Kutta methods. For such IBVPs, geometric structures arise that do not have an analog in ODE IVPs: boundary layers appear, induced by a mismatch between the approximation error in the interior and at the boundaries. To understand those boundary layers, an analysis of the modes of the numerical scheme is conducted, which explains under which circumstances boundary layers persist over many time steps. Based on this, two remedies to order reduction are studied: first, a new condition on the Butcher tableau, called weak stage order, that is compatible with diagonally implicit Runge-Kutta schemes; and second, the impact of modified boundary conditions on the boundary layer theory is analyzed.Comment: 41 pages, 9 figure