具有对称结构的稀疏非对称矩阵的迭代求解器的选择

计算科学 线性代数 线性求解器 稀疏矩阵 对称
2021-12-06 12:17:02

我有一个稀疏nxn矩阵 A 的结构非常有趣。它具有对称结构但不对称块的块结构。以数学方式表示块Ajk=AkjAjkAjkT,所以我的矩阵实际上不是对称的。有没有可以利用这种系统的算法?

1个回答

首先 - 你的数学描述似乎是错误的:Ajk=Akj, k,j{1,,n}意味着矩阵是对称的。所以你的意思并不完全清楚。还,AjkT精确定义为AjkT=Akj(参见例如维基百科)。这会以某种方式与您的方程式发生冲突。

至于实际问题:虽然我不确定您的结构实际上是什么 - 我也不知道任何可能利用类似于您的结构的方法。对于迭代方法的一个很好的概述,我推荐模板书但是最“结构利用”是关于对称性的。我希望这可以帮助一些..