diff --git a/OrbitDecomposition.jl b/OrbitDecomposition.jl index c9265d0..073451a 100644 --- a/OrbitDecomposition.jl +++ b/OrbitDecomposition.jl @@ -200,10 +200,8 @@ end function orthSVD(M::AbstractMatrix) M = full(M) -# matrixRank = rank(M) fact = svdfact(M) sings = fact[:S] -# @show sings[sings.>1e-14] M_rank = sum(fact[:S] .> maximum(size(M))*eps(eltype(fact[:S]))) Ufactor = fact[:U] return Ufactor[:,1:M_rank]