Name: Kyrylo Bondarenko Date: 02/28/15-12:07:32 PM Z


Another thing about DoPolarizationSums. Let
bla = LC[a, b, c, d] FV[w, c] PolarizationVector[w, d] LC[a, b, e, f] FV[w, e] Conjugate[PolarizationVector[w, f]]
tmp = Contract[bla]

Then DoPolarizationSums[bla, w, n] gives
6w^2 (as I expected)
but DoPolarizationSums[tmp, w, n] gives
4w^2 + 2(w^2 - n^2w^2/(nw)^2)
which is very different answer. For example, for n=w the last answer is 4w^2.