[论文解读] On the Equivalence of Primal and Dual Substructuring Preconditioners

Bedřich Sousedík, Jan Mandel|ArXiv.org|Feb 29, 2008

Matrix Theory and Algorithms参考文献 32被引用 13

一句话总结

本文建立了原始与对偶子结构化预条件子之间的数学等价性，证明了 P-FETI-DP 与 BDDC 为同一方法，并表明 BDD 与 FETI-1 的特定变体具有相同的特征值。该工作在统一框架下整合了关键的域分解方法，展示了其谱等价性，并为理解其收敛行为提供了新见解。

ABSTRACT

After a short historical review, we present four popular substructuring methods: FETI-1, BDD, FETI-DP, BDDC, and derive the primal versions to the two FETI methods, called P-FETI-1 and P-FETI-DP, as proposed by Fragakis and Papadrakakis. The formulation of the BDDC method shows that it is the same as P-FETI-DP and the same as a preconditioner introduced by Cros. We prove the equality of eigenvalues of a particular case of the FETI-1 method and of the BDD method by applying a recent abstract result by Fragakis.

研究动机与目标

• 建立原始与对偶子结构化预条件子之间的理论等价性，特别是 P-FETI-DP 与 BDDC 之间的等价性。
• 通过证明在特定变体下 FETI-1 与 BDD 方法的特征值谱完全相同，解决长期存在的关于二者关系的疑问。
• 通过推导 FETI 方法的原始形式并将其置于统一的数学框架中，统一分析域分解方法。
• 利用 Fragakis 的抽象结果，严格证明 BDD 与 FETI-1 的一个变体之间的谱等价性。
• 澄清 FETI、BDD、BDDC 及其原始对应方法之间的联系，深化对它们收敛性与预条件行为的理解。

提出的方法

• 通过从对偶形式到原始形式的系统性变换，推导出 FETI-1 和 FETI-DP 的原始形式，分别命名为 P-FETI-1 和 P-FETI-DP。
• 应用 Fragakis 提出的抽象谱等价性结果，证明 BDD 与 FETI-1 的特定变体具有相同的特征值。
• 利用涉及舒尔补、拉格朗日乘子和投影算子的矩阵变换，建立不同方法预条件算子之间的关系。
• 引入算子 $ T_D $ 和 $ T_P $，在 BDD 与 FETI-1 预条件系统之间映射特征向量，证明非平凡特征值的一一对应关系。
• 使用恒等式 $ B_D S \widetilde{S}^{+} B^T $ 作为投影算子，分析特征值重数，表明仅当映射非单射时，可能的特征值为 0 和 1。
• 证明预条件算子 $ M_{BDD} \widehat{S} $ 与 $ M_{FETI-1} \mathcal{F} $ 的谱相同，仅在 0 和 1 的特征值上存在差异，且对于 $ \lambda \neq 0,1 $，其重数完全一致。

实验结果

研究问题

• RQ1 FETI 方法的原始与对偶形式在谱性质上是否等价？
• RQ2 P-FETI-DP 方法在数学上是否等价于 BDDC 方法？
• RQ3 在 FETI-1 的特定变体下，BDD 与 FETI-1 方法是否具有相同的特征值谱？
• RQ4 能否将抽象谱等价性结果具体应用于域分解方法，以证明特征值的恒等性？
• RQ5 BDD 与 FETI-1 预条件算子的特征空间之间存在何种关系？

主要发现

• P-FETI-DP 方法在数学上等价于 BDDC 方法，确认二者为同一算法。
• Cros 提出的预条件子与 P-FETI-DP 和 BDDC 均等价，建立了现有方法之间的新联系。
• BDD 与 FETI-1 的特定变体的特征值完全相同，仅排除特征值 0 和 1 的比较。
• 对于每个共同的特征值 $ \lambda \neq 0,1 $，BDD 与 FETI-1 预条件算子的重数完全一致。
• BDD 与 FETI-1 的非平凡特征值对应的特征向量之间存在一一对应关系，从而证明了谱等价性。
• 该证明依赖于算子 $ B_D S \widetilde{S}^{+} B^T $ 是投影算子的性质，该性质限制了当映射非单射时，可能的特征值仅为 0 和 1。