GPUArrays.jl
GPUArrays.jl copied to clipboard
Linear solvers
A_ldiv_B!
and \
support for GPUArrays is pretty essential and I keep assuming it already exists.
Would it be appropriate to add some overloads for Diagonal{AbstractGPUArray}
analogous to those for cholesky
I added here?
Can you implement them using GPUArrays' abstractions, or a GPUArrays kernel? If so, then yes this would be a good fit here.
Actually, for diagonal matrices we can just use broadcasting syntax. I'll open a PR.