From 3ff649683a031f5cfab525acdf2f83b10e9e80cd Mon Sep 17 00:00:00 2001 From: kalmar Date: Wed, 26 Jul 2017 12:12:58 +0200 Subject: [PATCH] FIXXX: clamp to 0 only those eps-close to 0, not < eps --- src/Orbit-wise.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/Orbit-wise.jl b/src/Orbit-wise.jl index afc37b8..3590b22 100644 --- a/src/Orbit-wise.jl +++ b/src/Orbit-wise.jl @@ -58,7 +58,7 @@ dens(M::AbstractArray) = sum(abs.(M) .!= 0)/length(M) function sparsify{T}(U::AbstractArray{T}, check=true) W = deepcopy(U) - W[W .< eps(T)] .= zero(T) + W[abs.(W) .< eps(T)] .= zero(T) if check && rank(W) != rank(U) info("Sparsification would decrease the rank!") W = U