---
title: Equivariant Operator Discovery (EqOD)
url: https://www.emergentmind.com/topics/equivariant-operator-discovery-eqod
type: topic
---

# Equivariant Operator Discovery (EqOD)

Equivariant Operator Discovery (EqOD) is the task of learning or enumerating operators whose input–output relationship respects symmetry, while inferring the relevant operator structure from data rather than assuming the symmetry group a priori. In the literature, the term covers several closely related programs: learning latent-space operators that canonicalize transformed observations, discovering equivariance through parameter-sharing schemes, identifying symmetry-respecting Koopman or neural operators, reducing PDE libraries by detected invariances, and constructing complete algebraic or geometric classes of equivariant operators [2204.03640][2602.18406][2605.11524]. The name is also used specifically by “EqOD: Symmetry-Informed Stability Selection for PDE Identification,” a PDE-identification pipeline that combines symmetry detection, library reduction, sparse regression, and a residual-based fallback [2605.11524].

## 1. Core definitions and mathematical formulations

The common starting point is equivariance. A map $f$ is equivariant to a group $G$ acting on inputs via $g \cdot x$ and on outputs via a representation $\rho(g)$ if
$$
f(g \cdot x) = \rho(g) f(x).
$$
EqOD asks how such structure can be inferred or instantiated from data. In the latent-operator formulation, an encoder $\phi$ maps inputs to latent codes $z=\phi(x)$, and one learns latent operators $L_g$ such that
$$
L_g \phi(x) \approx \phi(g \cdot x),
$$
or, in canonicalization form,
$$
L_g^{-1}\phi(g \cdot x) \approx \phi(x).
$$
This places the discovery problem at the level of the representation $\phi$ and the operator family $\{L_g\}$ rather than the raw input architecture [2602.18406].

A second formulation represents equivariance structurally through parameter sharing. For a linear map $f(x)=Wx$ under input–output permutations $(\pi_N,\pi_M)$, equivariance is equivalent to the weight-tying constraint
$$
W_{m,n}=W_{\pi_M(m),\pi_N(n)}.
$$
Learning the sharing pattern is therefore a way of learning the discrete group action itself. In “Equivariance Discovery by Learned Parameter-Sharing,” the sharing configuration is encoded by an assignment matrix $A$ and free parameters $\psi$, with tied parameter vector $\theta=A\psi$; the discovery problem is posed as a bilevel optimization over $A$ [2204.03640].

A third formulation appears in operator-theoretic and Koopman settings. There, EqOD learns a finite-dimensional surrogate operator $K$ in the commutant of a representation $\rho(G)$, enforcing
$$
K \rho(g)=\rho(g) K.
$$
Under suitable assumptions on system equivariance, observable equivariance, and sampling-distribution invariance, the optimal EDMD matrix is itself equivariant, which permits group-convolutional parameterizations and block diagonalization by the generalized Fourier transform [2411.00905].

These formulations are not identical, but they share a common logic: symmetry is encoded as an operator relation, and the discovery problem consists in recovering the operator, the operator class, the canonical frame, or the symmetry-restricted hypothesis space from transformed data.

## 2. Principal methodological families

EqOD is not a single algorithmic template. The current literature contains several distinct families that differ in what is assumed known, what is learned, and where equivariance is enforced.

| Family | Learned or fixed object | Representative paper |
|---|---|---|
| Latent canonicalization | $\phi$, $L_g$, or Lie-algebra coordinates for input/output transforms | [2602.18406], [2605.18606] |
| Parameter-sharing discovery | Sharing matrix $A$ and induced equivariant operator class | [2204.03640] |
| Koopman / EDMD / OpInf | Equivariant $K$, stitched local operators, or reduced operators with symmetry coupling | [2106.15678], [2411.00905], [2507.18780] |
| Sparse PDE identification | Symmetry-reduced library and sparse coefficients | [2605.11524] |
| Algebraic or geometric classification | $\star_G$-algebra operators, multigraph-indexed operator bases, isotropic kernels | [2605.20440], [2605.16847], [2108.09541] |

One methodological axis concerns whether the group is known or discovered. Latent-operator approaches are described as learning latent operators “directly from data, using pairs or sets of transformed views, without requiring the group to be explicitly encoded in the network architecture,” although the fixed-operator variant still uses a generic cyclic shift template and the learned-operator variant uses a periodicity prior [2602.18406]. By contrast, PACE-FNO assumes known continuous symmetries of evolution equations on periodic domains and uses a learned Lie-algebra coordinate estimator only for the frame parameters, while the canonical operator itself is learned in the reduced frame [2605.18606].

Another axis concerns whether discovery is architectural, statistical, or algebraic. Learned parameter-sharing treats equivariance as a partition over parameters and evaluates recovery via the partition distance, which lower-bounds the discrepancy between equivariance groups [2204.03640]. Sparse PDE identification treats symmetry as a mechanism for pruning candidate libraries, for example by detecting Galilean invariance and removing terms ruled out by an exclusion theorem [2605.11524]. Algebraic approaches treat equivariance as intrinsic to the multiplication law or the combinatorics of contractions: in the $\star_G$ tensor algebra, equivariance is built into the product, while on constant-curvature spaces nonlinear polynomial differential operators are classified by equivalence classes of multigraphs [2605.20440][2605.16847].

A plausible implication is that “EqOD” is best understood as a research program organized around the discovery of symmetry-compatible operator structure, not as a single canonical architecture.

## 3. Latent operators, canonicalization, and robust recognition

The latent-operator approach in “Latent Equivariant Operators for Robust Object Recognition: Promise and Challenges” provides a minimal EqOD instantiation in which a linear encoder $f_E$ maps a flattened $3\times 28\times 28$ input to a $d$-dimensional latent with $d=70$, and a two-layer MLP classifier $f_{\mathrm{CLF}}$ operates on canonicalized latents [2602.18406]. The paper studies discrete rotations, with 10 steps of $36^\circ$, and toroidal translations on a $28\times 28$ grid with stride 2 and order 14 per axis. The method learns $\phi$ jointly with either a pre-defined block-circulant shift operator or a learned $d\times d$ operator with orthogonal initialization and a periodicity loss $L_{\mathrm{op}}=\|L^N-I\|_2$, with $N=d=70$.

The training objective combines task supervision with canonicalization consistency. For transformed views $x_1=T^{k_1}(x)$ and $x_2=T^{k_2}(x)$, the model forms
$$
Z_1=L^{-k_1}\phi(x_1), \qquad Z_2=L^{-k_2}\phi(x_2),
$$
and minimizes
$$
L_{\mathrm{reg}}=\|Z_1-Z_2\|^2,
$$
together with cross-entropy on $f_{\mathrm{CLF}}(Z_1)$. The total loss is
$$
L=L_{\mathrm{CE}}+\lambda L_{\mathrm{reg}} \;(+ \alpha L_{\mathrm{op}} \text{ if using learned operator}).
$$
This makes canonicalization, rather than architectural group convolution, the central mechanism [2602.18406].

The experiments use “noisy MNIST” with distractor backgrounds. Digits are thresholded, colored blue, then pasted over random black–white checkerboards; class 9 is excluded to avoid confusion with rotated 6. Training sees only rotations in $\{-72^\circ,-36^\circ,0^\circ,36^\circ,72^\circ\}$ and translations in $\{-4,-2,0,2,4\}$ pixels, while OOD consists of the remaining rotations and larger translations. Pose selection at inference uses a canonical reference bank and k-NN over candidate latent actions, with default main-plot settings $N=2000$ and $K=1$ [2602.18406].

The reported trends are sharply separated. The baseline using the same encoder and classifier but no operator peaks in-domain and collapses OOD: for rotation it reaches approximately $78$–$79\%$ in-domain and drops to approximately $25$–$31\%$ at $\pm 144^\circ$ and $180^\circ$; for vertical translation it reaches approximately $78.5$–$83.3\%$ in-domain and drops to $13.6$–$23.4\%$ at far shifts. Operator models are flat across the full range. With ground-truth pose, both fixed and learned operators achieve approximately $95$–$96\%$ across all rotations and translations; with automatic k-NN pose inference, rotation is approximately $85$–$89\%$ and translation approximately $91$–$94\%$ across the range. For compound translations, training uses only single-axis transforms, yet pre-defined and learned operators maintain high accuracy across most unseen $(x,y)$ combinations, and the learned operator is sometimes slightly better in corner regions [2602.18406].

The paper positions this as an alternative to classic group-equivariant neural networks. G-CNNs and related architectures build equivariance by design and need the group a priori; latent EqOD instead learns usable operators in a learned latent space from transformed examples. The trade-off is explicit: latent EqOD offers flexibility and data-driven discovery of operators, but the paper emphasizes the lack of formal guarantees on OOD equivariance, open questions about operator placement in deep stacks, inference overhead from k-NN pose selection, and added difficulty for continuous, non-commutative, or 3D transformations [2602.18406].

## 4. Dynamical systems, neural operators, and PDE identification

EqOD has developed particularly rapidly in dynamical systems and PDE learning. In the Koopman setting, “Data-Driven Operator Theoretic Methods for Phase Space Learning and Analysis” gives a symmetry-based transport theorem for local Koopman operators on invariant subspaces. If $M_j=g\cdot M_i$ and the dictionary is consistent, then the local operators satisfy
$$
K_j=\gamma^{-1}K_i\gamma,
$$
and a global operator is assembled as a block-diagonal stitched operator
$$
K_S=\operatorname{diag}(K_1,\dots,K_v),
$$
with stitched observables weighted by characteristic functions of the invariant subspaces. The same paper gives a residual criterion for when new data require a new dictionary and shows that the multiplicity of the unit eigenvalue increases as larger invariant subspaces are discovered [2106.15678].

“Group-Convolutional Extended Dynamic Mode Decomposition” extends this line by showing that, under equivariance assumptions on the system, the observables, and the sampling distribution, the optimal EDMD matrix is equivariant and can be represented as a group convolution. In Fourier space, the operator decomposes into independent blocks indexed by irreducible representations, reducing both learning and spectral analysis to smaller problems. On the 2D Kuramoto–Sivashinsky PDE with $G=\mathbb{Z}_{16}\times\mathbb{Z}_{16}$, GC-EDMD accurately learns dynamics using only 20 input-output pairs, while unconstrained EDMD fails; on the $\lambda$–$\omega$ spiraling wave system with $G=\mathbb{Z}_{48}\times\mathbb{Z}_{48}$, it computes approximately $69\mathrm{k}$ eigenvalues by solving 2304 independent $30\times 30$ problems rather than one $69{,}120\times 69{,}120$ problem [2411.00905].

A related shift-equivariant reduced-order setting appears in “Symmetry-reduced model reduction of shift-equivariant systems via operator inference.” There, the solution is rewritten in a traveling frame,
$$
x' = x-\xi(t), \qquad \hat u(x',t)=u(x'+\xi(t),t),
$$
so that
$$
\partial_t \hat u = f(\hat u) + \dot{\xi}(t)\,\partial_x \hat u.
$$
After POD on aligned snapshots, the reduced dynamics acquire an explicit symmetry-coupling term,
$$
\dot a_i = d_i + A_{ij} a_j + B_{ijk} a_j a_k + \dot c(a)\,(b_i + C_{ij}a_j),
$$
with a rational speed law
$$
\dot c(a) = -\,\frac{e + p_j a_j + Q_{jk} a_j a_k}{w + s_j a_j}.
$$
For the Kuramoto–Sivashinsky equation, SR-OpInf with re-projection attains effectively zero training loss, matches intrusive SR-Galerkin, and yields relative snapshot errors of approximately $2.743\%$ for $r=4$ and approximately $0.319\%$ for $r=7$, while the non-reprojected model is unstable [2507.18780].

PACE-FNO takes a canonicalization route closer in spirit to latent EqOD, but in PDE operator learning. It estimates a Lie-algebra coordinate $\xi_0$, maps the input to a canonical frame, applies a standard Fourier Neural Operator $G_\theta$, and restores the output to the target frame:
$$
\hat u_{\mathrm{pred}} = C_{\exp(\xi_T)} \circ G_\theta \circ C_{\exp(-\xi_0)}(a_{\mathrm{OOD}}).
$$
The FNO core is unchanged; equivariance is enforced by the input and output transformations. On 1-D and 2-D Burgers, shallow-water, and Navier–Stokes equations on periodic domains, PACE-FNO matches the in-distribution accuracy of standard neural operators and reduces out-of-distribution relative error by up to $12\times$ over FNO with symmetry augmentation under translations and Galilean shifts, with smaller gains for coupled rotation–translation shifts. In 1-D Burgers, the OOD relative error drops from $0.3222$ for FNO+Aug to $0.0267$ for PACE-FNO (one-shot) and $0.0259$ with test-time adaptation [2605.18606].

The 2026 PDE-identification method titled “EqOD: Symmetry-Informed Stability Selection for PDE Identification” uses the acronym in a narrower sense. It begins from a standard weak-form library
$$
L_{\mathrm{std}}=\{u,\,u^2,\,u^3,\,u_x,\,u_{xx},\,u_{xxx},\,u_{xxxx},\,u\,u_x,\,u\,u_{xx},\,u^2u_x\},
$$
and combines two library-reduction mechanisms. If Galilean invariance is detected from trajectory data by a six-term weak-form structural test, the method uses a reduced library; otherwise it applies randomized LASSO stability selection. The sparse regression stage solves weak-form LASSO, chooses $\lambda$ by 5-fold cross-validation, applies adaptive thresholding and OLS debiasing, and then compares the reduced-library residual to the full-library residual through a fallback rule with $\gamma=1.5$ on the symmetry path and $\gamma=1.2$ on the stability path [2605.11524].

This PDE-identification EqOD comes with explicit benchmark results. On 8 PDEs at 4 noise levels, it attains $F_1 = 1.000 \pm 0.000$ on Heat at $20\%$ noise, where WF-LASSO obtains $0.475 \pm 0.181$, official PySINDy 2.0 obtains $0.000$, and the WSINDy reimplementation obtains $0.789$. Under the strict criterion that the mean F1 difference exceeds the larger of the two standard deviations, EqOD wins 7 of 32 cells, WF-LASSO wins none, and the remaining 25 cells are ties. External validation on WeakIdent and PINN-SR datasets gives $F_1=1.000$ on all 5 clean benchmarks. At the same time, the paper is explicit that the Galilean library reduction is proved only under autonomy and library assumptions, and that formal guarantees for correlated PDE design matrices remain open [2605.11524].

## 5. Algebraic, representation-theoretic, and geometric viewpoints

A different EqOD strand treats equivariance as an algebraic object that can itself be scanned, decomposed, or classified. In “Group-Algebraic Tensors: Provably-optimal Equivariant Learning and Physical Symmetry Discovery,” the $\star_G$ tensor algebra makes equivariance intrinsic to multiplication. For tensors $A$ and $B$ with a group index, the product is
$$
(A\star_G B)_{ij}(c)=\sum_{k=1}^m \sum_{a\in G} A_{ik}(a) B_{kj}(a^{-1}c),
$$
and in the Fourier domain it decouples per irrep. The paper proves an Eckart–Young optimality theorem for the $\star_G$-SVD, shows that product groups factor through Kronecker products of generalized Fourier transforms, and uses a group scan over candidate $G$ to score symmetry hypotheses. On full QM9, $\star_G$-SVD with ridge regression provides closed-form predictions at approximately $50$–$90\times$ fewer parameters than parameter-matched MLPs, and per-irrep decompositions recover selection-rule structure such as scalar properties being $A_1$-dominated and dipole components being $T_1$-dominated [2605.20440].

On geometric PDE operator spaces, “Equivariant nonlinear partial differential operators on constant curvature spaces” gives a classifying space for nonlinear differential operators on simply connected constant-curvature manifolds that are equivariant under the full isometry group. The operator class consists of nonlinear operators that can be written as polynomials in linear operators, and the paper shows that the classifying space is spanned by equivalence classes of multigraphs. For a multigraph $\gamma$,
$$
N_\gamma f = \operatorname{tr}_{E[\gamma]}\bigl(\nabla^{V[\gamma]}f\bigr),
$$
where vertices encode symmetric covariant derivatives and edges encode contractions by the metric. The map from graph classes to operators is surjective, and for graph classes with at most $p$ edges and $p\le d$, the restriction is bijective. The same framework exposes dimension-dependent linear relations, such as the order-6 Hessian identities that vanish in $d=2$ [2605.16847].

A more classical linear-operator perspective appears in “Rotation Equivariant Operators for Machine Learning on Scalar and Vector Fields.” There, linear, translation-invariant, rotation-equivariant operators on scalar, vector, and tensor fields are characterized as tensor-field convolutions with radially symmetric kernel fields. For vectors, isotropy forces the matrix kernel to take the form
$$
K(r)=a(\|r\|)I+b(\|r\|)\hat r \hat r^\top
      = k_T(\|r\|)P_T(\hat r)+k_L(\|r\|)P_L(\hat r),
$$
and differential operators such as gradient, divergence, curl, and Laplacian appear as distributional-kernel instances of the same theory. This viewpoint places EqOD in an explicitly interpretable function class: one discovers radial functions rather than unrestricted filters [2108.09541].

These algebraic and geometric formulations broaden EqOD beyond data-driven canonicalization. They frame symmetry discovery as basis discovery, per-irrep decomposition, group scanning, or operator-library generation. A plausible implication is that, in some domains, the primary EqOD object is not a neural architecture or latent code but a structured operator algebra.

## 6. Guarantees, trade-offs, and open questions

A recurrent issue in EqOD is the tension between flexibility and guarantees. Classic group-equivariant neural networks provide stronger theoretical guarantees but require the group a priori; latent EqOD explicitly trades those guarantees for flexibility and data-driven discovery of operators [2602.18406]. PACE-FNO, by contrast, restores a degree of analytical control by assuming known continuous symmetry groups and proving an OOD error decomposition
$$
\sup_{a \in A_{\mathrm{OOD}}}
\| C_{g_T} G_\theta C_{g_0}^{-1}(a) - G_T(a) \|
\le \epsilon + M_0 \epsilon_{\mathrm{geo},0} + \epsilon_{\mathrm{geo},T},
$$
which separates canonical-set approximation from alignment error [2605.18606].

Another central issue is identifiability. In learned parameter-sharing, the sharing matrix $A$ is identifiable only up to permutation of cluster labels, and finite data can make spurious partitions locally optimal. The paper’s theoretical analysis is for Gaussian shared-mean estimation, where the mean-squared error decomposes into a bias term and a variance term depending on $\operatorname{rank}(A)$, and the finite-sample gap bound implies that, for strong sharing, relatively larger validation splits can be optimal—often much larger than the common 80/20 split [2204.03640]. This is a precise theoretical result, but it is tied to a restricted data model.

Several misconceptions are corrected by the existing literature. One is that EqOD necessarily removes the need for prior symmetry information. This is not uniformly true. The latent-operator formulation avoids explicit group-engineered layers but still uses cyclic templates or periodicity priors; PACE-FNO requires the continuous symmetry group and its field action; the constant-curvature classification presupposes the full isometry group; and GC-EDMD assumes an equivariant dictionary and group action [2602.18406][2605.18606][2411.00905][2605.16847]. Another misconception is that equivariance discovery is inherently neural. The literature includes bilevel partition learning, block-diagonal Koopman constructions, sparse regression over symmetry-pruned libraries, exact algebraic SVDs, and multigraph-indexed differential-operator bases [2204.03640][2106.15678][2605.11524][2605.20440].

Open problems are correspondingly diverse. Latent canonicalization faces open questions about operator placement in deep stacks, stability of learned operators, scaling to continuous, non-commutative, or 3D groups, and the computational cost of k-NN pose inference [2602.18406]. PDE-identification EqOD has a Galilean exclusion theorem under explicit assumptions, but formal guarantees for correlated weak-form PDE design matrices remain open [2605.11524]. Constant-curvature operator classification is restricted to nonlinear operators that are polynomials in linear operators and to simply connected model spaces [2605.16847]. The $\star_G$ algebra is finite-group based, so continuous symmetries are handled through finite subgroups or discretizations [2605.20440].

Taken together, these works suggest that EqOD is best regarded as a unifying perspective on symmetry-aware operator learning rather than a settled framework. Its central promise is to convert symmetry from a fixed architectural prior into a discoverable, testable, or classifiable property of operators. Its central challenge is that the meaning of “discovery” varies sharply across settings: discovering a latent operator, a parameter-sharing partition, a canonical frame, a sparse PDE library, a Koopman commutant element, or a multigraph basis are mathematically related but methodologically distinct enterprises.

Source: https://www.emergentmind.com/topics/equivariant-operator-discovery-eqod