x^12-1

Utilisation de Xcas

Modérateur : xcasadmin

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

x^12-1

Message par frederic han » sam. août 05, 2017 8:19 am

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)

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

Re: x^12-1

Message par parisse » sam. août 05, 2017 11:12 am

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!

compsystems
Messages : 518
Inscription : sam. févr. 04, 2017 11:34 pm
Localisation : Colombia
Contact :

Re: x^12-1

Message par compsystems » ven. août 25, 2017 12:43 pm

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])

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

Re: x^12-1

Message par parisse » ven. août 25, 2017 1:43 pm

The English documentation is unfortunately not as complete as the French documentation...

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

Re: x^12-1

Message par frederic han » mer. sept. 13, 2017 10:34 am

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

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

Re: x^12-1

Message par parisse » mer. sept. 13, 2017 11:05 am

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.

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

Re: x^12-1

Message par parisse » jeu. sept. 14, 2017 9:25 am

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

Répondre