If we have a vector (u), we can easily convert it into a sparse diagonal matrix using the command spdiagm(u).
I want to do the same operation using CUDA where u is defined as CuArray(u). The command spdiagm(u) gives the scalar indexing error and therefore it does not work.
I tried using CuSparseMatrixCSC(Diagonal(u)), it works but it is very slow. Is there a way to perform this task efficiently?
I am relatively new to CUDA so any help is appreciated.