Shearlet Neural Operators for Anisotropic-Shock-Dominated and Multi-scale parametric partial differential equations
For researchers using neural operators to solve parametric PDEs with shocks or multiscale features, SNO offers a more effective architecture than FNOs, though the improvement is incremental (replacing one spectral transform with another).
Shearlet Neural Operators (SNOs) replace Fourier transforms with shearlet-based representations to better resolve anisotropic structures, sharp gradients, and discontinuities in parametric PDEs. Across seven benchmark families, SNO consistently improves predictive accuracy and feature fidelity over Fourier Neural Operators, with largest gains in anisotropic and discontinuity-dominated settings.
Neural operators have emerged as powerful data-driven surrogates for learning solution operators of parametric partial differential equations (PDEs). However, widely used Fourier Neural Operators (FNOs) rely on global Fourier representations, which can be inefficient for resolving anisotropic structures, sharp gradients, and spatially localized discontinuities that arise in shock-dominated and multiscale regimes. To address these limitations, we introduce the Shearlet Neural Operator (SNO), a neural operator architecture that replaces the Fourier transform with a shearlet-based representation. Shearlets offer directional, multiscale, and spatially localized atoms with near-optimal sparse approximation of anisotropic features, providing an inductive bias aligned with PDE solutions containing edges, fronts, and shocks. SNO learns in the shearlet domain and reconstructs predictions via the inverse transform, retaining efficient spectral computation while improving locality and directional selectivity. Across seven benchmark PDE families, including strongly anisotropic advection, anisotropic diffusion, and nonlinear conservation laws with straight, curved, interacting, spiral, and polygonal shock structures, SNO consistently improves predictive accuracy and feature fidelity over FNO baselines, with the largest gains observed in anisotropic and discontinuity-dominated settings.