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

tidy ValidatedNumerics and move functions

Browse Source
This commit is contained in:
kalmar 2017-03-14 23:34:47 +01:00
parent dc528cc97c
commit 52a68e8429
1 changed files with 22 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

36
src/checksolution.jl
View File

@ -1,7 +1,10 @@
using ProgressMeter using ProgressMeter
using ValidatedNumerics
import Base: rationalize import Base: rationalize
using ValidatedNumerics
setrounding(Interval, :narrow)
setdisplay(:standard)
function EOI{T<:Number}(Δ::GroupAlgebraElement{T}, κ::T) function EOI{T<:Number}(Δ::GroupAlgebraElement{T}, κ::T)
return Δ*Δ - κ*Δ return Δ*Δ - κ*Δ
end end
@ -41,6 +44,22 @@ function correct_to_augmentation_ideal{T<:Rational}(sqrt_matrix::Array{T,2})
return sqrt_corrected return sqrt_corrected
end end
import ValidatedNumerics.±
function (±){T<:Number}(X::AbstractArray{T}, tol::Real)
r{T}(x::T) = (x == zero(T)? @biginterval(0) : x ± tol)
return r.(X)
end
(±)(X::GroupAlgebraElement, tol::Real) = GroupAlgebraElement(X.coefficients ± tol, X.product_matrix)
function Base.rationalize{T<:Integer, S<:Real}(::Type{T},
X::AbstractArray{S}; tol::Real=eps(eltype(X)))
r(x) = rationalize(T, x, tol=tol)
return r.(X)
end
ℚ(x, tol::Real) = rationalize(BigInt, x, tol=tol)
function check_solution{T<:Number}(κ::T, sqrt_matrix::Array{T,2}, Δ::GroupAlgebraElement{T}; verbose=true, augmented=false) function check_solution{T<:Number}(κ::T, sqrt_matrix::Array{T,2}, Δ::GroupAlgebraElement{T}; verbose=true, augmented=false)
result = compute_SOS(sqrt_matrix, Δ) result = compute_SOS(sqrt_matrix, Δ)
if augmented if augmented
@ -67,21 +86,10 @@ function check_solution{T<:Number}(κ::T, sqrt_matrix::Array{T,2}, Δ::GroupAlge
return distance_to_cone return distance_to_cone
end end
import ValidatedNumerics.±
function (±)(X::AbstractArray, tol::Real)
r{T}(x::T) = ( x==zero(T) ? @interval(x) : x ± tol)
return r.(X)
end end
(±)(X::GroupAlgebraElement, tol::Real) = GroupAlgebraElement(X.coefficients ± tol, X.product_matrix)
function Base.rationalize{T<:Integer, S<:Real}(::Type{T},
X::AbstractArray{S}; tol::Real=eps(eltype(X)))
r(x) = rationalize(T, x, tol=tol)
return r.(X)
end
ℚ(x, tol::Real) = rationalize(BigInt, x, tol=tol)
function ℚ_distance_to_positive_cone(Δ::GroupAlgebraElement, κ, A; function ℚ_distance_to_positive_cone(Δ::GroupAlgebraElement, κ, A;
tol=1e-7, verbose=true, rational=false) tol=1e-7, verbose=true, rational=false)