feat: support GCD for non-concrete polynomial types

Author: AayushSabharwal

src/gcd.jl

@@ -306,11 +306,9 @@ function MA.operate(
     defl,
 )
     return polynomial(
-        map(terms(p)) do t
-            return term(
-                coefficient(t),
-                MA.operate(op, monomial(t), shift, defl),
-            )
+        coefficients(p),
+        map(monomials(p)) do m
+            return MA.operate(op, m, shift, defl)
         end,
     )
 end

test/commutative/gcd.jl

@@ -291,4 +291,11 @@ end
             end
         end
     end
+    @testset "Non-concrete polynomial GCD" begin
+        Mod.@polyvar p q
+        p1 = MP.polynomial(p^2 + q^2, Number)
+        p2 = MP.polynomial(p * q, Number)
+        g = gcd(p1, p2)
+        @test isone(g)
+    end
 end