Krylov subspace residual and restarting for certain second order differential equations

2022/06/14 by M. A. Botchev, Botchev, M. A., Leonid Knizhnerman +3 · 1 citation
Mathematics · Engineering · Computer Science · #Numerical methods for differential equations #Electromagnetic Simulation and Numerical Methods #Digital Filter Design and Implementation

paper · pdf · doi:10.48550/arxiv.2206.06909

Abstract

We propose algorithms for efficient time integration of large systems of oscillatory second order ordinary differential equations (ODEs) whose solution can be expressed in terms of trigonometric matrix functions. Our algorithms are based on a residual notion for second order ODEs, which allows to extend the ``residual-time restarting'' Krylov subspace framework -- which was recently introduced for exponential and φ-functions occurring in time integration of first order ODEs -- to our setting. We then show that the computational cost can be further reduced in many cases by using our restarting in the Gautschi cosine scheme. We analyze residual convergence in terms of Faber and Chebyshev series and supplement these theoretical results by numerical experiments illustrating the efficiency of the proposed methods.

Cited by

• A general framework for Krylov ODE residuals with applications to randomized Krylov methods

Related

• Block Krylov subspace exact time integration of linear ODE systems. Part 1: algorithm description
• A residual concept for Krylov subspace evaluation of the φ matrix function
• A general framework for Krylov ODE residuals with applications to randomized Krylov methods
• A flexible short recurrence Krylov subspace method for matrices arising in the time integration of port Hamiltonian systems and ODEs/DAEs with a dissipative Hamiltonian
• On convergence of waveform relaxation for nonlinear systems of ordinary differential equations