Diagnosing stiffness for Runge-Kutta methods
SIAM Journal on Scientific and Statistical Computing
Numerical methods for ordinary differential systems: the initial value problem
Numerical methods for ordinary differential systems: the initial value problem
Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems
Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems
Calculation of pseudospectra by the Arnoldi iteration
SIAM Journal on Scientific Computing - Special issue on iterative methods in numerical linear algebra; selected papers from the Colorado conference
ACM Transactions on Mathematical Software (TOMS)
Hi-index | 0.01 |
A new stiffness detection scheme based on explicit Runge-Kutta methods is proposed. It uses a Krylov subspace approximation to estimate the eigenvalues of the Jacobian of the differential system. The numerical examples indicate that this technique is a worthwhile alternative to other known stiffness detection schemes, especially when the systems are large and when it is desirable to know more about the spectrum of the Jacobian than just the spectral radius.