Use mkl_sparse_d_mv() from oneAPI MKL for sparse matrix-vector multiplication. Ensure the matrix is in CSR format and handle descriptors correctly for accurate scalar product results.
mkl_sparse_d_mv()