SumOfSquares.jl
SumOfSquares.jl copied to clipboard
A trivial example that returns infeasible
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)