If I understand correctly, operations such as `v = A*v`, where A is a dense matrix and v is a vector, create temporaries and allocate memory; if one wishes to do the computation in place without allocation, the only solution is using directly `gemv!`.
What if `A` is a sparse matrix? Is there an in-place alternative to `v = A*v`? Thanks, --federico
