Community
cancel
Showing results for 
Search instead for 
Did you mean: 
damienhocking
New Contributor I
77 Views

gemv for sparse, rectangular matrices?

Jump to solution

Is there a recommended routine in the MKL that can do a GEMV with the transpose of a rectangular sparse matrix?  I know I can pad it to square but I'd prefer not to have to do that.  These sparse matrices could be large, 10 million+ nonzeroes and a million rows or columns. 

Damien

0 Kudos
1 Solution
Zhen_Z_Intel
Employee
77 Views

Hi,

mkl_?<format>gemv is designed for square, but mkl_?<format>mv can be used for rectangular matrix (m * k). Or, you could use Inspector-execute sparse function mkl_sparse_?_mv, normally would provide a better performance. 

View solution in original post

4 Replies
Zhen_Z_Intel
Employee
77 Views

Hi,

You provided BLAS level 2 function for sparse matrix. For sparse matrix, we normally use compressed data structure that Intel MKL supports CSR, CSC, COO, diagonal and skyline. You could learn more from developer reference about usage with different matrix structure& data compression format:

https://software.intel.com/en-us/node/520802

damienhocking
New Contributor I
77 Views

Yes, but the docs are all for a square m x m matrix.  I have rectangular matrices.  

Zhen_Z_Intel
Employee
78 Views

Hi,

mkl_?<format>gemv is designed for square, but mkl_?<format>mv can be used for rectangular matrix (m * k). Or, you could use Inspector-execute sparse function mkl_sparse_?_mv, normally would provide a better performance. 

View solution in original post

damienhocking
New Contributor I
77 Views

Thank you Fiona, that's exactly what I wanted.

Reply