Skip to content

算法Algorithm

加性 Schwarz 子域预条件器

Additive Schwarz preconditioner

把重叠子域的同一份残差校正相加,证明所得算子正定,说明稳定分解和粗空间为何控制长程误差,并实算九节点链。

形式陈述 ​

一个全局大系统可以拆成许多容易求解的局部系统。但重叠区域会收到多份校正,怎样组合才得到可靠的预条件器?

设 A∈Rn×n、n≥1 且 A=AT≻0,非空子域索引集合 Ω1,…,Ωs 覆盖全部未知量。矩阵 Ri 从全局向量提取 Ωi 分量,RiT 把局部向量以零延伸回全局。定义

Ai=RiARiT,B=∑i=1sRiTAi−1Ri.

Ai 是正定主子矩阵,局部逆均有意义。加性 Schwarz 预条件器应用 B:对同一输入残差 r,独立求解 Aizi=Rir,再返回 z=∑iRiTzi。重叠位置真正相加,不能一会儿覆盖、一会儿平均。

B 对称,而且对非零 r,覆盖性保证至少一个 Rir≠0,所以

rTBr=∑i(Rir)TAi−1(Rir)>0.

因此 B 是固定线性正定算子,可作 M−1。若各子域并未覆盖全部自由度,存在非零向量在所有子域限制都为零,B 就可能奇异。

直觉

所有子域都拿到同一时刻的残差,各自提出一份局部修正,然后一起合并。这让局部求解能够并行。若求完一个子域立即更新全局残差,再把新残差给下一个子域,得到的则是乘性 Schwarz,更新算子的顺序已经改变。

局部解相加与全局粗方向

BA 还能写成一组子空间投影之和:

BA=∑iQi,Qi=RiTAi−1RiA.

每个 Qi 都是到子域支撑空间的 A-正交投影,因为 Qi2=Qi,且 AQi=QiTA。因此 ‖Qi‖A=1,有 λmax(BA)≤s。这个粗界足以提醒我们:正定预条件器不等于步长为一的定常迭代一定收敛。做 x←x+τB(b−Ax) 时,还需 0<τ<2/λmax(BA)。

粗空间修补的是稳定分解 ​

若每个向量都能分解为 u=∑iui,其中 ui 支撑在子域内,并满足

∑i‖ui‖A2≤K0‖u‖A2,

则 λmin(BA)≥1/K0。证明使用 ui=Qiui,再在各子域空间的直和上应用 Cauchy–Schwarz 不等式:

‖u‖A2=∑i(ui,Qiu)A≤(∑i‖ui‖A2)1/2(∑i‖Qiu‖A2)1/2.

又有 ∑i‖Qiu‖A2=(BAu,u)A,平方整理便得到下界。若一个缓慢变化的全局向量被硬切成局部片段,片段边缘可能产生较大能量,K0 因而变坏。

加入满列秩粗基矩阵 Z,令

B2=B+Z(ZTAZ)−1ZT,

便允许分解额外使用一个跨域分量。它仍正定,但良好尺度无关的条件数还需要证明合适的稳定分解与子空间相互作用界,不能只凭“加了一层”宣布获得。

例子与边界

取九节点链 A=tridiag(−1,2,−1),两个子域为 Ω1={0,1,2,3,4,5}、Ω2={3,4,5,6,7,8}。每个局部矩阵都是六阶链矩阵。

对 r=(1,…,1)T,两份局部解都为 (3,5,6,6,5,3)T。延伸相加得到

Br=(3,5,6,9,10,9,6,5,3)T.

例如全局节点 4 在第一子域收到 5,在第二子域也收到 5,合并为 10。这份结果是预条件校正,不是宣称已经精确解出全局方程。

再取一个跨域帐篷方向

p=(1,2,3,4,5,4,3,2,1)T.

有 Ap=(0,0,0,0,2,0,0,0,0)T,且 pTAp=10。令 Z=p,粗校正为 B0=ppT/10,所以 B0Ap=p。

从误差 e=p 出发,以相同步长 τ=1/2 做一次校正。单层结果为

e1=17(5,10,15,15,15,15,15,10,5)T;

两层加性结果为

e2=114(3,6,9,2,−5,2,9,6,3)T.

能量平方分别为

‖p‖A2=10,‖e1‖A2=150/49,‖e2‖A2=125/98.

粗方向在这个例子中明显改善一步效果。不过两层加性校正并不把 p 一步归零:虽然 B0Ap=p,局部项 BAp 也同时参加,不能只看粗项的精确性就忽略其余校正。步长一还可能过度修正;本例的三个投影给出 λmax(B2A)≤3,所以 1/2<2/3 是有依据的安全选择。

推论与应用

局部矩阵来自稀疏矩阵的索引限制,可以预先分解并复用。若每个子域又有内部与界面变量,可用静态凝聚减少局部解的内部工作;这不改变加性组合要求所有子域使用同一输入残差。

一次应用成本是各局部求解成本之和,加上限制与重叠求和的索引工作;并行墙钟时间还受最大局部问题、通信和归约影响。增大重叠通常改善局部信息传播,却增加重复计算、存储与通信,应按总成本取舍。

只把局部解的一部分延伸回全局的 restricted 变体,或在输出端单独加权,可能失去对称性。它们有各自的算法用途,但不自动继承本页的正定证明。要用于标准 PCG,设置和每次应用都须保持同一个线性、自伴、正定算子。

参考资料
  • Yousef Saad, Iterative Methods for Sparse Linear Systems, 2nd ed., 2003, §§14.3.3–14.3.4, Algorithm 14.6 and Theorem 14.7:作者教材。加性应用、子空间投影与稳定分解下界。
关系图谱10 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系