Pb depassement?

Bugs

Modérateur : xcasadmin

frederic han
Messages : 1113
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Pb depassement?

Message par frederic han » mer. janv. 21, 2009 9:21 am

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

parisse
Messages : 5114
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Message par parisse » mer. janv. 21, 2009 1:38 pm

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.

Répondre