kalmar/PropertyT.jl
1
0
Fork 0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-12-26 02:30:29 +01:00
Code Issues Releases Wiki Activity

@info "checking in (...) arithmetic" before the actual check

Browse Source
This commit is contained in:
Marek Kaluba 2024-02-19 18:48:28 +01:00
parent 0f3e4311b3
commit 083e1b4735
No known key found for this signature in database
GPG Key ID: 8BF1A3855328FC15
1 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/certify.jl
View File

@ -118,10 +118,9 @@ function certify_solution(
!augmented && SA.aug(elt) == SA.aug(orderunit) == 0 !augmented && SA.aug(elt) == SA.aug(orderunit) == 0
Q = should_we_augment ? augment_columns!(Q) : Q Q = should_we_augment ? augment_columns!(Q) : Q
@info "Checking in $(eltype(Q)) arithmetic with" λ
@time sos = compute_sos(parent(elt), Q; augmented = augmented) @time sos = compute_sos(parent(elt), Q; augmented = augmented)
@info "Checking in $(eltype(sos)) arithmetic with" λ
λ_flpoint = sufficient_λ(elt, orderunit, λ, sos; halfradius = halfradius) λ_flpoint = sufficient_λ(elt, orderunit, λ, sos; halfradius = halfradius)
if λ_flpoint ≤ 0 if λ_flpoint ≤ 0
@ -130,6 +129,7 @@ function certify_solution(
λ_int = IntervalArithmetic.interval(λ) λ_int = IntervalArithmetic.interval(λ)
Q_int = IntervalMatrices.IntervalMatrix(IntervalArithmetic.interval.(Q)) Q_int = IntervalMatrices.IntervalMatrix(IntervalArithmetic.interval.(Q))
@info "Checking in $(eltype(Q_int)) arithmetic with" λ_int
check, sos_int = @time if should_we_augment check, sos_int = @time if should_we_augment
@info("Projecting columns of Q to the augmentation ideal...") @info("Projecting columns of Q to the augmentation ideal...")
@ -145,8 +145,6 @@ function certify_solution(
true, compute_sos(parent(elt), Q_int; augmented = augmented) true, compute_sos(parent(elt), Q_int; augmented = augmented)
end end
@info "Checking in $(eltype(sos_int)) arithmetic with" λ_int
λ_certified = λ_certified =
sufficient_λ(elt, orderunit, λ_int, sos_int; halfradius = halfradius) sufficient_λ(elt, orderunit, λ_int, sos_int; halfradius = halfradius)