Abstract:We consider commutator-free exponential integrators as put forward in [Alverman, A., Fehske, H.: High-order commutator-free exponential time-propagation of driven quantum systems. J. Comput. Phys. 230, 5930-5956 (2011)]. For parabolic problems, it is important for the well-definedness that such an integrator satisfies a positivity condition such that essentially it only proceeds forward in time. We prove that this requirement implies maximal convergence order of four for real coefficients, which has been conje… Show more
“…with coefficients a j,k given in Table 3. For these coefficients the positivity condition is not satisfied, in agreement with the fact that this cannot be the case for real coefficients, see [10]. For some applications, however, it is essential for stability reasons that this condition is satisfied.…”
Section: Th Order Commutator-free Magnus-type Integratorsmentioning
This paper discusses an efficient implementation of the generation of order conditions for the construction of exponential integrators like exponential splitting and Magnus-type methods in the computer algebra system Maple. At the core of this implementation is a new algorithm for the computation of coefficients of words in the formal expansion of the local error of the integrator. The underlying theoretical background including an analysis of the structure of the local error is briefly reviewed. As an application the coefficients of all 8th order self-adjoint commutator-free Magnus-type integrators involving the minimum number of 8 exponentials are computed.Keywords: Splitting methods • Magnus-type integrators • Local error • Order conditions • Computer algebra. 4 For our considerations, it is sufficient to discuss linear problems. The algebraic structure underlying method construction is the same for nonlinear problems due to the calculus of Lie derivatives [8, Section III.5.1]. 5 We have used Maple 18, Maple is a trademark of Waterloo Maple Inc.
“…with coefficients a j,k given in Table 3. For these coefficients the positivity condition is not satisfied, in agreement with the fact that this cannot be the case for real coefficients, see [10]. For some applications, however, it is essential for stability reasons that this condition is satisfied.…”
Section: Th Order Commutator-free Magnus-type Integratorsmentioning
This paper discusses an efficient implementation of the generation of order conditions for the construction of exponential integrators like exponential splitting and Magnus-type methods in the computer algebra system Maple. At the core of this implementation is a new algorithm for the computation of coefficients of words in the formal expansion of the local error of the integrator. The underlying theoretical background including an analysis of the structure of the local error is briefly reviewed. As an application the coefficients of all 8th order self-adjoint commutator-free Magnus-type integrators involving the minimum number of 8 exponentials are computed.Keywords: Splitting methods • Magnus-type integrators • Local error • Order conditions • Computer algebra. 4 For our considerations, it is sufficient to discuss linear problems. The algebraic structure underlying method construction is the same for nonlinear problems due to the calculus of Lie derivatives [8, Section III.5.1]. 5 We have used Maple 18, Maple is a trademark of Waterloo Maple Inc.
“…A proof of Theorem 2 was proposed in [7]. In Section 2 we will give a new independent proof by showing that Theorem 3 (and thus also Theorem 2) is an easy consequence of a recent result proved by the authors in [8].…”
Section: Theoremmentioning
confidence: 94%
“…2 Note that we have changed some denotations: H(t), H 0 , H 1 , s, a j , c j correspond respectively to the denotations A(t), A 0 , A 1 , J, b j , y j of [8].…”
Section: Proof Of Theoremmentioning
confidence: 99%
“…The essential step leading to the main result of [8] is comprised by the following proposition. 2 Proposition 1.…”
Section: Proof Of Theoremmentioning
confidence: 99%
“…Remark 2. Strictly speaking, only a version of Proposition 1 with the weaker conclusion that at least one of the coefficients a j is non-positive has been proved in [8]. Since we may assume form the outset that a j = 0 for j = 2, .…”
We prove that generalized exponential splitting methods making explicit use of commutators of the vector fields are limited to order four when only real coefficients are admitted. This generalizes the restriction to order two for classical splitting methods with only positive coefficients.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.