We discuss the structure of the local error of exponential operator splitting methods. In particular, it is shown that the leading error term is a Lie element, i.e., a linear combination of higher-degree commutators of the given operators. This structural assertion can be used to formulate a simple algorithm for the automatic generation of a minimal set of polynomial equations representing the order conditions, for the general case as well as in symmetric settings.
The concepts of stability regions, A- and A(α)-stability - albeit based on scalar models - turned out to be essential for the identification of implicit methods suitable for the integration of stiff ODEs. However, for multistep methods, knowledge of the stability region provides no information on the quantitative stability behavior of the scheme. In this paper we fill this gap for the important class of Backward Differentiation Formulas (BDF). Quantitative stability bounds are derived which are uniformly valid in the stability region of the method. Our analysis is based on a study of the separation of the characteristic roots and a special similarity decomposition of the associated companion matrix.
JavaScript jest wyłączony w Twojej przeglądarce internetowej. Włącz go, a następnie odśwież stronę, aby móc w pełni z niej korzystać.