形式陈述
一个全局大系统可以拆成许多容易求解的局部系统。但重叠区域会收到多份校正,怎样组合才得到可靠的预条件器公理库预条件Preconditioning · Preconditioner用易应用的近似逆改变等价线性系统的尺度与 Krylov 几何,并计入构造、存储和每步应用成本。?
设 、 且 ,非空子域索引集合 覆盖全部未知量。矩阵 从全局向量提取 分量, 把局部向量以零延伸回全局。定义
是正定主子矩阵,局部逆均有意义。加性 Schwarz 预条件器应用 :对同一输入残差 ,独立求解 ,再返回 。重叠位置真正相加,不能一会儿覆盖、一会儿平均。
对称,而且对非零 ,覆盖性保证至少一个 ,所以
因此 是固定线性正定公理库正定与半正定矩阵Positive definite matrix · Positive semidefinite matrix · PSD matrix由二次能量严格为正或非负定义的实对称与复 Hermitian 矩阵。算子,可作 。若各子域并未覆盖全部自由度,存在非零向量在所有子域限制都为零, 就可能奇异。
直觉
所有子域都拿到同一时刻的残差,各自提出一份局部修正,然后一起合并。这让局部求解能够并行。若求完一个子域立即更新全局残差,再把新残差给下一个子域,得到的则是乘性 Schwarz,更新算子的顺序已经改变。
局部解相加与全局粗方向 还能写成一组子空间投影之和:
每个 都是到子域支撑空间的 -正交投影,因为 ,且 。因此 ,有 。这个粗界足以提醒我们:正定预条件器不等于步长为一的定常迭代一定收敛。做 时,还需 。
粗空间修补的是稳定分解
若每个向量都能分解为 ,其中 支撑在子域内,并满足
则 。证明使用 ,再在各子域空间的直和上应用 Cauchy–Schwarz 不等式公理库Cauchy–Schwarz 不等式Cauchy–Schwarz inequality · 柯西–施瓦茨不等式内积的绝对值不超过两向量范数之积,且等号精确刻画线性相关。:
又有 ,平方整理便得到下界。若一个缓慢变化的全局向量被硬切成局部片段,片段边缘可能产生较大能量, 因而变坏。
加入满列秩粗基矩阵 ,令
便允许分解额外使用一个跨域分量。它仍正定,但良好尺度无关的条件数还需要证明合适的稳定分解与子空间相互作用界,不能只凭“加了一层”宣布获得。
例子与边界
取九节点链 ,两个子域为 、。每个局部矩阵都是六阶链矩阵。
对 ,两份局部解都为 。延伸相加得到
例如全局节点 在第一子域收到 ,在第二子域也收到 ,合并为 。这份结果是预条件校正,不是宣称已经精确解出全局方程。
再取一个跨域帐篷方向
有 ,且 。令 ,粗校正为 ,所以 。
从误差 出发,以相同步长 做一次校正。单层结果为
两层加性结果为
能量平方分别为
粗方向在这个例子中明显改善一步效果。不过两层加性校正并不把 一步归零:虽然 ,局部项 也同时参加,不能只看粗项的精确性就忽略其余校正。步长一还可能过度修正;本例的三个投影给出 ,所以 是有依据的安全选择。
推论与应用
局部矩阵来自稀疏矩阵公理库稀疏矩阵表示与运算Sparse matrix只存非零项及其位置,以稀疏模式组织矩阵运算、图结构、重排序与填充成本。的索引限制,可以预先分解并复用。若每个子域又有内部与界面变量,可用静态凝聚公理库Schur 补与静态凝聚Schur complement · Static condensation把内部变量的精确响应折入界面矩阵与右端,推导 Schur 补、恢复公式及最小能量性质,并逐项凝聚五节点链。减少局部解的内部工作;这不改变加性组合要求所有子域使用同一输入残差。
一次应用成本是各局部求解成本之和,加上限制与重叠求和的索引工作;并行墙钟时间还受最大局部问题、通信和归约影响。增大重叠通常改善局部信息传播,却增加重复计算、存储与通信,应按总成本取舍。
只把局部解的一部分延伸回全局的 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:作者教材。加性应用、子空间投影与稳定分解下界。