Fast Least-Squares Padé approximation of problems with normal operators and meromorphic structure

In this work, we consider the approximation of Hilbert space-valued meromorphic functions that arise as solution maps of parametric PDEs whose operator is the shift of an operator with normal and compact resolvent, e.g. the Helmholtz equation. In this restrictive setting, we propose a simplified version of the Least Squares Padé approximation technique introduced in [Bonizzoni eal, M2AN 2017]. In particular, the estimation of the poles of the target function reduces to a low-dimensional eigenproblem for a Gramian matrix, allowing for a robust and efficient numerical implementation (hence the "fast" in the name). Moreover, we prove several theoretical results that improve and extend those in [Bonizzoni eal, M2AN 2017], including the exponential decay of the error in the approximation of the poles, and the convergence in measure of the approximant to the target function. The latter result extends the classical one for scalar Padé approximation to our functional framework. We provide numerical results that confirm the improved accuracy of the proposed method with respect to the one introduced in [Bonizzoni eal, M2AN 2017] for differential operators with normal and compact resolvent.

Published in:
Other identifiers:

 Record created 2018-09-13, last modified 2019-06-19

Rate this document:

Rate this document:
(Not yet reviewed)