A Method for Dimensionally Adaptive Sparse Trigonometric Interpolation of Periodic Functions

08/28/2019
by   Zack Morrow, et al.
0

We present a method for dimensionally adaptive sparse trigonometric interpolation of multidimensional periodic functions belonging to a smoothness class of finite order. This method is useful for problems where periodicity must be preserved and the precise anisotropy is not known a priori. The motivating application behind this work is the adaptive approximation of a multi-input model for a molecular potential energy surface (PES) where each input represents an angle of rotation. Our method is based on an anisotropic quasi-optimal estimate for the decay rate of the Fourier coefficients of the model; a least-squares fit to the coefficients of the interpolant is used to estimate the anisotropy. Thus, our adaptive approximation strategy begins with a coarse isotropic interpolant, which is gradually refined using the estimated anisotropic rates. The procedure takes several iterations where ever-more accurate interpolants are used to generate ever-improving anisotropy rates. We present several numerical examples of our algorithm where the adaptives procedure successfully recovers the theoretical "best" convergence rate, including an application to a periodic PES approximation. An open-source implementation of our algorithm resides in the Tasmanian UQ library developed at Oak Ridge National Laboratory.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
09/06/2019

An analog of the Sinc approximation for periodic functions

In this paper, we propose an interpolation formula for periodic function...
research
02/02/2020

Approximation by periodic multivariate quasi-projection operators

Approximation properties of periodic quasi-projection operators with mat...
research
11/26/2020

Exact Reconstruction of Sparse Non-Harmonic Signals from Fourier Coefficients

In this paper, we derive a new reconstruction method for real non-harmon...
research
08/05/2022

The Object Oriented c++ library QIBSH++ for Hermite spline Quasi Interpolation

The library QIBSH++ is a C++ object oriented library for the solution of...
research
10/20/2020

Gradient flow finite element discretisations with energy-based adaptivity for excited states of Schrödingers equation

We present an effective numerical procedure, which is based on the compu...
research
08/09/2022

Decay of coefficients and approximation rates in Gabor Gaussian frames

Gabor frames are a standard tool to decompose functions into a discrete ...

Please sign up or login with your details

Forgot password? Click here to reset