factor -> dic euclidiene non exacte

Utilisation de Xcas

Modérateur : xcasadmin

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

factor -> dic euclidiene non exacte

Message par frederic han » jeu. août 07, 2014 6:21 pm

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)

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

Re: factor -> dic euclidiene non exacte

Message par parisse » ven. août 08, 2014 7:15 am

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...

Répondre