greduce

Bugs

Modérateur : xcasadmin

frederic han
Messages : 1113
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

greduce

Message par frederic han » mer. sept. 23, 2015 8:26 am

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
Messages : 5109
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: greduce

Message par parisse » mer. sept. 23, 2015 6:52 pm

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

frederic han
Messages : 1113
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: greduce

Message par frederic han » mer. sept. 23, 2015 8:29 pm

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
Messages : 1113
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: greduce

Message par frederic han » jeu. sept. 24, 2015 7:06 am

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
Messages : 5109
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: greduce

Message par parisse » jeu. sept. 24, 2015 8:06 am

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.

Répondre