{"title":"Stability of convergence rates: kernel interpolation on non-Lipschitz domains","authors":"Tizian Wenzel, Gabriele Santin, Bernard Haasdonk","doi":"10.1093/imanum/drae014","DOIUrl":null,"url":null,"abstract":"Error estimates for kernel interpolation in Reproducing Kernel Hilbert Spaces usually assume quite restrictive properties on the shape of the domain, especially in the case of infinitely smooth kernels like the popular Gaussian kernel. In this paper we prove that it is possible to obtain convergence results (in the number of interpolation points) for kernel interpolation for arbitrary domains $\\varOmega \\subset{\\mathbb{R}} ^{d}$, thus allowing for non-Lipschitz domains including e.g., cusps and irregular boundaries. Especially we show that, when going to a smaller domain $\\tilde{\\varOmega } \\subset \\varOmega \\subset{\\mathbb{R}} ^{d}$, the convergence rate does not deteriorate—i.e., the convergence rates are stable with respect to going to a subset. We obtain this by leveraging an analysis of greedy kernel algorithms. The impact of this result is explained on the examples of kernels of finite as well as infinite smoothness. A comparison to approximation in Sobolev spaces is drawn, where the shape of the domain $\\varOmega $ has an impact on the approximation properties. Numerical experiments illustrate and confirm the analysis.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.3000,"publicationDate":"2024-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"IMA Journal of Numerical Analysis","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1093/imanum/drae014","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
Error estimates for kernel interpolation in Reproducing Kernel Hilbert Spaces usually assume quite restrictive properties on the shape of the domain, especially in the case of infinitely smooth kernels like the popular Gaussian kernel. In this paper we prove that it is possible to obtain convergence results (in the number of interpolation points) for kernel interpolation for arbitrary domains $\varOmega \subset{\mathbb{R}} ^{d}$, thus allowing for non-Lipschitz domains including e.g., cusps and irregular boundaries. Especially we show that, when going to a smaller domain $\tilde{\varOmega } \subset \varOmega \subset{\mathbb{R}} ^{d}$, the convergence rate does not deteriorate—i.e., the convergence rates are stable with respect to going to a subset. We obtain this by leveraging an analysis of greedy kernel algorithms. The impact of this result is explained on the examples of kernels of finite as well as infinite smoothness. A comparison to approximation in Sobolev spaces is drawn, where the shape of the domain $\varOmega $ has an impact on the approximation properties. Numerical experiments illustrate and confirm the analysis.
期刊介绍:
The IMA Journal of Numerical Analysis (IMAJNA) publishes original contributions to all fields of numerical analysis; articles will be accepted which treat the theory, development or use of practical algorithms and interactions between these aspects. Occasional survey articles are also published.