From 7f09e7c8b5199c213aaf59d29d9994e06fe1fbb4 Mon Sep 17 00:00:00 2001 From: kalmarek Date: Mon, 20 Aug 2018 03:51:51 +0200 Subject: [PATCH] instead of recomputing Delta, just pass it --- src/PropertyT.jl | 8 +++----- 1 file changed, 3 insertions(+), 5 deletions(-) diff --git a/src/PropertyT.jl b/src/PropertyT.jl index 004ea56..cc62b3f 100644 --- a/src/PropertyT.jl +++ b/src/PropertyT.jl @@ -144,19 +144,17 @@ function check_property_T(name::String, S, solver, upper_bound, tol, radius, war isapprox(eigvals(P), abs.(eigvals(P)), atol=tol) || warn("The solution matrix doesn't seem to be positive definite!") - return interpret_results(name, S, radius, λ, P) + return interpret_results(name, Δ, radius,length(S), λ, P) end -function interpret_results(name, S, radius, λ, P) +function interpret_results(name::String, Δ::GroupRingElem, radius::Integer, length_S::Integer, λ::AbstractFloat, P) - RG = GroupRing(parent(first(S)), load(filename(name, :pm), "pm")) - Δ = GroupRingElem(load(filename(name, :Δ), "Δ")[:, 1], RG) @time Q = real(sqrtm(Symmetric(P))) sgap = distance_to_cone(Δ, λ, Q, wlen=2*radius) if sgap > 0 - Kazhdan_κ = Kazhdan(sgap, length(S)) + Kazhdan_κ = Kazhdan(sgap, length_S) if Kazhdan_κ > 0 info("κ($name, S) ≥ $Kazhdan_κ: Group HAS property (T)!") return true