SumOfSquares.jl icon indicating copy to clipboard operation
SumOfSquares.jl copied to clipboard

A trivial example that returns infeasible

Open AlexandreAmice opened this issue 1 year ago • 2 comments

The following univariate example returns infeasible even though it should be a very simple problem. For the purposes of debugging, is it possible to have the domain of the SOS constraint printed out?

using DynamicPolynomials
using SumOfSquares
using MosekTools
@polyvar x
S = @set -1 <= x && x <= 1

import Mosek
model = SOSModel(Mosek.Optimizer)
set_silent(model)
con_ref = @constraint(model, x^2 <= 1.1, domain = S)
optimize!(model)
solution_summary(model)

AlexandreAmice avatar Mar 29 '23 16:03 AlexandreAmice