modulo

[frederic han]

Salut, j'ai des Pb avec modulo selon les variables.
(en mode maple ca peut etre utile d'utiliser le Mod de pari et son lift pour palier au manque du %0)


Mini Pb:
pari()
a:=Mod(1,7);
P:=normal((x^2+5)*a);
lift(a)
lift(P)
marchent bien, mais si la variable est autre que x, il bascule en mode polynome lors du lift.
et s'il y a x^2+5*y il ne bascule en mode poly que pour y et pas pour x c'est un peu perturbant, sans etre faux.

En revanche:
en mode xcas ou maple, factor a plusieurs variables ne marche pas:

P:=normal((x^2+y^2)*Mod(1,2))

ou bien (x^2+y^2)%2 en xcas

a+

Fred

factor(P) donne 1.

[parisse]

Apparamment, c'est la fonction C de pari varn qui a un comportement different selon que x est variable ou pas, ou alors je n'ai pas bien compris. Je rajoute un test qui devrait renvoyer un polynome symbolique. Pour la factorisation modulaire multi-variables, je ne l'ai pas implemetee pour le moment (car mon code de factorisation sur Z s'en passe), je modifie le code pour renvoyer une erreur. En attendant, tu sais si on peut le faire avec PARI? Apparamment, ca ne marche pas avec factor ni avec factormod.

[frederic han]

Ah OK, je ne vois pas trop comment le faire avec pari, mais je connais tres peu.
il me dit:
*** factor: sorry, factor for general polynomials is not yet implemented.


Sous Macaulay c'est assez transparent (pour un corps):
ZZ/2[x,y]
factor(x^4+y^2)
fonctionne

Fred

[parisse]

Ca doit aussi surement exister avec CoCoA, il faut juste que je trouve la bonne routine et la conversion. Interfacer avec macaulay demanderait surement plus de boulot... ca peut valoir le coup si on récupère d'autres trucs intéressants, mais j'avais plutot en vue d'écrire une interface avec singular, tu sais si ça le fait?

[frederic han]

Oui singular le fait, je viens d'essayer:


> ring r=5,(x,y,z),dp; factorize((x^5+y^25+z^5),2);
// ** redefining r **
[1]:
_[1]=y5+x+z
[2]:
5

a+

Frederic