Page 1 sur 1

Pb depassement?

Publié : mer. janv. 21, 2009 9:21 am
par frederic han
J'ai un pb probablement de depassement?
(en mode maple)
P:=x^(2^30)-x
Q:=x^2-x
Gcd(P,Q) mod 2

fait planter xcas avec pleins de
*** segmentation fault: bug in PARI or calling program.
*** segmentation fault: bug in PARI or calling program.
*** segmentation fault: bug in PARI or calling program.
*** segmentation fault: bug in PARI or calling program.
*** segmentation fault: bug in PARI or calling program.

et
P:=x^(2^32)-x
Q:=x^2-x
Gcd(P,Q) mod 2

ne fait plus planter xcas mais retourne 1 immediatement.

a la limite, c'est plutot le 1 qui est plus genant que l'arret brutal?

Fred

Publié : mer. janv. 21, 2009 1:38 pm
par parisse
Oui, c'est un probleme de depassement. Les exposants des polynomes sont stockes sur des entiers signes 16 ou 32 bits selon la version de xcas et il n'y a effectivement a l'heure actuelle pas de detection d'overflow. Il faut donc imperativement construire x^(2^n) mod x^2-x,2 avec powermod.
De toutes facons, meme en l'absence de depassement, le cout de calcul serait prohibitif, car la 1ere division modulaire se ferait avec un polynome dense de 2^n+1 coefficients.