title: What sparse solver supports diagonal storage format
body: <p>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. </p>
  
  <p>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 =(</p>
  
  <p>Thank you.</p>