Re: giacpy et sage
Publié : sam. oct. 05, 2013 9:46 am
Je viens d'essayer cyclic9 mod en 32 bits, ca passe (1700 secondes au lieu de 650 sur la meme machine mais en 64 bits).
Tu peux avoir des infos avec export GIAC_DEBUG=2
Voila mon cyclic9mod (ou je lui ai fait calculer le nombre de monomes de la reponse, c'est un peu plus de 1 millions, et la taille des coefficients est de 1600 chiffres environ en base 10 pour cyclic9 pas mod).
Tu peux avoir des infos avec export GIAC_DEBUG=2
Voila mon cyclic9mod (ou je lui ai fait calculer le nombre de monomes de la reponse, c'est un peu plus de 1 millions, et la taille des coefficients est de 1600 chiffres environ en base 10 pour cyclic9 pas mod).
Code : Tout sélectionner
cyclic9:=[x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9, x1*x2 + x2*x3 + x3*x4 + x4*x5 + x5*x6 + x6*x7 + x7*x8 + x1*x9 + x8*x9, x1*x2*x3 + x2*x3*x4 + x3*x4*x5 + x4*x5*x6 + x5*x6*x7 + x6*x7*x8 + x1*x2*x9 + x1*x8*x9 + x7*x8*x9, x1*x2*x3*x4 + x2*x3*x4*x5 + x3*x4*x5*x6 + x4*x5*x6*x7 + x5*x6*x7*x8 + x1*x2*x3*x9 + x1*x2*x8*x9 + x1*x7*x8*x9 + x6*x7*x8*x9, x1*x2*x3*x4*x5 + x2*x3*x4*x5*x6 + x3*x4*x5*x6*x7 + x4*x5*x6*x7*x8 + x1*x2*x3*x4*x9 + x1*x2*x3*x8*x9 + x1*x2*x7*x8*x9 + x1*x6*x7*x8*x9 + x5*x6*x7*x8*x9, x1*x2*x3*x4*x5*x6 + x2*x3*x4*x5*x6*x7 + x3*x4*x5*x6*x7*x8 + x1*x2*x3*x4*x5*x9 + x1*x2*x3*x4*x8*x9 + x1*x2*x3*x7*x8*x9 + x1*x2*x6*x7*x8*x9 + x1*x5*x6*x7*x8*x9 + x4*x5*x6*x7*x8*x9, x1*x2*x3*x4*x5*x6*x7 + x2*x3*x4*x5*x6*x7*x8 + x1*x2*x3*x4*x5*x6*x9 + x1*x2*x3*x4*x5*x8*x9 + x1*x2*x3*x4*x7*x8*x9 + x1*x2*x3*x6*x7*x8*x9 + x1*x2*x5*x6*x7*x8*x9 + x1*x4*x5*x6*x7*x8*x9 + x3*x4*x5*x6*x7*x8*x9, x1*x2*x3*x4*x5*x6*x7*x8 + x1*x2*x3*x4*x5*x6*x7*x9 + x1*x2*x3*x4*x5*x6*x8*x9 + x1*x2*x3*x4*x5*x7*x8*x9 + x1*x2*x3*x4*x6*x7*x8*x9 + x1*x2*x3*x5*x6*x7*x8*x9 + x1*x2*x4*x5*x6*x7*x8*x9 + x1*x3*x4*x5*x6*x7*x8*x9 + x2*x3*x4*x5*x6*x7*x8*x9, x1*x2*x3*x4*x5*x6*x7*x8*x9 - 1]:;
p:=prevprime(45000);
//p:=prevprime(2^25); // 33554393
time(H:=gbasis(cyclic9 % p,indets(cyclic9),revlex,with_cocoa=false));
size(eval(H,1));
b:=0:; for j from 0 to size(eval(H,1))-1 do b+=size(eval(H,1)[j]); od:; b;
write("Hcyclic9",H);