diff --git a/runcomputations.jl b/runcomputations.jl index ee8fb1a..8be492b 100644 --- a/runcomputations.jl +++ b/runcomputations.jl @@ -53,11 +53,11 @@ let GROUP = ARGS[1] global_logger(logger) @info "" group_name - certified_λ, λ = check_propertyT(groups[GROUP], group_name, + spectral_gap, λ = check_propertyT(groups[GROUP], group_name, HALFRADIUS, Inf, AutomaticStructure, with_SCS(100_000, 50)) - if certified_λ < 0.0 && λ > 0.01 # there is still a chance to detect spectral gap - new_λ = round(0.8λ, sigdigits=2) + if spectral_gap < 0.0 && λ > 0.01 # there is still a chance to detect spectral gap + new_λ = round(0.8λ, sigdigits=3) new_dir = joinpath(group_name, "$new_λ") isdir(new_dir) || mkpath(new_dir) cp(joinpath(group_name, "Inf", "warmstart.jld"), joinpath(new_dir, "warmstart.jld"), force=true) diff --git a/src/utils.jl b/src/utils.jl index d5d17ee..e237598 100644 --- a/src/utils.jl +++ b/src/utils.jl @@ -38,6 +38,7 @@ function check_propertyT(sett::PropertyT.Settings) S = gens(RG.group) S = unique!([S; inv.(S)]) + @assert length(S) == 6 Δ = RG(length(S)) - sum(RG(s) for s in S) # small Laplacian @info "Checking the positivity of" Δ @@ -61,7 +62,7 @@ function check_propertyT(sett::PropertyT.Settings) PropertyT.interpret_results(sett, certified_λ/100) - return certified_λ/100, λ/100 + return certified_λ/100, λ end function load_basis!(RG::GroupRing, sett::PropertyT.Settings)