---
title: 'GeoFAE: Geometric Autoencoder for Inverse Problems'
url: https://www.emergentmind.com/topics/geometric-function-autoencoder-geofae
type: topic
---

# GeoFAE: Geometric Autoencoder for Inverse Problems

Geometric Function Autoencoder (GeoFAE) denotes a geometric function autoencoder introduced as the representation layer of GeoFunFlow for inverse problems on complex geometries. For a geometry \(\Omega\subset\mathbb{R}^d\) and a physical field \(u:\Omega\to\mathbb{R}^p\), it learns an encoder–decoder pair
\[
\mathcal{E}_\theta:\mathbf{c}\mapsto \mathbf{z}\in Z=\mathbb{R}^\ell,\qquad
\mathcal{D}_\phi:(\mathbf{z},q)\mapsto \hat u(q),
\]
where \(\mathbf{c}\) encodes the geometry as a point cloud plus sparse measurements, \(\mathbf{z}\) is a geometry-aware latent code, and \(q\in\Omega\) is any query point. The defining property of GeoFAE is that it encodes an entire field on an irregular domain into a fixed-length latent vector and decodes that vector into a continuous neural field evaluable at arbitrary coordinates; within GeoFunFlow, rectified flow then acts only in latent space rather than directly on discretized meshes [2509.24117].

## 1. Conceptual scope and problem setting

GeoFAE is formulated for geometric functions \(u\in\mathcal{U}(\Omega)=H^s(\Omega;\mathbb{R}^p)\), \(s\ge 1\), defined over sample-dependent geometries. The conditioning object is \(\mathbf{c}=(\Omega,X,\mathbf{y})\), where \(\Omega\) is the geometry, \(X\) is the sensor set, and \(\mathbf{y}\) are noisy readings. The global target in GeoFunFlow is a conditional probabilistic operator
\[
\mathcal{G}^\star:\mathcal{C}\to\bigsqcup_{\Omega\in\Lambda}\mathcal{P}_2(\mathcal{U}(\Omega)),
\qquad
\mathcal{G}^\star(\mathbf{c})=p^\ast(\cdot\mid \mathbf{c}),
\]
mapping conditioning instances to posteriors over fields on \(\Omega\). GeoFAE supplies the fixed-dimensional representation required to make this operator tractable on irregular domains [2509.24117].

The motivation is operational rather than purely notional. Direct diffusion or flow modeling on discretized fields over unstructured meshes is problematic because mesh sizes vary by orders of magnitude, mesh connectivity and geometry vary across samples, and a fixed tensor layout does not exist. GeoFAE addresses this by compressing each field and its geometry into a latent vector \(\mathbf{z}\in\mathbb{R}^\ell\) whose size is independent of mesh size and layout, while preserving continuous decoding back to the original domain [2509.24117].

In a broader geometric reading, GeoFAE belongs to the class of autoencoders understood as jointly learning a manifold and a coordinate chart. That viewpoint emphasizes that autoencoders do not merely compress data: they choose an embedded manifold through the decoder image and a coordinate system through the encoder–decoder pair, which is precisely the interpretive frame within which geometric regularization, latent metrics, and latent manifolds become meaningful design objects [2309.08247].

## 2. Representation of geometry, observations, and fields

Each geometry \(\Omega\) is discretized into a point cloud
\[
V_\Omega=\{x_i\}_{i=1}^m\subset\Omega,
\]
with \(m\) varying across samples and ranging from \(\mathcal{O}(10^3)\) to \(\mathcal{O}(10^5)\) in the reported experiments. Sensors form a subset \(X\subset V_\Omega\), and the observation pattern is encoded by the node-wise mask
\[
M_X(x)=\mathbf{1}_{\{x\in X\}}.
\]
At each node \(x_i\in V_\Omega\), the feature is
\[
\big(M_X(x_i),\;u(x_i)\cdot M_X(x_i)\big),
\]
so values appear only where \(M_X=1\), and the conditioning instance presented to the encoder is
\[
\mathbf{c}=\big(x_i,\;M_X(x_i),\;u(x_i)M_X(x_i)\big)_{i=1}^m
\]
[2509.24117].

GeoFAE is explicitly mesh-agnostic in the sense used by the paper. It does not assume any grid structure; geometric information is conveyed through coordinates \(x_i\), embedded with random Fourier features, and there is no dependence on graph connectivity or regular lattice. This is a specific design choice: the model treats the domain as a point cloud with coordinate embeddings, and its geometric expressivity comes from set processing and coordinate-conditional decoding rather than from explicit mesh operators [2509.24117].

A common misconception is to equate GeoFAE with a discrete mesh autoencoder that reconstructs nodal arrays at a fixed resolution. The decoder instead defines a continuous reconstruction. For any query set \(\{q_j\}_{j=1}^M\subset\Omega\), the latent representation is evaluated at those coordinates, so the same latent code can be queried on the training mesh, on a refined mesh, or at sensor locations chosen after inference. This continuous-query formulation is central to its role in inverse problems, where sensor configurations and target resolutions vary across instances [2509.24117].

## 3. Encoder–decoder architecture

GeoFAE uses a Perceiver-style encoder followed by a CViT-style decoder. At the encoder input, each node \(x_i\in V_\Omega\) receives a coordinate embedding
\[
e_{\rm coord}(x_i)=\big[\cos(2\pi Bx_i),\;\sin(2\pi Bx_i)\big],
\]
where \(B\) has entries sampled from \(\mathcal{N}(0,10)\), together with an MLP embedding of \((M_X(x_i),u(x_i)M_X(x_i))\). These are concatenated into node tokens \(\mathbf{z}_i\in\mathbb{R}^D\), stacked as \(\mathbf{z}\in\mathbb{R}^{m\times D}\) [2509.24117].

Compression from a variable-size point set to a fixed representation is handled by \(P\) learnable latent queries \(\mathbf{z}_q\in\mathbb{R}^{P\times D}\), which attend to the node tokens through Perceiver cross-attention:
\[
\mathbf{z}'=\mathbf{z}_q+\mathrm{MHA}\big(\mathrm{LN}(\mathbf{z}_q),\,\mathrm{LN}(\mathbf{z}),\,\mathrm{LN}(\mathbf{z})\big),
\]
followed by
\[
\mathbf{z}_{\rm agg}=\mathbf{z}'+\mathrm{MLP}(\mathrm{LN}(\mathbf{z}')).
\]
The aggregated tokens are then refined by \(L\) pre-norm Transformer blocks with multi-head self-attention over the \(P\) latent tokens, producing \(\mathbf{z}_L\in\mathbb{R}^{P\times D}\). This token set serves two roles: it is the encoder feature set supplied to the decoder by cross-attention, and it is aggregated into the latent vector \(\mathbf{z}\in\mathbb{R}^\ell\) used by the latent flow model [2509.24117].

The decoder is a coordinate-based neural field. Query coordinates \(\{q_j\}_{j=1}^M\) are embedded by the same type of Fourier mapping into query tokens \(\mathbf{x}_0\in\mathbb{R}^{M\times D}\). These are updated through \(K\) cross-attention blocks with keys and values from \(\mathbf{z}_L\):
\[
\mathbf{x}^{\prime}_{k}
=
\mathbf{x}_{k-1}
+
\mathrm{MHA}\big(\mathrm{LN}(\mathbf{x}_{k-1}),\,\mathrm{LN}(\mathbf{z}_L),\,\mathrm{LN}(\mathbf{z}_L)\big),
\]
\[
\mathbf{x}_{k}
=
\mathbf{x}^{\prime}_{k}
+
\mathrm{MLP}\big(\mathrm{LN}(\mathbf{x}^{\prime}_{k})\big),
\qquad k=1,\dots,K,
\]
followed by a final MLP \(f_{\rm out}\) giving
\[
\hat u(q_j)=f_{\rm out}(\mathbf{x}_{K,j}).
\]
Equivalently, the decoder defines an operator
\[
T_\Omega: Z\to L^2(\Omega;\mathbb{R}^p),\qquad
T_\Omega(\mathbf{z})(q)=\mathcal{D}_\phi(\mathbf{z})(q),
\]
so continuous evaluation is part of the model definition rather than a post-processing step [2509.24117].

Training of the autoencoder stage uses random query points \(\{q_{i,j}\}_{j=1}^M\subset\Omega_i\) and the reconstruction loss
\[
\mathcal{L}_{\mathrm{FAE}}
=
\frac{1}{NM}\sum_{i=1}^N\sum_{j=1}^M
\big|u_i(q_{i,j})-\hat u_i(q_{i,j})\big|^2.
\]
Instances and query points are resampled at each iteration. The paper states that no explicit PDE residual or boundary-term is used; GeoFAE is therefore a purely data-driven geometric function autoencoder, even though it is trained on PDE solution families [2509.24117].

## 4. Latent space, rectified flow, and posterior reconstruction

GeoFAE is the representation engine of GeoFunFlow, not the diffusion model itself. After GeoFAE is trained, it is frozen, and latent rectified flow is trained on its codes. The encoder produces \(\mathbf{z}_1=\mathcal{E}_\theta(\mathbf{c}_{\rm ref})\), the latent encoding of the full field on the full point cloud, and \(\mathbf{z}_c=\mathcal{E}_\theta(\mathbf{c})\), the latent encoding of partial noisy observations. With \(\mathbf{z}_0\sim\mathcal{N}(0,I)\) and \(t\sim\mathcal{U}[0,1]\), the interpolation
\[
\mathbf{z}_t=(1-t)\mathbf{z}_1+t\mathbf{z}_0
\]
is used to train a conditional velocity field \(\mathbf{g}_\psi\) to satisfy
\[
\mathbf{g}_\psi(\mathbf{z}_t,t,\mathbf{z}_c)\approx \mathbf{z}_1-\mathbf{z}_0,
\]
under the latent CRF objective
\[
\mathcal{L}_{\mathrm{CRF}}
=
\mathbb{E}
\left[
\left\|(\mathbf{z}_1-\mathbf{z}_0)-\mathbf{g}_\psi(\mathbf{z}_t,t,\mathbf{z}_c)\right\|_2^2
\right]
\]
[2509.24117].

Inference proceeds in three stages. First, the conditioning instance is encoded to \(\mathbf{z}_c\). Second, latent sampling is performed by solving the ODE
\[
\frac{d\mathbf{z}(t)}{dt}
=
\mathbf{g}_\psi(\mathbf{z}(t),t,\mathbf{z}_c),
\qquad
\mathbf{z}(0)\sim\mathcal{N}(0,I),
\]
whose terminal state \(\mathbf{z}_1=\mathbf{z}(1)\) is a sample from the model posterior over latents. Third, the sampled latent is decoded at any desired query set:
\[
\hat u(q_j)=\mathcal{D}_\phi(\mathbf{z}_1)(q_j).
\]
Repeated latent sampling yields posterior samples and uncertainty quantification in function space [2509.24117].

The paper also provides a function-space error bound that makes the decoder geometry explicit. If the decoder is Lipschitz in latent space,
\[
\|T_\Omega(z)-T_\Omega(z')\|_{L^2(\Omega)}\le L_D\|z-z'\|_2,
\]
if the latent flow approximation error satisfies \(W_2(Q_Z,P_Z^\ast)\le \epsilon_{\rm flow}\), and if the autoencoder reconstruction error is \(\epsilon_{\rm rec}(\mathbf{c})\), then
\[
W_2^{\Omega}\big(\widehat P(\cdot\mid\mathbf{c}),P^\ast(\cdot\mid\mathbf{c})\big)
\le
L_D\,\epsilon_{\rm flow}+\epsilon_{\rm rec}(\mathbf{c}).
\]
With additional Sobolev regularity \(s>d/2\) and quasi-uniform sensors, the paper further states
\[
W_2^{\Omega}(\widehat P,P^\ast)
\le
C\Big(L_{D,s}\epsilon_{\rm flow}
+
h_X^{\,s}\big(L_{D,s}\epsilon_{\rm flow}+\epsilon_{\rm rec,s}(\mathbf{c})\big)\Big),
\]
where \(h_X\sim m^{-1/d}\) is the fill distance of sensors. This makes decoder smoothness and autoencoder fidelity structurally central to posterior accuracy [2509.24117].

A second misconception is that GeoFAE is already physics-informed because it operates on PDE solution families. The paper is explicit that GeoFAE is purely data-driven and that future work could integrate PDE constraints or physics-informed losses into GeoFAE training [2509.24117].

## 5. Relation to the broader literature on geometric autoencoders

GeoFAE sits within a broader body of work that treats autoencoders geometrically, but it occupies a specific niche: continuous function encoding and decoding over irregular geometries. The general geometric perspective states that an autoencoder learns both a manifold and a coordinate chart, and it also highlights the non-uniqueness of both manifold and chart under reconstruction loss alone. That perspective motivates explicit geometric regularization and helps explain why GeoFAE’s architecture treats geometry as a first-class modeling object rather than as an incidental property of the latent code [2309.08247].

Several nearby lines of work illustrate alternative mechanisms for geometry preservation. Geometry-Regularized Autoencoders regularize the bottleneck of an autoencoder toward a precomputed PHATE embedding so that the learned latent representation follows intrinsic data geometry while retaining out-of-sample extension and reconstruction [2007.07142]. Multi-Scale Geometric Autoencoder imposes global distance constraints on the encoder and local Jacobian constraints on the decoder, explicitly separating global and local geometry in an asymmetric design [2509.24168]. Convergent low-bending, low-distortion autoencoder embeddings regularize an encoder using local geodesic distances and Riemannian midpoints so that the latent embedding approaches an isometric, extrinsically flat manifold embedding [2208.10193]. These works suggest a shared theme—control of geometry through distances, Jacobians, or manifold structure—even though GeoFAE itself realizes that theme through point-cloud set encoding and coordinate-conditional neural fields rather than through a direct latent isometry loss.

A second family of related works concerns latent manifold structure rather than encoder–decoder geometry alone. Geometry-Aware Hamiltonian VAE treats the latent space as a Riemannian manifold endowed with a learned metric tensor and uses Riemannian Hamiltonian dynamics for inference [2010.11518]. Geodesic Calculus on Latent Spaces represents latent manifolds as implicit submanifolds, learns an approximate projection onto the manifold, and then computes discrete geodesics and exponential maps [2510.09468]. Lie Group Auto-Encoder represents Gaussian latent distributions as points on a Lie group and uses Lie algebra coordinates, exponential maps, and intrinsic losses [1901.09970]. Relative to these approaches, GeoFAE does not endow latent space with an explicit learned Riemannian metric in the same manner; its principal geometric act is to map between irregular domains and continuous fields through a fixed-dimensional latent representation.

A third, distinct notion of geometric autoencoding appears in sparse coding and optimization-based encoding. “Sparse, Geometric Autoencoder Models of V1” introduces a weighted-\(\ell_1\) regularizer that favors basis functions close to the input in ambient space and yields latent representations “implicitly, locally organized for spectral clustering” [2302.11162]. Gradient Flow Encoding replaces the encoder by latent-space optimization driven only by the decoder and a reconstruction loss, using gradient-flow ODEs and adaptive step size [2105.05031]. These works use “geometric” in senses tied respectively to structured sparsity and optimization geometry, which clarifies that GeoFAE in GeoFunFlow is not a generic label for all geometry-aware autoencoders but a specific architecture for geometry-conditioned function representation.

## 6. Empirical behavior, limitations, and directions

On the five PDE benchmarks reported in the GeoFunFlow paper—Darcy flow on domains with notches, Cylinder flow, Plasticity, Airfoil flow, and Ahmed-body aerodynamics—GeoFAE alone, used as a deterministic reconstructor, achieves state-of-the-art or best deterministic accuracy: it has the lowest average relative \(L^2\) error on Darcy, Cylinder, Plasticity, and Airfoil, and is second-best on Ahmed body. The paper further states that it outperforms Geo-FNO and Transolver significantly on several tasks, despite having fewer parameters, with GeoFAE at approximately \(7\)M parameters versus approximately \(18\)M for those baselines [2509.24117].

The integration with latent rectified flow changes the objective from deterministic regression to posterior modeling. GeoFunFlow therefore gives slightly higher reconstruction error than GeoFAE alone, which the paper explicitly identifies as expected because the model is optimized for posterior modelling rather than pure regression. The gain is that posterior sampling and calibrated uncertainty quantification become available while maintaining continuous decoding over complex geometries [2509.24117].

Ablations sharpen the interpretation of the encoder and latent space. On the Cylinder dataset, increasing the number of latent queries \(P\) yields only small gains, and a single Perceiver cross-attention block is often sufficient, indicating that the Perceiver acts robustly as a set-to-fixed-latent compressor. Latent sampling is also efficient: even a single integration step for the rectified flow produces reconstructions with very low error, and additional ODE steps yield diminishing returns. The paper interprets this as evidence that GeoFAE’s latent space is well-structured so that the learned flow is nearly straight and short [2509.24117].

The stated limitations are equally specific. GeoFAE treats the domain as a point cloud with coordinate embeddings and does not exploit mesh connectivity or explicit geometric priors such as boundary normals or topology features beyond coordinates. Posterior variance is sometimes relatively high, which the authors take as evidence that richer geometric encoders, including graph-based or equivariant architectures, could improve representation quality. The framework handles up to \(\mathcal{O}(10^5)\) nodes in the reported experiments, but scaling to millions of elements would require additional architectural and computational optimization. Extensions to more complex multi-physics operators, higher-dimensional domains, and direct incorporation of physics-informed constraints are identified as natural next steps [2509.24117].

Taken together, these properties define GeoFAE as a geometry-aware, mesh-size-agnostic, continuous neural-field autoencoder specialized to inverse operator learning on complex geometries. A plausible implication is that its main contribution is not merely compression, but the establishment of a common latent interface between irregular geometric domains and latent generative modeling. Within the current formulation, however, that interface remains coordinate-based and data-driven rather than connectivity-aware or physics-constrained.

Source: https://www.emergentmind.com/topics/geometric-function-autoencoder-geofae