Add Support for evluation of system of polynomials and better homogen…

…ization
JuliaAlgebra · Nov 14, 2017 · 9d0809f · 9d0809f
1 parent 337278f
commit 9d0809f
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/src/config.jl b/src/config.jl
@@ -5,7 +5,7 @@ const Index = UInt16
 const Delimiter = UInt16
 const Exponent = UInt16
 
-struct PolyConfig{T}
+mutable struct PolyConfig{T}
     monomials_delimiters::Vector{Delimiter}
     monomials::Vector{Index}
     grad_monomials_delimiters::Vector{Vector{Tuple{Delimiter, Exponent}}}

diff --git a/src/poly.jl b/src/poly.jl
@@ -221,6 +221,7 @@ end
 
 (p::Polynomial)(x) = evaluate(p, x)
 
+evaluate(F::Vector{<:Polynomial}, x::AbstractVector) = map(f -> evaluate(f, x), F)
 """
     substitute(p::Polynomial, i, x)
 
@@ -349,8 +350,8 @@ end
 Makes `p` homogenous, if `ishomogenized(p)` is `true` this is just the identity.
 The homogenization variable will always be considered as the first variable of the polynomial.
 """
-function homogenize(p::Polynomial, variable::Symbol=:x0)
-    if p.homogenized
+function homogenize(p::Polynomial, variable::Symbol=:x0; respect_homogenous=true)
+    if p.homogenized || (respect_homogenous && ishomogenous(p))
         p
     else
         monomials_degree = sum(exponents(p), 1)