Hi Christian, thanks for your prompt answer from last Friday. Unfortunately I have another little problem: The following lines symbol eps("eps"); idx a(symbol("a"),8) ,b(symbol("b"),8), c(symbol("c"),8),d(symbol("d"),8); ex g = (1+eps*color_T(a)+pow(eps,2)*color_T(b))*(1-eps*color_T(c)+pow(eps,2)*color_T(d)); cout << expand(g) << endl; prints out eps^4*(T.b*T.d)+(T.a*T.d)*eps^3-eps^3*(T.b*T.c)-eps^2*(T.a*T.c) e.g. the terms which would have a "1" from either one of the two factors do not show up. Stefan
Hi! On Mon, Jun 25, 2001 at 04:47:28PM +0200, Stefan Weinzierl wrote:
ex g = (1+eps*color_T(a)+pow(eps,2)*color_T(b))*(1-eps*color_T(c)+pow(eps,2)*color_T(d)); cout << expand(g) << endl;
You have to use color_ONE() instead of just "1" where appropriate (terms of a sum must all be of the same type): ex g = (color_ONE()+eps*color_T(a)+pow(eps,2)*color_T(b)) * (color_ONE()-eps*color_T(c)+pow(eps,2)*color_T(d)); This will produce the right expression on expand(). Bye, Christian -- / Coding on PowerPC and proud of it \/ http://www.uni-mainz.de/~bauec002/
participants (2)
-
Christian Bauer
-
Stefan Weinzierl