From 1468986b4a46c3f25f985d1bc8bd39aa255d399c Mon Sep 17 00:00:00 2001 From: kalmar Date: Thu, 8 Jun 2017 20:07:34 +0200 Subject: [PATCH] remove Group specific code --- Orb_AutFN.jl | 43 ++++--------------------------------------- 1 file changed, 4 insertions(+), 39 deletions(-) diff --git a/Orb_AutFN.jl b/Orb_AutFN.jl index 91ed676..774bb2d 100644 --- a/Orb_AutFN.jl +++ b/Orb_AutFN.jl @@ -181,18 +181,6 @@ function λandP(m::JuMP.Model, data::ProblemData) return λ, recP end -function SOutFN_generating_set(logger, name::String, N::Int) - - SOutFN = AutGroup(FreeGroup(N), special=true, outer=true) - info(logger, SOutFN) - S = generators(SOutFN); - S = [S; [inv(s) for s in S]] - info(logger, "Symmetric generating set of size $(length(S))") - info(logger, S) - - return SOutFN, S -end - function init_orbit_data{T<:Nemo.GroupElem}(logger, name::String, G::Nemo.Group, S::Vector{T}, AutS::Nemo.Group; radius=2) @@ -207,35 +195,14 @@ function init_orbit_data{T<:Nemo.GroupElem}(logger, name::String, G::Nemo.Group, return 0 end -function main() - parsed_args = parse_commandline() +function orbit_check_propertyT(logger, dirname, G, S, AutS, solver, upper_bound) - N = parsed_args["N"] - radius = parsed_args["radius"] - dirname = "SOutF$(N)_E_r=$radius" + init_orbit_data(logger, dirname, G, S, AutS; radius=2) - isdir(dirname) || mkdir(dirname) - - logger = PropertyT.setup_logging(dirname) - - if parsed_args["cpus"] != nothing - if parsed_args["cpus"] > cpuinfo_physicalcores() - warn("Number of specified cores exceeds the physical core cound. Performance will suffer.") - end - Blas.set_num_threads(parsed_args["cpus"]) - end - - G, S = SOutFN_generating_set(logger, dirname, N) - - init_orbit_data(logger, dirname, G, S; radius=2) - - tol = parsed_args["tol"] - iterations = parsed_args["iterations"] - upper_bound = parsed_args["upper-bound"] - - solver = SCS.SCSSolver(eps=tol, max_iters=iterations, verbose=true, linearsolver=SCS.Indirect) + Δ = PropertyT.ΔandSDPconstraints(name, G)[1] fnames = PropertyT.λSDPfilenames(dirname) + if all(isfile.(fnames)) λ, P = PropertyT.λandP(dirname) else @@ -270,5 +237,3 @@ function main() info(logger, "κ($name, S) ≥ $λ < 0: Tells us nothing about property (T)") return false end - -main()