mirror of
https://github.com/kalmarek/PropertyT.jl.git
synced 2024-11-26 17:05:27 +01:00
Constraint sum of A to minimize |b|_1
This commit is contained in:
parent
913962d9e6
commit
26f4900e82
@ -54,10 +54,10 @@ function create_SDP_problem(matrix_constraints, Δ::GroupAlgebraElement)
|
|||||||
m = JuMP.Model();
|
m = JuMP.Model();
|
||||||
JuMP.@variable(m, A[1:N, 1:N], SDP)
|
JuMP.@variable(m, A[1:N, 1:N], SDP)
|
||||||
JuMP.@SDconstraint(m, A >= 0)
|
JuMP.@SDconstraint(m, A >= 0)
|
||||||
|
JuMP.@constraint(m, sum(A[i] for i in eachindex(A)) == 0)
|
||||||
JuMP.@variable(m, κ >= 0.0)
|
JuMP.@variable(m, κ >= 0.0)
|
||||||
JuMP.@constraint(m, κ <= 0.26)
|
JuMP.@constraint(m, κ <= 0.26)
|
||||||
JuMP.@objective(m, Max, κ)
|
JuMP.@objective(m, Max, κ)
|
||||||
JuMP.@constraint(m, sum(A[i] for i in eachindex(A)) == 0)
|
|
||||||
|
|
||||||
for (pairs, δ², δ) in zip(matrix_constraints, Δ².coefficients, Δ.coefficients)
|
for (pairs, δ², δ) in zip(matrix_constraints, Δ².coefficients, Δ.coefficients)
|
||||||
JuMP.@constraint(m, sum(A[i,j] for (i,j) in pairs) == δ² - κ*δ)
|
JuMP.@constraint(m, sum(A[i,j] for (i,j) in pairs) == δ² - κ*δ)
|
||||||
|
Loading…
Reference in New Issue
Block a user