factor with ** and extensions

[frederic han]

il y a un problème avec factor pour les polynomes en notation 'printpow(-1)' j'ai l'impression qu'avec des vieilles versions de pari cela fonctionne mais avec pari 2.13 sous fedora j'ai ceci.

Code : Tout sélectionner

13>> P:=x**4-x**2+1
x^4-x^2+1
// Time 0
14>> printpow(1)
"^"
// Time 0
15>> factor(P,[i,sqrt(3)])
(x+((-i)*3+3*sqrt(3))/6)*(x+(i*3-3*sqrt(3))/6)*(x+((-i)*3-3*sqrt(3))/6)*(x+(i*3+3*sqrt(3))/6)
// Time 0.01
16>> printpow(-1)
"**"
// Time 0
17>> factor(P,[i,sqrt(3)])

  ***   syntax error, unexpected '*': [x,y,64*x**4-64*x**2+64]
  ***                                                        ^-"Erreur dans PARI Erreur: Valeur Argument Incorrecte"
// Time 0.01

[parisse]

Ca devrait pouvoir se corriger en remplacant dans pari.cc

Code : Tout sélectionner

  GEN gen2GEN(const gen & e,const vecteur & vars,GIAC_CONTEXT){
    const char * ptr=__pow.s;
    __pow.s="^";
    cb_pari_err_recover=gp_err_recover;
    if (setjmp(env)){
      __pow.s=ptr;
      setsizeerr(gettext("Error in PARI subsystem"));
    }
    GEN G=ingen2GEN(e,vars,contextptr);
    __pow.s=ptr;
    return G;
  }