While using csrcoo to convert an unsorted (rows and colums) complex16 (double precision complex) matrix from COO to CSR using the option to sort the columns in the output, I discovered the array of matrix values gets garbled during the sort. The same operation appears to work if the input is already sorted, or using any other data type.
I'm guessing there's something wrong in how it handles swapping/gathering the data in the value arrray. Can anyone confirm?
Would you please provide a simple test case, including a sparse matrix and a code snippet showing how you call mkl_zcsrcoo?
And I also need information on the version of MKL and your OS.
Thanks.
This does look like a bug. This issue is now escalated to the MKL engineering team. Thanks for catching it!
For more complete information about compiler optimizations, see our Optimization Notice.