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

logging changes

Browse Source
This commit is contained in:
kalmarek 2019-01-14 17:46:13 +01:00
parent af009547b5
commit 80ac135f59
2 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/1712.07167.jl
View File

@ -138,6 +138,7 @@ function distance_to_positive_cone(Δ::GroupRingElem, λ, Q; R::Int=2)
@info("λ = $λ")
eoi = Δ^2-λ*Δ
@info("Computing sum of squares decomposition...")
@time residual = eoi - compute_SOS(parent(eoi), augIdproj(Q))
@info("ɛ(Δ² - λΔ - ∑ξᵢ*ξᵢ) ≈ $(@sprintf("%.10f", aug(residual)))")
L1_norm = norm(residual,1)
@ -163,7 +164,8 @@ function distance_to_positive_cone(Δ::GroupRingElem, λ, Q; R::Int=2)
@info("Checking that sum of every column contains 0.0... ")
@info((check ? "They do." : "FAILED!"))
check || @warn("The following numbers are meaningless!")
@info("Computing sum of squares decomposition...")
@time residual = eoi - compute_SOS(parent(eoi), Q)
@info("ɛ(Δ² - λΔ - ∑ξᵢ*ξᵢ) ∈ $(aug(residual))")
L1_norm = norm(residual,1)

3
src/sos_sdps.jl
View File

@ -136,8 +136,7 @@ function reconstruct(Ps::Vector{M},
transfP = [dims[π].*Uπs[π]*Ps[π]*Uπs[π]' for π in 1:lU]
tmp = [zeros(Float64, size(first(transfP))) for _ in 1:lU]
@time Threads.@threads for π in 1:lU
Threads.@threads for π in 1:lU
tmp[π] = perm_avg(tmp[π], transfP[π], values(preps))
end