Schur Decomposition for Stiff Differential Equations

Thoma Zoto, John C. Bowman

A quantitative definition of numerical stiffness for initial value problems is proposed. Exponential integrators can effectively integrate linearly stiff systems, but they become expensive when the linear coefficient is a matrix, especially when the time step is adapted to maintain a prescribed local error. Schur decomposition is shown to avoid the need for computing matrix exponentials in such simulations, while still circumventing linear stiffness.

Knowledge Graph



Sign up or login to leave a comment