Le forum de XCAS

Aie, J'ai aussi remarque cela:
gp:=2%7;
gpz:=1965198393271355449089580772376483480%1073741824000000000000000000061203284109000000000000000000000000008409;

pari_znorder(gp) donne bien 3 mais:

pari_znorder(gpz) donne

[New Thread 0xb56efb70 (LWP 31092)]
*** incorrect type in generic discrete logarithm (order factorization).

Y-a-t-il eu des changements dans flisexpr (les Mod ne sont pas traduits directement mais via l'interpreteur)?

bon apparamment c'est utoi qui ne fonctionne plus comme avant ou un changement pour signaler l'utilisation d'un argument par defaut.
Tu peux essayer ca ligne 736 environ et voir si ca pose probleme? j'ai aussi l'impression que l'erreur trapping deconne.

bon pour l'error trapping il fallait aussi definir cb_pari_err_recover

Oui ca corrige le probleme

a+
Fred

Salut,

Avec la version 0.9.3 du 1/9 je n'ai plus de Pb.

Le bug du n:=2^128+1;pari("factor",n) signale dans un autre post est meme repare alors qu'avec giac version pari2.3 pari("factor",n) etait faux mais pari(factor,n) correct.

A part la precision des flottants, je ne remarque rien.

Bravo, voici en attache un petit fichier de tests

a+

Fred
[code]
maple_mode(0):;
p:=nextprime(nextprime(10^30)*nextprime(20^30)) ;
g:=123456789101112 % p ;
pari() :;
("version de pari:",pari(version));//affiche la version de pari. Ca rale avec giac 093 et pari23
pari_znorder(g) ;
ifactor(p-1) ;
pari_factorint(p-1) ;
pari_isprime(p,1);//c'est normal qu'il y ait un terme de moins ([8097977,2,1])avec pari25 qu'avec pari23, le test a ete ameliore
znstar(8) ;
znstar(p) ;
tt1:=znstar(p-1) ;
dim(tt1[1]) ;
(dim(tt1[1])) ==6;// doit faire 1
v1:=[17,-60,-120,0,0] ;
EC:=ellinit(v1) ;
elltors(EC),elltors(ellinit([17,-60,-120,0,0])) ;
ellpow(EC,[30,-90],8);//doit etre nul
v2:=[17/3,-60*sqrt(2),-120,0,0] ;
ecnonrat:=ellinit(v2) ;
pari(sqrt(2)) ;
n:=2^128+1 ;
ifactors(n) ;
pari("factor",n);// test bug giac 093 et pari23. repare avec pari25
pari(factor,n);
normal(poly2symb(pari(x^2+y-99),x,y)-x^2-3*y);
pari(sqrt(2)) ;
normal(poly2symb(pari(x^2+y-99),x,y)-x^2-3*y) ;
nodisp(Digits := 50) ;
pari(1.12345678910111213141516171819202122232425262728293) ;//ici la precision fonctionne.
s2:=approx(sqrt(2)) ;
pari_sqrt(s2) ;//teste perte de precision (reste en 14 chiffres)
pari("",s2);pari(s2);//c'est en fait logique,le premier evalue s2 dans giac le second dans pari donc c'est le polynome s2
pari(sqrt(2)) ;
pari(c=sqrt(3));
pari(c*c);
M:=[[-73,68,-95,56,-41,0,-50,62,1,-23,-3,55,87,-12,95,-39,87,59,99,-63,75,24,65,-83,-25,-48,-38,39,-31,-87,0,79,-26,-7,87,-37,92,44,-99,85],[-73,89,60,-76,60,60,-47,27,19,-17,44,-74,-13,-87,32,10,-3,-96,-96,-80,97,-71,19,-81,-26,40,-87,-35,-51,73,-89,-24,46,-1,68,-11,32,10,-67,-91],[-92,66,5,-42,-8,-33,52,-89,55,74,-96,31,34,-8,92,0,8,59,4,-70,-54,23,43,-13,-80,-74,-76,93,96,-50,-27,34,-15,-91,-40,47,60,70,-60,20],[57,37,12,-35,25,-40,-90,-99,61,-2,89,24,3,-71,62,71,2,-8,23,-21,48,87,-88,-68,-17,4,-17,-86,62,9,6,44,9,63,-19,12,12,-6,55,-98],[12,-97,-6,-84,-57,-76,31,85,51,21,71,85,-17,-66,87,30,-67,-45,48,-28,-32,69,67,-13,19,-72,58,-35,42,-45,20,-19,44,-78,42,-31,80,44,-87,-18],[-76,-63,67,-43,94,98,-61,59,-32,54,35,-63,44,35,35,-49,93,99,-3,-86,58,1,-6,17,39,-97,61,50,39,26,-29,14,-23,-39,-87,31,76,92,4,48],[-81,-96,-81,26,63,-7,-76,-59,-60,-45,47,50,-94,-4,60,-10,-19,-7,-23,-72,-63,-95,-32,-62,67,14,5,53,-75,-73,-78,-75,-80,-1,-27,-42,93,-74,44,59],[-36,-35,-39,67,-71,95,12,-22,-51,40,13,66,-8,13,-33,4,-47,-77,50,62,-66,-91,28,-82,0,-63,-25,-54,-95,-79,12,-87,73,-98,-13,-45,-28,-36,-64,-41],[-74,-6,98,33,-74,-10,-68,-31,88,6,45,0,96,51,-35,-42,93,15,-1,-91,-34,-39,-59,2,76,38,76,3,-21,-42,18,13,-6,98,57,94,-43,-33,-74,85],[34,89,93,-90,-90,62,67,71,32,33,22,60,47,-66,41,-65,78,-51,81,6,-4,-60,-49,-80,-64,-95,14,16,5,-19,-87,-91,45,94,-6,-70,-34,-45,-24,-10],[-52,88,-78,-81,17,22,29,11,87,-96,-48,54,-34,68,-81,-43,32,58,-2,66,90,-33,-45,91,79,43,46,-96,-39,36,26,-60,71,-22,69,-99,47,-86,-82,38],[10,76,-61,-60,-65,5,12,-97,66,-92,-68,9,-68,-41,15,-12,-7,58,54,-28,88,-81,59,51,-63,-1,-54,-78,5,-8,17,-82,54,-59,-8,-48,45,-89,19,37],[-79,99,42,-87,-51,93,-40,51,-56,-36,-83,-24,44,-89,12,-46,-25,-50,-77,-1,61,-58,-54,54,-74,67,-36,83,-56,42,27,-60,-24,-17,55,-58,61,5,-51,36],[-12,80,35,57,21,23,-76,-24,-82,-72,77,47,85,-50,11,-48,31,-69,3,25,-25,-37,26,87,-85,4,-3,-35,88,29,46,-2,82,32,-89,57,-83,-33,-68,-54],[-2,-70,-13,84,45,99,31,13,79,-98,-36,-30,-19,-68,-97,61,-81,-87,-8,-54,-49,2,16,87,85,-57,-57,-81,-78,-49,73,-7,-53,54,-75,27,13,-22,-82,-55],[77,54,-6,7,-91,-91,45,-60,70,-46,66,-47,65,56,-78,-39,-46,32,33,80,-67,85,89,-54,-3,99,68,91,-43,67,83,-51,-70,-38,-55,-70,1,-46,-56,99],[-48,77,-52,-47,-45,29,-74,-29,-57,24,-17,86,-56,-18,15,-66,79,-36,92,-7,-54,-43,41,-71,45,95,-32,28,4,28,-16,22,44,-70,59,10,36,55,91,-28],[11,-57,91,-23,49,-51,-7,-14,27,41,-26,-51,-49,-71,5,89,-68,28,-86,94,-38,-9,55,85,-69,58,-71,93,30,-30,-17,-82,16,-91,-38,92,43,-2,56,30],[-26,33,-12,74,39,-31,89,-65,-5,32,-68,20,17,2,-53,-30,93,-87,-19,67,-43,-93,62,62,-85,-55,91,75,22,37,-52,96,-52,-64,-90,-11,-17,-38,-94,68],[-58,81,74,-3,17,-76,-53,51,-87,37,-58,99,-65,-97,93,-39,-37,-25,84,66,-20,79,-97,-86,83,46,15,4,49,93,-32,-6,-20,48,92,-85,-94,-85,41,10],[72,83,36,54,-16,92,-59,64,-83,21,40,-91,9,32,59,-89,12,-36,81,-58,10,49,89,-84,15,-95,36,-34,-74,33,-78,26,-46,-3,-7,-98,23,-88,-35,-70],[35,5,-40,-57,-93,-91,58,15,66,-15,-17,-85,-66,14,-16,-63,56,-76,56,41,23,43,-26,-22,38,24,33,-95,72,93,31,8,76,-47,-57,8,-34,16,36,57],[77,-93,95,-60,-10,73,63,62,10,77,3,-66,42,2,86,42,-90,-85,87,31,51,75,-97,56,94,-17,-8,49,-18,-33,32,-34,-21,27,96,-65,47,12,41,-4],[83,40,28,-67,-71,71,49,-98,48,-84,-21,-52,-19,-8,41,-73,48,-80,6,-94,47,69,-16,-30,0,-54,-40,-85,76,-85,-67,-59,-18,-17,-16,31,-75,-89,-50,-68],[36,-24,32,20,-22,76,12,-27,-5,85,-12,-81,-40,-65,-70,-3,76,-42,-38,-3,-53,-33,49,1,-85,-21,16,-80,-24,79,-24,13,-42,-60,76,-84,62,66,78,48],[89,-96,36,-69,11,72,7,76,37,62,78,-3,26,53,4,-18,-13,0,-19,-9,-73,39,15,-75,-8,-41,35,0,-43,-3,2,-90,-9,96,-46,82,-1,94,56,-60],[15,-43,-14,31,-12,-39,-87,-14,6,-80,26,-65,-77,-50,57,-54,22,77,-65,-60,76,1,47,7,55,86,39,15,13,-91,14,11,43,37,58,-83,-2,-27,-14,15],[80,-89,3,-57,6,69,-9,52,30,45,-99,-99,-55,-83,-23,96,-59,-2,95,71,52,-74,-80,24,52,3,-26,-29,-72,74,-62,-21,15,45,63,-93,-75,-38,2,-72],[98,-90,-41,-50,-22,76,54,32,3,49,22,91,10,56,-89,5,75,52,-94,7,-54,-48,-18,75,-80,96,57,17,-80,84,-18,-52,9,17,34,53,97,-82,42,-32],[-99,-32,20,17,80,-88,25,84,44,82,-67,-55,-78,-11,30,11,0,-62,-65,-31,-62,-1,-83,-8,44,17,-45,78,-77,-26,6,-43,-10,-25,-9,90,-42,0,1,89],[-9,89,27,-1,25,72,98,-8,-38,-60,-52,-51,22,-36,94,39,-15,56,44,-8,-27,80,-93,-47,-73,-9,-65,51,36,-68,78,-50,46,-27,68,86,-73,-13,98,-82],[12,-23,46,50,-81,-59,-94,-30,-86,-17,-70,57,67,-26,8,9,68,40,-48,-58,11,17,-5,94,43,-71,29,22,-5,16,42,-5,60,20,82,65,71,-11,36,-39],[12,-56,-63,43,-74,87,74,82,73,46,42,-95,77,68,-12,51,-85,55,-12,-74,98,-51,46,-73,-87,-56,4,84,24,68,75,91,-36,-13,89,16,-32,-42,-31,23],[98,30,-90,-21,-71,-7,-4,50,91,-48,-82,10,8,65,-79,-83,74,-80,65,99,22,40,-88,-58,-50,59,-84,-57,-87,18,52,-11,0,-83,-47,-92,-40,2,57,21],[-93,-86,96,-5,-6,-45,-19,69,-45,20,15,2,10,64,-66,-10,51,10,40,88,-19,6,-8,30,-64,-74,72,-68,64,6,-89,-8,-9,-60,98,93,23,-99,38,-11],[58,18,-94,-16,32,0,28,14,50,-14,-10,19,90,-44,-52,-97,-39,-18,69,-8,39,14,-5,90,76,-73,0,23,-31,87,8,-22,-4,-12,86,-65,32,9,-15,94],[31,34,-11,3,62,2,-16,60,34,-93,45,-69,-22,-61,45,-53,-3,73,-21,1,-2,1,-97,41,70,99,7,46,56,-49,91,9,-25,-96,-36,5,15,-19,64,56],[12,23,0,89,73,-12,-95,89,82,77,-92,-55,-20,27,35,83,90,27,21,11,48,-15,-63,37,-38,41,-59,92,-30,26,52,-55,52,-45,35,-90,21,71,-92,-51],[-33,-64,-36,67,84,51,19,-92,-35,-56,-21,21,51,-28,7,-41,-24,-5,44,31,43,-26,-89,79,11,23,-48,-68,16,1,94,51,-19,-54,22,88,-17,57,-31,16],[97,-79,28,-33,12,71,94,50,97,11,-97,87,56,-40,-97,-84,46,51,49,-2,-46,86,2,26,-43,-94,59,22,41,-69,94,16,-2,-69,23,38,-25,-55,83,92]]:;
det(M)==(pari_matdet(M));
M:=[[1,2,3/7],[5,1/3,-6],[-3,0,1/9]];
pari_matdet(M)==det(M);

[/code]

Bonsoir,

Je suis tombé sur le même problème et j'avais cru comprendre à la lecture des réponses qu'il existait une version qui compilait avec pari 2.5. Je viens de télécharger la version instable du 09/09/2011 dont la compilation échoue sur

../../../../../cvs/rpl/rplcas/giac-0.9.3/src/pari.cc:267:43: error: ‘lgeflist’ was not declared in this scope
...
../../../../../cvs/rpl/rplcas/giac-0.9.3/src/pari.cc:468:9: error: ‘struct gp_data’ has no member named ‘env’
...
../../../../../cvs/rpl/rplcas/giac-0.9.3/src/pari.cc:582:69: error: ‘combine_factors’ was not declared in this scope
...
../../../../../cvs/rpl/rplcas/giac-0.9.3/src/pari.cc:681:25: error: invalid conversion from ‘const char*’ to ‘char*’ [-fpermissive]

La dernière erreur est triviale à corriger. En revanche, les trois premières me posent problème...

Cordialement,

JB

ca ne devrait pas au sens ou dans pari.cc les lignes suivantes
#if PARI_VERSION_CODE<PARI_VERSION(2,4,0) // 132096
#define PARI23
#else
...
definissent PARI23 si la version de pari est < a 2.4.0.
IL doit y avoir un pb avec la version de giac desarchivee.

Oups, au temps pour moi. Mon makefile s'évertuait à recompiler l'ancienne version. La compilation passe correctement. Désolé pour le bruit.

Et merci naturellement pour la réponse rapide que je n'avais pas encore vue lorsque j'écrivais la mienne...