From fb3b51fd6e81f563a01d17f8632dcffd5a1977b8 Mon Sep 17 00:00:00 2001 From: Marek Kaluba Date: Sun, 19 Mar 2023 23:26:44 +0100 Subject: [PATCH] fix sos_problem_dual --- src/sos_sdps.jl | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/src/sos_sdps.jl b/src/sos_sdps.jl index 4eecd09..a91b240 100644 --- a/src/sos_sdps.jl +++ b/src/sos_sdps.jl @@ -12,7 +12,7 @@ function sos_problem_dual( @assert parent(elt) == parent(order_unit) algebra = parent(elt) moment_matrix = let m = algebra.mstructure - [m[-i, j] for i in axes(m, 1) for j in axes(m, 2)] + [m[-i, j] for i in axes(m, 1), j in axes(m, 2)] end # 1 variable for every primal constraint @@ -20,16 +20,16 @@ function sos_problem_dual( # Symmetrized: # 1 dual variable for every orbit of G acting on basis model = Model() - @variable model y[1:length(basis(algebra))] - @constraint model λ_dual dot(order_unit, y) == 1 - @constraint(model, psd, y[moment_matrix] in PSDCone()) + JuMP.@variable(model, y[1:length(basis(algebra))]) + JuMP.@constraint(model, λ_dual, dot(order_unit, y) == 1) + JuMP.@constraint(model, psd, y[moment_matrix] in PSDCone()) if !isinf(lower_bound) throw("Not Implemented yet") - @variable model λ_ub_dual - @objective model Min dot(elt, y) + lower_bound * λ_ub_dual + JuMP.@variable(model, λ_ub_dual) + JuMP.@objective(model, Min, dot(elt, y) + lower_bound * λ_ub_dual) else - @objective model Min dot(elt, y) + JuMP.@objective(model, Min, dot(elt, y)) end return model