title: Fortran 90/95 library for sparse matrices?
body: <p>I am looking for a library for dealing with sparse matrices in fortran 90/95. I only need very basic operations like matrix-vector multiplication. What do you suggest I use?</p>
  
  <p>I have searched around and an extension(?) to BLAS called "sparse blas", documented in chapter 3 of the blast technical forum specification:</p>
  
  <p><a href="http://www.netlib.org/blas/blast-forum/" rel="noreferrer">http://www.netlib.org/blas/blast-forum/</a>,</p>
  
  <p>seems ideal. According to that document, there should be a fortran 95 interface to the library. However, I haven't been able to find an actual implementation of this specification anywhere that I was able to download. I am kind of getting the impression that the fortran 95 sparse blas specification isn't really implemented anywhere? Not as open source anyway.</p>
  
  <p>I also found a library called sparsekit:</p>
  
  <p><a href="http://people.sc.fsu.edu/~jburkardt/f_src/sparsekit/sparsekit.html" rel="noreferrer">http://people.sc.fsu.edu/~jburkardt/f_src/sparsekit/sparsekit.html</a></p>
  
  <p>Anyone got any experience with either one of these two, or any other sparse matrix library for fortran 90/95? I would like advice on which one to use, and also where to get it (in the case of sparsekit the source code is available through the above link). Is there anything like a "standard" for sparse matrices, like BLAS is for dense ones?</p>
  
  <p>Cheers,</p>
  
  <p>Arne</p>