Domain Decomposition计算公式备忘

东北豆子哥

已于 2022-06-17 07:31:20 修改

阅读量506

点赞数

分类专栏：数值计算/数值优化文章标签： latex 算法

于 2020-09-25 08:31:13 首次发布

本文链接：https://blog.csdn.net/weixin_42849849/article/details/108788578

版权

138 篇文章

订阅专栏

本文介绍Domain Decomposition方法的并行计算流程与核心计算公式。该方法通过将问题分解到不同子域中，利用并行计算提高效率。文中详细介绍了如何进行并行分解、并行求解关键矩阵运算，并更新H矩阵。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

在这里插入图片描述

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

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