GiNaC/ginac clifford.cpp

[Do Hoang Son]

Hi,

There still has a bug somewhere in simplify_index(). Look at this case:

ex dt = dirac_slash(p1, D) + m1* dirac_ONE();
ex db = dirac_slash(p2, D) + m2* dirac_ONE();
cout << (dt*dirac_gamma5()*db*dirac_gamma5()).simplify_indexed()<< endl;

It returns:
-m2*p1\-p1\*p2\-m2*ONE*m1-p2\*m1 

The correct one should be:
m2*p1\-p1\*p2\+m2*ONE*m1+p2\*m1 

It seems that GiNaC does:
dirac_gamma5()*dirac_ONE() = -dirac_ONE()*dirac_gamma5()

Cheers,
Son