Symmetrized local error estimators for time-reversible one-step methods in nonlinear evolution equations

Abstract

Prior work on computable defect-based local error estimators for (linear) time-reversible integrators is extended to nonlinear and nonautonomous evolution equations. We prove that the asymptotic results from the linear case [W. Auzinger and O. Koch, An improved local error estimator for symmetric time-stepping schemes, Appl.Math.Lett. 82 (2018), pp. 106-110] remain valid, i.e., the modified estimators yield an improved asymptotic order as the step size goes to zero. Typically, the computational effort is only slightly higher than for conventional defect-based estimators, and it may even be lower in some cases. We illustrate this by some examples and present numerical results for evolution equations of Schrödinger type, solved by either time-splitting or Magnus-type integrators. Finally, we demonstrate that adaptive time-stepping schemes can be successfully based on our local error estimators.