From 66f860bac10bc5887cf0efa04cab7934d0bd0e53 Mon Sep 17 00:00:00 2001
From: kalmar <kalmar@amu.edu.pl>
Date: Wed, 26 Jul 2017 12:15:31 +0200
Subject: [PATCH] sparsify both on the forward and backward transform

---
 src/Orbit-wise.jl         | 2 +-
 src/OrbitDecomposition.jl | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/Orbit-wise.jl b/src/Orbit-wise.jl
index 3590b22..bf85296 100644
--- a/src/Orbit-wise.jl
+++ b/src/Orbit-wise.jl
@@ -83,7 +83,7 @@ function init_orbit_data(logger, sett::Settings; radius=2)
    return 0
 end
 
-function transform(U::AbstractArray, V::AbstractArray; sparse=false)
+function transform(U::AbstractArray, V::AbstractArray; sparse=true)
    if sparse
       return sparsify(U'*V*U)
    else
diff --git a/src/OrbitDecomposition.jl b/src/OrbitDecomposition.jl
index d2963b9..30b7ec3 100644
--- a/src/OrbitDecomposition.jl
+++ b/src/OrbitDecomposition.jl
@@ -123,7 +123,7 @@ 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 .+= dims[π].* (A * Us[π]*Ps[π]*Ust[π] * B)
+            recP .+= sparsify(dims[π].* (A * Us[π]*Ps[π]*Ust[π] * B))
         end
     end
     recP .*= 1/length(keys(mreps))