[GiNaC-devel] Vladimirs patch has been applied.

Vladimir Kisil kisilv at maths.leeds.ac.uk
Fri Aug 18 13:40:37 CEST 2006


		Dear Chris,

>>>>> "CD" == Chris Dams <Chris.Dams at mi.infn.it> writes:
    CD> I have another question for Vladimir: at some places in
    CD> clifford.cpp locally within some functions a varidx is
    CD> created. But now that clifford objects can also carry indices
    CD> without variance, is this still correct?  Maybe it should depend
    CD> on the kind of index that wat supplied to such a function?

	Thanks for pointing this, I have overseen it. Now I try to look on
  it, but....

    CD> I have improved the dummy index renaming. Now dummy indices are (if
    CD> necessary) also renamed after mul::map is called. I made the code a bit
    CD> nicer.

	... it seems to be a problem with dummy renaming in combination of
  indexes without variance and non-commutative products. Here is the
  example: 


#include <stdexcept>
#include <sstream>
using namespace std;
 
#include <ginac/ginac.h>
using namespace GiNaC;


int main() {
	try {
		idx i(symbol("i"),4);
		realsymbol t("t"), x("x"), y("y"), z("z");
		matrix M(4,1), M1(4,1);
		M= t, x,y,z;
		M1= 1,2,3,4;
		ex e=clifford_unit(i,diag_matrix(lst(1,-1,1,-1)));
		ex a=indexed(M,i)*e;
		ex a1=indexed(M1,i)*e;
		cout << (a+dirac_ONE())*a1 << endl;
		//-> [[1],[2],[3],[4]].i*((ONE+[[t],[x],[y],[z]].symbol8*e.symbol8)*e.symbol8)

	} catch  (exception &p) {
			cerr << "Got problem: " <<  p.what() << endl;
	}
}

		Have you got any idea how to fix it?

		Best wishes,
		Vladimir

-- 
Vladimir V. Kisil     email: kisilv at maths.leeds.ac.uk
--                      www: http://maths.leeds.ac.uk/~kisilv/


More information about the GiNaC-devel mailing list