我正在研究解决一个 PDE 问题,并且 3D 离散域可以在 6 个边界中的每一个上具有不同的边界条件(或全部相同)。
将此稀疏矩阵转换为压缩格式的最佳方法是什么?企业社会责任会是我在这里唯一的选择吗?我考虑过使用 ellpack,但我不确定它是如何改变边界条件的。
考虑 3d 空间的 2D 矩阵表示......它将主要以对角线为主,有 7 个对角线,但这些对角线可能会沿着边界发生变化。似乎我不能使用存储值的格式,以及每次都相同的对角线偏移量。
显然,我试图将这个问题设置为对我的 CG 求解器更加缓存友好,它正在做很多向量矩阵乘法