# Fowler Angles

Original by Rob Fowler
Edited by Paul Bourke
February 1997

If you want a monotonic function of the angle, rather than needing the value of the angle itself, you can avoid the arcsine and arctangent computation by using Fowler angles. These use dx and dy to return a value that can be used for angular comparison routines.

```/*
This function is due to Rob Fowler.  Given dy and dx between 2 points
A and B, we calculate a number in [0.0, 8.0) which is a monotonic
function of the direction from A to B.

(0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0) correspond to
(  0,  45,  90, 135, 180, 225, 270, 315, 360) degrees, measured
counter-clockwise from the positive x axis.
*/
double FowlerAngle(dy,dx)
double dy,dx;
{
int    code;        /* Angular Region Classification Code */

adx = (dx < 0) ? -dx : dx;  /* Compute the absolute values. */
ady = (dy < 0) ? -dy : dy;

if (dx < 0)  code += 2;
if (dy < 0)  code += 4;

switch (code)
{
case 0: return (dx==0) ? 0 : ady/adx;  /* [  0, 45] */