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

09/14/2020
by   Manuel Tiglio, et al.
0

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.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 14

01/27/2021

Reduced Order and Surrogate Models for Gravitational Waves

We present an introduction to some of the state of the art in reduced or...
10/17/2021

Gravitational wave surrogates through automated machine learning

We analyze a prospect for predicting gravitational waveforms from compac...
08/03/2021

Arby - Fast data-driven surrogates

The availability of fast to evaluate and reliable predictive models is h...
01/10/2022

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

This work investigates the use of sparse polynomial interpolation as a m...
09/02/2020

Filtered interpolation for solving Prandtl's integro-differential equations

In order to solve Prandtl-type equations we propose a collocation-quadra...
02/03/2020

Semi-active H_∞ damping optimization by adaptive interpolation

In this work we consider the problem of semi-active damping optimization...
11/13/2018

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.