Benoît Legat

It's not implemented for DynamicPolynomials yet, see the discussion in https://github.com/JuliaAlgebra/MultivariatePolynomials.jl/pull/230 We would need to do something very similar to https://github.com/JuliaAlgebra/DynamicPolynomials.jl/blob/005621b454915428fa1a4c6ad72f90dc42b47f6f/src/diff.jl#L1-L10

> It is a bit weird that this causes a `StackOverflowError` as opposed to some error signaling that the function is not implemented, which was what confused me. Yes, it's...

Agreed, this would be a nice feature to have

> In my mind, a @polyvar defined with extra dimensions can be treated just like a tensor, to be indexed at will, but this seems not the case. Yes, but...

Thanks for reporting these. I would find the root cause of what outputed `Polynomial{true,Any}`. The coefficient type of a polynomial is not supposed to be `Any`, e.g., we should be...

You need `reduce` ```julia julia> a = [[1, 2], [3, 4]] 2-element Vector{Vector{Int64}}: [1, 2] [3, 4] julia> hcat(a) 2×1 Matrix{Vector{Int64}}: [1, 2] [3, 4] julia> reduce(hcat, a) 2×2 Matrix{Int64}:...

The first matrix has `eltype` of `Any` which is making the multiplication get confused. Try converting it to an `Array{polynomial_type(rot_mat_flat_var_original[1], Float64)}` before you do the multiplication.

Yes, I would definitely not incorporate whether it's a real, imaginary part etc... as type as it would disallow mixing different ones in a polynomial etc... About whether there should...

@projekter I forgot about this PR, we should merge it. I'll try to take a look next week, remind me if I don't

I'm tempted to say that we don't want to define `in` for the same reason we do not define `isless`. It does not make sense to ask whether a symbolic...