Benoît Legat

Do you get the error simply by doing `using SymbolicWedderburn` ?

I don't know if it's worth debugging nightly if it's not yet a release candidate

The issue seems that the resulting gram matrix is a block diagonal gram matrix with no block hence inference fails to infer the element type.

Thanks for the suggestion but that I'll have to close :laughing:

Good point, I never thought of this. The polynomials `σ1` indeed won't have the right form and hence when the monomials of the polynomial `σ0` are computed using newton polytope,...

> is it possible to have the domain of the SOS constraint printed out? Yes, that shouldn't be too hard Maybe try increasing the `maxdegree` in the constraint keyword argument...

I'm not sure to understand, `hcat` does not turn sparse matrices into dense ones: ```julia julia> using SparseArrays julia> A = spzeros(2, 2) 2×2 SparseMatrixCSC{Float64, Int64} with 0 stored entries:...

I see, I guess the use of this `hcat` is to turn vectors to matrices so we can just create a `vector_to_matrix` function that calls `hcat` on `AbstractVector` and does...

This is an old syntax from a previous release I think, try `hrep(A, b)` instead. See https://juliapolyhedra.github.io/Polyhedra.jl/stable/representation/#Representation

Thanks for reporting this, the noncommutative part does not have enough tests unfortunately, especially with constraints