I have a serial implementation, an openMP implementation, and mkl_dcsrmv. I'm computing the average and minimum time over a number of runs, say, 10. Strangely, though, the openMP implementation beats MKL always. For the matrix sizes in the code, openMP has a min time of...
I think I can perform part of this operation using mkl_dcsrmv, which computes y := alpha*A*x + beta*y where x and y are the values array of the identity matrix, alpha is -1.0 and beta is 1.0. According to the manual, I can operate on the diagonal by...
Bug in 64-bit mkl_dcsrmv Subscribe More actions Lars_K_1 Beginner 05-16-2016 09:57 AM 513 Views Hello, I need sparse matrix-vector multiplication in my code and want to use the MKL, as I want to also use the Xeon Phi cards in my system. However, the 64-bit version (...