Page 1 sur 1

factor -> dic euclidiene non exacte

Publié : jeu. août 07, 2014 6:21 pm
par frederic han
J'ai une reponse bizarre:

Code : Tout sélectionner

p:=8/3*a^2*b^3+1/3*a^2*b^2*d+4/3*a^2*b^2-8*a^2*b*c^2-8/3*a^2*b*c*d+16/21*a^2*b*d-a^2*c^2*d-4*a^2*c^2-1/3*a^2*c*d^2-4/3*a^2*c*d+2/21*a^2*d^2+8/21*a^2*d+120*a*b^4*d+8*a*b^4+15*a*b^3*d^2-67*a*b^3*d+180/37*a*b^3-360*a*b^2*c^2*d-18*a*b^2*c^2-120*a*b^2*c*d^2-8*a*b^2*c*d+128/7*a*b^2*d^2-15956/259*a*b^2*d+16/37*a*b^2-45*a*b*c^2*d^2+201*a*b*c^2*d-540/37*a*b*c^2-15*a*b*c*d^3+67*a*b*c*d^2-180/37*a*b*c*d+30/7*a*b*d^3-134/7*a*b*d^2+360/259*a*b*d-18*a*c^4-6*a*c^3*d+48*a*c^2*d^2+50088/259*a*c^2*d-48/37*a*c^2+16*a*c*d^3+2364/37*a*c*d^2-16/37*a*c*d-32/7*a*d^3-4728/259*a*d^2+32/259*a*d+270*b^3*c^2*d+18*b^3*c^2-288*b^2*c^2*d+72/37*b^2*c^2-810*b*c^4*d-54*b*c^4-270*b*c^3*d^2-18*b*c^3*d+540/7*b*c^2*d^2+36/7*b*c^2*d+864*c^4*d-216/37*c^4+288*c^3*d^2-72/37*c^3*d-576/7*c^2*d^2+144/259*c^2*d:;
factor(p)

donne:
"Division euclidienne non exacte, peut-être des erreurs d'arrondis...

(avec 0.9.9 .... 1.1.2)

Re: factor -> dic euclidiene non exacte

Publié : ven. août 08, 2014 7:15 am
par parisse
Il y a un probleme dans l'algo de factorisation multivariables, on peut y remedier avec ca

Code : Tout sélectionner

diff ezgcd.cc ezgcd.cc~
418,419d417
< 	if (lgcd(pcurx1x2).coord.size()>1)
< 	  return false;
mais il faudrait voir si ca n'elimine pas trop de cas (qui devviendraient infaisables ou trop longs). Il faudrait alors que je me plonge plus precisement dans cet exemple pour voir pourquoi ca ne marche pas...