Which version are you using?
The latest version does not exhibit this bug.
I shortened your code a bit, and I hope I did not dot introduce
errors.
Regards,
Rolf Mertig
[rolf@colossos rolf]$ math
Mathematica 4.2 for Linux
Copyright 1988-2002 Wolfram Research, Inc.
-- Motif graphics initialized --
In[1]:= !!fct
<<HighEnergyPhysics`fc`;
ScalarProduct[p, p]=0; ScalarProduct[l, l]=0; ScalarProduct[o, o]=0;
BB = GAD[rho] . GSD[o] . GAD[rho] . (GSD[l] + GSD[o]) .
GAD[mu] . (1 - GA[5]) . GSD[p] . GAD[nu] . (1 - GA[5]) . (GSD[l] + GSD[o]);
trBB = Tr[BB];
contrazioneBB = Contract[trBB*I*LCD[mu,nu,al,bet] FVD[p,al] FVD[o+l-p, bet]];
contrazioneBB = ChangeDimension[contrazioneBB, 4];
contrazioneBB=-contrazioneBB/4/2/SP[l,o]/2/SP[l,o];
res1 = Expand[contrazioneBB /. D -> 4 + 2*eps]//FCI;
CC = GSD[o] . GAD[rho] . (GSD[l] + GSD[o]) . GAD[mu] . (1 - GA[5]) .
GSD[p] . GAD[nu] . (1 - GA[5]) . (GSD[l] + GSD[o]) . GAD[rho];
trCC = Tr[CC];
contrazioneCC = Contract[trCC*I*LCD[mu,nu,al,bet] FVD[p,al] FVD[o+l-p, bet]];
contrazioneCC = ChangeDimension[contrazioneCC, 4];
contrazioneCC=-contrazioneCC/4/2/SP[l,o]/2/SP[l,o];
res2 = Expand[contrazioneCC /. D -> 4 + 2*eps]//FCI;
Expand[res1 - res2]
In[1]:= <<fct
FeynCalc4.1.2 Type ?FeynCalc for help or visit http://www.feyncalc.org
$PrePrint is set to FeynCalcForm. Use FI and FC to change the display format.
Out[1]= 0
This archive was generated by hypermail 2b29 : 09/04/20-12:55:05 AM Z CEST