From 2e2add905b7b8860e8fa76e7fdf854de3620e93f Mon Sep 17 00:00:00 2001 From: kalmarek Date: Wed, 5 Dec 2018 10:25:43 +0100 Subject: [PATCH] save and use only Q while checking --- paper_data/check_positivity.jl | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) 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)