From 98cdeedc7aad62eaae5289c8d06c2b821eb9a5e3 Mon Sep 17 00:00:00 2001 From: kalmar Date: Fri, 31 Mar 2017 22:36:18 +0200 Subject: [PATCH] make A_sqrt_Q_aug \pm \delta a matrix of Float64 intervals --- src/checksolution.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/checksolution.jl b/src/checksolution.jl index ef8680e..3577e96 100644 --- a/src/checksolution.jl +++ b/src/checksolution.jl @@ -135,8 +135,8 @@ function check_distance_to_positive_cone(Δ::GroupAlgebraElement, κ, A; Δ_ℚ = ℚ(Δ, δ) info(logger, "Checking in interval arithmetic") - A_sqrt_ℚ_augᴵ = A_sqrt_ℚ_aug ± δ - t = @timed Interval_dist_to_Σ² = distance_to_cone(κ_ℚ, A_sqrt_ℚ_augᴵ, Δ_ℚ) + A_sqrt_ℚ_augⁱⁿᵗ = Float64.(A_sqrt_ℚ_aug) ± δ + t = @timed Interval_dist_to_Σ² = distance_to_cone(κ_ℚ, A_sqrt_ℚ_augⁱⁿᵗ, Δ_ℚ) info(logger, timed_msg(t)) info(logger, "The Augmentation-projected actual distance (to positive cone) ≥ $(@sprintf("%.10f", Interval_dist_to_Σ².lo))") info(logger, "------------------------------------------------------------")