Name: Vladyslav Shtabovenko Date: 10/21/17-07:10:22 AM Z
> I’ve stumbled upon a frustrating bug in the FeynCalcExternal form
of
> D-dimensional Eps expressions. Essentially, the expressions
behave
> differently depending on how many of their indices are contracted
with
> four vectors (and which indices are contracted!).
> Attached is a short notebook file that illustrates the point
> concisely. Has anyone else run into this?
>
> Thanks,
>
> Alex Williamson
> University of Toronto
With the current stable and development version this works as expected:
ex1 = Table[
IntegerDigits[n, 2, 4] /. {0 -> LorentzIndex, 1 ->
Momentum}, {n,
0, 15}]
ex2 = (Eps[#[[1]][mu, D],
#[[2]][nu, D], #[[3]][al,
D], #[[4]][be,
D]] &) /@ ex1
FCE /@ ex2
giving
{LCD[al, be, mu, nu], LCD[al, mu, nu][be],
-LCD[be, mu, nu][al],
LCD[mu, nu][al, be], LCD[al, be,
mu][nu], -LCD[al, mu][be, nu],
LCD[be, mu][al, nu], LCD[mu][al, be,
nu], -LCD[al, be, nu][mu],
LCD[al, nu][be, mu], -LCD[be, nu][al,
mu], -LCD[nu][al, be, mu],
LCD[al, be][mu, nu], LCD[al][be, mu,
nu], -LCD[be][al, mu, nu],
LCD[][al, be, mu, nu]}
I’ll nevertheless add this example to our testsuite.