kalmar
/
PropertyT.jl
mirror of https://github.com/kalmarek/PropertyT.jl.git
1
0
Fork 0
Code Issues Releases Wiki Activity

rework moment_matrix formulation

This commit is contained in:
Marek Kaluba 2023-10-06 14:16:49 +02:00
parent 759d4c9ed7
commit 5f54ecb515
No known key found for this signature in database
GPG Key ID: 8BF1A3855328FC15
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/sos_sdps.jl
View File

@ -8,21 +8,22 @@ function sos_problem_dual(
elt::StarAlgebras.AlgebraElement,
order_unit::StarAlgebras.AlgebraElement = zero(elt);
lower_bound = -Inf,
supp::AbstractVector{<:Integer} = axes(parent(elt).mstructure, 1),
)
@assert parent(elt) == parent(order_unit)
algebra = parent(elt)
moment_matrix = let m = algebra.mstructure
[m[-i, j] for i in axes(m, 1), j in axes(m, 2)]
end
A = parent(elt)
# 1 variable for every primal constraint
# 1 dual variable for every element of basis
# Symmetrized:
# 1 dual variable for every orbit of G acting on basis
model = Model()
JuMP.@variable(model, y[1:length(basis(algebra))])
JuMP.@variable(model, y[1:length(basis(A))])
JuMP.@constraint(model, λ_dual, dot(order_unit, y) == 1)
JuMP.@constraint(model, psd, y[moment_matrix] in PSDCone())
let mstr = A.mstructure
moment_matrix = [mstr[-i, j] for i in supp, j in supp]
JuMP.@constraint(model, psd, y[moment_matrix] in PSDCone())
end
if !isinf(lower_bound)
throw("Not Implemented yet")