I am wondering whether there is a matrix version of ?spmv, e.g. ?spmm. I am dealing with huge symmetric matrices (+50,000 row dimension) which I want to multiply with other matrices. Currently I store them in general format and call gemm. I am aware of calling ?spmv inside a loop but I am looking for a more elegant way. I found gemm_compute in the leatest manual, but it occurs to me that it is not doing what I want.
thanks for the quick response.
However, you have pointed me to the sparse routines, but non of the matrices in my example is actually sparse.
My current understanding is:
Let me know what I am getting wrong.