Multivariate Data-dependent Partition of Unity based on Moving Least Squares method

Data approximation is essential in fields such as geometric design, numerical PDEs, and curve modeling. Moving Least Squares (MLS) is a widely used method for data fitting; however, its accuracy degrades in the presence of discontinuities, often resulting in spurious oscillations similar to those associated with the Gibbs phenomenon. This work extends the integration of MLS with the Weighted Essentially Non-Oscillatory (WENO) method and with an innovative partition of unity approach to higher dimensions. We propose a data-dependent operator using the novel Non-Linear Partition of Unity based on Moving Least Squares method in $\mathbb{R}^n$, which improves accuracy near discontinuities and maintains high-order accuracy in smooth regions. We demonstrate some theoretical properties of the method and perform numerical experiments to validate its effectiveness.