From 2e13bfeb92ce55b3ee7e24cc6041f70c7a0eff4e Mon Sep 17 00:00:00 2001 From: Christian Bauer Date: Fri, 20 Apr 2001 23:03:43 +0000 Subject: [PATCH] - dirac_one() -> dirac_ONE() - more clifford checks --- check/exam_archive.cpp | 2 +- check/exam_clifford.cpp | 34 ++++++++++++++++++++++++++++++++-- 2 files changed, 33 insertions(+), 3 deletions(-) diff --git a/check/exam_archive.cpp b/check/exam_archive.cpp index af8f4905..2b8862cc 100644 --- a/check/exam_archive.cpp +++ b/check/exam_archive.cpp @@ -36,7 +36,7 @@ unsigned exam_archive(void) // This expression is complete nonsense but it contains every type of // GiNaC object - e = -42 * x * pow(y, sin(y*Catalan)) * dirac_one() + e = -42 * x * pow(y, sin(y*Catalan)) * dirac_ONE() * epsilon_tensor(idx(fail(), 3), idx(0, 3), idx(y/2, 3)) + lorentz_g( varidx(lst(x, -11*y, acos(2*x).series(x==3-5*I, 3)) * color_ONE() diff --git a/check/exam_clifford.cpp b/check/exam_clifford.cpp index d4491191..2b28be4c 100644 --- a/check/exam_clifford.cpp +++ b/check/exam_clifford.cpp @@ -54,13 +54,42 @@ static unsigned clifford_check1(void) varidx mu(symbol("mu"), dim), nu(symbol("nu"), dim); ex e; + e = dirac_ONE() * dirac_ONE(); + result += check_equal(e, dirac_ONE()); + + e = dirac_ONE() * dirac_gamma(mu) * dirac_ONE(); + result += check_equal(e, dirac_gamma(mu)); + + e = dirac_gamma(varidx(2, dim)) * dirac_gamma(varidx(1, dim)) * + dirac_gamma(varidx(1, dim)) * dirac_gamma(varidx(2, dim)); + result += check_equal(e, dirac_ONE()); + e = dirac_gamma(mu) * dirac_gamma(nu) * dirac_gamma(nu.toggle_variance()) * dirac_gamma(mu.toggle_variance()); - result += check_equal_simplify(e, pow(dim, 2) * dirac_one()); + result += check_equal_simplify(e, pow(dim, 2) * dirac_ONE()); e = dirac_gamma(mu) * dirac_gamma(nu) * dirac_gamma(mu.toggle_variance()) * dirac_gamma(nu.toggle_variance()); - result += check_equal_simplify(e, 2*dim*dirac_one()-pow(dim, 2)*dirac_one()); + result += check_equal_simplify(e, 2*dim*dirac_ONE()-pow(dim, 2)*dirac_ONE()); + + return result; +} + +static unsigned clifford_check2(void) +{ + // checks identities relating to gamma5 + + unsigned result = 0; + + symbol dim("D"); + varidx mu(symbol("mu"), dim), nu(symbol("nu"), dim); + ex e; + + e = dirac_gamma(mu) * dirac_gamma5() + dirac_gamma5() * dirac_gamma(mu); + result += check_equal(e, 0); + + e = dirac_gamma5() * dirac_gamma(mu) * dirac_gamma5() + dirac_gamma(mu); + result += check_equal(e, 0); return result; } @@ -73,6 +102,7 @@ unsigned exam_clifford(void) clog << "----------clifford objects:" << endl; result += clifford_check1(); cout << '.' << flush; + result += clifford_check2(); cout << '.' << flush; if (!result) { cout << " passed " << endl; -- 2.44.0