Hi, when I read reference manual of mkl, the mkl_?csrgemv routine calculates y=Ax (or y=A'x) where A is m-by-m matrix in 3-array CSR format.
However, when I tested, I figured out that this routine also works well in m-by-k matrix (both k<m and k>m).
Also, proper results are obtained even when I set ja and a not in order.
ex) for 2-by-2 matrix, [(1, 2), (3, 4)], I set ja(1)=2, a(1)=2, ja(2)=1, a(2)=1, ...
Does this routine updated in recent version?