kalmar/GroupsWithPropertyT
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

add normalizing 36 in comments

Browse Source
This commit is contained in:
kalmarek 2018-12-04 04:50:24 +01:00
parent c4dafc635e
commit 6f82292b52
1 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
paper_data/check_positivity.jl
View File

@ -30,7 +30,7 @@ if K == nothing || LAMBDA == nothing
throw(invalid_use_message) throw(invalid_use_message)
end end
info("Running checks for Adj₅ + $K·Op₅ - $LAMBDA·Δ₅") info("Running checks for 36(Adj₅ + $K·Op₅) - $LAMBDA·Δ₅")
G = AutGroup(FreeGroup(5), special=true) G = AutGroup(FreeGroup(5), special=true)
S = generating_set(G) S = generating_set(G)
@ -115,7 +115,7 @@ if !isfile(SOLUTION_FILE)
save(SOLUTION_FILE, "λ", λ, "P", P) save(SOLUTION_FILE, "λ", λ, "P", P)
end end
info("Checking the sum of squares solution for Adj₅ + $K·Op₅ - $LAMBDA·Δ₅") info("Checking the sum of squares solution for 36(Adj₅ + $K·Op₅) - $LAMBDA·Δ₅")
P, λ = load(SOLUTION_FILE, "P", "λ") P, λ = load(SOLUTION_FILE, "P", "λ")
info("Computing Q = √P") info("Computing Q = √P")
@ -139,4 +139,7 @@ Q_int = PropertyT.augIdproj(Q);
b_int = SOS_residual(EOI_int, Q_int) b_int = SOS_residual(EOI_int, Q_int)
@show norm(b_int, 1); @show norm(b_int, 1);
info("λ is certified to be > ", (@interval(λ) - 2^2*norm(b_int,1))).lo λ_cert = @interval(λ) - 2^2*norm(b_int,1)
info("λ is certified to be > ", λ_cert.lo)
info("i.e Adj₅ + $K·Op₅ - $((λ_cert/36).lo)·Δ₅ ∈ Σ²₂ ISAut(F₅)")