Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- gemv for sparse, rectangular matrices?

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

damienhocking

New Contributor I

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-08-2017
09:13 PM

77 Views

gemv for sparse, rectangular matrices?

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

1 Solution

Zhen_Z_Intel

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-09-2017
12:14 AM

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.

Link Copied

4 Replies

Zhen_Z_Intel

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-08-2017
11:21 PM

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:

damienhocking

New Contributor I

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-08-2017
11:35 PM

77 Views

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

Zhen_Z_Intel

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-09-2017
12:14 AM

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.

damienhocking

New Contributor I

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-09-2017
08:20 AM

77 Views

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

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

For more complete information about compiler optimizations, see our Optimization Notice.