GPU-friendly SVD + correct gaugefix

Author: lkdvos

src/implementations/svd.jl

@@ -255,11 +255,8 @@ function svd_full!(A::AbstractMatrix, USVᴴ, alg::DiagonalAlgorithm)
     T = eltype(Vᴴ)
     zero!(U)
     zero!(Vᴴ)
-    @inbounds for (i, pi) in enumerate(p)
-        s = Ad[pi]
-        U[pi, i] = sign_safe(s)
-        Vᴴ[i, pi] = one(T)
-    end
+    U[CartesianIndex.(enumerate(p))] .= Ref(one(T))
+    Vᴴ[CartesianIndex.(reverse.(enumerate(p)))] .= sign_safe.(view(Ad, p))
     return U, S, Vᴴ
 end
 function svd_compact!(A::AbstractMatrix, USVᴴ, alg::DiagonalAlgorithm)