0

I'm writing finite-difference method program using C. The stiffness matrix is symmetrical and band. For its storage I'd like to use Sparse Diagonal Storage format.

Could someone tell please, what solvers can use diagonal storage format? Intel MKL's Paradiso uses only CRS-format, SparseLib doesn't have, Spooles seems also =(

Thank you.

Ivan Kush
  • 2,958
  • 2
  • 23
  • 39

1 Answers1

0

The sparse libraries concentrate on general purpose sparse storage. For banded matrices you can use any of the standard linear algebra libraries, for example LAPACK.

David Heffernan
  • 601,492
  • 42
  • 1,072
  • 1,490