double FowlerAngle(double dy, double dx) {
  double adx, ady;
  int code;
  adx = (dx < 0) ? -dx : dx;
  ady = (dy < 0) ? -dy : dy;
  code = (adx < ady) ? 1 : 0;
  if (dx < 0)
    code += 2;
  if (dy < 0)
    code += 4;
  switch (code) {
  case 0:
    return (dx == 0) ? 0 : ady / adx;
  case 1:
    return (2.0 - (adx / ady));
  case 3:
    return (2.0 + (adx / ady));
  case 2:
    return (4.0 - (ady / adx));
  case 6:
    return (4.0 + (ady / adx));
  case 7:
    return (6.0 - (adx / ady));
  case 5:
    return (6.0 + (adx / ady));
  case 4:
    return (8.0 - (ady / adx));
  }
  return (0.0);
}


x86 assembly (clang, -O0):
<FowlerAngle>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    1142 <FowlerAngle+0x42>
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    114c <FowlerAngle+0x4c>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    118a <FowlerAngle+0x8a>
movsd  -0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1194 <FowlerAngle+0x94>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
xor    %eax,%eax
mov    $0x1,%ecx
ucomisd %xmm1,%xmm0
cmova  %ecx,%eax
mov    %eax,-0x2c(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    11d0 <FowlerAngle+0xd0>
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    11e7 <FowlerAngle+0xe7>
mov    -0x2c(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %rax,-0x48(%rbp)
sub    $0x7,%rax
ja     132f <FowlerAngle+0x22f>
mov    -0x48(%rbp),%rax
lea    0xe1d(%rip),%rcx
movslq (%rcx,%rax,4),%rax
add    %rcx,%rax
jmp    *%rax
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1231 <FowlerAngle+0x131>
jp     1231 <FowlerAngle+0x131>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1240 <FowlerAngle+0x140>
movsd  -0x28(%rbp),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
movsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  0xdb7(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
movsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  0xd97(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
movsd  -0x28(%rbp),%xmm1
divsd  -0x20(%rbp),%xmm1
movsd  0xd6f(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
movsd  -0x28(%rbp),%xmm1
divsd  -0x20(%rbp),%xmm1
movsd  0xd4f(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
movsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  0xd27(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
movsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  0xd07(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
movsd  -0x28(%rbp),%xmm1
divsd  -0x20(%rbp),%xmm1
movsd  0xcdf(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1337 <FowlerAngle+0x237>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<FowlerAngle>:
movapd %xmm0,%xmm2
movapd 0xef4(%rip),%xmm5
movapd %xmm1,%xmm4
xorpd  %xmm5,%xmm4
xorpd  %xmm0,%xmm0
movapd %xmm1,%xmm3
cmpltsd %xmm0,%xmm3
andpd  %xmm3,%xmm4
andnpd %xmm1,%xmm3
orpd   %xmm4,%xmm3
xorpd  %xmm2,%xmm5
movapd %xmm2,%xmm4
cmpltsd %xmm0,%xmm4
andpd  %xmm4,%xmm5
andnpd %xmm2,%xmm4
orpd   %xmm5,%xmm4
xor    %eax,%eax
ucomisd %xmm3,%xmm4
seta   %al
ucomisd %xmm1,%xmm0
lea    0x2(%rax),%ecx
cmovbe %eax,%ecx
ucomisd %xmm2,%xmm0
lea    0x4(%rcx),%eax
cmovbe %ecx,%eax
lea    0xec6(%rip),%rcx
movslq (%rcx,%rax,4),%rax
add    %rcx,%rax
jmp    *%rax
ucomisd %xmm0,%xmm1
jne    11f3 <FowlerAngle+0xf3>
jp     11f3 <FowlerAngle+0xf3>
ret
divsd  %xmm4,%xmm3
movsd  0xea0(%rip),%xmm0
subsd  %xmm3,%xmm0
ret
divsd  %xmm3,%xmm4
movsd  0xe87(%rip),%xmm0
subsd  %xmm4,%xmm0
ret
divsd  %xmm4,%xmm3
addsd  0xe7e(%rip),%xmm3
movapd %xmm3,%xmm0
ret
divsd  %xmm3,%xmm4
movsd  0xe55(%rip),%xmm0
subsd  %xmm4,%xmm0
ret
divsd  %xmm4,%xmm3
addsd  0xe4c(%rip),%xmm3
movapd %xmm3,%xmm0
ret
divsd  %xmm3,%xmm4
addsd  0xe43(%rip),%xmm4
movapd %xmm4,%xmm0
ret
divsd  %xmm4,%xmm3
movsd  0xe2a(%rip),%xmm0
subsd  %xmm3,%xmm0
ret
divsd  %xmm3,%xmm4
movapd %xmm4,%xmm0
ret

LLVM IR (-O0):
define double @FowlerAngle(double %0, double %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca i32
store double %0, double* %4
store double %1, double* %5
%9 = load double, double* %5
%10 = fcmp olt double %9, 0.000000e+00
br i1 %10, label %11, label %14
11: ; preds = %2
%12 = load double, double* %5
%13 = fneg double %12
br label %16
14: ; preds = %2
%15 = load double, double* %5
br label %16
16: ; preds = %14, %11
%17 = phi double [ %13, %11 ], [ %15, %14 ]
store double %17, double* %6
%18 = load double, double* %4
%19 = fcmp olt double %18, 0.000000e+00
br i1 %19, label %20, label %23
20: ; preds = %16
%21 = load double, double* %4
%22 = fneg double %21
br label %25
23: ; preds = %16
%24 = load double, double* %4
br label %25
25: ; preds = %23, %20
%26 = phi double [ %22, %20 ], [ %24, %23 ]
store double %26, double* %7
%27 = load double, double* %6
%28 = load double, double* %7
%29 = fcmp olt double %27, %28
%30 = zext i1 %29 to i64
%31 = select i1 %29, i32 1, i32 0
store i32 %31, i32* %8
%32 = load double, double* %5
%33 = fcmp olt double %32, 0.000000e+00
br i1 %33, label %34, label %37
34: ; preds = %25
%35 = load i32, i32* %8
%36 = add nsw i32 %35, 2
store i32 %36, i32* %8
br label %37
37: ; preds = %34, %25
%38 = load double, double* %4
%39 = fcmp olt double %38, 0.000000e+00
br i1 %39, label %40, label %43
40: ; preds = %37
%41 = load i32, i32* %8
%42 = add nsw i32 %41, 4
store i32 %42, i32* %8
br label %43
43: ; preds = %40, %37
%44 = load i32, i32* %8
switch i32 %44, label %90 [
i32 0, label %45
i32 1, label %55
i32 3, label %60
i32 2, label %65
i32 6, label %70
i32 7, label %75
i32 5, label %80
i32 4, label %85
]
45: ; preds = %43
%46 = load double, double* %5
%47 = fcmp oeq double %46, 0.000000e+00
br i1 %47, label %48, label %49
48: ; preds = %45
br label %53
49: ; preds = %45
%50 = load double, double* %7
%51 = load double, double* %6
%52 = fdiv double %50, %51
br label %53
53: ; preds = %49, %48
%54 = phi double [ 0.000000e+00, %48 ], [ %52, %49 ]
store double %54, double* %3
br label %91
55: ; preds = %43
%56 = load double, double* %6
%57 = load double, double* %7
%58 = fdiv double %56, %57
%59 = fsub double 2.000000e+00, %58
store double %59, double* %3
br label %91
60: ; preds = %43
%61 = load double, double* %6
%62 = load double, double* %7
%63 = fdiv double %61, %62
%64 = fadd double 2.000000e+00, %63
store double %64, double* %3
br label %91
65: ; preds = %43
%66 = load double, double* %7
%67 = load double, double* %6
%68 = fdiv double %66, %67
%69 = fsub double 4.000000e+00, %68
store double %69, double* %3
br label %91
70: ; preds = %43
%71 = load double, double* %7
%72 = load double, double* %6
%73 = fdiv double %71, %72
%74 = fadd double 4.000000e+00, %73
store double %74, double* %3
br label %91
75: ; preds = %43
%76 = load double, double* %6
%77 = load double, double* %7
%78 = fdiv double %76, %77
%79 = fsub double 6.000000e+00, %78
store double %79, double* %3
br label %91
80: ; preds = %43
%81 = load double, double* %6
%82 = load double, double* %7
%83 = fdiv double %81, %82
%84 = fadd double 6.000000e+00, %83
store double %84, double* %3
br label %91
85: ; preds = %43
%86 = load double, double* %7
%87 = load double, double* %6
%88 = fdiv double %86, %87
%89 = fsub double 8.000000e+00, %88
store double %89, double* %3
br label %91
90: ; preds = %43
store double 0.000000e+00, double* %3
br label %91
91: ; preds = %90, %85, %80, %75, %70, %65, %60, %55, %53
%92 = load double, double* %3
ret double %92
}

LLVM IR (-O2):
define double @FowlerAngle(double %0, double %1) {
%3 = fcmp olt double %1, 0.000000e+00
%4 = fneg double %1
%5 = select i1 %3, double %4, double %1
%6 = fcmp olt double %0, 0.000000e+00
%7 = fneg double %0
%8 = select i1 %6, double %7, double %0
%9 = fcmp olt double %5, %8
%10 = zext i1 %9 to i32
%11 = select i1 %9, i32 3, i32 2
%12 = select i1 %3, i32 %11, i32 %10
%13 = or i32 %12, 4
%14 = select i1 %6, i32 %13, i32 %12
switch i32 %14, label %40 [
i32 0, label %15
i32 1, label %19
i32 3, label %22
i32 2, label %25
i32 6, label %28
i32 7, label %31
i32 5, label %34
i32 4, label %37
]
15: ; preds = %2
%16 = fcmp oeq double %1, 0.000000e+00
%17 = fdiv double %8, %5
%18 = select i1 %16, double 0.000000e+00, double %17
br label %41
19: ; preds = %2
%20 = fdiv double %5, %8
%21 = fsub double 2.000000e+00, %20
br label %41
22: ; preds = %2
%23 = fdiv double %5, %8
%24 = fadd double %23, 2.000000e+00
br label %41
25: ; preds = %2
%26 = fdiv double %8, %5
%27 = fsub double 4.000000e+00, %26
br label %41
28: ; preds = %2
%29 = fdiv double %8, %5
%30 = fadd double %29, 4.000000e+00
br label %41
31: ; preds = %2
%32 = fdiv double %5, %8
%33 = fsub double 6.000000e+00, %32
br label %41
34: ; preds = %2
%35 = fdiv double %5, %8
%36 = fadd double %35, 6.000000e+00
br label %41
37: ; preds = %2
%38 = fdiv double %8, %5
%39 = fsub double 8.000000e+00, %38
br label %41
40: ; preds = %2
unreachable
41: ; preds = %37, %34, %31, %28, %25, %22, %19, %15
%42 = phi double [ %39, %37 ], [ %36, %34 ], [ %33, %31 ], [ %30, %28 ], [ %27, %25 ], [ %24, %22 ], [ %21, %19 ], [ %18, %15 ]
ret double %42
}