Atan2: Difference between revisions

Content deleted Content added

Inline

Latest revision as of 20:37, 15 July 2024

Graph of $\operatorname {atan2} (y,x)$ over $y/x$

In computing and mathematics, the function atan2 is the 2-argument arctangent. By definition, $\theta =\operatorname {atan2} (y,x)$ is the angle measure (in radians, with $-\pi <\theta \leq \pi$ ) between the positive $x$ -axis and the ray from the origin to the point $(x,\,y)$ in the Cartesian plane. Equivalently, $\operatorname {atan2} (y,x)$ is the argument (also called phase oder angle) of the complex number $x+iy.$

The $\operatorname {atan2}$ function first appeared in the programming language Fortran in 1961. It was originally intended to return a correct and unambiguous value for the angle $θ$ in converting from Cartesian coordinates $(x, y)$ to polar coordinates $(r, θ)$ . If $\theta =\operatorname {atan2} (y,x)$ and ${\textstyle r={\sqrt {x^{2}+y^{2}}}}$ , then $x=r\cos \theta$ and $y=r\sin \theta .$

If $x > 0$ , the desired angle measure is ${\textstyle \theta =\operatorname {atan2} (y,x)=\arctan \left(y/x\right).}$ However, when $x < 0$ , the angle $\arctan(y/x)$ is diametrically opposite the desired angle, and ± $π$ (a half turn) must be added to place the point in the correct quadrant.^[1] Using the $\operatorname {atan2}$ function does away with this correction, simplifying code and mathematical formulas.

Motivation

The ordinary single-argument arctangent function only returns angle measures in the interval ${\textstyle {\left(-{\frac {\pi }{2}},+{\frac {\pi }{2}}\right)},}$ and when invoking it to find the angle measure between the $x$ -axis and an arbitrary vector in the Cartesian plane, there is no simple way to indicate a direction in the left half-plane (that is, a point $(x,\,y)$ with $x<0$ ). Diametrically opposite angle measures have the same tangent because $y/x=(-y)/(-x),$ so the tangent $y/x$ is not in itself sufficient to uniquely specify an angle.

To determine an angle measure using the arctangent function given a point or vector $(x,y),$ mathematical formulas or computer code must handle multiple cases; at least one for positive values of $x$ and one for negative values of $x,$ and sometimes additional cases when $y$ is negative or one coordinate is zero. Finding angle measures and converting Cartesian to polar coordinates are common in scientific computing, and this code is redundant and error-prone.

To remedy this, computer programming languages introduced the $atan2$ function, at least as early as the Fortran IV language of the 1960s.^[2] The quantity $atan2(y, x)$ is the angle measure between the $x$ -axis and a ray from the origin to a point $(x, y)$ anywhere in the Cartesian plane. The signs of $x$ and $y$ are used to determine the quadrant of the result and select the correct branch of the multivalued function $Arctan(y / x)$ .

The $atan2$ function is useful in many applications involving Euclidean vectors such as finding the direction from one point to another or converting a rotation matrix to Euler angles.

The $atan2$ function is now included in many other programming languages, and is also commonly found in mathematical formulas throughout science and engineering.

Argument order

In 1961, Fortran introduced the $atan2$ function with argument order $(y,x)$ so that the argument (phase angle) of a complex number is $\operatorname {arg} z=\operatorname {atan2} (\operatorname {Im} z,\operatorname {Re} z).$ This follows the left-to-right order of a fraction written $y/x,$ so that $\operatorname {atan2} (y,x)=\operatorname {atan} (y/x)$ for positive values of $x.$ However, this is the opposite of the conventional component order for complex numbers, $z=x+iy,$ or as coordinates $(\operatorname {Re} z,\operatorname {Im} z).$ See section Definition and computation.

Some other programming languages (see § Realizations of the function in common computer languages) picked the opposite order instead. For example Microsoft Excel uses $\operatorname {Atan2} (x,y),$ OpenOffice Calc uses $\operatorname {arctan2} (x,y),$ and Mathematica uses $\operatorname {ArcTan} [x,y],$ defaulting to one-argument arctangent if called with one argument.

Definition and computation

The function $atan2$ computes the principal value of the argument function applied to the complex number $x + i y$ . That is, $atan2(y, x) = Pr arg(x + i y) = Arg(x + i y)$ . The argument could be changed by an arbitrary multiple of $2π$ (corresponding to a complete turn around the origin) without making any difference to the angle, but to define $atan2$ uniquely one uses the principal value in the range $(-\pi ,\pi ]$ , that is, $- π < atan2(y, x) \leq π$ .

In terms of the standard $arctan$ function, whose range is $(-π/2, π/2)$ , it can be expressed as follows to define a surface that has no discontinuities except along the semi-infinite line x<0 y=0:

$\operatorname {atan2} (y,x)={\begin{cases}\arctan \left({\frac {y}{x}}\right)&{\text{if }}x>0,\\[5mu]\arctan \left({\frac {y}{x}}\right)+\pi &{\text{if }}x<0{\text{ and }}y\geq 0,\\[5mu]\arctan \left({\frac {y}{x}}\right)-\pi &{\text{if }}x<0{\text{ and }}y<0,\\[5mu]+{\frac {\pi }{2}}&{\text{if }}x=0{\text{ and }}y>0,\\[5mu]-{\frac {\pi }{2}}&{\text{if }}x=0{\text{ and }}y<0,\\[5mu]{\text{undefined}}&{\text{if }}x=0{\text{ and }}y=0.\end{cases}}$

A compact expression with four overlapping half-planes is

$\operatorname {atan2} (y,x)={\begin{cases}\arctan \left({\frac {y}{x}}\right)&{\text{if }}x>0,\\[5mu]{\frac {\pi }{2}}-{\arctan }{\bigl (}{\frac {x}{y}}{\bigr )}&{\text{if }}y>0,\\[5mu]-{\frac {\pi }{2}}-{\arctan }{\bigl (}{\frac {x}{y}}{\bigr )}&{\text{if }}y<0,\\[5mu]\arctan \left({\frac {y}{x}}\right)\pm \pi &{\text{if }}x<0,\\[5mu]{\text{undefined}}&{\text{if }}x=0{\text{ and }}y=0.\end{cases}}$

The Iverson bracket notation allows for an even more compact expression:^{[note 1]}

${\begin{aligned}\operatorname {atan2} (y,x)&=\arctan \left({\frac {y}{x}}\right)[x\neq 0]\\[5mu]&\qquad +{\bigl (}1-2[y<0]{\bigr )}\left(\pi [x<0]+{\tfrac {1}{2}}\pi [x=0]\right)\\[5mu]&\qquad +{\text{undefined}}\;\![x=0\wedge y=0]\end{aligned}}$

Formula without apparent conditional construct: $\operatorname {atan2} (y,x)=\lim _{z\to x^{+}}\arctan \left({\frac {y}{z}}\right)+{\frac {\pi }{2}}\operatorname {sgn}(y)\operatorname {sgn}(x)\left(\operatorname {sgn}(x)-1\right)$

The following expression derived from the tangent half-angle formula can also be used to define $atan2$ : $\operatorname {atan2} (y,x)={\begin{cases}2\arctan \left({\frac {y}{{\sqrt {x^{2}+y^{2}}}+x}}\right)&{\text{if }}x>0{\text{ or }}y\neq 0,\\\pi &{\text{if }}x<0{\text{ and }}y=0,\\{\text{undefined}}&{\text{if }}x=0{\text{ and }}y=0.\end{cases}}$ This expression may be more suited for symbolic use than the definition above. However it is unsuitable for general floating-point computational use, as the effect of rounding errors in ${\textstyle {\sqrt {x^{2}+y^{2}}}}$ expand near the region $x < 0, y = 0$ (this may even lead to a division of y by zero).

A variant of the last formula that avoids these inflated rounding errors: $\operatorname {atan2} (y,x)={\begin{cases}2\arctan \left({\frac {y}{{\sqrt {x^{2}+y^{2}}}+x}}\right)&{\text{if }}x>0,\\2\arctan \left({\frac {{\sqrt {x^{2}+y^{2}}}-x}{y}}\right)&{\text{if }}x\leq 0{\text{ and }}y\neq 0,\\\pi &{\text{if }}x<0{\text{ and }}y=0,\\{\text{undefined}}&{\text{if }}x=0{\text{ and }}y=0.\end{cases}}$

Notes:

This produces results in the range $(-π, π]$ .^{[note 2]}
As mentioned above, the principal value of the argument $atan2(y, x)$ can be related to $arctan(y / x)$ by trigonometry. The derivation goes as follows:
If $(x, y) = (r cos θ, r sin θ)$ , then $tan(θ /2) = y / (r + x)$ . It follows that $\operatorname {atan2} (y,x)=\theta =2\,\theta /2=2\arctan {\frac {y}{{\sqrt {x^{2}+y^{2}}}+x}}.$ Note that $\sqrt x 2 + y 2 + x \neq 0$ in the domain in question.

Derivative

As the function $atan2$ is a function of two variables, it has two partial derivatives. At points where these derivatives exist, $atan2$ is, except for a constant, equal to $arctan(y / x)$ . Hence for $x > 0$ or $y \neq 0$ ,

{\begin{aligned}&{\frac {\partial }{\partial x}}\operatorname {atan2} (y,\,x)={\frac {\partial }{\partial x}}\arctan \left({\frac {y}{x}}\right)=-{\frac {y}{x^{2}+y^{2}}},\\[5pt]&{\frac {\partial }{\partial y}}\operatorname {atan2} (y,\,x)={\frac {\partial }{\partial y}}\arctan \left({\frac {y}{x}}\right)={\frac {x}{x^{2}+y^{2}}}.\end{aligned}}

Thus the gradient of atan2 is given by

\nabla {\text{atan2}}(y,x)=\left({-y \over x^{2}+y^{2}},\ {x \over x^{2}+y^{2}}\right).

Informally representing the function $atan2$ as the angle function $θ (x, y) = atan2(y, x)$ (which is only defined up to a constant) yields the following formula for the total differential:

{\begin{aligned}\mathrm {d} \theta &={\frac {\partial }{\partial x}}\operatorname {atan2} (y,\,x)\,\mathrm {d} x+{\frac {\partial }{\partial y}}\operatorname {atan2} (y,\,x)\,\mathrm {d} y\\[5pt]&=-{\frac {y}{x^{2}+y^{2}}}\,\mathrm {d} x+{\frac {x}{x^{2}+y^{2}}}\,\mathrm {d} y.\end{aligned}}

While the function $atan2$ is discontinuous along the negative $x$ -axis, reflecting the fact that angle cannot be continuously defined, this derivative is continuously defined except at the origin, reflecting the fact that infinitesimal (and indeed local) changes in angle can be defined everywhere except the origin. Integrating this derivative along a path gives the total change in angle over the path, and integrating over a closed loop gives the winding number.

In the language of differential geometry, this derivative is a one-form, and it is closed (its derivative is zero) but not exact (it is not the derivative of a 0-form, i.e., a function), and in fact it generates the first de Rham cohomology of the punctured plane. This is the most basic example of such a form, and it is fundamental in differential geometry.

The partial derivatives of $atan2$ do not contain trigonometric functions, making it particularly useful in many applications (e.g. embedded systems) where trigonometric functions can be expensive to evaluate.

Illustrations

This figure shows values of atan2 along selected rays from the origin, labelled at the unit circle. The values, in radians, are shown inside the circle. The diagram uses the standard mathematical convention that angles increase counterclockwise from zero along the ray to the right. Note that the order of arguments is reversed; the function $atan2(y, x)$ computes the angle corresponding to the point $(x, y)$ .

Comparison of arctan and atan2 functions

This figure shows the values of $\arctan(\tan(\theta ))$ along with $\operatorname {atan2} (\sin(\theta ),\cos(\theta ))$ for $0\leq \theta \leq 2\pi$ . Both functions are odd and periodic with periods $\pi$ and $2\pi$ , respectively, and thus can easily be supplemented to any region of real values of $\theta$ . One can clearly see the branch cuts of the $\operatorname {atan2}$ -function at $\theta =\pi$ , and of the $\arctan$ -function at $\theta \in \{{\tfrac {\pi }{2}},\;{\tfrac {3\pi }{2}}\}$ .^[3]

The two figures below show 3D views of respectively $atan2(y, x)$ and $arctan(.mw-parser-output .sfrac{white-space:nowrap}.mw-parser-output .sfrac.tion,.mw-parser-output .sfrac .tion{display:inline-block;vertical-align:-0.5em;font-size:85%;text-align:center}.mw-parser-output .sfrac .num{display:block;line-height:1em;margin:0.0em 0.1em;border-bottom:1px solid}.mw-parser-output .sfrac .den{display:block;line-height:1em;margin:0.1em 0.1em}.mw-parser-output .sr-only{border:0;clip:rect(0,0,0,0);clip-path:polygon(0px 0px,0px 0px,0px 0px);height:1px;margin:-1px;overflow:hidden;padding:0;position:absolute;width:1px}⁠y/x⁠)$ over a region of the plane. Note that for $atan2(y, x)$ , rays in the X/Y-plane emanating from the origin have constant values, but for $arctan(⁠ y / x ⁠)$ lines in the X/Y-plane passing through the origin have constant values. For $x > 0$ , the two diagrams give identical values.

Angle sum and difference identity

The sum or difference of multiple angles to be computed by ⁠ $\operatorname {atan2}$ ⁠ can alternately be computed by composing them as complex numbers. Given two coordinate pairs ⁠ $(x_{1},y_{1})$ ⁠ and ⁠ $(x_{2},y_{2})$ ⁠, their angles from the positive ⁠ $x$ ⁠ axis will be composed (and lengths multiplied) if they are treated as complex numbers and then multiplied together, ⁠ $(x_{1}+iy_{1})(x_{2}+iy_{2})={}$ ⁠⁠ $(x_{1}x_{2}-y_{1}y_{2})+i(y_{1}x_{2}+x_{1}y_{2})$ ⁠. The resulting angle can be found using a single ⁠ $\operatorname {atan2}$ ⁠ operation, so long as the as long as the resulting angle lies in ⁠ $(-\pi ,\pi ]$ ⁠:

\operatorname {atan2} (y_{1},x_{1})\pm \operatorname {atan2} (y_{2},x_{2})=\operatorname {atan2} (y_{1}x_{2}\pm x_{1}y_{2},x_{1}x_{2}\mp y_{1}y_{2}),

and likewise for more than two coordinate pairs. If the composed angle crosses the negative ⁠ $x$ ⁠-axis (i.e. exceeds the range ⁠ $(-\pi ,\pi ]$ ⁠), then the crossings can be counted and the appropriate integer multiple of ⁠ $2\pi$ ⁠ added to the final result to correct it.

This difference formula is frequently used in practice to compute the angle between two planar vectors, since the resulting angle is always in the range ⁠ $(-\pi ,\pi ]$ ⁠.

East-counterclockwise, north-clockwise and south-clockwise conventions, etc.

The $\mathrm {atan2}$ function was originally designed for the convention in pure mathematics that can be termed east-counterclockwise. In practical applications, however, the north-clockwise and south-clockwise conventions are often the norm. In atmospheric sciences, for instance, the wind direction can be calculated using the $\mathrm {atan2}$ function with the east- and north-components of the wind vector as its arguments;^[4] the solar azimuth angle can be calculated similarly with the east- and north-components of the solar vector as its arguments. The wind direction is normally defined in the north-clockwise sense, and the solar azimuth angle uses both the north-clockwise and south-clockwise conventions widely.^[5] These different conventions can be realized by swapping the positions and changing the signs of the x- and y-arguments as follows:

$\mathrm {atan2} (y,x),\;\;\;\;\;$ (East-Counterclockwise Convention)
$\mathrm {atan2} (x,y),\;\;\;\;\;$ (North-Clockwise Convention)
$\mathrm {atan2} (-x,-y)$ . (South-Clockwise Convention)

As an example, let $x_{0}={\frac {\sqrt {3}}{2}}$ and $y_{0}={\frac {1}{2}}$ , then the east-counterclockwise format gives $\mathrm {atan2} (y_{0},x_{0})\cdot {\frac {180}{\pi }}=30^{\circ }$ , the north-clockwise format gives $\mathrm {atan2} (x_{0},y_{0})\cdot {\frac {180}{\pi }}=60^{\circ }$ , and the south-clockwise format gives $\mathrm {atan2} (-x_{0},-y_{0})\cdot {\frac {180}{\pi }}=-120^{\circ }$ .

Changing the sign of the x- and/or y-arguments and/or swapping their positions can create 8 possible variations of the $\mathrm {atan2}$ function and they, interestingly, correspond to 8 possible definitions of the angle, namely, clockwise or counterclockwise starting from each of the 4 cardinal directions, north, east, south and west.

Realizations of the function in common computer languages

The realization of the function differs from one computer language to another:

In Microsoft Excel,^[6] OpenOffice.org Calc, LibreOffice Calc,^[7] Google Spreadsheets,^[8] and iWork Numbers,^[9] the 2-argument arctangent function has the two arguments in the standard sequence $(\operatorname {Re} ,\operatorname {Im} )$ (reversed relative to the convention used in the discussion above).
In Mathematica, the form ArcTan[x, y] is used where the one parameter form supplies the normal arctangent. Mathematica classifies ArcTan[0, 0] as an indeterminate expression.
On most TI graphing calculators (excluding the TI-85 and TI-86), the equivalent function is called R►Pθ and has the arguments $(\operatorname {Re} ,\operatorname {Im} )$ .
On TI-85 the $arg$ function is called angle(x,y) and although it appears to take two arguments, it really only has one complex argument which is denoted by a pair of numbers: $x + i y = (x, y)$ .

The $(\operatorname {Im} ,\operatorname {Re} )$ convention is used by:

The C function atan2, and most other computer implementations, are designed to reduce the effort of transforming cartesian to polar coordinates and so always define atan2(0, 0). On implementations without signed zero, or when given positive zero arguments, it is normally defined as 0. It will always return a value in the range $[-π, π]$ rather than raising an error or returning a NaN (Not a Number).
In Common Lisp, where optional arguments exist, the atan function allows one to optionally supply the x coordinate: (atan y x).^[10]
In Julia, the situation is similar to Common Lisp: instead of atan2, the language has a one-parameter and a two-parameter form for atan.^[11] However, it has many more than two methods, to allow for aggressive optimisation at compile time (see the section "Why don't you compile Matlab/Python/R/… code to Julia?" ^[12]).
For systems implementing signed zero, infinities, or Not a Number (for example, IEEE floating point), it is common to implement reasonable extensions which may extend the range of values produced to include − $π$ and −0 when $y$ = −0. These also may return NaN or raise an exception when given a NaN argument.
In the Intel x86 Architecture assembler code, atan2 is known as the FPATAN (floating-point partial arctangent) instruction.^[13] It can deal with infinities and results lie in the closed interval $[-π, π]$ , e.g. atan2(∞, x) = + $π$ /2 for finite x. Particularly, FPATAN is defined when both arguments are zero:
atan2(+0, +0) = +0;

atan2(+0, −0) = + $π$ ;

atan2(−0, +0) = −0;

atan2(−0, −0) = − $π$ .

This definition is related to the concept of signed zero.

In mathematical writings other than source code, such as in books and articles, the notations Arctan^[14] and Tan⁻¹^[15] have been utilized; these are capitalized variants of the regular arctan and tan⁻¹. This usage is consistent with the complex argument notation, such that $Atan(y, x) = Arg(x + i y)$ .
On HP calculators, treat the coordinates as a complex number and then take the ARG. Or << C->R ARG >> 'ATAN2' STO.
On scientific calculators the function can often be calculated as the angle given when $(x, y)$ is converted from rectangular coordinates to polar coordinates.
Systems supporting symbolic mathematics normally return an undefined value for $atan2(0, 0)$ or otherwise signal that an abnormal condition has arisen.
The free math library FDLIBM (Freely Distributable LIBM) available from netlib has source code showing how it implements atan2 including handling the various IEEE exceptional values.
For systems without a hardware multiplier the function $atan2$ can be implemented in a numerically reliable manner by the CORDIC method. Thus implementations of $atan(y)$ will probably choose to compute $atan2(y, 1)$ .

References

^ "The argument of a complex number" (PDF). Santa Cruz Institute for Particle Physics. Winter 2011.
^ Organick, Elliott I. (1966). A FORTRAN IV Primer. Addison-Wesley. p. 42. Some processors also offer the library function called ATAN2, a function of two arguments (opposite and adjacent).
^ "Wolf Jung: Mandel, software for complex dynamics". www.mndynamics.com. Retrieved 20 April 2018.
^ "Wind Direction Quick Reference". NCAR UCAR Earth Observing Laboratory.
^ Zhang, Taiping; Stackhouse, Paul W.; MacPherson, Bradley; Mikovitz, J. Colleen (2021). "A solar azimuth formula that renders circumstantial treatment unnecessary without compromising mathematical rigor: Mathematical setup, application and extension of a formula based on the subsolar point and atan2 function". Renewable Energy. 172: 1333–1340. Bibcode:2021REne..172.1333Z. doi:10.1016/j.renene.2021.03.047. S2CID 233631040.
^ "Microsoft Excel Atan2 Method". Microsoft. 14 June 2014.
^ "LibreOffice Calc ATAN2". Libreoffice.org.
^ "Functions and formulas – Docs Editors Help". support.google.com.
^ "Numbers' Trigonometric Function List". Apple.
^ "CLHS: Function ASIN, ACOS, ATAN". LispWorks.
^ "Mathematics · The Julia Language". docs.julialang.org.
^ "Frequently Asked Questions · The Julia Language". docs.julialang.org.
^ IA-32 Intel Architecture Software Developer’s Manual. Volume 2A: Instruction Set Reference, A-M, 2004.
^ Burger, Wilhelm; Burge, Mark J. (7 July 2010). Principles of Digital Image Processing: Fundamental Techniques. Springer Science & Business Media. ISBN 978-1-84800-191-6. Retrieved 20 April 2018 – via Google Books.
^ Glisson, Tildon H. (18 February 2011). Introduction to Circuit Analysis and Design. Springer Science & Business Media. ISBN 9789048194438. Retrieved 20 April 2018 – via Google Books.

External links

ATAN2 Online calculator
Java 1.6 SE JavaDoc
atan2 at Everything2
PicBasic Pro solution atan2 for a PIC18F

Other implementations/code for atan2

"Bearing Between Two Points". Archived from the original on 18 November 2020. Retrieved 21 February 2022.

"Arctan and Polar Coordinates". Archived from the original on 18 October 2018. Retrieved 21 February 2022.

"What's 'Arccos'?". Archived from the original on 6 September 2017. Retrieved 21 February 2022.

Notes

^ Assuming the definitions ${\text{undefined}}\cdot 0=0,$ ${\text{undefined}}\cdot 1={\text{undefined}},$ and $z+{\text{undefined}}={\text{undefined}}$ for any $z.$
^ One can apply the periodicity of the result to map to another desired range, e.g. mapping to $[0, 2π)$ by adding $2π$ to the negative results.

[1] "The argument of a complex number" (PDF). Santa Cruz Institute for Particle Physics. Winter 2011.

[2] Organick, Elliott I. (1966). A FORTRAN IV Primer. Addison-Wesley. p. 42. Some processors also offer the library function called ATAN2, a function of two arguments (opposite and adjacent).

[5] "Wolf Jung: Mandel, software for complex dynamics". www.mndynamics.com. Retrieved 20 April 2018.

[6] "Wind Direction Quick Reference". NCAR UCAR Earth Observing Laboratory.

[7] Zhang, Taiping; Stackhouse, Paul W.; MacPherson, Bradley; Mikovitz, J. Colleen (2021). "A solar azimuth formula that renders circumstantial treatment unnecessary without compromising mathematical rigor: Mathematical setup, application and extension of a formula based on the subsolar point and atan2 function". Renewable Energy. 172: 1333–1340. Bibcode:2021REne..172.1333Z. doi:10.1016/j.renene.2021.03.047. S2CID 233631040.

[8] "Microsoft Excel Atan2 Method". Microsoft. 14 June 2014.

[9] "LibreOffice Calc ATAN2". Libreoffice.org.

[10] "Functions and formulas – Docs Editors Help". support.google.com.

[11] "Numbers' Trigonometric Function List". Apple.

[12] "CLHS: Function ASIN, ACOS, ATAN". LispWorks.

[13] "Mathematics · The Julia Language". docs.julialang.org.

[14] "Frequently Asked Questions · The Julia Language". docs.julialang.org.

[15] IA-32 Intel Architecture Software Developer’s Manual. Volume 2A: Instruction Set Reference, A-M, 2004.

[16] Burger, Wilhelm; Burge, Mark J. (7 July 2010). Principles of Digital Image Processing: Fundamental Techniques. Springer Science & Business Media. ISBN 978-1-84800-191-6. Retrieved 20 April 2018 – via Google Books.

[17] Glisson, Tildon H. (18 February 2011). Introduction to Circuit Analysis and Design. Springer Science & Business Media. ISBN 9789048194438. Retrieved 20 April 2018 – via Google Books.

[3] Assuming the definitions ${\text{undefined}}\cdot 0=0,$ ${\text{undefined}}\cdot 1={\text{undefined}},$ and $z+{\text{undefined}}={\text{undefined}}$ for any $z.$

[4] One can apply the periodicity of the result to map to another desired range, e.g. mapping to $[0, 2π)$ by adding $2π$ to the negative results.

[1]

[2]

[note 1]

[note 2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

@@ Line 1: / Line 1: @@
 {{Short description|Arctangent function with two arguments}}
 {{lowercase title}}
-[[File:Atan2definition.svg|thumb|{{math|atan2(''y'', ''x'')}} returns the angle {{mvar|θ}} between the [[ray (geometry)|ray]] to the point {{math|(''x'', ''y'')}} and the positive ''x''-axis, confined to {{open-closed|−''π'', ''π''}}.]]
+[[File:Atan2 differs from arctan.png|thumb|{{math|atan2(''y'', ''x'')}} returns the angle {{mvar|θ}} between the positive ''x''-axis and the [[ray (geometry)|ray]] from the origin to the point {{math|(''x'', ''y'')}}, confined to {{open-closed|−''π'', ''π''}}.]]
 [[File:Arctangent2.svg|thumb|Graph of <math>\operatorname{atan2}(y, x)</math> over <math>y / x</math>]]
-In [[computing]] and [[mathematics]], the [[function (mathematics)|function]] '''atan2''' is the 2-[[Argument of a function|argument]] [[arctangent]]. By definition, <math>\theta = \operatorname{atan2}(y, x)</math> is the [[angle measure]] (in [[radian]]s, with <math>-\pi < \theta \leq \pi</math>) between the positive [[axis (mathematics)|<math>x</math>-axis]] and the [[ray (geometry)|ray]] from the [[origin (mathematics)|origin]] to the point <math>(x,\,y)</math> in the [[Cartesian coordinate system|Cartesian plane]]. Equivalently, <math>\operatorname{atan2}(y, x)</math> is the [[argument (complex analysis)|argument]] (also called ''phase'' or ''angle'') of the [[complex number]] <math>x + iy.</math>
+In [[computing]] and [[mathematics]], the [[function (mathematics)|function]] '''atan2''' is the 2-[[Argument of a function|argument]] [[arctangent]]. By definition, <math>\theta = \operatorname{atan2}(y, x)</math> is the [[angle measure]] (in [[radian]]s, with <math>-\pi < \theta \leq \pi</math>) between the positive [[Cartesian coordinate system|<math>x</math>-axis]] and the [[ray (geometry)|ray]] from the [[origin (mathematics)|origin]] to the point <math>(x,\,y)</math> in the [[Cartesian coordinate system|Cartesian plane]]. Equivalently, <math>\operatorname{atan2}(y, x)</math> is the [[argument (complex analysis)|argument]] (also called ''phase'' or ''angle'') of the [[complex number]] <math>x + iy.</math>
 The <math>\operatorname{atan2}</math> function first appeared in the programming language [[Fortran]] in 1961. It was originally intended to return a correct and unambiguous value for the angle {{mvar|θ}} in converting from Cartesian coordinates {{math|(''x'', ''y'')}} to [[polar coordinates]] {{math|(''r'', ''θ'')}}. If <math>\theta = \operatorname{atan2}(y, x)</math> and <math display=inline>r = \sqrt{x^2 + y^2}</math>, then <math>x = r \cos \theta</math> and <math>y = r \sin \theta.</math>
-If {{math|''x'' > 0}}, the desired angle measure is <math display=inline>\theta = \operatorname{atan2}(y,x) = \arctan\left( y / x \right).</math> However, when {{math|''x'' < 0}}, the angle <math>\arctan(y / x)</math> is [[antipodal point|diametrically opposite]] the desired angle, and ±{{pi}} (a half [[turn (angle)|turn]]) must be added to place the point in the correct [[quadrant (plane geometry)|quadrant]].<ref>http://scipp.ucsc.edu/~haber/ph116A/arg_11.pdf {{Bare URL PDF|date=March 2022}}</ref> Using the <math>\operatorname{atan2}</math> function does away with this correction, simplifying code and mathematical formulas.
+If {{math|''x'' > 0}}, the desired angle measure is <math display=inline>\theta = \operatorname{atan2}(y,x) = \arctan\left( y / x \right).</math> However, when {{math|''x'' < 0}}, the angle <math>\arctan(y / x)</math> is [[antipodal point|diametrically opposite]] the desired angle, and ±{{pi}} (a half [[turn (angle)|turn]]) must be added to place the point in the correct [[quadrant (plane geometry)|quadrant]].<ref>{{cite web|url=http://scipp.ucsc.edu/~haber/ph116A/arg_11.pdf|title=The argument of a complex number|date=Winter 2011|publisher=Santa Cruz Institute for Particle Physics}}</ref> Using the <math>\operatorname{atan2}</math> function does away with this correction, simplifying code and mathematical formulas.
 == Motivation ==
 [[File:Atan2 argument sign graph.svg|thumb|Graph of the tangent function from −{{pi}} to +{{pi}} with the corresponding signs of ''y''/''x''. The green arrows point to the results of atan2(−1,&nbsp;−1) and atan2(1,&nbsp;1).]]
-The ordinary single-argument [[arctangent]] function only returns angle measures in the interval <math>{\left[-\tfrac12\pi, +\tfrac12\pi\right]},</math> and when invoking it to find the angle measure between the {{mvar|x}}-axis and an arbitrary vector in the [[Cartesian coordinate system|Cartesian]] plane, there is no simple way to indicate a direction in the left half-plane (that is, a point <math>(x,\,y)</math> with <math>x < 0</math>). [[Antipodal point|Diametrically opposite]] angle measures have the same tangent because <math>y/x = (-y) / (-x),</math> so the tangent <math>y/x</math> is not in itself sufficient to uniquely specify an angle.
+The ordinary single-argument [[arctangent]] function only returns angle measures in the interval <math display="inline">{\left(-\frac{\pi}{2}, +\frac{\pi}{2}\right)},</math> and when invoking it to find the angle measure between the {{mvar|x}}-axis and an arbitrary vector in the [[Cartesian coordinate system|Cartesian]] plane, there is no simple way to indicate a direction in the left half-plane (that is, a point <math>(x,\,y)</math> with <math>x < 0</math>). [[Antipodal point|Diametrically opposite]] angle measures have the same tangent because <math>y/x = (-y) / (-x),</math> so the tangent <math>y/x</math> is not in itself sufficient to uniquely specify an angle.
 To determine an angle measure using the arctangent function given a point or vector <math>(x, y),</math> mathematical formulas or computer code must handle multiple cases; at least one for positive values of <math>x</math> and one for negative values of <math>x,</math> and sometimes additional cases when <math>y</math> is negative or one coordinate is zero. Finding angle measures and converting Cartesian to [[Polar coordinate system|polar coordinates]] are common in scientific computing, and this code is redundant and error-prone.
@@ Line 40: / Line 40: @@
 The function {{math|atan2}} computes the [[principal value]] of the [[argument (complex analysis)|argument]] function applied to the [[complex number]] {{math|''x'' + ''i''&hairsp;''y''}}. That is, {{math|1=atan2(''y'', ''x'') = Pr arg(''x'' + ''i''&hairsp;''y'') = Arg(''x'' + ''i''&hairsp;''y'')}}. The argument could be changed by an arbitrary multiple of {{math|2π}} (corresponding to a complete turn around the origin) without making any difference to the angle, but to define {{math|atan2}} uniquely one uses the principal value in the [[interval (mathematics)|range]] <math>( -\pi, \pi ]</math>, that is, {{math|−''π'' < atan2(''y'', ''x'') ≤ ''π''}}.
-In terms of the standard {{math|arctan}} function, whose range is {{open-closed|−π/2, π/2}}, it can be expressed as follows to define a surface that has no discontinuities except along the semi-infinite line x<0 y=0:
+In terms of the standard {{math|arctan}} function, whose range is {{open-open|−π/2, π/2}}, it can be expressed as follows to define a surface that has no discontinuities except along the semi-infinite line x<0 y=0:
 <math display="block"> \operatorname{atan2}(y, x) =
@@ Line 133: / Line 133: @@
 The two figures below show 3D views of respectively {{math|atan2(''y'', ''x'')}} and {{math|arctan({{sfrac|''y''|''x''}})}} over a region of the plane. Note that for {{math|atan2(''y'', ''x'')}}, ''rays'' in the ''X''/''Y''-plane emanating from the origin have constant values, but for {{math|arctan({{sfrac|''y''|''x''}})}} ''lines'' in the ''X''/''Y''-plane passing through the origin have constant values. For {{math|''x'' > 0}}, the two diagrams give identical values.
-{| class="wikitable" style="text-align: center;
+{| class="wikitable" style="text-align: center"
 | [[File:Atan2 diagram.svg|300px]]
 | [[File:Atan diagram.svg|300px]]
@@ Line 141: / Line 141: @@
 {{Main|List of trigonometric identities#Angle sum and difference identities}}
+The sum or difference of multiple angles to be computed by {{tmath|\operatorname{atan2} }} can alternately be computed by composing them as [[complex number]]s. Given two coordinate pairs {{tmath|(x_1, y_1)}} and {{tmath|(x_2, y_2)}}, their angles from the positive {{tmath|x}} axis will be composed (and lengths multiplied) if they are treated as complex numbers and then multiplied together, {{tmath|1=(x_1 + iy_1)(x_2 + iy_2)={} }}{{tmath|(x_1x_2 - y_1y_2) + i(y_1x_2 + x_1y_2)}}. The resulting angle can be found using a single {{tmath|\operatorname{atan2} }} operation, so long as the as long as the resulting angle lies in {{tmath|(-\pi, \pi]}}:
-Sums of <math>\operatorname{atan2}</math> may be collapsed into a single operation according to the following identity
-:<math>\operatorname{atan2} (y_1, x_1) \pm \operatorname{atan2} (y_2, x_2) = \operatorname{atan2} (y_1 x_2 \pm y_2 x_1, x_1 x_2 \mp y_1 y_2)</math>
+:<math>\operatorname{atan2} (y_1, x_1) \pm \operatorname{atan2} (y_2, x_2) = \operatorname{atan2} (y_1 x_2 \pm x_1 y_2, x_1 x_2 \mp y_1 y_2),</math>
+and likewise for more than two coordinate pairs. If the composed angle crosses the negative {{tmath|x}}-axis (i.e. exceeds the range {{tmath|(-\pi, \pi]}}), then the crossings can be counted and the appropriate integer multiple of {{tmath|2\pi}} added to the final result to correct it.
-...provided that <math>\operatorname{atan2} (y_1, x_1) \pm \operatorname{atan2} (y_2, x_2) \in (-\pi, \pi]</math>.
+This difference formula is frequently used in practice to compute the angle between two planar [[Euclidean vector|vector]]s, since the resulting angle is always in the range {{tmath|(-\pi, \pi]}}.
-The proof involves considering two cases, one where <math>y_2 \neq 0</math> or <math>x_2 > 0</math> and one where <math>y_2 = 0</math> and <math>x_2 < 0</math>.
-We only consider the case where <math>y_2 \neq 0</math> or <math>x_2 > 0</math>. To start, we make the following observations:
-# <math>-\operatorname{atan2}(y,x) = \operatorname{atan2}(-y,x)</math> provided that <math>y \neq 0</math> or <math>x > 0</math>.
-# <math>\operatorname{Arg} (x + i y) = \operatorname{atan2} (y, x)</math>, where <math>\operatorname{Arg}</math> is the [[Argument (complex analysis)#Computation|complex argument function]].
-# <math>\theta = \operatorname{Arg} e^{i \theta}</math> whenever <math>\theta \in (-\pi, \pi]</math>, a consequence of [[Euler's formula]].
-# <math>\operatorname{Arg} (e^{i \operatorname{Arg} \zeta_1} e^{i \operatorname{Arg} \zeta_2}) = \operatorname{Arg} (\zeta_1 \zeta_2)</math>.
-To see (4), we have the [[Argument (complex analysis)#Identities|identity]] <math>e^{i \operatorname{Arg} \zeta} = \bar{\zeta}</math> where <math>\bar{\zeta} = \zeta / \left|\zeta\right|</math>, hence <math>\operatorname{Arg} (e^{i \operatorname{Arg} \zeta_1} e^{i \operatorname{Arg} \zeta_2}) = \operatorname{Arg} (\bar{\zeta_1} \bar{\zeta_2})</math>.  Furthermore, since <math>\operatorname{Arg} \zeta = \operatorname{Arg} a \zeta</math> for any positive real value <math>a</math>, then if we let <math>\zeta = \zeta_1 \zeta_2</math> and <math>a = \frac{1}{\left|\zeta_1\right|\left|\zeta_2\right|}</math> then we have <math>\operatorname{Arg} (\bar{\zeta_1} \bar{\zeta_2}) = \operatorname{Arg} (\zeta_1 \zeta_2)</math>.
-From these observations have following equivalences:
-:<math>\begin{align}
-\operatorname{atan2} (y_1, x_1) \pm \operatorname{atan2} (y_2, x_2)
-&{} = \operatorname{atan2} (y_1, x_1) + \operatorname{atan2} (\pm y_2, x_2) &  \text{by (1)} \\
-&{} = \operatorname{Arg} (x_1 + i y_1) + \operatorname{Arg} (x_2 \pm i y_2) &  \text{by (2)}   \\
-&{} = \operatorname{Arg} e^{i (\operatorname{Arg} (x_1 + i y_1) + \operatorname{Arg} (x_2 \pm i y_2))} & \text{by (3)} \\
-&{} = \operatorname{Arg} (e^{i \operatorname{Arg} (x_1 + i y_1)} e^{i \operatorname{Arg} (x_2 \pm i y_2)}) \\
-&{} = \operatorname{Arg} ((x_1 + i y_1) ( x_2 \pm i y_2)) & \text{by (4)} \\
-&{} = \operatorname{Arg} (x_1 x_2 \mp y_1 y_2 + i (y_1 x_2 \pm y_2 x_1)) \\
-&{} = \operatorname{atan2} (y_1 x_2 \pm y_2 x_1, x_1 x_2 \mp y_1 y_2) & \text{by (2)}
-\end{align}</math>
-'''Corollary''': if <math>(y_1, x_1)</math> and <math>(y_2, x_2)</math> are 2-dimensional vectors, the difference formula is frequently used in practice to compute the angle between those vectors with the help of <math>\operatorname{atan2}</math>, since the resulting computation behaves benign in the range <math>(-\pi, \pi]</math> and can thus be used without range checks in many practical situations.
 == East-counterclockwise, north-clockwise and south-clockwise conventions, etc. ==
-The <math>\mathrm{atan2}</math> function was originally designed for the convention in pure mathematics that can be termed ''east-counterclockwise''. In practical applications, however, the ''north-clockwise'' and ''south-clockwise'' conventions are often the norm. In atmospheric sciences, for instance, the [[wind direction]] can be calculated using the <math>\mathrm{atan2}</math> function with the east- and north-components of the wind vector as its arguments;<ref>Wind Direction Quick Reference, NCAR UCAR Earth Observing Laboratory. https://www.eol.ucar.edu/content/wind-direction-quick-reference</ref> the [[solar azimuth angle]] can be calculated similarly with the east- and north-components of the solar vector as its arguments. The wind direction is normally defined in the north-clockwise sense, and the solar azimuth angle uses both the north-clockwise and south-clockwise conventions widely.<ref>{{cite journal|doi=10.1016/j.renene.2021.03.047|title=A solar azimuth formula that renders circumstantial treatment unnecessary without compromising mathematical rigor: Mathematical setup, application and extension of a formula based on the subsolar point and atan2 function|year=2021|last1=Zhang|first1=Taiping|last2=Stackhouse|first2=Paul W.|last3=MacPherson|first3=Bradley|last4=Mikovitz|first4=J. Colleen|journal=Renewable Energy|volume=172|pages=1333–1340|s2cid=233631040}}</ref> These different conventions can be realized by swapping the positions and changing the signs of the x- and y-arguments as follows:
+The <math>\mathrm{atan2}</math> function was originally designed for the convention in pure mathematics that can be termed ''east-counterclockwise''. In practical applications, however, the ''north-clockwise'' and ''south-clockwise'' conventions are often the norm. In atmospheric sciences, for instance, the [[wind direction]] can be calculated using the <math>\mathrm{atan2}</math> function with the east- and north-components of the wind vector as its arguments;<ref>{{cite web|url=https://www.eol.ucar.edu/content/wind-direction-quick-reference|title=Wind Direction Quick Reference|publisher=NCAR UCAR Earth Observing Laboratory}}</ref> the [[solar azimuth angle]] can be calculated similarly with the east- and north-components of the solar vector as its arguments. The wind direction is normally defined in the north-clockwise sense, and the solar azimuth angle uses both the north-clockwise and south-clockwise conventions widely.<ref>{{cite journal|doi=10.1016/j.renene.2021.03.047|title=A solar azimuth formula that renders circumstantial treatment unnecessary without compromising mathematical rigor: Mathematical setup, application and extension of a formula based on the subsolar point and atan2 function|year=2021|last1=Zhang|first1=Taiping|last2=Stackhouse|first2=Paul W.|last3=MacPherson|first3=Bradley|last4=Mikovitz|first4=J. Colleen|journal=Renewable Energy|volume=172|pages=1333–1340|s2cid=233631040|doi-access=free|bibcode=2021REne..172.1333Z }}</ref> These different conventions can be realized by swapping the positions and changing the signs of the x- and y-arguments as follows:
 * <math>\mathrm{atan2}(y, x),\;\;\;\;\;</math> (East-Counterclockwise Convention)
 * <math>\mathrm{atan2}(x, y),\;\;\;\;\;</math> (North-Clockwise Convention)
@@ Line 182: / Line 158: @@
 As an example, let <math>x_{0}=\frac{\sqrt{3}}{2}</math> and <math>y_{0}=\frac{1}{2}</math>, then the east-counterclockwise format gives <math>\mathrm{atan2}(y_{0}, x_{0})\cdot\frac{180}{\pi}=30^{\circ}</math>, the north-clockwise format gives <math>\mathrm{atan2}(x_{0}, y_{0})\cdot\frac{180}{\pi}=60^{\circ}</math>, and the south-clockwise format gives <math>\mathrm{atan2}(-x_{0}, -y_{0})\cdot\frac{180}{\pi}=-120^{\circ}</math>.
-Apparently, changing the sign of the x- and/or y-arguments and swapping their positions can create 8 possible variations of the <math>\mathrm{atan2}</math> function and they, interestingly, correspond to 8 possible definitions of the angle, namely, clockwise or counterclockwise starting from each of the 4 [[cardinal direction]]s, north, east, south and west.
+Changing the sign of the x- and/or y-arguments and/or swapping their positions can create 8 possible variations of the <math>\mathrm{atan2}</math> function and they, interestingly, correspond to 8 possible definitions of the angle, namely, clockwise or counterclockwise starting from each of the 4 [[cardinal direction]]s, north, east, south and west.
 == Realizations of the function in common computer languages ==
 The realization of the function differs from one computer language to another:
-* In [[Microsoft Excel]],<ref>{{cite web|url=http://msdn2.microsoft.com/en-us/library/bb238940.aspx|title=Microsoft Excel Atan2 Method|publisher=Microsoft}}</ref> [[OpenOffice.org Calc]], [[LibreOffice Calc]],<ref>{{cite web|url=https://help.libreoffice.org/Calc/Mathematical_Functions#ATAN2|title=LibreOffice Calc ATAN2|publisher=Libreoffice.org}}</ref> [[Google Docs|Google Spreadsheets]],<ref>{{Cite web|url=https://support.google.com/docs/topic/1361471 |title=Functions and formulas – Docs Editors Help|website=support.google.com}}</ref> [[Numbers (spreadsheet)|iWork Numbers]],<ref>{{cite web|url=https://www.apple.com/mac/numbers/compatibility/functions.html#trigonometric |title=Numbers' Trigonometric Function List|publisher=Apple }}</ref> and [[SQL:2008|ANSI SQL:2008 standard]],<ref>{{cite web|url=http://www.info.teradata.com/HTMLPubs/DB_TTU_15_00/index.html#page/SQL_Reference/B035_1145_015K/Arithmetic.062.225.html#ww15697556 |title=ANSI SQL:2008 standard |publisher=Teradata |url-status=dead |archive-url=https://web.archive.org/web/20150820094929/http://www.info.teradata.com/HTMLPubs/DB_TTU_15_00/index.html |archive-date=2015-08-20 }}</ref> the 2-argument arctangent function has the two arguments in the standard sequence <math>(\operatorname{Re}, \operatorname{Im})</math> (reversed relative to the convention used in the discussion above).
+* In [[Microsoft Excel]],<ref>{{cite web|url=http://msdn2.microsoft.com/en-us/library/bb238940.aspx|title=Microsoft Excel Atan2 Method|date=14 June 2014 |publisher=Microsoft}}</ref> [[OpenOffice.org Calc]], [[LibreOffice Calc]],<ref>{{cite web|url=https://help.libreoffice.org/Calc/Mathematical_Functions#ATAN2|title=LibreOffice Calc ATAN2|publisher=Libreoffice.org}}</ref> [[Google Docs|Google Spreadsheets]],<ref>{{Cite web|url=https://support.google.com/docs/topic/1361471 |title=Functions and formulas – Docs Editors Help|website=support.google.com}}</ref> and [[Numbers (spreadsheet)|iWork Numbers]],<ref>{{cite web|url=https://www.apple.com/mac/numbers/compatibility/functions.html#trigonometric |title=Numbers' Trigonometric Function List|publisher=Apple }}</ref> the 2-argument arctangent function has the two arguments in the standard sequence <math>(\operatorname{Re}, \operatorname{Im})</math> (reversed relative to the convention used in the discussion above).
 * In [[Mathematica]], the form <code>ArcTan[''x'',&hairsp;''y'']</code> is used where the one parameter form supplies the normal arctangent. Mathematica classifies <code>ArcTan[0,&hairsp;0]</code> as an indeterminate expression.
 * On most TI graphing calculators (excluding the [[TI-85]] and [[TI-86]]), the equivalent function is called '''R►Pθ''' and has the arguments <math>(\operatorname{Re}, \operatorname{Im})</math>.

v t e Trigonometric and hyperbolic functions
Groups	Trigonometric Sine and cosine Inverse trigonometric Hyperbolic Inverse hyperbolic
Other	Versine Exsecant Jyā, koti-jyā and utkrama-jyā atan2