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

@logtime Q = sqrtm(P)

Browse Source
This commit is contained in:
kalmarek 2017-11-05 19:30:11 +01:00
parent c3ea80e6a9
commit 36e5979774
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/PropertyT.jl
View File

@ -218,7 +218,7 @@ function check_property_T(name::String, S, Id, solver, upper_bound, tol, radius)
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))) @logtime logger Q = real(sqrtm(Symmetric(P)))
sgap = distance_to_positive_cone(Δ, λ, Q, 2*radius) sgap = distance_to_positive_cone(Δ, λ, Q, 2*radius)
if isa(sgap, Interval) if isa(sgap, Interval)