diff --git a/Orb_AutFN.jl b/Orb_AutFN.jl
index 5ab2bb5..e14bad2 100644
--- a/Orb_AutFN.jl
+++ b/Orb_AutFN.jl
@@ -130,7 +130,9 @@ end
 
 function reconstructP(m::JuMP.Model, data::OrbitData)
     computedPs = [getvalue(P) for P in data.Ps]
-    return sum(data.dims[π]*data.Us[π]*computedPs[π]*data.Us[π]' for π in 1:endof(data.Ps))
+    return sum(
+      data.dims[π]*transform(data.Us[π]',computedPs[π], sparse=false)
+      for π in 1:endof(data.Ps))
 end
 
 function create_SDP_problem(name::String; upper_bound=Inf)
diff --git a/Projections.jl b/Projections.jl
index 2a8d946..4374b35 100644
--- a/Projections.jl
+++ b/Projections.jl
@@ -125,8 +125,8 @@ function rankOne_projections(BN::WreathProducts.WreathProduct, T::Type=Rational{
    for i in 1:N-1
         Sk_first = [RBN(p, g->incl(0,g)) for p in SNprojs_nc[i]]
         Sk_last  = [RBN(p, g->incl(i,g)) for p in SNprojs_nc[N-i]]
-        append!(all_projs, [Qs[i+1]*p1*p2
-                for (p1,p2) in Base.product(Sk_first,Sk_last)])
+        append!(all_projs,
+            [Qs[i+1]*p1*p2 for (p1,p2) in Base.product(Sk_first,Sk_last)])
    end
 
    append!(all_projs, [Qs[N+1]*RBN(p, g-> incl(0,g)) for p in SNprojs_nc[N]])