当边界条件发生变化时，如何以稀疏格式存储离散化的3D域（用于求解PDE）？

Question

我正在寻找解决PDE问题的问题，并且3D离散化域可以在6个边界中的每一个上具有不同的边界条件（或者全部相同）。

将此稀疏矩阵放入压缩格式的最佳方法是什么？ CSR会成为我唯一的选择吗？我想过使用ellpack，但我不确定它是如何改变边界条件的。

考虑三维空间的二维矩阵表示......它主要是对角线占优势的​​7对角线，但这些对角线可能会沿着边界发生变化。似乎我不能使用我存储值的格式，以及每次都与对角线相同的偏移量。

显然我正试图将这个问题设置为对我的CG求解器更加缓存友好，这对于做很多向量矩阵乘法

Answer 1

坐标减少行格式..有效。