Remove buffer in gemv

Author: michel2323

lib/mkl/wrappers_sparse.jl

@@ -158,13 +158,14 @@ for (fname, elty) in ((:onemklSsparse_gemv, :Float32),
                     #   conj(y_new) = conj(alpha) * (A^T) * conj(x) + conj(beta) * conj(y)
                     # Since S=A^T and op='N' computes S*x = A^T*x, we can realize this with one call.
                     y .= conj.(y)
-                    x_conj = similar(x)
-                    x_conj .= conj.(x)
+                    x .= conj.(x)
 
                     queue = global_queue(context(x), device())
-                    $fname(sycl_queue(queue), 'N', conj(alpha), A.handle, x_conj, conj(beta), y)
+                    $fname(sycl_queue(queue), 'N', conj(alpha), A.handle, x, conj(beta), y)
 
                     y .= conj.(y)
+                    # Restore x
+                    x .= conj.(x)
                     return y
                 else
                     # real eltype: A^H == A^T