Make realdiag! work for non-square matrices

mohamed82008 · haampie · commit c078e9977b34 · 2018-07-11T22:59:50.000+03:00
diff --git a/src/lobpcg.jl b/src/lobpcg.jl
@@ -288,7 +288,7 @@ end
 
 realdiag!(M) = nothing
 function realdiag!(M::AbstractMatrix{TC}) where TC <: Complex
-    @inbounds for i in 1:size(M, 1)
+    @inbounds for i in 1:minimum(size(M))
         M[i,i] = real(M[i,i])
     end
     return M
