Compare commits
19 Commits
enh/global
...
master
Author | SHA1 | Date | |
---|---|---|---|
fa326a147f | |||
f1d0ff9c14 | |||
8666cc6b2c | |||
9fc23c43d7 | |||
c075b93313 | |||
bf8fc67bb4 | |||
523d783614 | |||
7b06dc1eb2 | |||
7dd061cdb7 | |||
b45adc2e00 | |||
8ee8eebd9e | |||
4cd2212554 | |||
b84b23d486 | |||
c7b8776e2d | |||
1b45625a71 | |||
1b6793f37c | |||
8a2aa6542e | |||
aaa36b2d5c | |||
88fa1ded31 |
2
.gitignore
vendored
2
.gitignore
vendored
@ -11,3 +11,5 @@ SL*_*
|
|||||||
*.gws
|
*.gws
|
||||||
.*
|
.*
|
||||||
tests*
|
tests*
|
||||||
|
*.py
|
||||||
|
*.pyc
|
||||||
|
65
AutFn.jl
65
AutFn.jl
@ -1,65 +0,0 @@
|
|||||||
using ArgParse
|
|
||||||
|
|
||||||
###############################################################################
|
|
||||||
#
|
|
||||||
# Parsing command line
|
|
||||||
#
|
|
||||||
###############################################################################
|
|
||||||
|
|
||||||
function parse_commandline()
|
|
||||||
s = ArgParseSettings()
|
|
||||||
|
|
||||||
@add_arg_table s begin
|
|
||||||
"--tol"
|
|
||||||
help = "set numerical tolerance for the SDP solver"
|
|
||||||
arg_type = Float64
|
|
||||||
default = 1e-6
|
|
||||||
"--iterations"
|
|
||||||
help = "set maximal number of iterations for the SDP solver"
|
|
||||||
arg_type = Int
|
|
||||||
default = 50000
|
|
||||||
"--upper-bound"
|
|
||||||
help = "Set an upper bound for the spectral gap"
|
|
||||||
arg_type = Float64
|
|
||||||
default = Inf
|
|
||||||
"--cpus"
|
|
||||||
help = "Set number of cpus used by solver (default: auto)"
|
|
||||||
arg_type = Int
|
|
||||||
required = false
|
|
||||||
"--radius"
|
|
||||||
help = "Radius of ball B_r(e,S) to find solution over"
|
|
||||||
arg_type = Int
|
|
||||||
default = 2
|
|
||||||
"--warmstart"
|
|
||||||
help = "Use warmstart.jld as the initial guess for SCS"
|
|
||||||
action = :store_true
|
|
||||||
"--nosymmetry"
|
|
||||||
help = "Don't use symmetries of the Laplacian"
|
|
||||||
action = :store_true
|
|
||||||
"N"
|
|
||||||
help = "Compute for the automorphisms group of the free group on N generators"
|
|
||||||
arg_type = Int
|
|
||||||
required = true
|
|
||||||
end
|
|
||||||
return parse_args(s)
|
|
||||||
end
|
|
||||||
const PARSEDARGS = parse_commandline()
|
|
||||||
|
|
||||||
#=
|
|
||||||
Note that the element
|
|
||||||
α(i,j,k) = ϱ(i,j)*ϱ(i,k)*inv(ϱ(i,j))*inv(ϱ(i,k)),
|
|
||||||
which surely belongs to ball of radius 4 in Aut(Fₙ) becomes trivial under the representation
|
|
||||||
Aut(Fₙ) → GLₙ(ℤ)⋉ℤⁿ → GL_(n+1)(ℂ).
|
|
||||||
Moreover, due to work of Potapchik and Rapinchuk [1] every real representation of Aut(Fₙ) into GLₘ(ℂ) (for m ≤ 2n-2) factors through GLₙ(ℤ)⋉ℤⁿ, so will have the same problem.
|
|
||||||
|
|
||||||
We need a different approach: Here we actually compute in (S)Aut(𝔽ₙ)
|
|
||||||
=#
|
|
||||||
|
|
||||||
include("CPUselect.jl")
|
|
||||||
set_parallel_mthread(PARSEDARGS, workers=true)
|
|
||||||
|
|
||||||
include("main.jl")
|
|
||||||
|
|
||||||
G = PropertyTGroups.SpecialAutomorphismGroup(PARSEDARGS)
|
|
||||||
|
|
||||||
main(G)
|
|
62
FPgroup.jl
62
FPgroup.jl
@ -1,62 +0,0 @@
|
|||||||
using ArgParse
|
|
||||||
|
|
||||||
function parse_commandline()
|
|
||||||
args = ArgParseSettings()
|
|
||||||
|
|
||||||
@add_arg_table args begin
|
|
||||||
"--tol"
|
|
||||||
help = "set numerical tolerance for the SDP solver"
|
|
||||||
arg_type = Float64
|
|
||||||
default = 1e-6
|
|
||||||
"--iterations"
|
|
||||||
help = "set maximal number of iterations for the SDP solver"
|
|
||||||
arg_type = Int
|
|
||||||
default = 50000
|
|
||||||
"--upper-bound"
|
|
||||||
help = "Set an upper bound for the spectral gap"
|
|
||||||
arg_type = Float64
|
|
||||||
default = Inf
|
|
||||||
"--cpus"
|
|
||||||
help = "Set number of cpus used by solver (default: auto)"
|
|
||||||
arg_type = Int
|
|
||||||
required = false
|
|
||||||
"--radius"
|
|
||||||
help = "Radius of ball B_r(e,S) to find solution over"
|
|
||||||
arg_type = Int
|
|
||||||
default = 2
|
|
||||||
"--warmstart"
|
|
||||||
help = "Use warmstart.jl as the initial guess for SCS"
|
|
||||||
action = :store_true
|
|
||||||
"--MCG"
|
|
||||||
help = "Compute for mapping class group of surface of genus N"
|
|
||||||
arg_type = Int
|
|
||||||
required = false
|
|
||||||
"--Higman"
|
|
||||||
help = "Compute for Higman Group"
|
|
||||||
action = :store_true
|
|
||||||
"--Caprace"
|
|
||||||
help = "Compute for Higman Group"
|
|
||||||
action = :store_true
|
|
||||||
end
|
|
||||||
|
|
||||||
return parse_args(args)
|
|
||||||
end
|
|
||||||
const PARSEDARGS = parse_commandline()
|
|
||||||
|
|
||||||
include("CPUselect.jl")
|
|
||||||
set_parallel_mthread(PARSEDARGS, workers=false)
|
|
||||||
|
|
||||||
include("main.jl")
|
|
||||||
include("FPGroups_GAP.jl")
|
|
||||||
|
|
||||||
if PARSEDARGS["Caprace"]
|
|
||||||
G = PropertyTGroups.CapraceGroup(PARSEDARGS)
|
|
||||||
elseif PARSEDARGS["Higman"]
|
|
||||||
G = PropertyTGroups.HigmanGroup(PARSEDARGS)
|
|
||||||
elseif PARSEDARGS["MCG"] != nothing
|
|
||||||
G = PropertyTGroups.MappingClassGroup(PARSEDARGS)
|
|
||||||
else
|
|
||||||
throw("You need to specify one of the --Higman, --Caprace, --MCG N")
|
|
||||||
end
|
|
||||||
|
|
||||||
main(G)
|
|
@ -1,4 +1,10 @@
|
|||||||
This repository contains code for computations in [Certifying Numerical Estimates of Spectral Gaps](https://arxiv.org/abs/1703.09680).
|
# DEPRECATED!
|
||||||
|
|
||||||
|
This repository has not been updated for a while!
|
||||||
|
If You are interested in replicating results for [1712.07167](https://arxiv.org/abs/1712.07167) please check [these instruction](https://kalmar.faculty.wmi.amu.edu.pl/post/1712.07176/)
|
||||||
|
Also [this notebook](https://nbviewer.jupyter.org/gist/kalmarek/03510181bc1e7c98615e86e1ec580b2a) could be of some help. If everything else fails the [zenodo dataset](https://zenodo.org/record/1133440) should contain the last-resort instructions.
|
||||||
|
|
||||||
|
This repository contains some legacy code for computations in [Certifying Numerical Estimates of Spectral Gaps](https://arxiv.org/abs/1703.09680).
|
||||||
|
|
||||||
# Installing
|
# Installing
|
||||||
To run the code You need `julia-v0.5` (should work on `v0.6`, but with warnings).
|
To run the code You need `julia-v0.5` (should work on `v0.6`, but with warnings).
|
||||||
|
62
SLn.jl
62
SLn.jl
@ -1,62 +0,0 @@
|
|||||||
using ArgParse
|
|
||||||
|
|
||||||
###############################################################################
|
|
||||||
#
|
|
||||||
# Parsing command line
|
|
||||||
#
|
|
||||||
###############################################################################
|
|
||||||
|
|
||||||
function parse_commandline()
|
|
||||||
settings = ArgParseSettings()
|
|
||||||
|
|
||||||
@add_arg_table settings begin
|
|
||||||
"--tol"
|
|
||||||
help = "set numerical tolerance for the SDP solver"
|
|
||||||
arg_type = Float64
|
|
||||||
default = 1e-6
|
|
||||||
"--iterations"
|
|
||||||
help = "set maximal number of iterations for the SDP solver"
|
|
||||||
arg_type = Int
|
|
||||||
default = 50000
|
|
||||||
"--upper-bound"
|
|
||||||
help = "Set an upper bound for the spectral gap"
|
|
||||||
arg_type = Float64
|
|
||||||
default = Inf
|
|
||||||
"--cpus"
|
|
||||||
help = "Set number of cpus used by solver"
|
|
||||||
arg_type = Int
|
|
||||||
required = false
|
|
||||||
"--radius"
|
|
||||||
help = "Radius of ball B_r(e,S) to find solution over"
|
|
||||||
arg_type = Int
|
|
||||||
default = 2
|
|
||||||
"--warmstart"
|
|
||||||
help = "Use warmstart.jld as the initial guess for SCS"
|
|
||||||
action = :store_true
|
|
||||||
"--nosymmetry"
|
|
||||||
help = "Don't use symmetries of the Laplacian"
|
|
||||||
action = :store_true
|
|
||||||
"-p"
|
|
||||||
help = "Matrices over field of p-elements (p=0 => over ZZ)"
|
|
||||||
arg_type = Int
|
|
||||||
default = 0
|
|
||||||
"-X"
|
|
||||||
help = "Consider EL(N, ZZ⟨X⟩)"
|
|
||||||
action = :store_true
|
|
||||||
"N"
|
|
||||||
help = "Compute with the group generated by elementary matrices of size n×n"
|
|
||||||
arg_type = Int
|
|
||||||
default = 2
|
|
||||||
end
|
|
||||||
return parse_args(settings)
|
|
||||||
end
|
|
||||||
const PARSEDARGS = parse_commandline()
|
|
||||||
|
|
||||||
include("CPUselect.jl")
|
|
||||||
set_parallel_mthread(PARSEDARGS, workers=true)
|
|
||||||
|
|
||||||
include("main.jl")
|
|
||||||
|
|
||||||
G = PropertyTGroups.SpecialLinearGroup(PARSEDARGS)
|
|
||||||
|
|
||||||
main(G)
|
|
@ -4,6 +4,7 @@ using PropertyT
|
|||||||
using AbstractAlgebra
|
using AbstractAlgebra
|
||||||
using Nemo
|
using Nemo
|
||||||
using Groups
|
using Groups
|
||||||
|
using GroupRings
|
||||||
|
|
||||||
export PropertyTGroup, SymmetrizedGroup, GAPGroup,
|
export PropertyTGroup, SymmetrizedGroup, GAPGroup,
|
||||||
SpecialLinearGroup,
|
SpecialLinearGroup,
|
||||||
@ -50,5 +51,6 @@ end
|
|||||||
include("mappingclassgroup.jl")
|
include("mappingclassgroup.jl")
|
||||||
include("higman.jl")
|
include("higman.jl")
|
||||||
include("caprace.jl")
|
include("caprace.jl")
|
||||||
|
include("actions.jl")
|
||||||
|
|
||||||
end # of module PropertyTGroups
|
end # of module PropertyTGroups
|
||||||
|
92
groups/actions.jl
Normal file
92
groups/actions.jl
Normal file
@ -0,0 +1,92 @@
|
|||||||
|
function (p::perm)(A::GroupRingElem)
|
||||||
|
RG = parent(A)
|
||||||
|
result = zero(RG, eltype(A.coeffs))
|
||||||
|
|
||||||
|
for (idx, c) in enumerate(A.coeffs)
|
||||||
|
if c!= zero(eltype(A.coeffs))
|
||||||
|
result[p(RG.basis[idx])] = c
|
||||||
|
end
|
||||||
|
end
|
||||||
|
return result
|
||||||
|
end
|
||||||
|
|
||||||
|
###############################################################################
|
||||||
|
#
|
||||||
|
# Action of WreathProductElems on Nemo.MatElem
|
||||||
|
#
|
||||||
|
###############################################################################
|
||||||
|
|
||||||
|
function matrix_emb(n::DirectProductGroupElem, p::perm)
|
||||||
|
Id = parent(n.elts[1])()
|
||||||
|
elt = diagm([(-1)^(el == Id ? 0 : 1) for el in n.elts])
|
||||||
|
return elt[:, p.d]
|
||||||
|
end
|
||||||
|
|
||||||
|
function (g::WreathProductElem)(A::MatElem)
|
||||||
|
g_inv = inv(g)
|
||||||
|
G = matrix_emb(g.n, g_inv.p)
|
||||||
|
G_inv = matrix_emb(g_inv.n, g.p)
|
||||||
|
M = parent(A)
|
||||||
|
return M(G)*A*M(G_inv)
|
||||||
|
end
|
||||||
|
|
||||||
|
import Base.*
|
||||||
|
|
||||||
|
doc"""
|
||||||
|
*(x::AbstractAlgebra.MatElem, P::Generic.perm)
|
||||||
|
> Apply the pemutation $P$ to the rows of the matrix $x$ and return the result.
|
||||||
|
"""
|
||||||
|
function *(x::AbstractAlgebra.MatElem, P::Generic.perm)
|
||||||
|
z = similar(x)
|
||||||
|
m = rows(x)
|
||||||
|
n = cols(x)
|
||||||
|
for i = 1:m
|
||||||
|
for j = 1:n
|
||||||
|
z[i, j] = x[i,P[j]]
|
||||||
|
end
|
||||||
|
end
|
||||||
|
return z
|
||||||
|
end
|
||||||
|
|
||||||
|
function (p::perm)(A::MatElem)
|
||||||
|
length(p.d) == A.r == A.c || throw("Can't act via $p on matrix of size ($(A.r), $(A.c))")
|
||||||
|
return p*A*inv(p)
|
||||||
|
end
|
||||||
|
|
||||||
|
###############################################################################
|
||||||
|
#
|
||||||
|
# Action of WreathProductElems on AutGroupElem
|
||||||
|
#
|
||||||
|
###############################################################################
|
||||||
|
|
||||||
|
function AutFG_emb(A::AutGroup, g::WreathProductElem)
|
||||||
|
isa(A.objectGroup, FreeGroup) || throw("Not an Aut(Fₙ)")
|
||||||
|
parent(g).P.n == length(A.objectGroup.gens) || throw("No natural embedding of $(parent(g)) into $A")
|
||||||
|
elt = A()
|
||||||
|
Id = parent(g.n.elts[1])()
|
||||||
|
flips = Groups.AutSymbol[Groups.flip_autsymbol(i) for i in 1:length(g.p.d) if g.n.elts[i] != Id]
|
||||||
|
Groups.r_multiply!(elt, flips, reduced=false)
|
||||||
|
Groups.r_multiply!(elt, [Groups.perm_autsymbol(g.p)])
|
||||||
|
return elt
|
||||||
|
end
|
||||||
|
|
||||||
|
function AutFG_emb(A::AutGroup, p::perm)
|
||||||
|
isa(A.objectGroup, FreeGroup) || throw("Not an Aut(Fₙ)")
|
||||||
|
parent(p).n == length(A.objectGroup.gens) || throw("No natural embedding of $(parent(p)) into $A")
|
||||||
|
return A(Groups.perm_autsymbol(p))
|
||||||
|
end
|
||||||
|
|
||||||
|
function (g::WreathProductElem)(a::Groups.Automorphism)
|
||||||
|
A = parent(a)
|
||||||
|
g = AutFG_emb(A,g)
|
||||||
|
res = A()
|
||||||
|
Groups.r_multiply!(res, g.symbols, reduced=false)
|
||||||
|
Groups.r_multiply!(res, a.symbols, reduced=false)
|
||||||
|
Groups.r_multiply!(res, [inv(s) for s in reverse!(g.symbols)])
|
||||||
|
return res
|
||||||
|
end
|
||||||
|
|
||||||
|
function (p::perm)(a::Groups.Automorphism)
|
||||||
|
g = AutFG_emb(parent(a),p)
|
||||||
|
return g*a*inv(g)
|
||||||
|
end
|
@ -1,21 +1,13 @@
|
|||||||
struct SpecialAutomorphismGroup <: SymmetrizedGroup
|
struct SpecialAutomorphismGroup{N} <: SymmetrizedGroup
|
||||||
args::Dict{String,Any}
|
|
||||||
group::AutGroup
|
group::AutGroup
|
||||||
N::Int
|
end
|
||||||
|
|
||||||
function SpecialAutomorphismGroup(args::Dict)
|
function SpecialAutomorphismGroup(args::Dict)
|
||||||
N = args["SAut"]
|
N = args["SAut"]
|
||||||
return new(args, AutGroup(FreeGroup(N), special=true), N)
|
return SpecialAutomorphismGroup{N}(AutGroup(FreeGroup(N), special=true))
|
||||||
end
|
|
||||||
end
|
end
|
||||||
|
|
||||||
function name(G::SpecialAutomorphismGroup)
|
name(G::SpecialAutomorphismGroup{N}) where N = "SAutF$(N)"
|
||||||
if G.args["nosymmetry"]
|
|
||||||
return "SAutF$(G.N)"
|
|
||||||
else
|
|
||||||
return "oSAutF$(G.N)"
|
|
||||||
end
|
|
||||||
end
|
|
||||||
|
|
||||||
group(G::SpecialAutomorphismGroup) = G.group
|
group(G::SpecialAutomorphismGroup) = G.group
|
||||||
|
|
||||||
@ -24,44 +16,6 @@ function generatingset(G::SpecialAutomorphismGroup)
|
|||||||
return unique([S; inv.(S)])
|
return unique([S; inv.(S)])
|
||||||
end
|
end
|
||||||
|
|
||||||
function autS(G::SpecialAutomorphismGroup)
|
function autS(G::SpecialAutomorphismGroup{N}) where N
|
||||||
return WreathProduct(PermutationGroup(2), PermutationGroup(G.N))
|
return WreathProduct(PermutationGroup(2), PermutationGroup(N))
|
||||||
end
|
|
||||||
|
|
||||||
###############################################################################
|
|
||||||
#
|
|
||||||
# Action of WreathProductElems on AutGroupElem
|
|
||||||
#
|
|
||||||
###############################################################################
|
|
||||||
|
|
||||||
function AutFG_emb(A::AutGroup, g::WreathProductElem)
|
|
||||||
isa(A.objectGroup, FreeGroup) || throw("Not an Aut(Fₙ)")
|
|
||||||
parent(g).P.n == length(A.objectGroup.gens) || throw("No natural embedding of $(parent(g)) into $A")
|
|
||||||
elt = A()
|
|
||||||
Id = parent(g.n.elts[1])()
|
|
||||||
flips = Groups.AutSymbol[Groups.flip_autsymbol(i) for i in 1:length(g.p.d) if g.n.elts[i] != Id]
|
|
||||||
Groups.r_multiply!(elt, flips, reduced=false)
|
|
||||||
Groups.r_multiply!(elt, [Groups.perm_autsymbol(g.p)])
|
|
||||||
return elt
|
|
||||||
end
|
|
||||||
|
|
||||||
function AutFG_emb(A::AutGroup, p::perm)
|
|
||||||
isa(A.objectGroup, FreeGroup) || throw("Not an Aut(Fₙ)")
|
|
||||||
parent(p).n == length(A.objectGroup.gens) || throw("No natural embedding of $(parent(g)) into $A")
|
|
||||||
return A(Groups.perm_autsymbol(p))
|
|
||||||
end
|
|
||||||
|
|
||||||
function (g::WreathProductElem)(a::Groups.Automorphism)
|
|
||||||
A = parent(a)
|
|
||||||
g = AutFG_emb(A,g)
|
|
||||||
res = A()
|
|
||||||
Groups.r_multiply!(res, g.symbols, reduced=false)
|
|
||||||
Groups.r_multiply!(res, a.symbols, reduced=false)
|
|
||||||
Groups.r_multiply!(res, [inv(s) for s in reverse!(g.symbols)])
|
|
||||||
return res
|
|
||||||
end
|
|
||||||
|
|
||||||
function (p::perm)(a::Groups.Automorphism)
|
|
||||||
g = AutFG_emb(parent(a),p)
|
|
||||||
return g*a*inv(g)
|
|
||||||
end
|
end
|
||||||
|
@ -1,6 +1,4 @@
|
|||||||
struct CapraceGroup <: GAPGroup
|
struct CapraceGroup <: GAPGroup end
|
||||||
args::Dict{String,Any}
|
|
||||||
end
|
|
||||||
|
|
||||||
name(G::CapraceGroup) = "CapraceGroup"
|
name(G::CapraceGroup) = "CapraceGroup"
|
||||||
|
|
||||||
|
@ -1,6 +1,4 @@
|
|||||||
struct HigmanGroup <: GAPGroup
|
struct HigmanGroup <: GAPGroup end
|
||||||
args::Dict{String,Any}
|
|
||||||
end
|
|
||||||
|
|
||||||
name(G::HigmanGroup) = "HigmanGroup"
|
name(G::HigmanGroup) = "HigmanGroup"
|
||||||
|
|
||||||
|
@ -1,17 +1,14 @@
|
|||||||
struct MappingClassGroup <: GAPGroup
|
struct MappingClassGroup{N} <: GAPGroup end
|
||||||
args::Dict{String,Any}
|
|
||||||
N::Int
|
|
||||||
|
|
||||||
MappingClassGroup(args) = MappingClassGroup(args, args["MCG"])
|
MappingClassGroup(args::Dict) = MappingClassGroup{args["MCG"]}()
|
||||||
end
|
|
||||||
|
|
||||||
name(G::MappingClassGroup) = "MCG($(G.N))"
|
name(G::MappingClassGroup{N}) where N = "MCG(N)"
|
||||||
|
|
||||||
function group(G::MappingClassGroup)
|
function group(G::MappingClassGroup{N}) where N
|
||||||
|
|
||||||
if G.N < 2
|
if N < 2
|
||||||
throw("Genus must be at least 2!")
|
throw("Genus must be at least 2!")
|
||||||
elseif G.N == 2
|
elseif N == 2
|
||||||
MCGroup = Groups.FPGroup(["a1","a2","a3","a4","a5"]);
|
MCGroup = Groups.FPGroup(["a1","a2","a3","a4","a5"]);
|
||||||
S = gens(MCGroup)
|
S = gens(MCGroup)
|
||||||
|
|
||||||
@ -31,7 +28,7 @@ function group(G::MappingClassGroup)
|
|||||||
return MCGroup
|
return MCGroup
|
||||||
|
|
||||||
else
|
else
|
||||||
MCGroup = Groups.FPGroup(["a$i" for i in 0:2G.N])
|
MCGroup = Groups.FPGroup(["a$i" for i in 0:2N])
|
||||||
S = gens(MCGroup)
|
S = gens(MCGroup)
|
||||||
|
|
||||||
a0 = S[1]
|
a0 = S[1]
|
||||||
@ -76,7 +73,7 @@ function group(G::MappingClassGroup)
|
|||||||
(A[2i+3]*A[2i+2]*A[2i+4]*A[2i+3])*( n(i+1)*A[2i+2]*A[2i+1]*A[2i] )
|
(A[2i+3]*A[2i+2]*A[2i+4]*A[2i+3])*( n(i+1)*A[2i+2]*A[2i+1]*A[2i] )
|
||||||
end
|
end
|
||||||
|
|
||||||
# push!(relations, X*n(G.N)*inv(n(G.N)*X))
|
# push!(relations, X*n(N)*inv(n(N)*X))
|
||||||
|
|
||||||
relations = [relations; [inv(rel) for rel in relations]]
|
relations = [relations; [inv(rel) for rel in relations]]
|
||||||
Groups.add_rels!(MCGroup, Dict(rel => MCGroup() for rel in relations))
|
Groups.add_rels!(MCGroup, Dict(rel => MCGroup() for rel in relations))
|
||||||
|
@ -1,60 +1,44 @@
|
|||||||
struct SpecialLinearGroup <: SymmetrizedGroup
|
struct SpecialLinearGroup{N} <: SymmetrizedGroup
|
||||||
args::Dict{String,Any}
|
|
||||||
group::AbstractAlgebra.Group
|
group::AbstractAlgebra.Group
|
||||||
N::Int
|
p::Int
|
||||||
|
X::Bool
|
||||||
|
end
|
||||||
|
|
||||||
function SpecialLinearGroup(args::Dict)
|
function SpecialLinearGroup(args::Dict)
|
||||||
n = args["SL"]
|
N = args["SL"]
|
||||||
p = args["p"]
|
p = args["p"]
|
||||||
X = args["X"]
|
X = args["X"]
|
||||||
|
|
||||||
if p == 0
|
if p == 0
|
||||||
G = MatrixSpace(Nemo.ZZ, n, n)
|
G = MatrixSpace(Nemo.ZZ, N, N)
|
||||||
else
|
else
|
||||||
R = Nemo.NmodRing(UInt(p))
|
R = Nemo.NmodRing(UInt(p))
|
||||||
G = MatrixSpace(R, n, n)
|
G = MatrixSpace(R, N, N)
|
||||||
end
|
|
||||||
return new(args, G, n)
|
|
||||||
end
|
end
|
||||||
|
return SpecialLinearGroup{N}(G, p, X)
|
||||||
end
|
end
|
||||||
|
|
||||||
function name(G::SpecialLinearGroup)
|
function name(G::SpecialLinearGroup{N}) where N
|
||||||
p = G.args["p"]
|
if G.p == 0
|
||||||
X = G.args["X"]
|
R = (G.X ? "Z[x]" : "Z")
|
||||||
|
|
||||||
if p == 0
|
|
||||||
R = (X ? "Z[x]" : "Z")
|
|
||||||
else
|
else
|
||||||
R = "F$p"
|
R = "F$(G.p)"
|
||||||
end
|
|
||||||
if G.args["nosymmetry"]
|
|
||||||
return "SL($(G.N),$R)"
|
|
||||||
else
|
|
||||||
return "oSL($(G.N),$R)"
|
|
||||||
end
|
end
|
||||||
|
return SL($(G.N),$R)
|
||||||
end
|
end
|
||||||
|
|
||||||
group(G::SpecialLinearGroup) = G.group
|
group(G::SpecialLinearGroup) = G.group
|
||||||
|
|
||||||
function E(i::Int, j::Int, M::MatSpace, val=one(M.base_ring))
|
function generatingset(G::SpecialLinearGroup{N}) where N
|
||||||
@assert i≠j
|
G.p > 0 && G.X && throw("SL(n, F_p[x]) not implemented")
|
||||||
m = one(M)
|
|
||||||
m[i,j] = val
|
|
||||||
return m
|
|
||||||
end
|
|
||||||
|
|
||||||
function generatingset(G::SpecialLinearGroup)
|
|
||||||
p = G.args["p"]
|
|
||||||
X = G.args["X"]
|
|
||||||
p > 0 && X && throw("SL(n, F_p[x]) not implemented")
|
|
||||||
SL = group(G)
|
SL = group(G)
|
||||||
r = G.args["radius"]
|
return generatingset(SL, G.X)
|
||||||
return generatingset(SL, r, X)
|
|
||||||
end
|
end
|
||||||
|
|
||||||
|
# r is the injectivity radius of
|
||||||
|
# SL(n, Z[X]) -> SL(n, Z) induced by X -> 100
|
||||||
|
|
||||||
function generatingset(SL::MatSpace, radius::Integer, X::Bool=false)
|
function generatingset(SL::MatSpace, X::Bool=false, r=5)
|
||||||
|
|
||||||
n = SL.cols
|
n = SL.cols
|
||||||
indexing = [(i,j) for i in 1:n for j in 1:n if i≠j]
|
indexing = [(i,j) for i in 1:n for j in 1:n if i≠j]
|
||||||
|
|
||||||
@ -66,49 +50,13 @@ function generatingset(SL::MatSpace, radius::Integer, X::Bool=false)
|
|||||||
return unique([S; inv.(S)])
|
return unique([S; inv.(S)])
|
||||||
end
|
end
|
||||||
|
|
||||||
function autS(G::SpecialLinearGroup)
|
function E(i::Int, j::Int, M::MatSpace, val=one(M.base_ring))
|
||||||
return WreathProduct(PermutationGroup(2), PermutationGroup(G.N))
|
@assert i≠j
|
||||||
|
m = one(M)
|
||||||
|
m[i,j] = val
|
||||||
|
return m
|
||||||
end
|
end
|
||||||
|
|
||||||
###############################################################################
|
function autS(G::SpecialLinearGroup{N}) where N
|
||||||
#
|
return WreathProduct(PermutationGroup(2), PermutationGroup(N))
|
||||||
# Action of WreathProductElems on Nemo.MatElem
|
|
||||||
#
|
|
||||||
###############################################################################
|
|
||||||
|
|
||||||
function matrix_emb(n::DirectProductGroupElem, p::perm)
|
|
||||||
Id = parent(n.elts[1])()
|
|
||||||
elt = diagm([(-1)^(el == Id ? 0 : 1) for el in n.elts])
|
|
||||||
return elt[:, p.d]
|
|
||||||
end
|
|
||||||
|
|
||||||
function (g::WreathProductElem)(A::MatElem)
|
|
||||||
g_inv = inv(g)
|
|
||||||
G = matrix_emb(g.n, g_inv.p)
|
|
||||||
G_inv = matrix_emb(g_inv.n, g.p)
|
|
||||||
M = parent(A)
|
|
||||||
return M(G)*A*M(G_inv)
|
|
||||||
end
|
|
||||||
|
|
||||||
import Base.*
|
|
||||||
|
|
||||||
doc"""
|
|
||||||
*(x::AbstractAlgebra.MatElem, P::Generic.perm)
|
|
||||||
> Apply the pemutation $P$ to the rows of the matrix $x$ and return the result.
|
|
||||||
"""
|
|
||||||
function *(x::AbstractAlgebra.MatElem, P::Generic.perm)
|
|
||||||
z = similar(x)
|
|
||||||
m = rows(x)
|
|
||||||
n = cols(x)
|
|
||||||
for i = 1:m
|
|
||||||
for j = 1:n
|
|
||||||
z[i, j] = x[i,P[j]]
|
|
||||||
end
|
|
||||||
end
|
|
||||||
return z
|
|
||||||
end
|
|
||||||
|
|
||||||
function (p::perm)(A::MatElem)
|
|
||||||
length(p.d) == A.r == A.c || throw("Can't act via $p on matrix of size ($(A.r), $(A.c))")
|
|
||||||
return p*A*inv(p)
|
|
||||||
end
|
end
|
||||||
|
94
main.jl
94
main.jl
@ -1,87 +1,61 @@
|
|||||||
using PropertyT
|
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("FPGroups_GAP.jl")
|
||||||
|
|
||||||
include("groups/Allgroups.jl")
|
include("groups/Allgroups.jl")
|
||||||
using PropertyTGroups
|
using PropertyTGroups
|
||||||
|
|
||||||
function summarize(groupdir, iterations, tol, upper_bound, radius, G, S)
|
import PropertyT.Settings
|
||||||
info("Group: $groupdir")
|
|
||||||
info("Iterations: $iterations")
|
function summarize(sett::PropertyT.Settings)
|
||||||
info("Precision: $tol")
|
|
||||||
info("Upper bound: $upper_bound")
|
|
||||||
info("Radius: $radius")
|
|
||||||
info("Threads: $(Threads.nthreads())")
|
info("Threads: $(Threads.nthreads())")
|
||||||
info("Workers: $(workers())")
|
info("Workers: $(workers())")
|
||||||
info(string(G))
|
info("GroupDir: $(PropertyT.prepath(sett))")
|
||||||
info("with generating set of size $(length(S))")
|
info(string(sett.G))
|
||||||
|
info("with generating set of size $(length(sett.S))")
|
||||||
|
|
||||||
|
info("Radius: $(sett.radius)")
|
||||||
|
info("Precision: $(sett.tol)")
|
||||||
|
info("Upper bound: $(sett.upper_bound)")
|
||||||
|
info("Solver: $(sett.solver)")
|
||||||
end
|
end
|
||||||
|
|
||||||
function params(Gr::PropertyTGroup)
|
function Settings(Gr::PropertyTGroup, args, solver)
|
||||||
radius = Gr.args["radius"]
|
r = get(args, "radius", 2)
|
||||||
tol = Gr.args["tol"]
|
gr_name = PropertyTGroups.name(Gr)*"_r$r"
|
||||||
iterations = Gr.args["iterations"]
|
|
||||||
upper_bound = Gr.args["upper-bound"]
|
|
||||||
warm = Gr.args["warmstart"]
|
|
||||||
return radius, tol, iterations, upper_bound, warm
|
|
||||||
end
|
|
||||||
|
|
||||||
function Settings(Gr::PropertyTGroup)
|
|
||||||
r = Gr.args["radius"]
|
|
||||||
ub = Gr.args["upper-bound"]
|
|
||||||
groupdir = "$(PropertyTGroups.name(Gr))_r$r"
|
|
||||||
|
|
||||||
radius, tol, iterations, upper_bound, warm = params(Gr)
|
|
||||||
|
|
||||||
G = PropertyTGroups.group(Gr)
|
G = PropertyTGroups.group(Gr)
|
||||||
S = PropertyTGroups.generatingset(Gr)
|
S = PropertyTGroups.generatingset(Gr)
|
||||||
|
|
||||||
summarize(groupdir, iterations, tol, upper_bound, radius, G, S)
|
sol = solver
|
||||||
|
ub = get(args,"upper-bound", Inf)
|
||||||
|
tol = get(args,"tol", 1e-10)
|
||||||
|
ws = get(args, "warmstart", false)
|
||||||
|
|
||||||
solver = scs_solver(tol, iterations)
|
if get(args, "nosymmetry", false)
|
||||||
|
return PropertyT.Settings(gr_name, G, S, r, sol, ub, tol, ws)
|
||||||
return PropertyT.Settings(groupdir, G, S, radius,
|
|
||||||
solver, upper_bound, tol, warm)
|
|
||||||
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
|
else
|
||||||
autS = PropertyTGroups.autS(Gr)
|
autS = PropertyTGroups.autS(Gr)
|
||||||
info("Symmetrising with $(autS)")
|
return PropertyT.Settings(gr_name, G, S, r, sol, ub, tol, ws, autS)
|
||||||
sett.autS = autS
|
|
||||||
return PropertyT.check_property_T(PropertyT.Symmetrize, sett)
|
|
||||||
end
|
end
|
||||||
end
|
end
|
||||||
|
|
||||||
function main(Gr::GAPGroup)
|
function main(::PropertyTGroup, sett::PropertyT.Settings)
|
||||||
sett = Settings(Gr)
|
isdir(PropertyT.fullpath(sett)) || mkpath(PropertyT.fullpath(sett))
|
||||||
|
|
||||||
|
summarize(sett)
|
||||||
|
|
||||||
|
return PropertyT.check_property_T(sett)
|
||||||
|
end
|
||||||
|
|
||||||
|
function main(::GAPGroup, sett::PropertyT.Settings)
|
||||||
|
isdir(PropertyT.fullpath(sett)) || mkpath(PropertyT.fullpath(sett))
|
||||||
|
|
||||||
|
summarize(sett)
|
||||||
|
|
||||||
S = [s for s in sett.S if s.symbols[1].pow == 1]
|
S = [s for s in sett.S if s.symbols[1].pow == 1]
|
||||||
relations = [k*inv(v) for (k,v) in sett.G.rels]
|
relations = [k*inv(v) for (k,v) in sett.G.rels]
|
||||||
|
|
||||||
prepare_pm_delta(PropertyT.prepath(sett), GAP_groupcode(S, relations), sett.radius)
|
prepare_pm_delta(PropertyT.prepath(sett), GAP_groupcode(S, relations), sett.radius)
|
||||||
|
|
||||||
return PropertyT.check_property_T(PropertyT.Naive, sett)
|
return PropertyT.check_property_T(sett)
|
||||||
end
|
end
|
||||||
|
197
positivity/check_positivity.jl
Normal file
197
positivity/check_positivity.jl
Normal file
@ -0,0 +1,197 @@
|
|||||||
|
using AbstractAlgebra
|
||||||
|
using Groups
|
||||||
|
using GroupRings
|
||||||
|
using PropertyT
|
||||||
|
|
||||||
|
using SCS
|
||||||
|
solver(tol, iterations) =
|
||||||
|
SCSSolver(linearsolver=SCS.Direct,
|
||||||
|
eps=tol, max_iters=iterations,
|
||||||
|
alpha=1.95, acceleration_lookback=1)
|
||||||
|
|
||||||
|
include("../main.jl")
|
||||||
|
|
||||||
|
using PropertyTGroups
|
||||||
|
|
||||||
|
args = Dict("SAut" => 5, "upper-bound" => 50.0, "radius" => 2, "nosymmetry"=>false, "tol"=>1e-12, "iterations" =>200000, "warmstart" => true)
|
||||||
|
|
||||||
|
Gr = PropertyTGroups.PropertyTGroup(args)
|
||||||
|
sett = PropertyT.Settings(Gr, args,
|
||||||
|
solver(args["tol"], args["iterations"]))
|
||||||
|
|
||||||
|
@show sett
|
||||||
|
|
||||||
|
fullpath = PropertyT.fullpath(sett)
|
||||||
|
isdir(fullpath) || mkpath(fullpath)
|
||||||
|
# setup_logging(PropertyT.filename(fullpath, :fulllog), :fulllog)
|
||||||
|
|
||||||
|
function small_generating_set(RG::GroupRing{AutGroup{N}}, n) where N
|
||||||
|
indexing = [(i,j) for i in 1:n for j in 1:n if i≠j]
|
||||||
|
|
||||||
|
rmuls = [Groups.rmul_autsymbol(i,j) for (i,j) in indexing]
|
||||||
|
lmuls = [Groups.lmul_autsymbol(i,j) for (i,j) in indexing]
|
||||||
|
gen_set = RG.group.([rmuls; lmuls])
|
||||||
|
|
||||||
|
return [gen_set; inv.(gen_set)]
|
||||||
|
end
|
||||||
|
|
||||||
|
function computeX(RG::GroupRing{AutGroup{N}}) where N
|
||||||
|
Tn = small_generating_set(RG, N-1)
|
||||||
|
|
||||||
|
ℤ = Int64
|
||||||
|
Δn = length(Tn)*one(RG, ℤ) - RG(Tn, ℤ);
|
||||||
|
|
||||||
|
Alt_N = [g for g in elements(PermutationGroup(N)) if parity(g) == 0]
|
||||||
|
|
||||||
|
@time X = sum(σ(Δn)*sum(τ(Δn) for τ ∈ Alt_N if τ ≠ σ) for σ in Alt_N);
|
||||||
|
return X
|
||||||
|
end
|
||||||
|
|
||||||
|
function Sq(RG::GroupRing{AutGroup{N}}) where N
|
||||||
|
T2 = small_generating_set(RG, 2)
|
||||||
|
ℤ = Int64
|
||||||
|
Δ₂ = length(T2)*one(RG, ℤ) - RG(T2, ℤ);
|
||||||
|
|
||||||
|
Alt_N = [g for g in elements(PermutationGroup(N)) if parity(g) == 0]
|
||||||
|
elt = sum(σ(Δ₂)^2 for σ in Alt_N)
|
||||||
|
return elt
|
||||||
|
end
|
||||||
|
|
||||||
|
function Adj(RG::GroupRing{AutGroup{N}}) where N
|
||||||
|
T2 = small_generating_set(RG, 2)
|
||||||
|
|
||||||
|
ℤ = Int64
|
||||||
|
Δ₂ = length(T2)*one(RG, ℤ) - RG(T2, ℤ);
|
||||||
|
|
||||||
|
Alt_N = [g for g in elements(PermutationGroup(N)) if parity(g) == 0]
|
||||||
|
|
||||||
|
adj(σ::perm, τ::perm, i=1, j=2) = Set([σ[i], σ[j]]) ∩ Set([τ[i], τ[j]])
|
||||||
|
adj(σ::perm) = [τ for τ in Alt_N if length(adj(σ, τ)) == 1]
|
||||||
|
|
||||||
|
@time elt = sum(σ(Δ₂)*sum(τ(Δ₂) for τ in adj(σ)) for σ in Alt_N);
|
||||||
|
# return RG(elt.coeffs÷factorial(N-2)^2)
|
||||||
|
return elt
|
||||||
|
end
|
||||||
|
|
||||||
|
function Op(RG::GroupRing{AutGroup{N}}) where N
|
||||||
|
T2 = small_generating_set(RG, 2)
|
||||||
|
|
||||||
|
ℤ = Int64
|
||||||
|
Δ₂ = length(T2)*one(RG, ℤ) - RG(T2, ℤ);
|
||||||
|
|
||||||
|
Alt_N = [g for g in elements(PermutationGroup(N)) if parity(g) == 0]
|
||||||
|
|
||||||
|
adj(σ::perm, τ::perm, i=1, j=2) = Set([σ[i], σ[j]]) ∩ Set([τ[i], τ[j]])
|
||||||
|
adj(σ::perm) = [τ for τ in Alt_N if length(adj(σ, τ)) == 0]
|
||||||
|
|
||||||
|
@time elt = sum(σ(Δ₂)*sum(τ(Δ₂) for τ in adj(σ)) for σ in Alt_N);
|
||||||
|
# return RG(elt.coeffs÷factorial(N-2)^2)
|
||||||
|
return elt
|
||||||
|
end
|
||||||
|
|
||||||
|
const ELT_FILE = joinpath(dirname(PropertyT.filename(sett, :Δ)), "SqAdjOp_coeffs.jld")
|
||||||
|
const WARMSTART_FILE = PropertyT.filename(sett, :warmstart)
|
||||||
|
|
||||||
|
if isfile(PropertyT.filename(sett,:Δ)) && isfile(ELT_FILE) &&
|
||||||
|
isfile(PropertyT.filename(sett, :OrbitData))
|
||||||
|
# cached
|
||||||
|
Δ = PropertyT.loadGRElem(PropertyT.filename(sett,:Δ), sett.G)
|
||||||
|
RG = parent(Δ)
|
||||||
|
orbit_data = load(PropertyT.filename(sett, :OrbitData), "OrbitData")
|
||||||
|
sq_c, adj_c, op_c = load(ELT_FILE, "Sq", "Adj", "Op")
|
||||||
|
# elt = ELT_FILE, sett.G)
|
||||||
|
sq = GroupRingElem(sq_c, RG)
|
||||||
|
adj = GroupRingElem(adj_c, RG)
|
||||||
|
op = GroupRingElem(op_c, RG);
|
||||||
|
else
|
||||||
|
info("Compute Laplacian")
|
||||||
|
Δ = PropertyT.Laplacian(sett.S, sett.radius)
|
||||||
|
RG = parent(Δ)
|
||||||
|
|
||||||
|
info("Compute Sq, Adj, Op")
|
||||||
|
@time sq, adj, op = Sq(RG), Adj(RG), Op(RG)
|
||||||
|
|
||||||
|
PropertyT.saveGRElem(PropertyT.filename(sett, :Δ), Δ)
|
||||||
|
save(ELT_FILE, "Sq", sq.coeffs, "Adj", adj.coeffs, "Op", op.coeffs)
|
||||||
|
|
||||||
|
info("Compute OrbitData")
|
||||||
|
if !isfile(PropertyT.filename(sett, :OrbitData))
|
||||||
|
orbit_data = PropertyT.OrbitData(parent(Y), sett.autS)
|
||||||
|
save(PropertyT.filename(sett, :OrbitData), "OrbitData", orbit_data)
|
||||||
|
else
|
||||||
|
orbit_data = load(PropertyT.filename(sett, :OrbitData), "OrbitData")
|
||||||
|
end
|
||||||
|
end;
|
||||||
|
|
||||||
|
orbit_data = PropertyT.decimate(orbit_data);
|
||||||
|
|
||||||
|
elt = adj+2op;
|
||||||
|
|
||||||
|
const SOLUTION_FILE = PropertyT.filename(sett, :solution)
|
||||||
|
|
||||||
|
if !isfile(SOLUTION_FILE)
|
||||||
|
|
||||||
|
SDP_problem, varλ, varP = PropertyT.SOS_problem(elt, Δ, orbit_data; upper_bound=sett.upper_bound)
|
||||||
|
|
||||||
|
begin
|
||||||
|
using SCS
|
||||||
|
scs_solver = SCS.SCSSolver(linearsolver=SCS.Direct,
|
||||||
|
eps=sett.tol,
|
||||||
|
max_iters=args["iterations"],
|
||||||
|
alpha=1.95,
|
||||||
|
acceleration_lookback=1)
|
||||||
|
|
||||||
|
JuMP.setsolver(SDP_problem, scs_solver)
|
||||||
|
end
|
||||||
|
|
||||||
|
λ = Ps = nothing
|
||||||
|
ws = PropertyT.warmstart(sett)
|
||||||
|
|
||||||
|
# using ProgressMeter
|
||||||
|
|
||||||
|
# @showprogress "Running SDP optimization step... " for i in 1:args["repetitions"]
|
||||||
|
while true
|
||||||
|
λ, Ps, ws = PropertyT.solve(PropertyT.filename(sett, :solverlog),
|
||||||
|
SDP_problem, varλ, varP, ws);
|
||||||
|
|
||||||
|
if all((!isnan).(ws[1]))
|
||||||
|
save(WARMSTART_FILE, "warmstart", ws, "λ", λ, "Ps", Ps)
|
||||||
|
save(WARMSTART_FILE[1:end-4]*"_$(now()).jld", "warmstart", ws, "λ", λ, "Ps", Ps)
|
||||||
|
else
|
||||||
|
warn("No valid solution was saved!")
|
||||||
|
end
|
||||||
|
end
|
||||||
|
|
||||||
|
info("Reconstructing P...")
|
||||||
|
@time P = PropertyT.reconstruct(Ps, orbit_data);
|
||||||
|
save(SOLUTION_FILE, "λ", λ, "P", P)
|
||||||
|
end
|
||||||
|
|
||||||
|
P, λ = load(SOLUTION_FILE, "P", "λ")
|
||||||
|
@show λ;
|
||||||
|
|
||||||
|
@time const Q = real(sqrtm(P));
|
||||||
|
|
||||||
|
function SOS_residual(eoi::GroupRingElem, Q::Matrix)
|
||||||
|
RG = parent(eoi)
|
||||||
|
@time sos = PropertyT.compute_SOS(RG, Q);
|
||||||
|
return eoi - sos
|
||||||
|
end
|
||||||
|
|
||||||
|
info("Floating Point arithmetic:")
|
||||||
|
EOI = elt - λ*Δ
|
||||||
|
b = SOS_residual(EOI, Q)
|
||||||
|
@show norm(b, 1);
|
||||||
|
|
||||||
|
info("Interval arithmetic:")
|
||||||
|
using IntervalArithmetic
|
||||||
|
Qint = PropertyT.augIdproj(Q);
|
||||||
|
@assert all([zero(eltype(Q)) in sum(view(Qint, :, i)) for i in 1:size(Qint, 2)])
|
||||||
|
|
||||||
|
EOI_int = elt - @interval(λ)*Δ;
|
||||||
|
Q_int = PropertyT.augIdproj(Q);
|
||||||
|
@assert all([zero(eltype(Q)) in sum(view(Q_int, :, i)) for i in 1:size(Q_int, 2)])
|
||||||
|
b_int = SOS_residual(EOI_int, Q_int)
|
||||||
|
@show norm(b_int, 1);
|
||||||
|
|
||||||
|
info("λ is certified to be > ", (@interval(λ) - 2^2*norm(b_int,1)).lo)
|
29
run.jl
29
run.jl
@ -78,10 +78,31 @@ include("CPUselect.jl")
|
|||||||
include("logging.jl")
|
include("logging.jl")
|
||||||
include("main.jl")
|
include("main.jl")
|
||||||
|
|
||||||
const G = PropertyTGroups.PropertyTGroup(PARSEDARGS)
|
using SCS.SCSSolver
|
||||||
|
# using Mosek
|
||||||
|
# using CSDP
|
||||||
|
# using SDPA
|
||||||
|
|
||||||
fullpath = joinpath(name(G), string(G.args["upper-bound"]))
|
solver(tol, iterations) =
|
||||||
|
SCSSolver(linearsolver=SCS.Direct,
|
||||||
|
eps=tol, max_iters=iterations,
|
||||||
|
alpha=1.95, acceleration_lookback=1)
|
||||||
|
|
||||||
|
# Mosek.MosekSolver(
|
||||||
|
# MSK_DPAR_INTPNT_CO_TOL_REL_GAP=tol,
|
||||||
|
# MSK_IPAR_INTPNT_MAX_ITERATIONS=iterations,
|
||||||
|
# QUIET=false)
|
||||||
|
|
||||||
|
# CSDP.CSDPSolver(axtol=tol, atytol=tol, objtol=tol, minstepp=tol*10.0^-1, minstepd=tol*10.0^-1)
|
||||||
|
|
||||||
|
# SDPA.SDPASolver(epsilonStar=tol, epsilonDash=tol)
|
||||||
|
|
||||||
|
const Gr = PropertyTGroups.PropertyTGroup(PARSEDARGS)
|
||||||
|
const sett = PropertyT.Settings(Gr, PARSEDARGS,
|
||||||
|
solver(PARSEDARGS["tol"], PARSEDARGS["iterations"]))
|
||||||
|
|
||||||
|
fullpath = PropertyT.fullpath(sett)
|
||||||
isdir(fullpath) || mkpath(fullpath)
|
isdir(fullpath) || mkpath(fullpath)
|
||||||
logger=setup_logging(PropertyT.filename(fullpath, :fulllog), :fulllog)
|
setup_logging(PropertyT.filename(fullpath, :fulllog), :fulllog)
|
||||||
|
|
||||||
main(G)
|
main(Gr, sett)
|
||||||
|
Loading…
Reference in New Issue
Block a user