A computer program for determining truncation error coefficients for Runge-Kutta methods

The basic structure of a program to generate the truncation error coefficients for Runge-Kutta (RK) methods is reformulated to reduce storage requirements significantly and to accommodate variable dimensioning. This FORTRAN program, SUBROUTINE RKEQ, determines truncation error coefficients for RK algorithms for orders 1 through 10 and extends the order of coefficients through 12 with the 11th- and 12th-order terms determined following the patterns used to establish the lower order coefficients. Both subroutines (the original and RKEQ) are also written to treat RK m-fold methods which utilize m known derivatives of f to increase the order of the algorithm. Setting m = 0 gives the classical RK algorithm

Reducing the error of geoid undulation computations by modifying Stokes' function

The truncation theory as it pertains to the calculation of geoid undulations based on Stokes' integral, but from limited gravity data, is reexamined. Specifically, the improved procedures of Molodenskii et al. are shown through numerical investigations to yield substantially smaller errors than the conventional method that is often applied in practice. In this improved method, as well as in a simpler alternative to the conventional approach, the Stokes' kernel is suitably modified in order to accelerate the rate of convergence of the error series. These modified methods, however, effect a reduction in the error only if a set of low-degree potential harmonic coefficients is utilized in the computation. Consider, for example, the situation in which gravity anomalies are given in a cap of radius 10 deg and the GEM 9 (20,20) potential field is used. Then, typically, the error in the computed undulation (aside from the spherical approximation and errors in the gravity anomaly data) according to the conventional truncation theory is 1.09 m; with Meissl's modification it reduces to 0.41m, while Molodenskii's improved method gives 0.45 m. A further alteration of Molodenskii's method is developed and yields an RMS error of 0.33 m. These values reflect the effect of the truncation, as well as the errors in the GEM 9 harmonic coefficients. The considerable improvement, suggested by these results, of the modified methods over the conventional procedure is verified with actual gravity anomaly data in two oceanic regions, where the GEOS-3 altimeter geoid serves as the basis for comparison. The optimal method of truncation, investigated by Colombo, is extremely ill-conditioned. It is shown that with no corresponding regularization, this procedure is inapplicable

Linear-implicit strong schemes for Itô-Galkerin approximations of stochastic PDEs

Linear-implicit versions of strong Taylor numerical schemes for finite dimensional Itô stochastic differential equations (SDEs) are shown to have the same order as the original scheme. The combined truncation and global discretization error of an gamma strong linear-implicit Taylor scheme with time-step delta applied to the N dimensional Itô-Galerkin SDE for a class of parabolic stochastic partial differential equation (SPDE) with a strongly monotone linear operator with eigenvalues lambda 1 <= lambda 2 <= ... in its drift term is then estimated by K(lambda N -½ + 1 + delta gamma) where the constant K depends on the initial value, bounds on the other coefficients in the SPDE and the length of the time interval under consideration. AMS subject classifications: 35R60, 60H15, 65M15, 65U05

Theoretical and Phenomenological Constraints on Form Factors for Radiative and Semi-Leptonic B-Meson Decays

We study transition form factors for radiative and rare semi-leptonic B-meson decays into light pseudoscalar or vector mesons, combining theoretical constraints and phenomenological information from Lattice QCD, light-cone sum rules, and dispersive bounds. We pay particular attention to form factor parameterisations which are based on the so-called series expansion, and study the related systematic uncertainties on a quantitative level. In this context, we also provide the NLO corrections to the correlation function between two flavour-changing tensor currents, which enters the unitarity constraints for the coefficients in the series expansion.Comment: 52 pages; v2: normalization error in (29ff.) corrected, conclusion about relevance of unitarity bounds modified; form factor fits unaffected; references added; v3: discussion on truncation of series expansion added, matches version to be published in JHEP; v4: corrected typos in Tables 5 and