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
- matrix version of ?spmv

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

may_ka

Beginner

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

11-05-2017
05:05 AM

52 Views

matrix version of ?spmv

Hi there,

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.

Any suggestions??

Thanks.

Link Copied

2 Replies

Gennady_F_Intel

Moderator

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

11-05-2017
06:24 AM

52 Views

may_ka

Beginner

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

11-05-2017
12:14 PM

52 Views

Hi Gennady,

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:

- the packed format used for the matrix in ?spmv is not equal to any sparse format because it lacks the integer index vectors.
- a dense symmetric matrix does not need any index vectors at all because the full upper or lower triangular is stored as a vector
- BLAS does not have a level 3 equivalent to ?spmv https://software.intel.com/en-us/mkl-developer-reference-fortran-blas-level-3-routines
- a symmetric dense matrix which can be stored packed but is stored sparse will be less efficient in terms of RAM because it has an unnesseary overhead for the index vectors. In fact it will have no storage advantage at all compared to storing the full symmetric matrix. It will also be less efficient in terms of speed during multiplication from going through these index vectors.

Let me know what I am getting wrong.

- 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.