๐ค AI Summary
To address the slow convergence and poor numerical stability of single-method solvers for partial differential equations (PDEs), this paper proposes a hybrid iterative solving framework based on approximate greedy routing. At each iteration, the framework dynamically selects the optimal solver from a heterogeneous ensemble comprising neural operators and classical numerical solversโwithout requiring ground-truth error feedback. By integrating an error estimation mechanism with a differentiable greedy selection strategy, it strikes a balance between combinatorial optimization complexity and theoretical interpretability. Theoretical analysis establishes a provable suboptimality bound for the routing policy. Experiments on Poisson and Helmholtz equations demonstrate that our method significantly outperforms both individual solvers and state-of-the-art hybrid approaches (e.g., HINTS) in terms of convergence rate and numerical stability.
๐ Abstract
When solving PDEs, classical numerical solvers are often computationally expensive, while machine learning methods can suffer from spectral bias, failing to capture high-frequency components. Designing an optimal hybrid iterative solver--where, at each iteration, a solver is selected from an ensemble of solvers to leverage their complementary strengths--poses a challenging combinatorial problem. While the greedy selection strategy is desirable for its constant-factor approximation guarantee to the optimal solution, it requires knowledge of the true error at each step, which is generally unavailable in practice. We address this by proposing an approximate greedy router that efficiently mimics a greedy approach to solver selection. Empirical results on the Poisson and Helmholtz equations demonstrate that our method outperforms single-solver baselines and existing hybrid solver approaches, such as HINTS, achieving faster and more stable convergence.