topic Re: Most efficient storage format for sparse blas matrix-matrix in Intel® oneAPI Math Kernel Library

Most efficient storage format for sparse blas matrix-matrix multiply

soliton — Thu, 11 Sep 2008 17:09:23 GMT

Dear All,

I am starting to write a new application that will have to do sparse*dense matrix multiplies. Before just coding all the alternatives and benchmarking them I'd like to know whether it is know if one of the alternative storage formats is more efficient than the others. I am just guessing here, but it seems to me that CSR or CSC might help with data locality more that the coordinate format since in this case the ordering of the matrix can be basically random. Any clues?

Cheers,

Cristian

Re: Most efficient storage format for sparse blas matrix-matrix

Sergey_K_Intel1 — Mon, 15 Sep 2008 10:45:05 GMT

Dear Christian,

First of all,I'd recommend to take a lookat the structure of matrices. For example if it is known that a matrix consists of a few diagonals, the diagonal format is the best storage scheme to gain efficiency. If your matrix consists of many square dense blocks, please use the block sparse row.

If the matrix is not regularly structured, MKL Sparse BLAS routines for the CSR and CSC formats perform significantlybetter than for the coordinate format. If you are going to do non-transposed sparse-dense matrix multiplies for general matrices like A x, Id recommend to use the CSR format. The CSC format is preferable for transposed case A^T x. Id alsolike toremind that the routines for the CSR and CSC formats are OpenMP threaded.

All the best

Sergey

Re: Most efficient storage format for sparse blas matrix-matrix

soliton — Tue, 30 Sep 2008 00:03:17 GMT

Dear Sergey,

Many thanks for your swift reply.

Cheers,

Cristian

Quoting - SERGEY KUZNETSOV (Intel)

Dear Christian,

First of all,I'd recommend to take a lookat the structure of matrices. For example if it is known that a matrix consists of a few diagonals, the diagonal format is the best storage scheme to gain efficiency. If your matrix consists of many square dense blocks, please use the block sparse row.

If the matrix is not regularly structured, MKL Sparse BLAS routines for the CSR and CSC formats perform significantlybetter than for the coordinate format. If you are going to do non-transposed sparse-dense matrix multiplies for general matrices like A x, Id recommend to use the CSR format. The CSC format is preferable for transposed case A^T x. Id alsolike toremind that the routines for the CSR and CSC formats are OpenMP threaded.

All the best

Sergey