We analyze a goal-oriented adaptive algorithm that aims to efficiently compute the quantity of interest $G(u^\star)$ with a linear goal functional $G$ and the solution $u^\star$ to a general second-order nonsymmetric linear elliptic partial differential equation. The current state of the analysis of iterative algebraic solvers for nonsymmetric systems lacks the contraction property in the norms that are prescribed by the functional analytic setting. This seemingly prevents their application in the optimality analysis of goal-oriented adaptivity. As a remedy, this paper proposes a goal-oriented adaptive iteratively symmetrized finite element method (GOAISFEM). It employs a nested loop with a contractive symmetrization procedure, e.g., the Zarantonello iteration, and a contractive algebraic solver, e.g., an optimal multigrid solver. The various iterative procedures require well-designed stopping criteria such that the adaptive algorithm can effectively steer the local mesh refinement and the computation of the inexact discrete approximations. The main results consist of full linear convergence of the proposed adaptive algorithm and the proof of optimal convergence rates with respect to both degrees of freedom and total computational cost (i.e., optimal complexity). Numerical experiments confirm the theoretical results and investigate the selection of the parameters.
翻译:本文分析了一种面向目标的自适应算法,该算法旨在高效计算量$G(u^\star)$,其中$G$为线性目标泛函,$u^\star$为一般二阶非对称线性椭圆型偏微分方程的解。当前非对称系统迭代代数求解器的分析缺乏泛函分析框架所规定的范数下的压缩性质,这看似阻碍了其在面向目标自适应最优性分析中的应用。为此,本文提出一种面向目标自适应迭代对称化有限元方法(GOAISFEM)。该方法采用嵌套循环结构,包含具压缩性的对称化过程(如Zarantonello迭代)和具压缩性的代数求解器(如最优多重网格求解器)。多种迭代过程需要精心设计的停止准则,以使自适应算法能有效引导局部网格细化与不精确离散近似计算。主要结果包括所提算法的完全线性收敛性,以及关于自由度和总计算成本(即最优复杂度)的最优收敛速率证明。数值实验验证了理论结果并探讨了参数选取。