mirror of
https://github.com/kalmarek/PropertyT.jl.git
synced 2024-12-25 18:25:30 +01:00
remove old functions
This commit is contained in:
parent
62ba245c89
commit
b0f273ce45
70
SL3Z.jl
70
SL3Z.jl
@ -96,84 +96,14 @@ function prepare_Δ_sdp_constraints(identity, S)
|
||||
return Δ, sdp_constraints
|
||||
end
|
||||
|
||||
function load_Δ_sdp_constraints(name::String;cached=true)
|
||||
pm_filename = "$name.product_matrix.jld"
|
||||
Δ_coeff_filename = "$name.delta.coefficients.jld"
|
||||
f₁ = isfile(pm_filename)
|
||||
f₂ = isfile(Δ_coeff_filename)
|
||||
if cached && f₁ && f₂
|
||||
println("Loading precomputed pm, Δ, sdp_constraints...")
|
||||
product_matrix = load(pm_filename, "pm")
|
||||
L = load(Δ_coeff_filename, "Δ")[:, 1]
|
||||
Δ = GroupAlgebraElement(L, Array{Int,2}(product_matrix))
|
||||
sdp_constraints = constraints_from_pm(product_matrix)
|
||||
else
|
||||
println("Computing pm, Δ, sdp_constraints...")
|
||||
ID = eye(Int, 3)
|
||||
S = SL₃ℤ_generatingset()
|
||||
Δ, sdp_constraints = prepare_Δ_sdp_constraints(ID, S)
|
||||
|
||||
save(pm_filename, "pm", Δ.product_matrix)
|
||||
save(Δ_coeff_filename, "Δ", Δ.coefficients)
|
||||
|
||||
end
|
||||
return Δ, sdp_constraints
|
||||
end
|
||||
|
||||
|
||||
function compute_κ_A(name::String, Δ, sdp_constraints;
|
||||
cached = true,
|
||||
tol = 1e-7,
|
||||
verbose = false,
|
||||
# solver = MosekSolver(INTPNT_CO_TOL_REL_GAP=tol, QUIET=!verbose))
|
||||
solver = SCSSolver(eps=tol, max_iters=20000, cg_rate=3, verbose=verbose))
|
||||
|
||||
f₁ = isfile("$name.kappa")
|
||||
f₂ = isfile("$name.SDPmatrixA")
|
||||
|
||||
if cached && f₁ && f₂
|
||||
println("Loading precomputed κ, A...")
|
||||
A = readdlm("$name.SDPmatrixA")
|
||||
κ = readdlm("$name.kappa")[1]
|
||||
else
|
||||
println("Solving SDP problem maximizing κ...")
|
||||
κ, A = solve_SDP(sdp_constraints, Δ, solver, verbose=verbose)
|
||||
# writedlm("$name.kappa", kappa)
|
||||
# writedlm("$name.SDPmatrixA", A)
|
||||
end
|
||||
return κ, A
|
||||
end
|
||||
|
||||
function main()
|
||||
const NAME = "SL3Z"
|
||||
const VERBOSE = true
|
||||
const TOL=1e-7
|
||||
const Δ, sdp_constraints = load_Δ_sdp_constraints(NAME)
|
||||
const κ, A = compute_κ_A(NAME, Δ, sdp_constraints, cached=false, verbose=VERBOSE)
|
||||
|
||||
if maximum(A) < 1e-2
|
||||
warn("Solver might not solved the problem successfully and the positive solution is due to floating-point error, proceeding anyway...")
|
||||
end
|
||||
|
||||
if κ > 0
|
||||
@assert A == Symmetric(A)
|
||||
const A_sqrt = real(sqrtm(A))
|
||||
|
||||
T = ℚ_distance_to_positive_cone(Δ, κ, A, tol=TOL, verbose=VERBOSE)
|
||||
|
||||
if T < 0
|
||||
println("$NAME HAS property (T)!")
|
||||
else
|
||||
println("$NAME may NOT HAVE property (T)!")
|
||||
end
|
||||
|
||||
else
|
||||
println("$κ < 0: $NAME may NOT HAVE property (T)!")
|
||||
end
|
||||
end
|
||||
|
||||
@everywhere push!(LOAD_PATH, "./")
|
||||
using GroupAlgebras
|
||||
include("property(T).jl")
|
||||
|
||||
main()
|
||||
|
Loading…
Reference in New Issue
Block a user