On the stability and accuracy of the Empirical Interpolation Method and Gravitational Wave Surrogates

by   Manuel Tiglio, et al.

The combination of the Reduced Basis (RB) and the Empirical Interpolation Method (EIM) approaches have produced outstanding results in gravitational wave (GW) science and in many other disciplines. In GW science in particular, these results range from building non-intrusive surrogate models for gravitational waves to fast parameter estimation adding the use of Reduced Order Quadratures. These surrogates have the salient feature of being essentially indistinguishable from or very close to supercomputer simulations of the Einstein equations but can be evaluated in less than a second on a laptop. In this note we analyze in detail how the EIM at each iteration attempts to choose the interpolation nodes so as to make the related Vandermonde-type matrix "as invertible as possible" as well as attempting to optimize the conditioning of its inversion and minimizing a Lebesgue-type constant for accuracy. We then compare through numerical experiments the EIM performance with fully optimized nested variations. We also discuss global optimal solutions through Fekete nodes. We find that in these experiments the EIM is actually close to optimal solutions but can be improved with small variations and relative low computational cost.



There are no comments yet.


page 14


Reduced Order and Surrogate Models for Gravitational Waves

We present an introduction to some of the state of the art in reduced or...

Gravitational wave surrogates through automated machine learning

We analyze a prospect for predicting gravitational waveforms from compac...

Arby - Fast data-driven surrogates

The availability of fast to evaluate and reliable predictive models is h...

Model Reduction Using Sparse Polynomial Interpolation for the Incompressible Navier-Stokes Equations

This work investigates the use of sparse polynomial interpolation as a m...

Filtered interpolation for solving Prandtl's integro-differential equations

In order to solve Prandtl-type equations we propose a collocation-quadra...

Semi-active H_∞ damping optimization by adaptive interpolation

In this work we consider the problem of semi-active damping optimization...

ROMAN: Reduced-Order Modeling with Artificial Neurons

Gravitational-wave data analysis is rapidly absorbing techniques from de...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.