Well, it seems to be an accuracy issue, so the first question would be
what the values on both sides of the equality are (just start up R and
enter the expressions on te command line). You might have
1) Completely wrong results in complex arithmetic
2) Slightly less than optimal accuracy
3) An underestimated .Machine$double.eps
Case 2) seems most likely, but case 3) has been observed with buggy
compilers that optimize calculations where they shouldn't.
Thanks for the quick response. The two sides of the equality are
definately different. Here's what I'm seeing
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
pi - log(exp(pi*1i)) / 1i
Mod(pi - log(exp(pi*1i)) / 1i)
Mod(pi - log(exp(pi*1i)) / 1i) < .Machine$double.eps
[1] FALSE
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
I get the same thing from R 1.6.2, which I compiled about six months
ago. Is there anything I can/should do to fix this?