correct geomean dual feasibility check

jump-dev · Jul 5, 2020 · 14907a2 · 14907a2
1 parent cd8ec3b
commit 14907a2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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