[论文解读] Termination criteria for inexact fixed point methods

Philipp Birken|arXiv (Cornell University)|Sep 30, 2013

Numerical methods for differential equations参考文献 8被引用 4

一句话总结

本文提出了一种非标准的相对终止准则，用于不精确的定点迭代，证明了无论内迭代求解器的容差如何，都能保证收敛到精确解。结果表明，与标准的相对或绝对准则相比，使用该准则可实现高效且粗粒度的内迭代，尤其适用于流固耦合及PDE中的Picard和Dirichlet-Neumann格式。

ABSTRACT

We analyze inexact fixed point iterations where the generating function contains an inexact solve of an equation system to answer the question of how tolerances for the inner solves influence the iteration error of the outer fixed point iteration. Important applications are the Picard iteration and partitioned fluid structure interaction. We prove that the iteration converges irrespective of how accurate the inner systems are solved, provided that a specific relative termination criterion is employed, whereas standard relative and absolute criteria do not have this property. For the analysis, the iteration is modelled as a perturbed fixed point iteration and existing analysis is extended to the nested case x=F(S(x)).

研究动机与目标

• 解决在内线性系统以不精确方式求解时，控制外定点迭代误差的挑战。
• 分析不同内求解器终止准则对定点外迭代收敛性与精度的影响。
• 将定点迭代的扰动理论扩展至形式为 x = F(S(x)) 的嵌套格式。
• 识别出能确保收敛到精确解、且与内求解容差无关的终止准则。
• 为实际应用（如Picard迭代和流固耦合）提供理论基础坚实且高效的实现策略。

提出的方法

• 将不精确定点迭代建模为由内求解不精确引起的扰动定点格式。
• 将现有扰动理论扩展至处理形式为 x = F(S(x)) 的嵌套定点方程，其中 S(x) 求解子问题。
• 提出一种非标准的相对终止准则：||r|| / ||f(x)|| ≤ τ_r，其中 r 为残差，f(x) 为函数值。
• 证明在该准则下，只要 τ_r < 1，外定点迭代将收敛到精确解，且与 τ_r 无关。
• 与标准相对准则（||r|| ≤ τ_r ||f(x)||）和绝对准则（||r|| ≤ τ_a）对比，后者无法保证收敛到精确解。
• 通过在传输问题和Picard迭代上的数值实验验证理论结果，比较不同准则下的总CG迭代次数与定点迭代次数。

实验结果

研究问题

• RQ1 能否设计一种内求解的终止准则，使得外定点迭代无论内求解容差如何，均能收敛到精确解？
• RQ2 标准相对与绝对终止准则对内求解的影响如何，其对定点外迭代的收敛性与精度有何影响？
• RQ3 在嵌套定点格式中，使用非标准相对终止准则的理论与数值影响为何？
• RQ4 为何标准准则在内求解较粗时无法保证收敛到精确解？
• RQ5 终止准则的选择如何影响流固耦合与Picard方法等应用中的计算效率？

主要发现

• 所提出的非标准相对终止准则（||r|| / ||f(x)|| ≤ τ_r）可确保外定点迭代收敛到精确解，且与 τ_r 无关，只要 τ_r < 1。
• 数值结果表明，使用 τ_r = 1e-1 时，该准则的解精度与内系统精确求解至机器精度相当，但CG迭代次数显著减少。
• 在传输问题中，非标准准则在 Δx = 1/80 时仅需 29,056 次CG迭代，而标准绝对准则在 TOL = 1e-4 时需 70,848 次迭代。
• 非标准准则在 TOL = 1e-1 时的解误差为 8.215e-1，而标准绝对准则的误差为 2.476e-0，表明尽管计算成本更高，精度反而更差。
• 标准相对准则（||r|| ≤ τ_r ||f(x)||）在 TOL 变小时需更精确地求解内系统，导致迭代次数增加——例如，TOL = 1e-4 时需 22,769 次CG迭代，而 TOL = 1e-1 时仅需 12,204 次。
• 非标准准则比标准绝对准则快至五倍，比标准相对准则快至两倍，同时提供更高的解精度。