---
title: Unstructured MLS for Manifolds & MPM
url: https://www.emergentmind.com/topics/unstructured-moving-least-squares-umls-mpm
type: topic
---

# Unstructured MLS for Manifolds & MPM

Unstructured Moving Least Squares for Manifolds and Material Point Methods (UMLS-MPM) encompasses a set of algorithms applying moving least squares (MLS) techniques to function approximation on sampled manifolds and Lagrangian–Eulerian discretizations in computational mechanics. The unifying characteristic is their ability to yield smooth, high-accuracy approximations and stable field reconstructions on data defined over non-uniform and unstructured domains—without reliance on global embeddings or structured grids. UMLS-MPM approaches are of critical importance in geometric data analysis, machine learning on manifolds, and large-deformation computational mechanics, particularly in settings where the underlying domain is known only through noisy samples or unstructured tessellation.

## 1. Manifold MLS: Chart Reconstruction and Atlas Generation

Let $M\subset\mathbb R^n$ be an unknown $d$-dimensional $C^2$ manifold. Given a set of noisy samples $R=\{r_i = \tilde r_i + n_i\}_1^I$ near $M$ ($\|\mathbf n_i\|\leq\sigma_M$), the first stage constructs a local coordinate chart around any query $r$ by solving
\[
J(r; q, H) = \sum_{i=1}^I d(r_i, H)^2\, \theta\left(\|r_i - q\|/h\right)
\]
where $H = q + \text{Span}\{e_1, ..., e_d\}$ is an affine $d$-subspace, $q\in\mathbb R^n$, $\theta$ is a compactly supported, smooth weight, and $h$ is the fill distance of clean samples. The minimization is subject to
- $r - q \perp H$,
- $q$ within a ball $B_\mu(r)$, $\mu < \text{rch}(M)/2$,
- $B_{h+\sigma_M}(q)$ contains at least one sample.

Under mild hypotheses, there exists a unique, $C^\infty$-smooth mapping $r\mapsto (q(r), H(r))$. The projection property $q(\tilde r) = q(r)$, $H(\tilde r) = H(r)$ for $\tilde r - q(r)\perp H(r)$ provides chart consistency over neighborhoods. The collection $\{(U_j, \varphi_j)\}$ with $\varphi_j(r) = \text{Proj}_{H_j}(r)$ defines a $C^\infty$ atlas for a manifold $\widetilde M$ close to $M$ in Hausdorff distance [1711.00765].

## 2. High-Order Polynomial MLS on Local Charts

With the local chart $(q(r), H(r))$, each sample $(r_i, \psi_i)$ is paired as $x_i = \text{Proj}_{H(r)}(r_i) \in \mathbb R^d$ and $f_i = \psi_i$. The moving least squares approximant solves for $p_r\in\Pi_m^d$, the space of polynomials of degree $\leq m$ in $d$ variables:
\[
p_r = \arg\min_{p\in\Pi_m^d} \sum_{i=1}^I w\left(\frac{\|r_i - q(r)\|}{h}\right) [p(x_i) - f_i]^2
\]
where $w$ is a $C^\infty$ weight of compact support. The approximation at $r$ is $\widetilde\psi(r) = p_r(0)$. The Backus–Gilbert representation gives
\[
p_r(0) = \sum_i a_i(r) f_i
\]
where $a_i(r)$ are smooth in $r$ [1711.00765].

## 3. Approximation Order, Smoothness, and Error Guarantees

Let $h$ denote the fill distance for $\{\tilde r_i\}$. Under $C^{m+1}$ smoothness of the target and sufficient density, projected samples in each chart are $O(h)$-dense. The resulting MLS fit achieves
\[
\max_{r\in M} |\widetilde\psi(r)-\psi(r)| < C h^{m+1}
\]
in the noiseless case, recovering the classical $O(h^{m+1})$ order [1711.00765]. The full map $r \mapsto \widetilde\psi(r)$ is $C^\infty$-smooth, guaranteed by the $C^\infty$ chart construction and smooth, compactly supported weights.

## 4. UMLS-MPM Algorithms for Function and Tensor Field Approximation

Two algorithmic forms are canonical:

**Manifold function regression [1711.00765]:**
- \( \textbf{ChartReconstruct}(\{r_i\}, r)\to(q, U): \) Iterative least-squares fits return an affine subspace around $r$.
- \( \textbf{Approximate}(\{(r_i,\psi_i)\}, r): \) Projects samples, builds the polynomial MLS normal matrix, and evaluates at the origin to yield $\widetilde\psi(r)$.

**Material Point Method on Tessellations [2312.10338]:**
- For a point $x$ in an unstructured mesh, neighbors $\{x_i\}$ are assigned raw MLS weights $W(\|x-x_i\|/h)$, then scaled by a “diminishing” factor $\eta_i(x)$ based on mesh adjacency and barycentric coordinates.
- The sample weight is $w_i(x) = W(\|x-x_i\|/h)\,\eta_i(x)$, ensuring $C^1$ continuity of the resulting MLS shape functions $N_i(x)$ across cell faces.
- Shape functions are constructed as
\[
N_i(x) = p^T(0) M(x)^{-1} \left[ w_i(x) p(x-x_i) \right]
\]
where $M(x)$ aggregates the weighted second moments. The field and its gradient are then
\[
u_h(x) = \sum_i N_i(x) u_i \qquad \nabla u_h(x) = \sum_i u_i \nabla N_i(x)
\]
Steps for explicit MPM include particle–grid transfers, mass/momentum assembly, and velocity updates using MLS kernels [2312.10338].

## 5. Stability, Mesh Independence, and $C^1$ Kernel Construction

The central innovation of UMLS-MPM for material point simulation is the construction of weights via a product of a smooth, compactly supported kernel $W$ and a diminishing factor $\eta_i(x)$. For each mesh cell, $\eta_i(x)$ is given by
\[
\eta_i(x) = \sum_{n\in V^0(x)} \lambda_n(x) A_{i,n}
\]
where $V^0(x)$ is the simplex containing $x$, $\lambda_n(x)$ are its barycentric coordinates, and $A_{i,n}$ encodes mesh adjacency. This guarantees:
- $\eta_i(x)=1$ for all vertices in the current cell; 
- $\eta_i(x)\to0$ as $x$ approaches a boundary losing/gaining the node.

When these conditions hold, both MLS moment matrices and hence shape functions $N_i(x)$ are at least $C^1$, which removes the cell-crossing artifacts characteristic of piecewise-linear interpolation on unstructured meshes [2312.10338].

## 6. Computational Complexity and Implementation Properties

For manifold regression, the dominant cost is in local chart construction: linear least-squares steps cost $O(nd^2)$ in ambient dimension $n$, with constant iterations; chart-based local MLS in $\mathbb R^d$ requires $O(d^3)$ per fit and $O(I d)$ for evaluation. The overall complexity is $O(n d^2 + I d + \text{poly}(d))$, yielding linear scaling in $n$. No global eigendecomposition, kernel matrix, or dimension reduction is performed [1711.00765].

In MPM, MLS assembly over a first-ring vertex neighborhood is local per timestep. Explicit UMLS-MPM remains subject to a standard CFL stability constraint, $\Delta t\leq C\,h/\sqrt{E/\rho}$ [2312.10338]. Mesh quality (simplex conditioning, support radius $h$) crucially affects invertibility and convergence; extremely skewed elements can lead to poor conditioning.

## 7. Numerical Experiments and Performance

Empirical evaluations demonstrate:

- **Function regression**: On smooth manifolds (helix in $\mathbb R^3$, $S^2\to(\phi,\theta)$), the error decays as $O(h^{m+1})$. High-dimensional tests (image regression in $\mathbb R^{448\times416}$) achieve RMSE $\approx$ 0.0066 with $n=72$ samples and $m=1$. Against regression baselines (MALLER, NEDE), UMLS-MPM attains lower RMSE and 2–3$\times$ faster CPU time for Klein-bottle tasks in $\mathbb R^4$ [1711.00765].
- **MPM on unstructured meshes**: In all cases, UMLS-MPM eliminates cell-crossing errors, matches or exceeds accuracy of structured-grid B-spline MPM, and delivers $2$nd-order spatial convergence (e.g., in the 1D vibrating bar, 2D cantilever, and 3D slope failure tests). Complex geometries and contact boundaries (3D sphere expansion) are resolved without producing stress discontinuities or artifacts [2312.10338].

## Table: Algorithmic Core Components

| Domain                | Charting/Neighborhood      | MLS Core              | Continuity      |
|-----------------------|---------------------------|-----------------------|-----------------|
| Manifold regression   | Affine tangent chart      | Weighted poly MLS     | $C^\infty$      |
| MPM on meshes         | Vertex/element adjacency  | $C^1$ diminished MLS  | $C^1$           |

The provided empirical and theoretical guarantees establish UMLS-MPM as an intrinsic, local, mesh-agnostic, and smooth high-order regression/interpolation framework for both geometric learning and continuum mechanics without global embeddings or kernel matrices [1711.00765], [2312.10338].

Source: https://www.emergentmind.com/topics/unstructured-moving-least-squares-umls-mpm