Benoît Legat

Sorry for the late reply, I missed the notification. You need to give a vector while here you are giving a tuple. ```julia monosY = [v^2 for v in y]...

All the names I know are "pseudo-expectations" and "Lasserre moment series"

Glad to hear that! I created the two milestones 0.0.1 and 0.1.0 to split issues that are a priority or easier with the other issues that might take more time...

It seems that the algebraic set is empty, it is defined by ``` -0.961538zr₁² - 0.961538zi₁² == 3.5 && -4.80769zr₁² - 4.80769zi₁² == -3.5 ``` which is rewritten into ```...

I agree that the fix would be to fix the newton polytope computation for empty monomial vector. We should both add a test of newton polytope for empty monomial vector...

It should return an empty monomial vector. Then it will create a 0x0 PSD matrix: https://github.com/JuliaOpt/SumOfSquares.jl/blob/fdb6adbd87e17e777b7f6abdb124bc8ddaef01fb/src/psd2x2.jl#L3 which is bridged into adding no constraint: https://github.com/JuliaOpt/SumOfSquares.jl/blob/fdb6adbd87e17e777b7f6abdb124bc8ddaef01fb/src/empty_bridge.jl#L1-L8 so it is equivalent to ignoring...

> That might be a better solution, but would add more code and overhead. Not if we use `MOI.Utilities.ModelFilter`

Indeed, maybe `MOI.set(model, MOI.Utilities.FeasibilityRelaxation(Dict(c => 2.0)))` could be `MOI.modify(model, c, MOI.Utilities.FeasibilityRelaxation(2.0))`

Good point, then https://github.com/jump-dev/MathOptInterface.jl/pull/2011 might be better

What do you mean exactly by theoretical basis ? It happens often in control that you get quadratic expressions (or bilinear if you have no square). When you are looking...