remove Group specific code

This commit is contained in:
kalmar 2017-06-08 20:07:34 +02:00
parent f9127c59b4
commit 1468986b4a

View File

@ -181,18 +181,6 @@ function λandP(m::JuMP.Model, data::ProblemData)
return λ, recP return λ, recP
end 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, function init_orbit_data{T<:Nemo.GroupElem}(logger, name::String, G::Nemo.Group,
S::Vector{T}, AutS::Nemo.Group; radius=2) 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 return 0
end end
function main() function orbit_check_propertyT(logger, dirname, G, S, AutS, solver, upper_bound)
parsed_args = parse_commandline()
N = parsed_args["N"] init_orbit_data(logger, dirname, G, S, AutS; radius=2)
radius = parsed_args["radius"]
dirname = "SOutF$(N)_E_r=$radius"
isdir(dirname) || mkdir(dirname) Δ = PropertyT.ΔandSDPconstraints(name, G)[1]
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)
fnames = PropertyT.λSDPfilenames(dirname) fnames = PropertyT.λSDPfilenames(dirname)
if all(isfile.(fnames)) if all(isfile.(fnames))
λ, P = PropertyT.λandP(dirname) λ, P = PropertyT.λandP(dirname)
else else
@ -270,5 +237,3 @@ function main()
info(logger, "κ($name, S) ≥ < 0: Tells us nothing about property (T)") info(logger, "κ($name, S) ≥ < 0: Tells us nothing about property (T)")
return false return false
end end
main()