Polynomial interpolation--regression on the sphere

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Dell'Accio, Francesco, Nudo, Federico, Pérez, Teresa E., Piñar, Miguel A.
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909039608201216
author Dell'Accio, Francesco
Nudo, Federico
Pérez, Teresa E.
Piñar, Miguel A.
author_facet Dell'Accio, Francesco
Nudo, Federico
Pérez, Teresa E.
Piñar, Miguel A.
contents We introduce an interpolation--regression operator for polynomial approximation on the unit sphere $\mathbb{S}^2$ from discrete samples. The approximant is a spherical polynomial of degree $r$ which interpolates the data on a prescribed subset of nodes and uses the remaining sampling nodes to minimize the residual in a least squares sense. Under natural rank assumptions on the associated Vandermonde matrices, the approximant is unique and is characterized by an orthogonality condition with respect to the discrete inner product on the sampling set. We then focus on the case in which the sampling and interpolation nodes are antipodally symmetric. In this setting, when the polynomial is expressed in real spherical harmonics, the constrained problem can be decomposed into independent even and odd components. In the same framework, we prove equivariance under the antipodal map and, more generally, under orthogonal transformations preserving the node sets. We also consider spherical designs. In this case, the normal matrix is a scalar matrix. Consequently, the spectral condition number of the associated KKT matrix can be written explicitly. Numerical experiments in both antipodal and non-antipodal settings illustrate the effectiveness of the proposed method.
format Preprint
id arxiv_https___arxiv_org_abs_2605_13124
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Polynomial interpolation--regression on the sphere
Dell'Accio, Francesco
Nudo, Federico
Pérez, Teresa E.
Piñar, Miguel A.
Numerical Analysis
We introduce an interpolation--regression operator for polynomial approximation on the unit sphere $\mathbb{S}^2$ from discrete samples. The approximant is a spherical polynomial of degree $r$ which interpolates the data on a prescribed subset of nodes and uses the remaining sampling nodes to minimize the residual in a least squares sense. Under natural rank assumptions on the associated Vandermonde matrices, the approximant is unique and is characterized by an orthogonality condition with respect to the discrete inner product on the sampling set. We then focus on the case in which the sampling and interpolation nodes are antipodally symmetric. In this setting, when the polynomial is expressed in real spherical harmonics, the constrained problem can be decomposed into independent even and odd components. In the same framework, we prove equivariance under the antipodal map and, more generally, under orthogonal transformations preserving the node sets. We also consider spherical designs. In this case, the normal matrix is a scalar matrix. Consequently, the spectral condition number of the associated KKT matrix can be written explicitly. Numerical experiments in both antipodal and non-antipodal settings illustrate the effectiveness of the proposed method.
title Polynomial interpolation--regression on the sphere
topic Numerical Analysis
url https://arxiv.org/abs/2605.13124