[GiNaC-list] BUG or incoherence: matrix A; A*pow(A,-1) =! 1

jros jros at unavarra.es
Wed Sep 16 16:46:34 CEST 2009


We now that for matrix class objects A, GiNaC simplifies 

matrix A;
cout << A-A << endl;

give 0

The same way that 

symbol a
cout << a-a << endl;

gives 0

But I think that for coherence

matrix A;
cout << A*pow(A,-1) << endl;

should give 1, (I mean reduced to numeric(1) regardless the particular
form of A, or at least if it is square)

the  same way as

symbol a
cout << a*pow(a,-1) << endl;

is 1

By the way it would be wonderful if A-A give matrix-es instead of
numeric-s, although I think you've already had long discussions about

Thank you,

PD: I wouldn't mind to derive a class from matrix to change this
behavior the easy way, but I don't know if this can be done this way.
Any light in this direction would be appreciated.

More information about the GiNaC-list mailing list