SAV Schemes with Decomposition-Induced Pullback Corrections for Gradient Flows
Abstract
The scalar auxiliary variable (SAV) method constructs linear, unconditionally energy-stable time discretizations of gradient flows. Eliminating the auxiliary variable in a first-order SAV step shows that the state equation is a semi-implicit update augmented by a rank-one positive semidefinite correction from the previous nonlinear force. The multiple-SAV (MSAV) method produces this correction componentwise, with rank up to the number of energy components. This separates two mechanisms usually coupled in MSAV: the number of scalar variables tracking the nonlinear energy and the rank of the correction. We introduce a pullback-corrected SAV (PB-SAV) family that keeps a single scalar auxiliary variable but replaces the rank-one SAV correction by the pullback correction induced by an admissible component decomposition. The correction remains positive semidefinite, has rank at most the number of components, and may change from step to step without changing the scalar auxiliary variable. We prove modified-energy dissipation laws for fixed and step-dependent decompositions, establish first-order convergence after a fixed spatial discretization under standard smoothness and positivity assumptions, derive a refinement identity whose gain is an explicit weighted variance, and give a Sherman--Morrison--Woodbury implementation of the low-rank perturbation of the standard semi-implicit solve. We also show, in finite dimensions, that the pullback correction is the Gauss--Newton matrix of a least-squares representation of the nonlinear energy. Numerical experiments on finite-dimensional gradient flows, Allen--Cahn dynamics, a finite-rank nonlocal gradient flow, and a linear nonlocal Cahn--Hilliard model show regimes in which PB-SAV and SAV differ only in the first-order error constant and regimes in which PB-SAV reduces the trajectory error by a large factor.
Related
- Maximal tails,character fibres and induced modules for pullback Kumjian-Pask algebras
- On Convergence of Gradient Temporal-Difference Learning with Regularized Corrections and Eligibility Traces
- Positivity-preserving scalar auxiliary variable schemes for gradient flows via a quadratic reformulation
- Low regularity maps on Lipschitz manifolds: lifting, pullback of forms, and approximation
- A note on the cohomology of pullback Lie algebroids
BibTeX
@misc{zhang2026schemes,
title = {{SAV Schemes with Decomposition-Induced Pullback Corrections for Gradient Flows}},
author = {Shiheng Zhang and Jie Shen},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2606.18551},
eprint = {2606.18551},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2606.18551},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - SAV Schemes with Decomposition-Induced Pullback Corrections for Gradient Flows
AU - Shiheng Zhang
AU - Jie Shen
PY - 2026
DA - 2026/09/12
PB - arXiv
DO - 10.48550/arxiv.2606.18551
UR - https://arxiv.org/abs/2606.18551
AB - The scalar auxiliary variable (SAV) method constructs linear, unconditionally energy-stable time discretizations of gradient flows. Eliminating the auxiliary variable in a first-order SAV step shows that the state equation is a semi-implicit update augmented by a rank-one positive semidefinite correction from the previous nonlinear force. The multiple-SAV (MSAV) method produces this correction componentwise, with rank up to the number of energy components. This separates two mechanisms usually coupled in MSAV: the number of scalar variables tracking the nonlinear energy and the rank of the correction. We introduce a pullback-corrected SAV (PB-SAV) family that keeps a single scalar auxiliary variable but replaces the rank-one SAV correction by the pullback correction induced by an admissible component decomposition. The correction remains positive semidefinite, has rank at most the number of components, and may change from step to step without changing the scalar auxiliary variable. We prove modified-energy dissipation laws for fixed and step-dependent decompositions, establish first-order convergence after a fixed spatial discretization under standard smoothness and positivity assumptions, derive a refinement identity whose gain is an explicit weighted variance, and give a Sherman--Morrison--Woodbury implementation of the low-rank perturbation of the standard semi-implicit solve. We also show, in finite dimensions, that the pullback correction is the Gauss--Newton matrix of a least-squares representation of the nonlinear energy. Numerical experiments on finite-dimensional gradient flows, Allen--Cahn dynamics, a finite-rank nonlocal gradient flow, and a linear nonlocal Cahn--Hilliard model show regimes in which PB-SAV and SAV differ only in the first-order error constant and regimes in which PB-SAV reduces the trajectory error by a large factor.
ER -