Le forum de XCAS

Xcas: un logiciel libre de calcul formel
Nous sommes actuellement le Mar Déc 12, 2017 6:23 am

Heures au format UTC




Publier un nouveau sujet Répondre au sujet  [ 7 messages ] 
Auteur Message
 Sujet du message: x^12-1
MessagePublié: Sam Aoû 05, 2017 8:19 am 
Hors-ligne

Inscrit le: Dim Mai 20, 2007 7:09 am
Messages: 1022
Localisation: Paris
c'est bizarre avec la 1.2.3.57 j'ai un resultat correct mais bien plus compliqué qu'avant:
Code:
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:
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)


Haut
 Profil  
 
 Sujet du message: Re: x^12-1
MessagePublié: Sam Aoû 05, 2017 11:12 am 
Hors-ligne

Inscrit le: Mar Déc 20, 2005 4:02 pm
Messages: 4153
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!


Haut
 Profil  
 
 Sujet du message: Re: x^12-1
MessagePublié: Ven Aoû 25, 2017 12:43 pm 
Hors-ligne

Inscrit le: Sam Fév 04, 2017 11:34 pm
Messages: 126
Localisation: Colombia
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])


Haut
 Profil  
 
 Sujet du message: Re: x^12-1
MessagePublié: Ven Aoû 25, 2017 1:43 pm 
Hors-ligne

Inscrit le: Mar Déc 20, 2005 4:02 pm
Messages: 4153
The English documentation is unfortunately not as complete as the French documentation...


Haut
 Profil  
 
 Sujet du message: Re: x^12-1
MessagePublié: Mer Sep 13, 2017 10:34 am 
Hors-ligne

Inscrit le: Dim Mai 20, 2007 7:09 am
Messages: 1022
Localisation: Paris
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:
// 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


Haut
 Profil  
 
 Sujet du message: Re: x^12-1
MessagePublié: Mer Sep 13, 2017 11:05 am 
Hors-ligne

Inscrit le: Mar Déc 20, 2005 4:02 pm
Messages: 4153
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.


Haut
 Profil  
 
 Sujet du message: Re: x^12-1
MessagePublié: Jeu Sep 14, 2017 9:25 am 
Hors-ligne

Inscrit le: Mar Déc 20, 2005 4:02 pm
Messages: 4153
j'ai trouve, c'est dans la remontee de l'appel a factorisation de pari que ca differe.
Code:
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).


Haut
 Profil  
 
Afficher les messages publiés depuis:  Trier par  
Publier un nouveau sujet Répondre au sujet  [ 7 messages ] 

Heures au format UTC


Qui est en ligne ?

Utilisateurs parcourant actuellement ce forum : Aucun utilisateur inscrit et 1 invité


Vous ne pouvez pas publier de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas insérer de pièces jointes dans ce forum

Rechercher pour:
Sauter vers:  
cron
Powered by phpBB® Forum Software © phpBB Group
Traduction réalisée par Maël Soucaze © 2009 phpBB.fr