Question with mkl_dcsrmultcsr

patrick_bao — Fri, 04 Nov 2011 03:51:20 GMT

Hi,

I have trouble using mkl_dcsrmultcsr when calculating A'A if A is not a square matrix. The problem is the size of array ia.
The doc says ia is an array of length m + 1 when trans = 'N' or 'n', or n + 1 otherwise. In my case, I was using 'T', so I need to allocate an array of length n + 1. How can a matrix having m rows, but rowIndices array only have length n? I'm really confusing.

Thanks,

Question with mkl_dcsrmultcsr

Sergey_K_Intel1 — Wed, 09 Nov 2011 13:13:55 GMT

Many thanks for catching this. It's a bug in documentation. The description of ia must be as follows ia

INTEGER. Array of length m + 1.

This array contains indices of elements in the array a, such that ia(I) is the index in the arraya of the first non-zero element from the row I. The value of the last element ia(m + 1) is equal to the number of non-zero elements of the matrixA plus one. Refer to rowIndex array description in Sparse Matrix Storage Formats for more details

but array ib should be described as array ia in the current version ib

INTEGER. Array of lengthn + 1 when trans = 'N' or 'n', orm + 1 otherwise.

This array contains indices of elements in the array b, such that ib(I) is the index in the arrayb of the first non-zero element from the row I. The value of the last element ib(m + 1) or ib(n + 1) is equal to the number of non-zero elements of the matrixB plus one. Refer to rowIndex array description in Sparse Matrix Storage Formats for more details.

Many thanks again All the best Sergey

topic Question with mkl_dcsrmultcsr in Intel® oneAPI Math Kernel Library

Question with mkl_dcsrmultcsr

Question with mkl_dcsrmultcsr