kalmar/SmallHyperbolic
1
0
Fork 0
mirror of https://github.com/kalmarek/SmallHyperbolic synced 2024-11-23 23:40:28 +01:00
Code Issues Releases Wiki Activity

to aid the numerical solver check the positivity of 100Δ²

Browse Source
This commit is contained in:
kalmarek 2020-01-20 03:18:51 +01:00
parent f984e91dd6
commit d0c2fcdb83
2 changed files with 9 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
runcomputations.jl
View File

@ -50,11 +50,12 @@ let GROUP = ARGS[1]
global_logger(logger) global_logger(logger)
@info "" group_name @info "" group_name
check_propertyT(groups[GROUP], group_name, λ = check_propertyT(groups[GROUP], group_name,
HALFRADIUS, Inf, AutomaticStructure, with_SCS(100_000, 50)) HALFRADIUS, Inf, AutomaticStructure, with_SCS(100_000, 50))
check_propertyT(groups[GROUP], group_name, if λ > 1e-2 # there is still a chance to detect spectral gap
HALFRADIUS, Inf, AutomaticStructure, with_SCS(1_000_000, 0)) check_propertyT(groups[GROUP], group_name,
HALFRADIUS, Inf, AutomaticStructure, with_SCS(500_000, 0))
end
end end
true # to keep make happy
end end

7
src/utils.jl
View File

@ -44,8 +44,8 @@ function check_propertyT(sett::PropertyT.Settings)
@assert iszero(aug(Δ)) @assert iszero(aug(Δ))
ELT = Δ^2; ELT = 100*(Δ^2);
ELT_NAME = "Δ²" ELT_NAME = "100Δ²"
λ, P = PropertyT.approximate_by_SOS(sett, ELT, Δ, λ, P = PropertyT.approximate_by_SOS(sett, ELT, Δ,
solverlog=PropertyT.filename(sett, :solverlog)) 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) 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 end
function load_basis!(RG::GroupRing, sett::PropertyT.Settings) function load_basis!(RG::GroupRing, sett::PropertyT.Settings)