A block Krylov subspace time-exact solution method for linear ODE systems

We propose a time-exact Krylov-subspace-based method for solving linear ODE (ordinary differential equation) systems of the form $y'=-Ay + g(t)$ and $y''=-Ay + g(t)$, where $y(t)$ is the unknown function. The method consists of two stages. The first stage is an accurate piecewise polynomial approximation of the source term $g(t)$, constructed with the help of the truncated SVD (singular value decomposition). The second stage is a special residual-based block Krylov subspace method. The accuracy of the method is only restricted by the accuracy of the piecewise polynomial approximation and by the error of the block Krylov process. Since both errors can, in principle, be made arbitrarily small, this yields, at some costs, a time-exact method. Numerical experiments are presented to demonstrate efficiency of the new method, as compared to an exponential time integrator with Krylov subspace matrix function evaluations

GMRES for oscillatory matrix-valued differential equations

We investigate the use of Krylov subspace methods to solve linear, oscillatory ODEs. When we apply a Krylov subspace method to a properly formulated equation, we retain the asymptotic accuracy of the asymptotic expansion whilst converging to the exact solution. We will demonstrate the effectiveness of this method by computing Error and Mathieu functions

A Class of Preconditioners for Large Indefinite Linear Systems, as by-product of Krylov subspace Methods: Part II

In this paper we consider the parameter dependent class of preconditioners M(a,d,D) defined in the companion paper The latter was constructed by using information from a Krylov subspace method, adopted to solve the large symmetric linear system Ax = b. We first estimate the condition number of the preconditioned matrix M(a,d,D). Then our preconditioners, which are independent of the choice of the Krylov subspace method adopted, proved to be effective also when solving sequences of slowly changing linear systems, in unconstrained optimization and linear algebra frameworks. A numerical experience is provided to give evidence of the performance of M(a,d,D).preconditioners; large indefinite linear systems; large scale nonconvex optimization; Krylov subspace methods