Stable interpolation with isotropic and anisotropic Gaussians using Hermite generating function

Abstract

Gaussian kernels can be an efficient and accurate tool for multivariate interpolation. In practice, high accuracies are often achieved in the flat limit where the interpolation matrix becomes increasingly ill-conditioned. Stable evaluation algorithms for isotropic Gaussians (Gaussian radial basis functions) have been proposed based on a Chebyshev expansion of the Gaussians by Fornberg, Larsson & Flyer and based on a Mercer expansion with Hermite polynomials by Fasshauer & McCourt. In this paper, we propose a new stabilization algorithm for the multivariate interpolation with isotropic or anisotropic Gaussians derived from the generating function of the Hermite polynomials. We also derive and analyse a new analytic cut-off criterion for the generating function expansion that allows to automatically adjust the number of the stabilizing basis functions.

References in corpus (3)

• Fully kinetic versus reduced-kinetic modelling of collisionless plasma turbulence
• Comparison of time splitting and backsubstitution methods for integrating Vlasov's equation with magnetic fields
• Verification of nonlinear particle simulation of radio frequency waves in tokamak

Cited by in corpus (11)

• A quantum-inspired method for solving the Vlasov-Poisson equations
• A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations
• Investigating the use of field solvers for simulating classical systems
• An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events
• Importance of accurate consideration of the electron inertia in hybrid-kinetic simulations of collisionless plasma turbulence: 1. The 2D limit
• High frequency non-gyrokinetic turbulence at tokamak edge parameter
• A mass-conserving sparse grid combination technique with biorthogonal hierarchical basis functions for kinetic simulations
• Convergence of splitting methods on rotating grids for the magnetized Vlasov equation
• Neural semi-Lagrangian method for high-dimensional advection-diffusion problems
• Local cubic spline interpolation for Vlasov-type equations on a multi-patch geometry
• An efficient stochastic particle method for high-dimensional nonlinear PDEs