* Here we test GiNaC's Clifford algebra objects. */
/*
- * GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2002 Johannes Gutenberg University Mainz, Germany
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
+ dim * dirac_gamma5() * dirac_gamma(nu) * dirac_gamma(rho) * dirac_gamma(sig) * dirac_gamma(kap);
e = dirac_trace(e).simplify_indexed();
e = (e / (dim - 4)).normal();
- result += check_equal(e, 8 * I * eps0123(nu, rho, sig, kap));
+ result += check_equal(e, 8 * I * lorentz_eps(nu.replace_dim(4), rho.replace_dim(4), sig.replace_dim(4), kap.replace_dim(4)));
// one-loop vacuum polarization in QED
e = dirac_gamma(mu) *
- (dirac_slash(l, dim) + dirac_slash(q, dim) + m * dirac_ONE()) *
+ (dirac_slash(l, dim) + dirac_slash(q, 4) + m * dirac_ONE()) *
dirac_gamma(mu.toggle_variance()) *
(dirac_slash(l, dim) + m * dirac_ONE());
e = dirac_trace(e).simplify_indexed(sp);
result += check_equal(e, 4*((2-dim)*l*l + (2-dim)*ldotq + dim*m*m).expand());
- e = dirac_slash(q, dim) *
- (dirac_slash(l, dim) + dirac_slash(q, dim) + m * dirac_ONE()) *
- dirac_slash(q, dim) *
+ e = dirac_slash(q, 4) *
+ (dirac_slash(l, dim) + dirac_slash(q, 4) + m * dirac_ONE()) *
+ dirac_slash(q, 4) *
(dirac_slash(l, dim) + m * dirac_ONE());
e = dirac_trace(e).simplify_indexed(sp);
result += check_equal(e, 4*(2*ldotq*ldotq + q*q*ldotq - q*q*l*l + q*q*m*m).expand());
+ // stuff that had problems in the past
+ ex prop = dirac_slash(q, dim) - m * dirac_ONE();
+ e = dirac_slash(l, dim) * dirac_gamma5() * dirac_slash(l, dim) * prop;
+ e = dirac_trace(dirac_slash(q, dim) * e) - dirac_trace(m * e)
+ - dirac_trace(prop * e);
+ result += check_equal(e, 0);
+
+ e = (dirac_gamma5() + dirac_ONE()) * dirac_gamma5();
+ e = dirac_trace(e);
+ result += check_equal(e, 4);
+
return result;
}
result += clifford_check3(); cout << '.' << flush;
result += clifford_check4(); cout << '.' << flush;
result += clifford_check5(); cout << '.' << flush;
-
+
if (!result) {
cout << " passed " << endl;
clog << "(no output)" << endl;