Normal() again

Do Hoang Son dhson at ThEP.Physik.Uni-Mainz.DE
Fri Nov 19 18:41:07 CET 1999


Hi all,

Consider this:

When one types in ginsh:

> normal(q10^2-q10^4*(q10^2-q21^2)^(-1)+q21^2*q10^2*(q10^2-q21^2)^(-1));
power::eval(): division by zero

The result, however, should be zero.

Maybe the next bug is related to the above one:

>
normal(q10^2*(I*Pi^2*(q10^2-q21^2)^(-1/2)*q10^2*(-(1/2)*(q10^2-q21^2)^(-1/2)*q10^2+(1/2)*q21^2*(q10^2-q21^2)^(-1/2))+1/2*I*Pi^2*q10^2));
q10^2*(1/2*I*Pi^2*q21^2*(q10^2-q21^2)^(-1)*q10^2-1/2*I*Pi^2*(q10^2-q21^2)^(-1)*q10^4+1/2*I*Pi^2*q10^2)
> normal(");
0

The result is correct, but it should already found it after the first
line.

Son.




More information about the GiNaC-devel mailing list