SymPlex: A Structure-Aware Transformer for Symbolic PDE Solving
Yesom Park, Annie C. Lu, Shao-Ching Huang, Qiyang Hu, Y. Sungtaek Ju, Stanley Osher
TL;DR
SymPlex formulates symbolic PDE solving as a tree-structured reinforcement learning problem and introduces SymFormer, a structure-aware Transformer that uses tree-relative self-attention and grammar-constrained autoregressive decoding to generate valid symbolic expressions. By optimizing a PDE-based reward and employing a curriculum learning strategy, it can recover exact analytical solutions, including non-smooth and parametric ones, while providing explicit parametric dependencies. Theoretical guarantees establish conditions under which exact symbolic recovery is possible, and experiments demonstrate superior symbolic recovery and interpretability compared with numerical solvers and neural baselines. This framework offers a principled, human-readable alternative to discretized or purely neural PDE solvers, with potential impact on inverse problems, control, and analytical understanding of PDEs, while leaving open avenues for scaling to high-dimensional problems and vocabulary design.
Abstract
We propose SymPlex, a reinforcement learning framework for discovering analytical symbolic solutions to partial differential equations (PDEs) without access to ground-truth expressions. SymPlex formulates symbolic PDE solving as tree-structured decision-making and optimizes candidate solutions using only the PDE and its boundary conditions. At its core is SymFormer, a structure-aware Transformer that models hierarchical symbolic dependencies via tree-relative self-attention and enforces syntactic validity through grammar-constrained autoregressive decoding, overcoming the limited expressivity of sequence-based generators. Unlike numerical and neural approaches that approximate solutions in discretized or implicit function spaces, SymPlex operates directly in symbolic expression space, enabling interpretable and human-readable solutions that naturally represent non-smooth behavior and explicit parametric dependence. Empirical results demonstrate exact recovery of non-smooth and parametric PDE solutions using deep learning-based symbolic methods.
