From 14907a2f77162160d5f4b0b25758c8511dbcbc1a Mon Sep 17 00:00:00 2001 From: lkapelevich Date: Sun, 5 Jul 2020 16:58:59 -0400 Subject: [PATCH] correct geomean dual feasibility check --- src/Cones/hypogeomean.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/Cones/hypogeomean.jl b/src/Cones/hypogeomean.jl index 7a7a91e6f..f8d59f000 100644 --- a/src/Cones/hypogeomean.jl +++ b/src/Cones/hypogeomean.jl @@ -113,7 +113,7 @@ function is_dual_feas(cone::HypoGeomean{T}) where {T} @views w = cone.dual_point[2:end] alpha = cone.alpha if u < -eps(T) && all(>(eps(T)), w) - @inbounds dual_wprodu = exp(sum(alpha[i] * log(w[i] / alpha[i]) for i in eachindex(alpha))) - u + @inbounds dual_wprodu = exp(sum(alpha[i] * log(w[i] / alpha[i]) for i in eachindex(alpha))) + u return (dual_wprodu > eps(T)) end return false