摘要 Abstract
本文提出了一种用于计算由非线性偏微分方程空间离散化产生的非线性矩阵方程低秩解的新框架:低秩Anderson加速(lrAA)。lrAA是对Anderson加速(AA)的一种适应,AA是一种已知的求解非线性不动点问题的方法,被调整为低秩格式。具体来说,lrAA在低秩形式下执行所有的线性和非线性运算,并通过自适应截断容差进行秩截断。我们提出了一种简单的调度策略,根据残差指示器在整个迭代过程中更新截断容差,从而有效控制中间秩和迭代次数。为了对非线性函数进行秩截断,我们提出了一个新的交叉逼近方法,称为基于离散经验插值方法(DEIM)的交叉逼近(Cross-DEIM),并采用自适应误差控制。Cross-DEIM在近似奇异值分解(SVD)和交叉逼近之间进行迭代更新,自然地为每次lrAA迭代引入了热启动策略。我们展示了lrAA在一系列线性和非线性问题上的优越性能,包括来自Laplace方程、Bratu问题、椭圆Monge-Ampère方程和Allen-Cahn方程有限差分离散化的例子。
This paper proposes a new framework for computing low-rank solutions to nonlinear matrix equations arising from spatial discretization of nonlinear partial differential equations: low-rank Anderson acceleration (lrAA). lrAA is an adaptation of Anderson acceleration (AA), a well-known approach for solving nonlinear fixed point problems, to the low-rank format. In particular, lrAA carries out all linear and nonlinear operations in low-rank form with rank truncation using an adaptive truncation tolerance. We propose a simple scheduling strategy to update the truncation tolerance throughout the iteration according to a residual indicator. This controls the intermediate rank and iteration number effectively. To perform rank truncation for nonlinear functions, we propose a new cross approximation, which we call Cross-DEIM, with adaptive error control that is based on the discrete empirical interpolation method (DEIM). Cross-DEIM employs an iterative update between the approximate singular value decomposition (SVD) and cross approximation. It naturally incorporates a warm-start strategy for each lrAA iterate. We demonstrate the superior performance of lrAA applied to a range of linear and nonlinear problems, including those arising from finite difference discretizations of Laplace's equation, the Bratu problem, the elliptic Monge-Amp\'ere equation and the Allen-Cahn equation.