Speed up (and test) sparse matrix-vector multiplication