Domain Decomposition计算公式备忘

Domain Decompostion 计算公式

示意图

推导

• 基本公式：$$\{\begin{array}{l}s={\mathbf{S}}^{-1}{F}_s-{\mathbf{S}}^{-1}{\mathbf{W}}_{R}h\\ h=(\mathbf{H}-{\mathbf{W}}_L{\mathbf{S}}^{-1}{\mathbf{W}}_R{)}^{-1}(F_h-{\mathbf{W}}_L{\mathbf{S}}^{-1}{F}_s)\end{array}$$对正定对称阵有$\mathbf{S}=\mathbf{L}{\mathbf{L}}^t$, ${\mathbf{W}}_L={\mathbf{W}}_R^t$, 因而上面公式可化为：
  $$\{\begin{array}{l}s={\mathbf{L}}^{-t}({\mathbf{L}}^{-1}{F}_s-{\mathbf{L}}^{-1}{\mathbf{W}}_{R}h)\\ h=(\mathbf{H}-({{\mathbf{L}}^{-1}\mathbf{W}}_R{)}^t({{\mathbf{L}}^{-1}\mathbf{W}}_R){)}^{-1}(F_h-({{\mathbf{L}}^{-1}\mathbf{W}}_R{)}^t({\mathbf{L}}^{-1}{F}_s))\end{array}$$所以仅需要计算 ${{\mathbf{L}}^{-1}\mathbf{W}}_R$ 和 ${\mathbf{L}}^{-1}{F}_s$。

计算步骤

• 并行分解$\mathbf{S}=\mathbf{L}{\mathbf{L}}^t$
• 并行求解${\mathbf{L}}^{-1}{F}_s$
• 并行求解${{\mathbf{L}}^{-1}\mathbf{W}}_R$, 同时更新$\mathbf{H}$矩阵。此处注意${\mathbf{W}}_R$的非零稀疏性，可以利用加速计算。更新$\mathbf{H}$矩阵需要考虑多个线程对$\mathbf{H}$相同块的更新，需要锁和同步机制
• 更新后的$\mathbf{H}$也是对称正定块状三对角阵，可以用块状分解计算LLT
• 先计算$h$,之后计算$s$