W-methods with automatic partitioning by Krylov techniques for large stiff systems

Solving the stage equations in W-methods approximately by a Krylov process may be interpreted as an automatic partitioning method, where some of the components are integrated by an implicit scheme whereas others are treated by an explicit Runge–Kutta method. The authors consider an implementation which uses only one family of Krylov spaces for all stages, introducing an additional error of the size of the discretization error. Two main results for the linear autonomous case show that (i) the method stays at the asymptotic limit solution under mild restrictions, which may be enforced in the numerical computation, and (ii) the dimensions of the Krylov spaces need only be slightly larger than the number of “fast” solution components with nonnegligible contribution to the local solution. Numerical examples of dimensions 20 to 302 indicate that the method may perform well even for Krylov spaces with low dimension.