# [Mathematica] ArcTan[] bug ?

• Mathematica
guerom00
Hello everyone,

I copy-paste here a little bit of code:
Code:
In:= x = (1.5 + 0. I);
y = (0. + 0. I);

(* 1 *) ArcTan[x, y]
(* 2 *) ArcTan[x // Chop, y // Chop]
(* 3 *) -I Log[(x + I y)/Sqrt[x^2 + y^2]]
{y == 0, Re[y] == 0, Im[y] == 0, Im[x] == 0}

Out= -1.5708 + 0. I

Out= 0

Out= 0. + 0. I

Out= {True, True, True, True}

The quantities x and y are calculated elsewhere.
The expression 1 for the arctangent is wrong since it gives me -pi/2 In expression 2, I chop the arguments and the result is right.
Expression 3 is used when x or y is complex and gives the right answer as well.
Finally, some tests…

Any thoughts ?
TIA

I agree it looks wrong, but I think ArcTan[x,y] is assuming x and y are both real. It must have a problem when thy are complex.

guerom00
The help page of ArcTan[x,y] indicates that it can handle complex arguments. It then uses the expression which involves the Log[] I gave.
Problem is: my quantities x and y are not complex according to the tests I make… Which is strange because I thought that the number 0.` is different from 0 (== wise)…

ArcTan[x] should handle complex arguments, but what does ArcTan[x,y] even mean when the arguments are complex?

guerom00

Result : 0
Alternative form : -1.5708 WTF ?! This is clearly broken, no ? Edit: Haha Compare http://www.wolframalpha.com/input/?i=ArcTan[1.+0+I,+0.] and http://www.wolframalpha.com/input/?i=ArcTan[1.+0.+I,+0.]
Difference is subtle…

Last edited:
Staff Emeritus
Gold Member
WTF ?! This is clearly broken, no ?
No and yes.

No, in the sense that you should expect numerical algorithms to misbehave at special values. Yes, in the sense that it's a little surprising these values are special.

What's most likely happening internally is that there is an algebraic identity being used and your point is right on a branch cut of the terms involved. And because your numbers are approximate, different sides of the branch cut wind up being used in different terms.

guerom00
I see now… :)
Well, that doesn't solve my problem: I now have a function of this angle (I use this ArcTan[] to calculate an angle, obviously) which is discontinuous at one point… I'll have to make a specific test to avoid this very point.

Thanks :)