29 Jul
2001
29 Jul
'01
9:19 p.m.
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