Greedy Adaptive Local Recovery of Functions in Sobolev Spaces

Fuente: arXiv
Saved in:
Bibliographic Details
Main Author: Schaback, Robert
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909272562991104
author Schaback, Robert
author_facet Schaback, Robert
contents There are many ways to upsample functions from multivariate scattered data locally, using only a few neighbouring data points of the evaluation point. The position and number of the actually used data points is not trivial, and many cases like Moving Least Squares require point selections that guarantee local recovery of polynomials up to a specified order. This paper suggests a kernel-based greedy local algorithm for point selection that has no such constraints. It realizes the optimal $L_\infty$ convergence rates in Sobolev spaces using the minimal number of points necessary for that purpose. On the downside, it does not care for smoothness, relying on fast $L_\infty$ convergence to a smooth function. The algorithm ignores near-duplicate points automatically and works for quite irregularly distributed point sets by proper selection of points. Its computational complexity is constant for each evaluation point, being dependent only on the Sobolev space parameters. Various numerical examples are provided. As a byproduct, it turns out that the well-known instability of global kernel-based interpolation in the standard basis of kernel translates arises already locally, independent of global kernel matrices and small separation distances.
format Preprint
id arxiv_https___arxiv_org_abs_2407_19864
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Greedy Adaptive Local Recovery of Functions in Sobolev Spaces
Schaback, Robert
Numerical Analysis
65D12, 65D05, 41A05, 65D25, 65D40
There are many ways to upsample functions from multivariate scattered data locally, using only a few neighbouring data points of the evaluation point. The position and number of the actually used data points is not trivial, and many cases like Moving Least Squares require point selections that guarantee local recovery of polynomials up to a specified order. This paper suggests a kernel-based greedy local algorithm for point selection that has no such constraints. It realizes the optimal $L_\infty$ convergence rates in Sobolev spaces using the minimal number of points necessary for that purpose. On the downside, it does not care for smoothness, relying on fast $L_\infty$ convergence to a smooth function. The algorithm ignores near-duplicate points automatically and works for quite irregularly distributed point sets by proper selection of points. Its computational complexity is constant for each evaluation point, being dependent only on the Sobolev space parameters. Various numerical examples are provided. As a byproduct, it turns out that the well-known instability of global kernel-based interpolation in the standard basis of kernel translates arises already locally, independent of global kernel matrices and small separation distances.
title Greedy Adaptive Local Recovery of Functions in Sobolev Spaces
topic Numerical Analysis
65D12, 65D05, 41A05, 65D25, 65D40
url https://arxiv.org/abs/2407.19864