PDE である問題を解決しようとしています.3D 離散化ドメインは、6 つの境界のそれぞれに異なる境界条件を持つことができます (またはすべて同じです)。
このスパース行列を圧縮形式にする最良の方法は何ですか? ここでは、CSR が唯一の選択肢になるのでしょうか? ellpack を使用することも考えましたが、境界条件の変更に対してどのように機能するかわかりません。
3D 空間の 2D マトリックス表現について考えると、主に 7 つの対角線で対角線が支配的になりますが、これらの対角線は境界に沿って変化する可能性があります。値を保存する形式と、毎回同じ対角線からのオフセットを使用できるようには思えません。
明らかに、私はこの問題を設定して、多くのベクトル行列乗算を行っている CG ソルバーのキャッシュ フレンドリーにしようとしています。