atan2 is not consistent with arctan2 of numpy

[LutzGross]

atan2(X,Y) should be X>0 Y>0 = 0 to 90 (in rad) (is correct) X>0 Y<0 = 90 to 180 (is incorrect) X <0, Y>0 = 0 to -90 (is correct) X<0 Y<0 = -90 to -180 (is incorrect)

[LutzGross]

This could provide a fix in util.atan2

s=atan(arg0/(arg1+m))(1-m)+(np.pi/2)(1-2m2) m s+=(wherePositive(arg0)whereNegative(arg1)-whereNegative(arg0)whereNegative(arg1))*np.pi

but maybe this can be implemented better.

[LutzGross]

A better version:

z1=whereZero(arg1, rtol=EPSILON) s=(atan(arg0/(arg1+z1)) + sign(arg0)whereNegative(arg1)np.pi )(1-z1)+sign(arg0)whereNegative(arg1)np.pi/2z1

[aellery]

I have fixed this.