sage: K.<a> = NumberField(x^2 + 1)
sage: R.<y, z> = PolynomialRing(K)
sage: f = 2*y^2 + 2*z^2
sage: F = f.factor(); F
2 * (y + (-a)*z) * (y + a*z)
sage: F.unit_part()
1
Shouldn't the unit part be 2? It seems to be listing 2 as a bona fide
factor.
(This was reported by Genya Zaytman.)
david
John