Benoît Legat

Fix rrule for * and add support for constant operations

2

I am a bit unsure what to do here. I think we should define `rrule` based on `LinearAlgebra.dot` so we should remove the definition of the rule for `differentiate` and...

https://github.com/JuliaAlgebra/MultivariatePolynomials.jl/pull/201 rebased on master

```julia julia> @polyvar x[1:3] (PolyVar{true}[x₁, x₂, x₃],) julia> 1.0x ≈ 1.0x ERROR: MethodError: no method matching rtoldefault(::Type{Term{true, Float64}}, ::Type{Term{true, Float64}}, ::Int64) Closest candidates are: rtoldefault(::Union{Type{T}, T}, ::Union{Type{S}, S}, ::Real) where...

- [x] https://github.com/JuliaAlgebra/DynamicPolynomials.jl/pull/105 - [x] Bump DynamicPolynomials See https://github.com/JuliaAlgebra/DynamicPolynomials.jl/issues/104

Fix differentiation of zero polynomial

1

Closes https://github.com/JuliaAlgebra/MultivariatePolynomials.jl/issues/143

Simple simple and inefficient implementation of the MultivariatePolynomials API I wrote some time ago. If we make this work, this could remove the dependency on DynamicPolynomials and TypedPolynomials for testing.

``` ERROR: StackOverflowError: Stacktrace: [1] -(p1::SOSDecomposition{Float64, Polynomial{true, Float64}, Float64}, p2::SOSDecomposition{Float64, Polynomial{true, Float64}, Float64}) (repeats 79984 times) @ MultivariatePolynomials ~/.julia/packages/MultivariatePolynomials/kvnmd/src/operators.jl:32 ```

`(x - x) / 0` does not produce any error and is not NaN as it tries to divide all term coefficients but there is no term.

Add the test added in https://github.com/JuliaAlgebra/TypedPolynomials.jl/pull/59 to the differentiate tests

Safe evaluation

4

The number 1 issue with DP is with the following: ```julia @polyvar x y p = x + y p(x => 1) # Garbage ``` What prevents us from showing...