Weighted Essentially Non-Oscillatory Shepard method

TL;DR

This paper introduces a nonlinear Shepard interpolation method inspired by WENO techniques, improving accuracy and stability near discontinuities in scattered data interpolation.

Contribution

It develops a novel nonlinear Shepard method that incorporates WENO-inspired adaptive weighting to better handle discontinuities and reduce oscillations.

Findings

Enhanced accuracy near discontinuities

Reduced oscillations in interpolations

Superior performance in complex datasets

Abstract

Shepard method is a fast algorithm that has been classically used to interpolate scattered data in several dimensions. This is an important and well-known technique in numerical analysis founded in the main idea that data that is far away from the approximation point should contribute less to the resulting approximation. Approximating piecewise smooth functions in $\mathbb{R}^n$ near discontinuities along a hypersurface in $\mathbb{R}^{n-1}$ is challenging for the Shepard method or any other linear technique for sparse data due to the inherent difficulty in accurately capturing sharp transitions and avoiding oscillations. This letter is devoted to constructing a non-linear Shepard method using the basic ideas that arise from the weighted essentially non-oscillatory interpolation method (WENO). The proposed method aims to enhance the accuracy and stability of the traditional Shepard method by incorporating WENO's adaptive and nonlinear weighting mechanism. To address this challenge, we will nonlinearly modify the weight function in a general Shepard method, considering any weight function, rather than relying solely on the inverse of the distance squared. This approach effectively reduces oscillations near discontinuities and improves the overall interpolation quality. Numerical experiments demonstrate the superior performance of the new method in handling complex datasets, making it a valuable tool for various applications in scientific computing and data analysis.