---
title: Point-based Neural Operators
url: https://www.emergentmind.com/topics/point-based-neural-operators
type: topic
---

# Point-based Neural Operators

Searching arXiv for recent papers on point-based neural operators and closely related operator-learning methods on point clouds, arbitrary point sets, and geometric domains.
Point-based neural operators are neural operator models that operate directly on point clouds or coordinate samples of a domain rather than on structured grids or purely graph-discretized meshes. In the PDE setting, they are used to approximate solution operators on arbitrary geometries, often through a continuous formulation that is later discretized on sampled points. A common abstraction is the solution map
\[
\mathcal{G}^{\dagger} : (a,\Omega)\mapsto u,
\]
where \(a\) denotes parameter fields or forcing data on a domain \(\Omega\) and \(u\) is the corresponding solution field; in discrete engineering settings, the queried output may be written as
\[
u(x_q) = M(x_q, F, G, S, p),
\]
with \(x_q\) a query coordinate, \(F\) sampled forcing data, \(G\) sampled boundary or initial data, \(S\) a geometry point cloud, and \(p\) an optional parametric design vector [2501.14475] [2510.05995].

## 1. Conceptual scope and operator-learning setting

The defining characteristic of point-based neural operators is that geometry enters through sampled coordinates and pointwise features rather than through a fixed Cartesian lattice. This is central when the domain is complex, variable, adaptively sampled, or topologically changing. In the engineering benchmark literature, point-based neural operators are described as models that operate directly on point clouds or coordinate samples of the PDE domain and learn mappings from geometric and physical inputs to field outputs on arbitrary \(3\)D geometries [2510.05995].

A recurrent formulation begins from an abstract PDE \(R(u,a)=0\), with neural operator training based on paired samples \(\{(a_i,\Omega_i,u_i)\}_{i=1}^n\) and an empirical objective
\[
\theta^*=\arg\min_\theta \frac1n\sum_{i=1}^n c\!\left(G_\theta(a_i,\Omega_i),u_i\right).
\]
This perspective emphasizes that the learned object is a map between function spaces, not merely a predictor on one fixed discretization. The stated motivation is robustness to arbitrary discretizations, irregular or adaptive meshes, variable geometries, and topology changes [2501.14475].

This functional viewpoint also explains why point-based neural operators appear in several distinct subfields. Some works address surrogate modeling for parametric PDEs on point clouds, some learn differential operators directly on point sets for geometry processing, some extend spectral neural operators to arbitrary non-equispaced samples, and some use pointwise information to calibrate spectral or physics-informed models. The common denominator is that the operator acts on sampled fields and geometry without requiring a rigid grid.

## 2. Continuous formulations and discretization on point sets

A major design principle is to define the operator in a continuous setting and discretize only afterward. In the analysis literature, neural operators are written as maps \(\mathcal{G}_\theta: U\to V\) between Banach or Sobolev spaces, for example \(U=H^s(D)\) and \(V=H^t(D)\), with stability expressed as
\[
\|\mathcal{G}_\theta(u)-\mathcal{G}_\theta(v)\|_{H^t(D)} \le L \|u-v\|_{H^s(D)}.
\]
The same work makes its universality argument by discretizing the domain into points \(\{x_i\}_{i=1}^N\) and approximating \(\mathcal{T}(u)(x_i)\) from the sampled values \(u(x_1),\dots,u(x_N)\), which places sampled-point representations directly inside the theoretical operator-learning framework [2410.21481].

An early operator-theoretic prototype is Point Convolutional Neural Networks, where a point cloud function is lifted to a volumetric function, processed by an ordinary Euclidean operator, and then sampled back on the point cloud. The resulting point-cloud operator is written as
\[
O_X = {}_X \circ O \circ {}_X,
\]
with extension and restriction providing the two outer maps. In this construction, point convolution is the pull-back of Euclidean volumetric convolution through an extension-restriction mechanism, and the paper explicitly interprets the framework as a prototype for point-based neural operators [1803.10091].

A different continuous-to-discrete strategy appears in the Resolution Independent Neural Operator. Instead of feeding raw sampled values to a branch net, it learns a continuous dictionary \(\{\psi_l(x)\}_{l=1}^Q\) parameterized as implicit neural representations and projects arbitrarily sampled input functions onto coefficient vectors,
\[
u^{(i)}(x) \approx \sum_{l=1}^Q \alpha^{(i)}_l \psi_l(x),
\]
with coefficients obtained by regularized projection. DeepONet then acts on the coefficient vector \(\boldsymbol{\alpha}(u)\) rather than on a resolution-dependent sensor vector. This directly removes the requirement that all input functions be observed at identical locations, while preserving the operator-learning structure of branch and trunk networks [2407.13010].

These constructions suggest a useful unifying interpretation: point-based neural operators differ less in the abstract learning target than in how they realize lifting, aggregation, and projection for irregular samples. Some volumetrize the data, some infer continuous basis coefficients, and some discretize integral or differential operators directly on the point set.

## 3. Point-set and token-based architectural families

Within this common framework, several architectural families recur. In the engineering benchmark literature, a generic point-based operator is described as predicting the field at each point \(x_i\) through a stacked point-operator composition,
\[
u(x_i) = \mathrm{Attn}^l \circ \mathrm{Attn}^{l-1} \circ \cdots \circ \mathrm{Attn}^1(v_i),
\]
where \(v_i\) denotes point features and the attention rule depends on the model class [2510.05995].

The simplest baseline is PointNet-style processing: shared MLPs followed by symmetric pooling on unordered point sets. This is easy to deploy, but the benchmark states that it is weaker at capturing local structure. More expressive point-based neural operators replace global pooling with nonlocal attention or learned tokenization. In GNOT, the update is globally nonlocal,
\[
v^{l+1}_i = \sum_j \frac{\exp(Q_i^l \cdot K_j^l / \tau)}{S_i} V_j^l,
\qquad
S_i = \sum_j \exp(Q_i^l \cdot K_j^l / \tau),
\]
with learned \(Q_i^l\), \(K_j^l\), and \(V_j^l\). In Transolver, points are first pooled into “physics-aware slice tokens,”
\[
w_j^i=\mathrm{Softmax}\!\left(\mathcal{U}_{\text{slice}^i}(v_j)\right),\qquad
s^i = \frac{\sum_j w^i_j v_j}{\sum_j w^i_j},
\]
then processed globally, and finally projected back to points through
\[
u(x_j)=\sum_i w_j^i \cdot s^{i,l}.
\]
This tokenization mechanism is a distinctive point-based aggregation-interpolation pattern: points are pooled into latent physics tokens, processed globally, and reconstructed pointwise [2510.05995].

Architectural reuse has also been studied at the operator-block level. DPNO introduces a dual-path composition with a ResNet-like stream,
\[
U_{k+1}(x)=G_k(U_k)(x)+U_k(x),
\]
and a DenseNet-like stream,
\[
V_{k+1}(x)=G_k([V_0,V_1,\ldots,V_k])(x).
\]
For DeepONet, this dual-path mechanism is applied to the trunk net rather than the branch net, with the stated interpretation that it enriches the output-space basis family while keeping the input-function encoding unchanged [2507.12719].

## 4. Spectral and basis-based methods on arbitrary point distributions

A second major line of work adapts spectral neural operators to non-equispaced point sets. Direct Spectral Evaluation replaces FFT-based spectral layers with direct evaluation of a truncated spectral basis on arbitrary samples. For point coordinates \(P\in[0,1]^{N\times D}\), the transform is written as
\[
\mathbf{X}=\mathbf{V}\mathbf{x},
\]
where \(\mathbf{V}\) is a Vandermonde-like matrix of basis evaluations on the sampled points. The same idea is used with Fourier bases in Euclidean domains and spherical harmonics on spherical point clouds. Since only a small truncated set of modes is retained, the paper gives the complexity of DSE as \(O(mN)\), contrasting it with \(O(N\log N)\) FFTs and \(O(N^2)\) brute-force DFTs [2305.19663].

This arbitrary-point spectralization is presented as a drop-in replacement for the spectral transform blocks in FNO, UFNO, FFNO, and SFNO. The stated practical consequence is that Fourier-based neural operators become compatible with irregular grids, lattices, point clouds, and spherical domains without interpolation onto a regular mesh [2305.19663].

A more adaptive variant appears in the Holistic Physics Solver, whose core block is described in the details as a Point-Calibrated Spectral Neural Operator. Here the spectral basis is not fixed globally. Instead, each point predicts frequency preferences through a gate
\[
g=[g_1,\ldots,g_{N^k}] = \text{Softmax}\!\left(\text{MLP}_{\text{gate}^{N^k}}(x)\right),
\]
and the truncated Laplace-Beltrami transform is replaced by a point-calibrated transform
\[
\mathcal{T}_{\text{PC-LBT}}(x)=\left[x^T(g_1\odot \phi_1),\,\ldots,\,x^T(g_{N^k}\odot \phi_{N^k})\right]^T.
\]
If \(g_i\equiv 1\), the construction reduces to the ordinary Laplace-Beltrami spectral transform, so classical spectral neural operators appear as a special case. The intended effect is to combine spectral continuity priors with point-level flexibility, especially on irregular geometries and under limited training data [2410.11382].

Taken together, these works show that “spectral” and “point-based” are not opposed categories. In current operator-learning practice, spectral processing can be executed directly on point sets, or the spectral basis itself can be made point-conditioned.

## 5. Geometry-native differential operators and learned geometric structure

Another large branch of point-based neural operators focuses on learning differential operators or geometric quantities directly on point clouds. The Neural Laplacian Operator is exemplary: it constructs a symmetric KNN graph, uses a GNN to produce per-point features, decodes edge weights \(w_{ij}\) and masses \(M_{ii}\), and trains the resulting operator not by matching matrix entries but by matching its action on probe functions,
\[
\mathcal{L}_{laplacian} = \sum_{f\in\mathcal{F}} w_f \left\|M^{-1}Lf - M_{gt}^{-1}L_{gt}f\right\|_2^2.
\]
This functional imitation objective is introduced specifically to bypass the connectivity mismatch between the KNN graph and the mesh on which the cotangent Laplacian is defined. The method is reported to reduce MSE by roughly an order of magnitude relative to the best prior method on ShapeNet testing data and to support downstream geometry processing tasks such as heat diffusion, geodesic distance, Laplacian smoothing, spectral filtering, and ARAP deformation directly on point clouds [2409.06506].

DeltaConv takes a more intrinsic differential-geometric route. It separates features into scalar and vector streams and connects them through discrete gradient, divergence, curl, co-gradient, and Hodge-Laplacian operators. The vector stream stores tangential vector coefficients at each point, allowing direction-dependent processing without a global tangent coordinate system. The core claim is that this enables anisotropic filtering directly on point clouds while remaining coordinate-independent and robust to irregular sampling [2111.08799].

Geometric Neural Operators push this direction toward reusable geometry models. They apply operator layers on local point-cloud neighborhoods \(\mathcal{N}_\epsilon(x)\), discretize integral operators by message passing,
\[
\tilde{\mathcal{K}[v_t](x_j)} = \frac{1}{N}\sum_{x_k\in B_r(x_j)} k(x_j,x_k)\,v_t(x_k),
\]
and train pretrained models to estimate Monge-gauge height functions, normals, inverse first fundamental forms, second fundamental forms, and Gaussian curvature. These pretrained GNPs are then reused in geometry estimation, mean-curvature flow, and Laplace-Beltrami PDE solves on noisy point clouds of arbitrary topology [2503.04649].

For parametric PDEs on complex domains, Point Cloud Neural Operator provides perhaps the most explicit PDE-native point-cloud formulation. Its input augments field values with coordinates and a density term,
\[
\tilde{a}(x)=
\begin{bmatrix}
a(x)\\
x\\
\rho(x;\Omega)
\end{bmatrix},
\]
and each layer combines a local linear map, a Fourier-kernel integral operator, and a local differential aggregation based on least-squares gradients:
\[
f_{\rm out}(x) = \sigma\Bigl( W^l f_{\rm in}(x)+b +\sum_k \int_\Omega e^{2\pi i \frac{k}{L}\cdot(x-y)}\,W_k^v f_{\rm in}(y)\rho(y;\Omega)\,dy +W^g \widetilde{\nabla} f_{\rm in}(x) \Bigr).
\]
The model is stated to be permutation invariant, to have layer complexity linear in the number of points when connectivity is sparse, and to admit a universal approximation result at the continuous level [2501.14475].

## 6. Physics-informed training and symmetry-aware loss design

A distinct, but related, use of pointwise structure appears in physics-informed neural operators. In this literature, “point symmetry” refers not to point clouds but to Lie point symmetries acting on the independent and dependent variables of a PDE,
\[
v=\sum_{i=1}^p \xi^i(x,u)\frac{\partial}{\partial x^i}+\sum_{\alpha=1}^q \phi_\alpha(x,u)\frac{\partial}{\partial u^\alpha}.
\]
The baseline PINO objective is the residual loss
\[
\min_\theta \|\Delta[u_\theta]\|_2^2.
\]
Earlier symmetry-aware loss augmentation added prolonged point-symmetry actions on the residual, whereas the generalized-symmetry formulation replaces standard generators by their evolutionary representatives,
\[
v_Q=\sum_{\alpha=1}^q Q_\alpha[u]\frac{\partial}{\partial u^\alpha},\qquad
Q_\alpha=\phi_\alpha-\sum_{i=1}^p \xi^i \partial u^\alpha/\partial x^i,
\]
and trains with
\[
\min_\theta\; \|\Delta[u_{\theta}]\|^2_2 + \gamma \sum_i\left\|\operatorname{pr} [v_i]_{Q}\left[\Delta\right][u_{\theta}]\right\|^2_2.
\]
The underlying proposition is that a generalized vector field is a symmetry of a differential system if and only if its evolutionary representative is [2502.00373].

The motivation for this replacement is that standard point symmetry losses can be uninformative. For Burgers’ equation,
\[
\Delta_B = u_t + u u_x - \nu u_{xx},
\]
the translation generator \(v_1=\partial_x\) yields
\[
\operatorname{pr} v_1[\Delta_B]=0,
\]
while its evolutionary representative satisfies
\[
\operatorname{pr}[v_1]_Q[\Delta_B]=-D_x[\Delta_B].
\]
The former contributes no training signal; the latter contributes a derivative of the residual. Analogous behavior is given for Darcy flow, where finite-dimensional point symmetries can collapse to zero while evolutionary representatives produce derivative terms such as \(D_x[\Delta_D]\) and \(D_y[\Delta_D]\). The stated practical effect is improved data efficiency, accuracy during training, and zero-shot resolution transfer on Darcy flow, together with gains on \(1\)D Burgers’ equation [2502.00373].

This line of work is conceptually important because it reframes symmetry information as operator regularization on sampled fields. It also clarifies that point-based inductive biases in neural operators are not limited to geometry representation; they also include pointwise physics priors and derivative-aware losses.

## 7. Applications, empirical behavior, and limitations

Empirically, point-based neural operators are especially prominent in settings where geometry is irregular, freeform, or naturally represented as a point cloud. The engineering benchmark on six industry-scale datasets concludes that point-based operators are the most natural family for freeform \(3\)D geometry and irregular discretizations, with transformer-based variants such as GNOT and especially Transolver substantially stronger than plain PointNet. The same benchmark also states that point-based methods are not universally best: branch-trunk models are usually better on Heat sink and Bracket, grid-based models outperform point-based models on Heat sink and JEB, and time-dependent Bracket-time benefits from explicit temporal branch models. Parametric input fusion by direct concatenation or branch-enhanced fusion improves point-based models substantially on parameterized engineering problems [2510.05995].

Beyond surrogate PDE prediction, learned point-based operators now support direct operator approximation and downstream numerical tasks. NeLo applies a learned Laplacian to point-cloud geometry processing, GNPs are used inside mean-curvature flow and Laplace-Beltrami solvers, PCNO is evaluated on boundary-layer problems, variable-domain Darcy flow, airfoil-plus-flap topology variation, vehicle pressure-load prediction, and parachute inflation dynamics, and RINO demonstrates operator learning from arbitrarily sampled inputs in antiderivative, nonlinear Darcy, and Burgers examples [2409.06506] [2503.04649] [2501.14475] [2407.13010].

The limitations reported across these works are equally consistent. NeLo states that there is no theoretical convergence guarantee as point density grows, that performance may degrade if test point density differs greatly from training, and that the current method is not rotation-invariant [2409.06506]. DSE assumes that only a small number of modes are needed, incurs per-sample transform-matrix construction when point clouds differ across samples, and notes orthogonality and inversion issues outside standard Fourier settings [2305.19663]. PCNO states that topology generalization is not guaranteed, since training on one topology in the airfoil experiment did not generalize well to the other, and it leaves time-dependent PDE operators as an open direction [2501.14475]. RINO emphasizes that resolution independence still requires sufficiently rich sampling and therefore does not evade Nyquist–Shannon-type limits [2407.13010].

A plausible implication is that the field is converging on a shared design doctrine rather than a single dominant architecture. Coordinates, pointwise features, quadrature or density information, spectral bases, tokenized nonlocal interactions, and local differential estimators are increasingly combined rather than treated as mutually exclusive choices. Point-based neural operators, in this sense, are best understood not as one architecture but as a family of operator-learning strategies for irregularly sampled function data and geometry.

Source: https://www.emergentmind.com/topics/point-based-neural-operators