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

arrow_drop_up

Comments

Sign up or login to leave a comment