diff --git a/paper_data/check_positivity.jl b/paper_data/check_positivity.jl
index 8c834a1..b812079 100644
--- a/paper_data/check_positivity.jl
+++ b/paper_data/check_positivity.jl
@@ -112,14 +112,14 @@ if !isfile(SOLUTION_FILE)
 
     info("Reconstructing P...")
     @time P = PropertyT.reconstruct(Ps, orbit_data);
-    save(SOLUTION_FILE, "λ", λ, "P", P)
+    info("Computing Q = √P")
+    @time const Q = real(sqrtm(P));
+    
+    save(SOLUTION_FILE, "λ", λ, "Q", Q)
 end
 
 info("Checking the sum of squares solution for 36(Adj₅ + $K·Op₅) - $LAMBDA·Δ₅") 
-P, λ = load(SOLUTION_FILE, "P", "λ")
-
-info("Computing Q = √P")
-@time const Q = real(sqrtm(P));
+Q, λ = load(SOLUTION_FILE, "Q", "λ")
 
 function SOS_residual(eoi::GroupRingElem, Q::Matrix)
     RG = parent(eoi)