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
Hi! On Sun, Jul 29, 2001 at 09:19:02PM +0200, Do Hoang Son wrote:
dirac_gamma5()*dirac_ONE() = -dirac_ONE()*dirac_gamma5()
No, but it rewrote (a+b)*gamma5 as -gamma5*(a+b). Anyway, this is fixed now. Bye, Christian -- / Coding on PowerPC and proud of it \/ http://www.uni-mainz.de/~bauec002/
On Mon, 30 Jul 2001, Christian Bauer wrote:
Hi!
On Sun, Jul 29, 2001 at 09:19:02PM +0200, Do Hoang Son wrote:
dirac_gamma5()*dirac_ONE() = -dirac_ONE()*dirac_gamma5()
No, but it rewrote (a+b)*gamma5 as -gamma5*(a+b). Anyway, this is fixed now.
It seems OK now, at least the t -> bH decay now can be completely writen in C++ with GiNaC and Xloops-GiNaC libs. Thank Christian and cheer up! Son
participants (2)
-
Christian Bauer
-
Do Hoang Son