Hard Thresholding Hyperinterpolation Over General Regions

• Published in: Journal of scientific computing 2025-02, Vol.102 (2), p.37, Article 37
• Main Authors: An, Congpei, Ran, Jiashu
• Format: Article
• Language: English
• Subjects: Algebra; Algorithms; Approximation; Computational Mathematics and Numerical Analysis; Error analysis; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematical functions; Mathematics; Mathematics and Statistics; Noise reduction; Norms; Operators (mathematics); Quadratures; Regularization; Theoretical; Upper bounds
• ISSN: 0885-7474; 1573-7691
• DOI: 10.1007/s10915-024-02754-4

This paper proposes a novel variant of hyperinterpolation, called hard thresholding hyperinterpolation. This approximation scheme of degree n leverages a hard thresholding operator to filter all hyperinterpolation coefficients, which approximate the Fourier coefficients of a continuous function through a quadrature rule with algebraic exactness 2 n . We prove that hard thresholding hyperinterpolation is the unique solution to an ℓ 0 -regularized weighted discrete least squares approximation problem. Hard thresholding hyperinterpolation is not only idempotent and commutative with hyperinterpolation, but also adheres to the Pythagorean theorem in terms of the discrete (semi) inner product. By the estimate of the reciprocal of Christoffel function, we present the upper bound of the uniform norm of hard thresholding hyperinterpolation operator. Additionally, hard thresholding hyperinterpolation possesses denoising and basis selection abilities akin to Lasso hyperinterpolation. To judge the L 2 errors of both hard thresholding and Lasso hyperinterpolations, we propose a criterion that integrates the regularization parameter with the product of noise coefficients and the signs of hyperinterpolation coefficients. Numerical examples on the sphere, the spherical triangle and the cube demonstrate the denoising ability of hard thresholding hyperinterpolation.