Какой разреженный решатель поддерживает диагональный формат хранения

Я пишу программу метода конечных разностей, используя C. Матрица жесткости является симметричной и ленточной. Для его хранения я бы хотел использовать формат Sparse Diagonal Storage.

Может кто-нибудь сказать, пожалуйста, какие решатели могут использовать диагональный формат хранения? Intel MKL's Paradiso использует только CRS-формат, SparseLib не имеет, Spooles вроде тоже =(

Спасибо.


person Ivan Kush    schedule 20.03.2013    source источник
comment
Я предложил ответ здесь на ваш пост SciComp SE.   -  person Michael Grant    schedule 21.03.2013


Ответы (1)


Разреженные библиотеки концентрируются на разреженном хранилище общего назначения. Для ленточных матриц вы можете использовать любую из стандартных библиотек линейной алгебры, например, LAPACK.

person David Heffernan    schedule 20.03.2013