Karim Belabas on Fri, 05 Dec 2014 17:17:31 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: Problem with Flx_ffisom function


Hi,

* Jean-Pierre Flori [2014-12-05 17:03]:
> The following programs errors out claiming that the produced
> irreducible polynomial is not anymore
> (***   not an irreducible polynomial in FpX_ffintersect):
> 
> #include <stdlib.h>
> #include <stdio.h>
> #include <time.h>
> 
> #include <pari/pari.h>
> 
> int main (int argc, char **argv)
> {
>         int i;
>         unsigned long l;
>         long n, v;
>         GEN p, f, r;
> 
>         pari_init(1<<30, 1UL<<20);
> 
>         v = 0L;
> 
>         l = unextprime(1UL<<20);
>         p = stoi(l);
>         n = 179L;
>         f = init_Fq(p, n, v);
>         i = Flx_is_irred(ZX_to_Flx(f, l), l);
                           ^^^^^^^^^^^^^^^
>         pari_printf("%d\n", i);
> 
>         t = clock();
          ^^^^^^^^^^^^  [ I'd declare t ... ]
> /*
>         The following works fine.
> */
>         r = FpX_ffisom(f, f, p);
>         pari_printf("%Ps\n", r);
> /*
>         This errors out
> */
>         r = Flx_ffisom(f, f, l);

I'd apply Flx_ffisom to a Flx, as correctly done above.
Not to the original FpX.

After

  f = ZX_to_Flx(f,l)

(and declaring clock_t t) the code above works fine.

> Thanks in advance for the support!

You're welcome :-)

    K.B.
--
Karim Belabas, IMB (UMR 5251)  Tel: (+33) (0)5 40 00 26 17
Universite de Bordeaux         Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux1.fr/~kbelabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux1.fr/  [PARI/GP]
`