fix rref for complex matrices (issue JuliaLang#1869)

stevengj · Jan 2, 2013 · c30a381 · c30a381
1 parent 464f407
commit c30a381
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/base/linalg_dense.jl b/base/linalg_dense.jl
@@ -203,7 +203,7 @@ end
 
 function rref{T}(A::Matrix{T})
     nr, nc = size(A)
-    U = copy_to(similar(A,Float64), A)
+    U = copy_to(similar(A, T <: Complex ? Complex128 : Float64), A)
     e = eps(norm(U,Inf))
     i = j = 1
     while i <= nr && j <= nc