Polynomials.jl
Polynomials.jl copied to clipboard
Polynomial manipulations in Julia
```julia julia> using SnoopCompileCore julia> invalidations = @snoopr using Polynomials; julia> using SnoopCompile julia> trees = invalidation_trees(invalidations); julia> length(trees) 6 julia> trees[end].method convert(::Type{T}, p::P) where {T, P trees[end-1].method promote_rule(::Type{
I want to write more algebraic style code , so do this here is the code for eigen polynomials with pluto note 1. get algebraic expression with symbolcs.jl ```julia using...
I tried to check my packages with [Aqua.jl](https://github.com/JuliaTesting/Aqua.jl) (just for fun). Performing the ambiguity test strictly for my package I obtained no ambiguities in my code ``` julia> Aqua.test_ambiguities(MatrixPencils,recursive=false) Test...
This polynomial shows that multiplication for different polynomial types can vary significantly: ``` julia> rts = [1.0, sqrt(2), sqrt(3)]; julia> ls = [2, 3, 4]; julia> x, xx = variable(Polynomial{Float64}),...
CUDA
In https://discourse.julialang.org/t/using-polynomials-jl-with-cuda-jl/80580 we find the following hangs: ``` using CUDA, Polynomials a = CuVector([1.0f0, 1.0f0, 1.0f0]) Polynomial(a) ```
The fallback is producing results that are not mathematically correct: ```julia julia> using Polynomials julia> using LinearAlgebra julia> p = Polynomial([0,1]) julia> p ⋅ p 1 ```
A bug?
I started some cleaning of my tests for rational transfer functions. The following fails and was commented out for `RationalTransferFunctions`. It seems however, it doesn't work for `RationalFunctions` too. ```...
This is possibly a silly question, but is related to the previous issue [#308](https://github.com/JuliaMath/Polynomials.jl/issues/308) . Admitedly we have an empty polynomial matrix, created for example with the command: ``` julia>...
I would appreciate very much if you could reveal how you enforce the polynomial type of elements in a polynomial matrix construct, as for example in the following vertical concatenation...
_Issue discovered by @AndiMD originally in [#295](https://github.com/JuliaMath/Polynomials.jl/issues/295#issuecomment-737641577):_ Possibly related, [Unitful.jl](https://github.com/PainterQubits/Unitful.jl) compatibility: ```julia julia> using Polynomials, Unitful julia> fit([1.0,2.0,3.0,4.0].*u"m", [6.0,4.0,7.0,7.5].*u"m",2) ERROR: DimensionError ... ``` --- Response from @jverzani on 5th Dec.,...