diff --git a/property(T).jl b/property(T).jl index 10b2cc3..9ff3cd3 100644 --- a/property(T).jl +++ b/property(T).jl @@ -54,10 +54,10 @@ function create_SDP_problem(matrix_constraints, Δ::GroupAlgebraElement) m = JuMP.Model(); JuMP.@variable(m, A[1:N, 1:N], SDP) JuMP.@SDconstraint(m, A >= 0) + JuMP.@constraint(m, sum(A[i] for i in eachindex(A)) == 0) JuMP.@variable(m, κ >= 0.0) JuMP.@constraint(m, κ <= 0.26) JuMP.@objective(m, Max, κ) - JuMP.@constraint(m, sum(A[i] for i in eachindex(A)) == 0) for (pairs, δ², δ) in zip(matrix_constraints, Δ².coefficients, Δ.coefficients) JuMP.@constraint(m, sum(A[i,j] for (i,j) in pairs) == δ² - κ*δ)