From d0c2fcdb830b6f14b0f649f89b8fe2418223c8ff Mon Sep 17 00:00:00 2001 From: kalmarek Date: Mon, 20 Jan 2020 03:18:51 +0100 Subject: [PATCH] =?UTF-8?q?to=20aid=20the=20numerical=20solver=20check=20t?= =?UTF-8?q?he=20positivity=20of=20100=CE=94=C2=B2?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- runcomputations.jl | 9 +++++---- src/utils.jl | 7 ++++--- 2 files changed, 9 insertions(+), 7 deletions(-) diff --git a/runcomputations.jl b/runcomputations.jl index d51c7a8..3fb154c 100644 --- a/runcomputations.jl +++ b/runcomputations.jl @@ -50,11 +50,12 @@ let GROUP = ARGS[1] global_logger(logger) @info "" group_name - check_propertyT(groups[GROUP], group_name, + λ = check_propertyT(groups[GROUP], group_name, HALFRADIUS, Inf, AutomaticStructure, with_SCS(100_000, 50)) - check_propertyT(groups[GROUP], group_name, - HALFRADIUS, Inf, AutomaticStructure, with_SCS(1_000_000, 0)) + if λ > 1e-2 # there is still a chance to detect spectral gap + check_propertyT(groups[GROUP], group_name, + HALFRADIUS, Inf, AutomaticStructure, with_SCS(500_000, 0)) + end end - true # to keep make happy end diff --git a/src/utils.jl b/src/utils.jl index e9bdc57..ed3753c 100644 --- a/src/utils.jl +++ b/src/utils.jl @@ -44,8 +44,8 @@ function check_propertyT(sett::PropertyT.Settings) @assert iszero(aug(Δ)) - ELT = Δ^2; - ELT_NAME = "Δ²" + ELT = 100*(Δ^2); + ELT_NAME = "100Δ²" λ, P = PropertyT.approximate_by_SOS(sett, ELT, Δ, solverlog=PropertyT.filename(sett, :solverlog)) @@ -59,7 +59,8 @@ function check_propertyT(sett::PropertyT.Settings) certified_λ = PropertyT.certify_SOS_decomposition(ELT, Δ, λ, Q, R=sett.halfradius) - return PropertyT.interpret_results(sett, certified_λ) + PropertyT.interpret_results(sett, certified_λ/100) + return certified_λ/100 end function load_basis!(RG::GroupRing, sett::PropertyT.Settings)