Page 1 sur 1
x^12-1
Publié : sam. août 05, 2017 8:19 am
par frederic han
c'est bizarre avec la 1.2.3.57 j'ai un resultat correct mais bien plus compliqué qu'avant:
Code : Tout sélectionner
3>> factor(x^12-1,[sqrt(3),i])
1/4194304*(x-1)*(x+1)*(2*x+sqrt(3)*(sqrt(3)*1/6*2*i+1))*(2*x-sqrt(3)*(sqrt(3)*1/6*2*i+1))*(4*x+sqrt(3)*2*i-2)*(4*x+sqrt(3)*2*i+2)*(4*x-sqrt(3)*2*i+2)*(4*x-sqrt(3)*2*i-2)*(8*x-sqrt(3)*4*(-sqrt(3)*1/96*32*i+1))*(8*x+8*i)*(8*x+sqrt(3)*4*(-sqrt(3)*1/96*32*i+1))*(8*x-8*i)
est ce que ca dépend de la version de pari ou NTL car sous mingw64 avec pari 2.9 j'ai:
Code : Tout sélectionner
0>> factor(x^12-1,[sqrt(3),i])
(x-1)*(x+1)*(x-i)*(x+((-i)*sqrt(3)-1)/2)*(x+((-i)*3+3*sqrt(3))/6)*(x+((-i)*sqrt(3)+1)/2)*(x+((-i)*3-3*sqrt(3))/6)*(x+(i*3-3*sqrt(3))/6)*(x+(i*sqrt(3)-1)/2)*(x+i)*(x+(i*sqrt(3)+1)/2)*(x+(i*3+3*sqrt(3))/6)
Re: x^12-1
Publié : sam. août 05, 2017 11:12 am
par parisse
Peut-etre car la version windows 64 bits de Xcas renvoie la factorisation la plus simple, et elle n'est pas linkee avec NTL. Mais je n'ai aucune idee du pourquoi!
Re: x^12-1
Publié : ven. août 25, 2017 12:43 pm
par compsystems
In the help or documentation I do not find that factor cmd can accept a list of parameters in the second argument, How many hidden things there must be, and we do not know
factor(x^12-1,[sqrt(3),i])
Re: x^12-1
Publié : ven. août 25, 2017 1:43 pm
par parisse
The English documentation is unfortunately not as complete as the French documentation...
Re: x^12-1
Publié : mer. sept. 13, 2017 10:34 am
par frederic han
parisse a écrit :Peut-etre car la version windows 64 bits de Xcas renvoie la factorisation la plus simple, et elle n'est pas linkee avec NTL. Mais je n'ai aucune idee du pourquoi!
Ca n'a pas l'air de venir de NTL mais plutot d'une optimisation qui ne serait faite que sour unix car avec la version windows que je compile moi même avec NTL et mingw64 j'ai le resultat le plus simple.
(J'ai mis cette version de icas que j'utiliser pour tester la giac.dll de giacpy ici:
http://www.imj-prg.fr/~frederic.han/xca ... 2.3.57.zip
)
avec cette version ca commence comme cela:
Code : Tout sélectionner
// Maximum number of parallel threads 1
Help file /usr/bin/doc/fr/aide_cas not found
Added 0 synonyms
Help file /usr/bin/aide_cas not found
47 icas_eval factor(x^12-1,[sqrt(3),i])
gcd modular algo begin 48
4.8e-005 NTL factor begin
square-free decomposition...0
factoring multiplicity 1, deg = 2
SFFactor -- deflation: [2]
begin - step 0, 0; deg = 1
*** SFFactor, trivial case 1.
end - step 0, 0; time = 0
Re: x^12-1
Publié : mer. sept. 13, 2017 11:05 am
par parisse
J'ai teste avec une 1.2.1 et la 1.2.3 de l'agreg, mais ils renvoient la meme chose que la version actuelle, il n'est donc pas possible de faire du debug en parallele pour voir ou ca diverge.
Re: x^12-1
Publié : jeu. sept. 14, 2017 9:25 am
par parisse
j'ai trouve, c'est dans la remontee de l'appel a factorisation de pari que ca differe.
Code : Tout sélectionner
diff gausspol.cc gausspol.cc~
5576d5575
< // p_y=p_y/Tcontent(p_y);
5727,5728d5725
< // ? unitarize res
< *res._POLYptr=*res._POLYptr/res._POLYptr->coord.front().value;
Du coup j'ai du changer quelques input/output dans 3 tests (chk_cas, chk_fhan0 et chk_fhan21).