Benoît Legat

It is using `default_solver` too https://github.com/JuliaPolyhedra/Polyhedra.jl/blob/485f7dfb2abb93b9032d711bc017c355796ad02e/src/linearity.jl#L35

Then that means that `p.solver` does not have presolve. There is no default solver in Polyhedra if `p.solver` is nothing so if a solver is used, the only possibility is...

Maybe the solver encounters numerical errors even if the presolve is on ? The warning is just a hint that might be useful, it does not mean that the presolve...

We are going to release JuMP v1.0 in the next few days. Any plan on doing the upgrade to MOI v0.10/JuMP v0.22 ? The upgrade to MOI v1/JuMP v0.23 is...

Thanks, could you add a test ?

Adding ```julia promote_type(Any, typeof(MP.changecoefficienttype(x + 1, Any))) ``` in https://github.com/JuliaAlgebra/MultivariatePolynomials.jl/blob/master/test/promote.jl should cover it.

Here is a one line: ```julia coefficient(leadingterm(isolate_variable(p, var))) ``` but it's not the fastest as it computes the coefficient for each degree of `var`, not just the maximal one. This...

No worry, I'll leave the issue open as it might be a nice additional feature of the package.

How do we need to treat those conjugates ? Simply as other variables or do we need to use identities like `x' * x` is `1` for othogonal matrices and...

I think it is doable. We will probably not create a type for every of the 16 combinations (commutativity, hermitianity, orthogonality and idempotency). I suggest to have 2 types (`PolyVar`...