Compatibility issue between SymbolicWedderburn with SumsOfSquares

[kneshat]

Symmetry reduction does not work properly while doing SOS optimization. Here is a minimal example, which is similar to an example in the SumsOfSquares document. Here, I used poly = sum(x) + sum(x.^4) instead of poly = sum(x) + sum(x.^2)

import MutableArithmetics as MA
using MultivariatePolynomials
using MultivariateBases
using DynamicPolynomials
@polyvar x[1:3]
poly = sum(x) + sum(x.^4)
using PermutationGroups
G = PermGroup([perm"(1,2,3)",perm"(1,2)"])
import CSDP
solver = CSDP.Optimizer
model = Model(solver)
@variable(model, t)
@objective(model, Max, t)
pattern = Symmetry.Pattern(G, Symmetry.VariablePermutation())
con_ref = @constraint model poly - t in SOSCone() symmetry = pattern
optimize!(model)
value(t)

@kalmarek has checked the problem and said "Indeed this seems to be the problem with SumOfSquares.jl formulation, as SymbolicWedderburn.jl correctly reduces the optimization problem"

For more info:
https://discourse.julialang.org/t/symmetry-reduction-in-sumsofsquares/115727

[blegat]

Thanks, once #355 is merged, I'll rework the symmetry reduction part and take a look that this example

[kalmarek]

@Khashayar-Neshat until then, if all you want is to solve the optimization problem and get the optimal value you may use sos_problem (or a modification of it). The actual model is created here

[kneshat]

@blegat is there any change regarding this issue? I saw #355 was merged.

[blegat]

Things were on hold while I was away for JuMP-dev and ISMP but we're working on it. We need to figure out what's the best way to do this now with the new StarAlgebras

[kneshat]

Do you happen to have any new updates? I am still not able to use symmetry reduction in SumOfSquares properly.

[blegat]

We're finishing StarAlgebras right now JuliaAlgebra/StarAlgebras.jl#63
Then we should be able to make a new release of SumOfSquares.
I will then rework the symmetry reduction in SumOfSquares