From f6674337eda7865b7a7a48e57616c4c2bea81956 Mon Sep 17 00:00:00 2001
From: kalmar <kalmar@amu.edu.pl>
Date: Fri, 28 Jul 2017 12:44:25 +0200
Subject: [PATCH] sparsify only once during the reconstruction

---
 src/OrbitDecomposition.jl | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/OrbitDecomposition.jl b/src/OrbitDecomposition.jl
index fed0591..28c9c95 100644
--- a/src/OrbitDecomposition.jl
+++ b/src/OrbitDecomposition.jl
@@ -123,10 +123,11 @@ function reconstruct_sol{T<:GroupElem, S<:AbstractArray}(mreps::Dict{T, S},
     for g in keys(mreps)
         A, B = mreps[g], mreps[inv(g)]
         for π in 1:length(Us)
-            recP .+= sparsify(dims[π].* (A * Us[π]*Ps[π]*Ust[π] * B))
+            recP .+= dims[π].* (A * Us[π]*Ps[π]*Ust[π] * B)
         end
     end
     recP .*= 1/length(keys(mreps))
+    recP[abs.(recP) .< eps(eltype(recP))] = zero(eltype(recP))
     return recP
 end