2025/08/09 by Francesco Dell’Accio, Francesco Larosa, Dell'Accio, Francesco +5 · 2 citations
Mathematics · Engineering · Computer Science · #Mathematical Approximation and Integration #Advanced Numerical Analysis Techniques #Polynomial and algebraic computation
paper · pdf · doi:10.48550/arxiv.2508.07070
Histopolation, or interpolation on segments, is a mathematical technique used to approximate a function f over a given interval I=[a,b] by exploiting integral information over a set of subintervals of I. Unlike classical polynomial interpolation, which is based on pointwise function evaluations, histopolation reconstructs a function using integral data. However, similar to classical polynomial interpolation, histopolation suffers from the well-known Runge phenomenon when integral data are based on a grid with many equispaced nodes, as well as the Gibbs phenomenon when approximating discontinuous functions. In contrast, quasi-histopolation is designed to relax the strict requirement of passing through all the given data points. This inherent flexibility can reduce the likelihood of oscillatory behavior using, for example, rational approximation operators. In this work, we introduce a C∞ rational quasi-histopolation operator, for bounded (integrable) functions, which reconstruct a function by defeating both the Runge and Gibbs phenomena. A key element of our approach is to blend local histopolation polynomials on a few nodes using multinode Shepard functions as blending functions. Several numerical experiments demonstrate the accuracy of our method.