Calculating ATAN2

[Albert Chan]

How about using acos for atan2(y,x) ?
This covered edges cases, except for (y,x) == (0,0)

Code:

function sign(y) return signbit(y) and -1 or 1 end
function myatan2(y,x) return acos(x/hypot(x,y)) * sign(y) end

lua> function test(y,x) return atan2(y,x), myatan2(y,x) end

lua> test(4,3) --> 0.9272952180016122      0.9272952180016123
lua> test(4,-3) --> 2.214297435588181      2.214297435588181
lua> test(-4,3) --> -0.9272952180016122      -0.9272952180016123
lua> test(-4,-3) --> -2.214297435588181      -2.214297435588181

X axis:

lua> test(0, 1) --> 0      0
lua> test(0,-1) --> 3.141592653589793      3.141592653589793

Y axis:

lua> test(1, 0) --> 1.5707963267948966      1.5707963267948968
lua> test(-1,0) --> -1.5707963267948966      -1.5707963267948968