Another GPU update

Author: kshyatt

test/testsuite/ad_utils.jl

@@ -41,7 +41,8 @@ function stabilize_eigvals!(D::AbstractVector)
     # rescale eigenvalues so that they lie on distinct radii in the complex plane
     # that are chosen randomly in non-overlapping intervals [k/n, (k+0.5)/n)] for k=1,...,n
     radii = ((1:n) .+ rand(real(eltype(D)), n) ./ 2) ./ n
-    D .= sign.(D) .* radii[p]
+    hD .= sign.(collect(D)) .* radii[p]
+    copyto!(D, hD)
     return D
 end
 function make_eig_matrix(T, sz)