greduce

[frederic han]

Salut,
pourquoi me faut il 2 greduce pour avoir 0? est ce un probleme du a la facon de nommer mes variables?

Code : Tout sélectionner

n:=2;
A:=matrix(n,n,(u,v)->a[u,v]);
B:=matrix(n,n,(u,v)->b[u,v]);
tmp:=A*B-A-B:;
I:=[seq(seq(tmp[u,v],u=0..n-1),v=0..n-1)]:;
X:=[seq(seq(a[u,v],u=0..n-1),v=0..n-1),seq(seq(b[u,v],u=0..n-1),v=0..n-1)]:;
gb:=gbasis(I,X):;
greduce(greduce((B*A-A*B)[1,1],gb),gb); // donne 0
greduce((B*A-A*B)[1,1],gb); // ne donne pas 0

[parisse]

oui. le code actuel ne gere pas les variables indiciees. Je regarderai vendredi si c'est facile a ameliorer...

[frederic han]

En 1.2.1-15 J'ai aussi le probleme avec des variables a,b,c,d,x,y,z,t et avec:

Code : Tout sélectionner

n:=2;
fa(u,v):=#("a"+string(u)+string(v));
fb(u,v):=#("b"+string(u)+string(v));
A:=matrix(n,n,fa);
B:=matrix(n,n,fb);
X:=simplify([seq(seq(A[u,v],u=0..n-1),v=0..n-1),seq(seq(B[u,v],u=0..n-1),v=0..n-1)]);
tmp:=A*B-A-B;
I:=[seq(seq(tmp[u,v],u=0..n-1),v=0..n-1)];
gb2:=gbasis(I,X);
greduce((B*A-A*B)[1,1],gb2);

mais pas en 1.2.0-13!


aussi j'ai remarque un petit plantage avec:

Code : Tout sélectionner

greduce(0,[x])

[frederic han]

Est ce que ca ne serait pas juste un effet de bord du changement d'ordre dans le retour de gbasis pour etre comme singular.

Code : Tout sélectionner

greduce((B*A-A*B)[1,1],revlist(gb2));

[parisse]

ah oui, bien vu, ca a l'air de passer en renversant l'ordre.
Mais il faudra quand meme que je regarde les variables indiciees ne sont pas traitees correctement.