DynamicPolynomials.jl
DynamicPolynomials.jl copied to clipboard
Multivariate polynomials implementation of commutative and non-commutative variables
This is in conjuction with https://github.com/JuliaAlgebra/MultivariatePolynomials.jl/pull/213 and provides the DynamicPolynomials implementation of the complex-valued variable functionality.
This addresses #20 coming from the recent observation that in distributed code (with `Distributed`) `gensym` doesn't guarantee unique variable ids anymore. I used the following script to benchmark: ```julia using...
I am using the following version `[7c1d4256] DynamicPolynomials v0.5.5` in my project. I wanted to compute the antiderivative of a polynomial. The `Examples` sections show the following: Examples –––––––– p...
I work with some massive polynomials (close to 4000 terms) and when calling `gcd` I get the following error ``` ERROR: LoadError: DivideError: integer division error Stacktrace: [1] div @...
Hi, Thanks for implementing monomial orderings, I think it is a nice feature. A question: is it possible to query the ordering on existing monomials? Currently, I have ```julia using...
```julia using DynamicPolynomials @polyvar x a = x^2 b = x # This runs without error randn(2, 2) * [a, b] # This throws StackOverflowError randn(2, 2) * [p for...
```julia using DynamicPolynomials @polyvar x # Passes @assert 1.2x ≈ 1.2x # StackOverflowError x ≈ x ``` The infinite recursion seems to be triggered by: ```julia promote(x, x) == (x,...
It would be nice if DP automatically promoted coefficients to big nums when they are too big to fit in Int. See ```julia julia> @polyvar x y (x, y) julia>...
Hi, Thanks for the great library! I'm wondering how I can operate on matrices or vectors of polynomials. I have tried this: `@polyvar rot_mat_flat_var_original[1:num_internal_bodies, 6] ` creating what (I think)...
``` (code) pkg> st Status `.../Project.toml` [7c1d4256] DynamicPolynomials v0.4.6 [102ac46a] MultivariatePolynomials v0.4.7 julia> using DynamicPolynomials julia> @ncpolyvar x[1:3] (PolyVar{false}[x₁, x₂, x₃],) julia> m1 = x[2]*x[1] x₂x₁ julia> m2 = monomials(x,...