Я смотрю на решение проблемы, которая представляет собой PDE, и трехмерный дискретизированный домен может иметь разные граничные условия на каждой из 6 границ (или все одинаковые).
Как лучше всегопоместить эту разреженную матрицу в сжатый формат?CSR будет моим единственным выбором здесь?Я думал об использовании ellpack, но я не уверен, как это работает для изменения граничных условий.
Думая о представлении 2D-матрицы в трехмерном пространстве ... оно будет главным образом по диагонали с 7 диагоналями, но этидиагонали могут измениться вдоль границ.Не похоже, что я могу использовать формат, в котором я храню значения и смещение от диагонали, которое каждый раз одно и то же.
Очевидно, я пытаюсь установить эту проблему, чтобы она была более дружественной к кэшу для моей компьютерной графики.решатель, который много умножает вектор-матрицу