oops

Author: jverzani

src/polynomials/standard-basis.jl

@@ -114,7 +114,7 @@ function Base.divrem(num::P, den::Q) where {T, P <: StandardBasisPolynomial{T},
             r_coeff[i - m + j + 1] -= elem
         end
     end
-    resize!(r_coeff, min(length(r_coeff), m+1))
+    resize!(r_coeff, min(length(r_coeff), m))
 
     return ⟒(P)(q_coeff, var), ⟒(P)(r_coeff, var)
 

test/StandardBasis.jl

@@ -180,7 +180,7 @@ end
         q, r = divrem(num,den)  # expected degrees: degree(q) = degree(num)-degree(den) = 1, degree(r) = degree(den)-1 = 1
         @test num ≈ den*q+r  # true
         @test degree(q) == 1 # true
-        degree(r) <= degree(den)
+        degree(r) < degree(den)
     end
 end