title: Need a sparse matrix library
body: <p>Is there a sparse matrix library that copies the functionality of dense BLAS?</p>
  
  <p>I'd want at least:</p>
  
  <ol>
  <li>efficient SYR and SYRK (rank-k update) with sparse input (and possibly dense output),</li>
  <li>option for sparse symmetric storage,</li>
  <li>DOT, NRM2, transpose and other standard goodies.</li>
  </ol>
  
  <p>The naming and API doesn't need to be BLAS-centric though, as long as it does the job. C or Fortran preferred, but C++ is also ok.</p>
  
  <p>I looked at <a href="http://software.intel.com/en-us/articles/intel-mkl/" rel="nofollow" title="MKL">MKL</a>, <a href="http://docs.scipy.org/doc/scipy/reference/sparse.html" rel="nofollow">scipy.sparse</a>, <a href="http://www.boost.org/doc/libs/1_47_0/libs/numeric/ublas/doc/index.htm" rel="nofollow">uBLAS</a> and <a href="http://eigen.tuxfamily.org/index.php?title=Main_Page" rel="nofollow">Eigen3</a>, but they either didn't support what I need or I couldn't figure out whether they do.</p>