Neither Precision Nor Architecture Alone: Controlled Tests of Failure Remedies for Physics-Informed Neural Networks
Organizations: Hubei University, Wuhan, China
Abstract
Physics-Informed Neural Networks (PINNs) frequently fail on stiff or advection-dominated PDEs, and two recent accounts offer competing remedies: switching from FP32 to FP64 to repair an L-BFGS stopping artifact, or replacing the MLP with a state-space-model (SSM) backbone plus sub-sequence alignment to counter architectural simplicity bias. We test both under matched, seed-paired controls in a pre-registered 144-run study spanning convection, reaction, and wave, plus an independent 85-run convection/wave study; success is relative error below . The two remedies act on disjoint regime-and-seed slices: neither substitutes for the other. On hard convection (), alignment recovers 2/5 seeds in FP32 and 3/5 in FP64, where the unaligned SSM succeeds on 0/5 seeds at either precision and the vanilla MLP moves only from 0/5 to 1/5 across the precision switch---the recoveries trace to the alignment objective, not the backbone. On reaction the backbone alone already succeeds on 3/5--4/5 seeds, so each remedy covers a regime the other does not. Responses are also seed-specific: the same precision switch flips individual seeds in opposite directions and, on wave, lowers median error with no statistically significant success gain. Tightening the inner L-BFGS tolerance in an independent repeated-step runner likewise lowers median error at a large runtime cost, with success counts unchanged. Precision, stopping, backbone, and alignment must therefore be evaluated jointly and reported per seed.
Figures & tables
| Slice | RQ | |
|---|---|---|
| Local study: five-seed conv/react slice | 1–2,4 | 60 |
| Local study: instrumented conv/react subset | 5 | |
| Local study: model-protocol comparison | 2 | 36 |
| Independent study: convection sweep | 2 | 25 |
| Independent study: convection tolerance sweep | 3 | 30 |
| Independent study: wave MLP precision slice | 4 | 20 |
| PDE | Regime | s=0 | s=1 | s=2 | s=3 | s=4 | Rate |
|---|---|---|---|---|---|---|---|
| Conv- 50 | mlp FP32 | 1.005 | 0.896 | 0.806 | 0.914 | 0.800 | 0/5 |
| mlp FP64 | 0.080 | 0.738 | 0.077 | 0.045 | 0.052 | 1/5 | |
| ssm_noalign FP32 | 2.027 | 1.361 | 1.188 | 1.327 | 1.372 | 0/5 | |
| ssm_noalign FP64 | 1.312 | 1.400 | 1.298 | 1.238 | 1.326 | 0/5 | |
| ssm_align FP32 | 0.011 | 0.0498 | 0.549 | 0.386 | 0.904 | 2/5 | |
| ssm_align FP64 | 0.006 | 1.287 | 0.005 | 0.168 | 0.026 | 3/5 |
Neural PDE Solvers
Neural PDE solvers use trainable models to approximate solutions or solution operators for partial differential equations, often with physical constraints.
See all 268 papers