82 lines
2.3 KiB
Julia
82 lines
2.3 KiB
Julia
using PropertyT
|
|
|
|
using SCS.SCSSolver
|
|
# using Mosek
|
|
# using CSDP
|
|
# using SDPA
|
|
|
|
scs_solver(tol, iterations) = SCSSolver(eps=tol, max_iters=iterations, linearsolver=SCS.Direct, alpha=1.95, acceleration_lookback=1)
|
|
|
|
# solver = Mosek.MosekSolver(
|
|
# MSK_DPAR_INTPNT_CO_TOL_REL_GAP=tol,
|
|
# MSK_IPAR_INTPNT_MAX_ITERATIONS=iterations,
|
|
# QUIET=false)
|
|
|
|
# solver = CSDP.CSDPSolver(axtol=tol, atytol=tol, objtol=tol, minstepp=tol*10.0^-1, minstepd=tol*10.0^-1)
|
|
|
|
# solver = SDPA.SDPASolver(epsilonStar=tol, epsilonDash=tol)
|
|
|
|
include("FPGroups_GAP.jl")
|
|
|
|
include("groups/Allgroups.jl")
|
|
using PropertyTGroups
|
|
|
|
function summarize(groupdir, iterations, tol, upper_bound, radius, G, S)
|
|
info("Group: $groupdir")
|
|
info("Iterations: $iterations")
|
|
info("Precision: $tol")
|
|
info("Upper bound: $upper_bound")
|
|
info("Radius: $radius")
|
|
import PropertyT.Settings
|
|
|
|
info("Threads: $(Threads.nthreads())")
|
|
info("Workers: $(workers())")
|
|
info(string(G))
|
|
info("with generating set of size $(length(S))")
|
|
end
|
|
|
|
function Settings(Gr::PropertyTGroup, args, solver):
|
|
r = get(args, "radius", 2)
|
|
gr_name = PropertyTGroups.name(Gr)*"_r$r"
|
|
G = PropertyTGroups.group(Gr)
|
|
S = PropertyTGroups.generatingset(Gr)
|
|
|
|
sol = solver
|
|
ub = get(args,"upper_bound", Inf)
|
|
tol = get(args,"tol", 1e-10)
|
|
ws = get(args, "warmstart", false)
|
|
|
|
if get(args, "nosymmetry", false)
|
|
autS = PropertyTGroups.autS(Gr)
|
|
return PropertyT.Settings(gr_name, G, S, r, sol, ub, tol, ws, autS)
|
|
else
|
|
return PropertyT.Settings(gr_name, G, S, r, sol, ub, tol, ws)
|
|
end
|
|
end
|
|
|
|
function main(Gr::SymmetrizedGroup)
|
|
sett = Settings(Gr)
|
|
|
|
isdir(PropertyT.fullpath(sett)) || mkpath(PropertyT.fullpath(sett))
|
|
|
|
if Gr.args["nosymmetry"]
|
|
return PropertyT.check_property_T(PropertyT.Naive, sett)
|
|
else
|
|
autS = PropertyTGroups.autS(Gr)
|
|
info("Symmetrising with $(autS)")
|
|
sett.autS = autS
|
|
return PropertyT.check_property_T(PropertyT.Symmetrize, sett)
|
|
end
|
|
end
|
|
|
|
function main(Gr::GAPGroup)
|
|
sett = Settings(Gr)
|
|
|
|
S = [s for s in sett.S if s.symbols[1].pow == 1]
|
|
relations = [k*inv(v) for (k,v) in sett.G.rels]
|
|
|
|
prepare_pm_delta(PropertyT.prepath(sett), GAP_groupcode(S, relations), sett.radius)
|
|
|
|
return PropertyT.check_property_T(PropertyT.Naive, sett)
|
|
end
|