用于最小化天际线矩阵的 ram 使用的重新排序算法

计算科学 有限元 线性求解器 稀疏矩阵 内存管理
2021-12-19 03:33:44

线性方程组的刚度矩阵是以对称矩阵,与 3-D 弹性固体的有限元模型相关联,当使用反向 Cuthill–McKee 算法不会产生具有有限 RAM 存储空间要求的天际线矩阵,而二维问题会发生这种情况。Ax=BAn×n

3-D 案例的最佳算法是什么?

1个回答

对于 3D 中的直接解决方案,您可能应该使用某种嵌套解剖 (ND) 或最小度数 (MD)。这些直接攻击 A=LL' 因式分解的存储需求,而不是带宽(对填充只有间接影响)。就其本身而言,带宽减少不足以使 3D 直接求解易于处理。

可以尝试的好的 ND 代码是 METIS 或 Scotch,而 AMD 最容易访问 MD 代码。许多稀疏直接求解器将在其输入参数中直接公开/引用这些重新排序包。

对于为您隐藏/处理所有这些细节的良好“交钥匙”包,我推荐英特尔的 MKL PARDISO。尽管它是封闭源代码,但在相当宽松的许可条款下可以使用二进制文件。