Page 1 sur 1

gcd

Publié : mer. nov. 20, 2013 1:25 pm
par frederic han
Salut,
avec la 1.1.1-26 amd 64 mais aussi les sources d'aujourd'hui j'ai le plantage dans le calcul du:

Code : Tout sélectionner

m:=((t*x^4-2*t*x^2*y^2-t*x^2*y*z+3*t*x^2*z^2+4*x^5+4*x^3*y*z-4*x^2*y^3+2*x^2*z^3+x*y^2*z^2)*(5*t*x^4-4*t*x^2*y^2+4*t*x^2*y*z+6*t*x^2*z^2-4*t*x*y^3+2*t*x*z^3+4*t*y^4-12*t*y^2*z^2+9*t*z^4+4*x^5-8*x^3*y^2+12*x^3*z^2-4*x^2*y^3+2*x^2*z^3-8*x*y^3*z-4*x*y^2*z^2+12*x*y*z^3+8*y^5+4*y^4*z-12*y^3*z^2-4*y^2*z^3-2*y*z^4+6*z^5)-(t*x^3*y-2*t*x*y^3-t*x*y^2*z+3*t*x*y*z^2+4*x^4*y+4*x^2*y^2*z-4*x*y^4+2*x*y*z^3+y^3*z^2)*(t*x^3*z-2*t*x*y^2*z-t*x*y*z^2+3*t*x*z^3+4*x^4*z+4*x^2*y*z^2-4*x*y^3*z+2*x*z^4+y^2*z^3))*(t*x^4-2*t*x^2*y^2-t*x^2*y*z+3*t*x^2*z^2+4*x^5+4*x^3*y*z-4*x^2*y^3+2*x^2*z^3+x*y^2*z^2),((t*x^4-2*t*x^2*y^2-t*x^2*y*z+3*t*x^2*z^2+4*x^5+4*x^3*y*z-4*x^2*y^3+2*x^2*z^3+x*y^2*z^2)*(5*t*x^4-4*t*x^2*y^2+4*t*x^2*y*z+6*t*x^2*z^2-4*t*x*y^3+2*t*x*z^3+4*t*y^4-12*t*y^2*z^2+9*t*z^4+4*x^5-8*x^3*y^2+12*x^3*z^2-4*x^2*y^3+2*x^2*z^3-8*x*y^3*z-4*x*y^2*z^2+12*x*y*z^3+8*y^5+4*y^4*z-12*y^3*z^2-4*y^2*z^3-2*y*z^4+6*z^5)-(t*x^3*y-2*t*x*y^3-t*x*y^2*z+3*t*x*y*z^2+4*x^4*y+4*x^2*y^2*z-4*x*y^4+2*x*y*z^3+y^3*z^2)*(t*x^3*z-2*t*x*y^2*z-t*x*y*z^2+3*t*x*z^3+4*x^4*z+4*x^2*y*z^2-4*x*y^3*z+2*x*z^4+y^2*z^3))*(t*x^3*y-2*t*x*y^3-t*x*y^2*z+3*t*x*y*z^2+4*x^4*y+4*x^2*y^2*z-4*x*y^4+2*x*y*z^3+y^3*z^2);
gcd(m);
idem avec modgcd

Bon courage!

Frederic

Re: gcd

Publié : jeu. nov. 21, 2013 9:24 am
par parisse
Ca semble marcher chez moi (mais j'ai fait hier apres-midi d'autres changements que juste le diff indique dans le post sur factor).

Re: gcd

Publié : jeu. nov. 21, 2013 5:53 pm
par frederic han
Salut,

Je l'ai encore avec le paquet 1.1.0-27 et aussi les sources du 21. (ca plante aussi avec un vieux 1.0.0 32 bits)

a suivre donc
...
Fred

Re: gcd

Publié : jeu. nov. 21, 2013 6:43 pm
par parisse
en effet, c'est sur mac que ca marchait (testé en 32 bits et en 64 mais compilé avec -g). Je regarde demain si j'ai le temps...

Re: gcd

Publié : ven. nov. 22, 2013 9:42 am
par parisse
Bon, une autre optimisation qui s'avere incorrecte, je la commente pour le moment.

Code : Tout sélectionner

diff poly.h poly.h~
325c325
<     if (false &&
---
>     if (//false &&