X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=check%2Fexam_clifford.cpp;h=b6ec71e455317bccb9b56eed2e96372f35ac47a3;hp=bea0a219f1d0cb045adf9737950bed8fa1a3905f;hb=6872b42cd09e72e7ca9906f4fc2fe775ace2f043;hpb=c029603c04ff63a6b171ff1f949d43c50413f49d diff --git a/check/exam_clifford.cpp b/check/exam_clifford.cpp index bea0a219..b6ec71e4 100644 --- a/check/exam_clifford.cpp +++ b/check/exam_clifford.cpp @@ -3,7 +3,7 @@ * Here we test GiNaC's Clifford algebra objects. */ /* - * GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2004 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 @@ -44,7 +44,7 @@ static unsigned check_equal_simplify(const ex &e1, const ex &e2) return 0; } -static unsigned clifford_check1(void) +static unsigned clifford_check1() { // checks general identities and contractions @@ -80,7 +80,7 @@ static unsigned clifford_check1(void) return result; } -static unsigned clifford_check2(void) +static unsigned clifford_check2() { // checks identities relating to gamma5 @@ -99,7 +99,7 @@ static unsigned clifford_check2(void) return result; } -static unsigned clifford_check3(void) +static unsigned clifford_check3() { // checks traces @@ -107,7 +107,7 @@ static unsigned clifford_check3(void) symbol dim("D"), m("m"), q("q"), l("l"), ldotq("ldotq"); varidx mu(symbol("mu"), dim), nu(symbol("nu"), dim), rho(symbol("rho"), dim), - sig(symbol("sig"), dim), kap(symbol("kap"), dim), lam(symbol("lam"), 4); + sig(symbol("sig"), dim), kap(symbol("kap"), dim), lam(symbol("lam"), dim); ex e; e = dirac_gamma(mu); @@ -156,27 +156,38 @@ static unsigned clifford_check3(void) + 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; } -static unsigned clifford_check4(void) +static unsigned clifford_check4() { // simplify_indexed()/dirac_trace() cross-checks @@ -210,7 +221,7 @@ static unsigned clifford_check4(void) return result; } -static unsigned clifford_check5(void) +static unsigned clifford_check5() { // canonicalize_clifford() checks @@ -221,7 +232,7 @@ static unsigned clifford_check5(void) ex e; e = dirac_gamma(mu) * dirac_gamma(nu) + dirac_gamma(nu) * dirac_gamma(mu); - result += check_equal(canonicalize_clifford(e), 2*lorentz_g(mu, nu)); + result += check_equal(canonicalize_clifford(e), 2*dirac_ONE()*lorentz_g(mu, nu)); e = (dirac_gamma(mu) * dirac_gamma(nu) * dirac_gamma(lam) + dirac_gamma(nu) * dirac_gamma(lam) * dirac_gamma(mu) @@ -238,7 +249,7 @@ static unsigned clifford_check5(void) return result; } -unsigned exam_clifford(void) +unsigned exam_clifford() { unsigned result = 0; @@ -250,7 +261,7 @@ unsigned exam_clifford(void) result += clifford_check3(); cout << '.' << flush; result += clifford_check4(); cout << '.' << flush; result += clifford_check5(); cout << '.' << flush; - + if (!result) { cout << " passed " << endl; clog << "(no output)" << endl;