Compact Rational Krylov for Parametrized Systems with Application to BEM Frequency Sweeping
For computational scientists solving parametrized systems (e.g., in acoustics), this method reduces the cost of multi-query simulations by reusing information across parameter values.
The paper adapts the compact rational Krylov (CORK) framework to efficiently solve parametrized linear systems for many parameter values at once, demonstrating speedups in frequency sweeping for Helmholtz BEM problems.
In parametrized linear systems $\mathsf{P}(μ)\mathsf{x}=\mathsf{b}$ the system matrix $\mathsf{P}$ depends nonlinearly on a parameter $μ$ and solutions are sought for many values of this parameter. We show that the compact rational Krylov (CORK) framework, originally introduced to solve nonlinear eigenvalue problems, can be used to efficiently produce approximate solutions to such a system for many values of the parameter at once. In this approach the parametrized system is first linearized, resulting in a large shifted linear system $(\boldsymbol{\mathsf{A}}-μ\boldsymbol{\mathsf{B}})\boldsymbol{\mathsf{y}}=\boldsymbol{\mathsf{d}}$. We formulate a left- and right-preconditioned rational Krylov GMRES method for shifted linear systems and show how the CORK framework can be used to speed up these methods in the setting of parametrized linear systems. Additionally, we show how to incorporate a right-hand side $\mathsf{b}(μ)$ that also depends on the parameter, how to choose the shifts to steer convergence and how to allow for inexact solves at these shifts throughout the iterations. As an application we consider the 'frequency sweeping' of Helmholtz scattering problems through the Boundary Element Method (BEM), enabled via an efficient representation of the dense but data-sparse wavenumber-dependent system matrix.