A New Optimality Property of Strang's Splitting

Abstract

For systems of the form q = M-1 p, p = -Aq+f(q), common in many applications, we analyze splitting integrators based on the (linear/nonlinear) split systems q = M-1 p, p = -Aq and q = 0, p = f(q). We show that the well-known Strang splitting is optimally stable in the sense that, when applied to a relevant model problem, it has a larger stability region than alternative integrators. This generalizes a well-known property of the common St\"ormer/Verlet/leapfrog algorithm, which of course arises from Strang splitting based on the (kinetic/potential) split systems q = M-1 p, p = 0 and q = 0, p = -Aq+f(q).

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…