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))