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