Bill Allombert on Tue, 12 Nov 2002 22:47:30 +0100 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: bug in polredabs? |
On Mon, Oct 14, 2002 at 10:32:25AM -0700, Carl R. Witty wrote: > Version: Debian packages version 2.1.4-1 (x86 linux). > > There is an example in the documentation for polcompositum that > demonstrates taking the result of polcompositum and calling polredabs > on it, and using subst on the result of that. After the last line of > the example, you can type "a^5 - 5" and get 0. > > I'm trying to do exactly the same thing, except with different > polynomials. Here's my example: > > ? z = polcompositum(x^2-3, 12*x^2-1, 1)[1]; > ? pol = z[1] > %2 = 12*x^2 - 49 > ? a = z[2]; a^2-3 > %3 = 0 > ? z = polredabs(pol, 1); > ? pol = z[1] > %5 = x^2 - 3 > ? a = subst(a.pol, x, z[2]) > %6 = Mod(-6*x, x^2 - 3) > ? a^2 - 3 > %7 = Mod(105, x^2 - 3) > > I was hoping, based on the polcompositum example, that the last line > would return 0; this would have worked if the second element returned > by polredabs were 7/6*x instead of 7*x. > > I don't know if this is a bug, because the polredabs documentation > says the second element returned by polredabs(..., 1) is > "an element expressed on a root of the polynomial $P$, > whose minimal polynomial is equal to $x$", and I don't know what > "minimal polynomial" means. Karim have fixed this bug in the stable branch of the CVS, so the fix will appear in GP 2.1.5. Thanks for your report ! Cheers, Bill.