Kernel interpolation on generalized sparse grids

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Griebel, Michael, Harbrecht, Helmut, Multerer, Michael
Format: Preprint
Veröffentlicht: 2025
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866911617458896896
author Griebel, Michael
Harbrecht, Helmut
Multerer, Michael
author_facet Griebel, Michael
Harbrecht, Helmut
Multerer, Michael
contents We consider scattered data approximation on product regions of equal and different dimensionality. On each of these regions, we assume quasi-uniform but unstructured data sites and construct optimal sparse grids for scattered data interpolation on the product region. For this, we derive new improved error estimates for the respective kernel interpolation error by invoking duality arguments. An efficient algorithm to solve the underlying linear system of equations is proposed. The algorithm is based on the sparse grid combination technique, where a sparse direct solver is used for the elementary anisotropic tensor product kernel interpolation problems. The application of the sparse direct solver is facilitated by applying a samplet matrix compression to each univariate kernel matrix, resulting in an essentially sparse representation of the latter. In this way, we obtain a method that is able to deal with large problems up to billions of interpolation points, especially in case of reproducing kernels of nonlocal nature. Numerical results are presented to qualify and quantify the approach.
format Preprint
id arxiv_https___arxiv_org_abs_2505_12282
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Kernel interpolation on generalized sparse grids
Griebel, Michael
Harbrecht, Helmut
Multerer, Michael
Numerical Analysis
41A46, 41A63, 46E35
We consider scattered data approximation on product regions of equal and different dimensionality. On each of these regions, we assume quasi-uniform but unstructured data sites and construct optimal sparse grids for scattered data interpolation on the product region. For this, we derive new improved error estimates for the respective kernel interpolation error by invoking duality arguments. An efficient algorithm to solve the underlying linear system of equations is proposed. The algorithm is based on the sparse grid combination technique, where a sparse direct solver is used for the elementary anisotropic tensor product kernel interpolation problems. The application of the sparse direct solver is facilitated by applying a samplet matrix compression to each univariate kernel matrix, resulting in an essentially sparse representation of the latter. In this way, we obtain a method that is able to deal with large problems up to billions of interpolation points, especially in case of reproducing kernels of nonlocal nature. Numerical results are presented to qualify and quantify the approach.
title Kernel interpolation on generalized sparse grids
topic Numerical Analysis
41A46, 41A63, 46E35
url https://arxiv.org/abs/2505.12282