Search papers, labs, and topics across Lattice.
This paper introduces the Hartley Neural Operator (HNO), which utilizes the real Discrete Hartley Transform to learn solution operators for partial differential equations (PDEs), contrasting with the Fourier Neural Operator (FNO) that employs complex FFTs. The authors demonstrate that HNO retains twice as many frequency corners while using fewer parameters, leading to improved performance for self-adjoint elliptic operators, while FNO excels in time-dependent scenarios due to its ability to represent phase. Their findings establish a predictive framework for selecting the appropriate spectral basis based on the symmetry of the solution operator, revealing a nuanced relationship between operator characteristics and spectral representation.
The Hartley Neural Operator outperforms its Fourier counterpart for elliptic PDEs by leveraging real-valued representations, challenging the dominance of complex methods in neural operator design.
Fourier Neural Operators (FNO) learn solution operators of partial differential equations by parameterizing global convolutions in the complex Fourier domain. For real-valued PDE solutions, the complex FFT carries representational redundancy through conjugate symmetry. We introduce the Hartley Neural Operator (HNO), the exact real-valued mirror of FNO: it replaces the FFT with the purely real Discrete Hartley Transform and learns a single real multiplier per retained spectral mode, with no complex arithmetic. Because the real Hartley spectrum is not halved by conjugate symmetry, HNO retains twice as many frequency corners as FNO but one real weight where FNO carries a complex pair, so the two operators are iso-parametric at equal width and differ only in spectral basis. Our central thesis is that the best basis is a property of the operator. Self-adjoint elliptic operators (Poisson, biharmonic) have real, symmetric Green's functions that the real Hartley multiplier diagonalizes exactly, and HNO is favored there. Time-dependent operators carry phase, from oscillation in the wave equation to transport in advection, Burgers, and Navier-Stokes, which a real diagonal multiplier cannot represent, so FNO is favored there, and increasingly so with the operator's phase content, leaving the phaseless heat equation as the borderline case. Training both operators identically and benchmarking across PDE classes, initial-condition families, and boundary conditions, we find an elliptic-versus-time-dependent split that is monotone in operator phase content and matches the Green's-function theory we develop. Rather than a universal winner, our findings give a predictive rule: match the spectral basis to the symmetry of the solution operator.
