kalmar/PropertyT.jl
1
0
Fork 0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-11-27 01:10:28 +01:00
Code Issues Releases Wiki Activity

indentation

Browse Source
This commit is contained in:
kalmar 2017-08-04 20:33:21 +02:00
parent 9ed8679c95
commit 845d7733a4
1 changed files with 11 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
src/PropertyT.jl
View File

@ -193,24 +193,24 @@ function check_property_T(name::String, S, Id, solver, upper_bound, tol, radius)
if λ > 0 if λ > 0
isapprox(eigvals(P), abs(eigvals(P)), atol=tol) || isapprox(eigvals(P), abs(eigvals(P)), atol=tol) ||
warn("The solution matrix doesn't seem to be positive definite!") warn("The solution matrix doesn't seem to be positive definite!")
# @assert P == Symmetric(P) # @assert P == Symmetric(P)
Q = real(sqrtm(Symmetric(P))) Q = real(sqrtm(Symmetric(P)))
sgap = check_distance_to_positive_cone(Δ, λ, Q, 2*radius, tol=tol, rational=false) sgap = check_distance_to_positive_cone(Δ, λ, Q, 2*radius, tol=tol)
if isa(sgap, Interval) if isa(sgap, Interval)
sgap = sgap.lo sgap = sgap.lo
end end
if sgap > 0 if sgap > 0
info(logger, "λ ≥ $(Float64(trunc(sgap,12)))") info(logger, "λ ≥ $(Float64(trunc(sgap,12)))")
Kazhdan_κ = Kazhdan_from_sgap(sgap, length(S)) Kazhdan_κ = Kazhdan_from_sgap(sgap, length(S))
Kazhdan_κ = Float64(trunc(Kazhdan_κ, 12)) Kazhdan_κ = Float64(trunc(Kazhdan_κ, 12))
info(logger, "κ($name, S) ≥ $Kazhdan_κ: Group HAS property (T)!") info(logger, "κ($name, S) ≥ $Kazhdan_κ: Group HAS property (T)!")
return true return true
else else
sgap = Float64(trunc(sgap, 12)) sgap = Float64(trunc(sgap, 12))
info(logger, "λ($name, S) ≥ $sgap: Group may NOT HAVE property (T)!") info(logger, "λ($name, S) ≥ $sgap: Group may NOT HAVE property (T)!")
return false return false
end end
end end
info(logger, "κ($name, S) ≥ $λ < 0: Tells us nothing about property (T)") info(logger, "κ($name, S) ≥ $λ < 0: Tells us nothing about property (T)")