#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double atan(double);
double md_get_pho(int x0, int y0, int x1, int y1) {
  if ((x1 - x0) == 0) {
    if ((y1 - y0) > 0) {
      return (3.141592653589793238462643383279502884 / 2);
    } else {
      return (-3.141592653589793238462643383279502884 / 2);
    }
  } else {
    if ((((x1 - x0) < 0.0) && ((y1 - y0) >= 0.0)) ||
        (((x1 - x0) < 0.0) && ((y1 - y0) < 0.0))) {
      return (3.141592653589793238462643383279502884 +
              atan((double)(y1 - y0) / (double)(x1 - x0)));
    } else {
      return (atan((double)(y1 - y0) / (double)(x1 - x0)));
    }
  }
}


x86 assembly (clang, -O0):
<md_get_pho>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x14(%rbp),%eax
sub    -0xc(%rbp),%eax
cmp    $0x0,%eax
jne    1166 <md_get_pho+0x56>
mov    -0x18(%rbp),%eax
sub    -0x10(%rbp),%eax
cmp    $0x0,%eax
jle    1154 <md_get_pho+0x44>
movsd  0xec6(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121a <md_get_pho+0x10a>
movsd  0xeac(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121a <md_get_pho+0x10a>
mov    -0x14(%rbp),%eax
sub    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    1194 <md_get_pho+0x84>
mov    -0x18(%rbp),%eax
sub    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jae    11c2 <md_get_pho+0xb2>
mov    -0x14(%rbp),%eax
sub    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    11f8 <md_get_pho+0xe8>
mov    -0x18(%rbp),%eax
sub    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    11f8 <md_get_pho+0xe8>
mov    -0x18(%rbp),%eax
sub    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x14(%rbp),%eax
sub    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1030 <atan@plt>
movaps %xmm0,%xmm1
movsd  0xe16(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121a <md_get_pho+0x10a>
mov    -0x18(%rbp),%eax
sub    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x14(%rbp),%eax
sub    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1030 <atan@plt>
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<md_get_pho>:
push   %rbx
mov    %edx,%ebx
sub    %edi,%ebx
je     1138 <md_get_pho+0x28>
cvtsi2sd %ebx,%xmm1
sub    %esi,%ecx
cvtsi2sd %ecx,%xmm0
divsd  %xmm1,%xmm0
call   1030 <atan@plt>
test   %ebx,%ebx
jns    1136 <md_get_pho+0x26>
addsd  0xeca(%rip),%xmm0
pop    %rbx
ret
xor    %eax,%eax
cmp    %esi,%ecx
setg   %al
lea    0xec2(%rip),%rcx
movsd  (%rcx,%rax,8),%xmm0
pop    %rbx
ret

LLVM IR (-O0):
define double @md_get_pho(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca double
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store i32 %3, i32* %9
%10 = load i32, i32* %8
%11 = load i32, i32* %6
%12 = sub nsw i32 %10, %11
%13 = icmp eq i32 %12, 0
br i1 %13, label %14, label %21
14: ; preds = %4
%15 = load i32, i32* %9
%16 = load i32, i32* %7
%17 = sub nsw i32 %15, %16
%18 = icmp sgt i32 %17, 0
br i1 %18, label %19, label %20
19: ; preds = %14
store double 0x3FF921FB54442D18, double* %5
br label %68
20: ; preds = %14
store double 0xBFF921FB54442D18, double* %5
br label %68
21: ; preds = %4
%22 = load i32, i32* %8
%23 = load i32, i32* %6
%24 = sub nsw i32 %22, %23
%25 = sitofp i32 %24 to double
%26 = fcmp olt double %25, 0.000000e+00
br i1 %26, label %27, label %33
27: ; preds = %21
%28 = load i32, i32* %9
%29 = load i32, i32* %7
%30 = sub nsw i32 %28, %29
%31 = sitofp i32 %30 to double
%32 = fcmp oge double %31, 0.000000e+00
br i1 %32, label %45, label %33
33: ; preds = %27, %21
%34 = load i32, i32* %8
%35 = load i32, i32* %6
%36 = sub nsw i32 %34, %35
%37 = sitofp i32 %36 to double
%38 = fcmp olt double %37, 0.000000e+00
br i1 %38, label %39, label %57
39: ; preds = %33
%40 = load i32, i32* %9
%41 = load i32, i32* %7
%42 = sub nsw i32 %40, %41
%43 = sitofp i32 %42 to double
%44 = fcmp olt double %43, 0.000000e+00
br i1 %44, label %45, label %57
45: ; preds = %39, %27
%46 = load i32, i32* %9
%47 = load i32, i32* %7
%48 = sub nsw i32 %46, %47
%49 = sitofp i32 %48 to double
%50 = load i32, i32* %8
%51 = load i32, i32* %6
%52 = sub nsw i32 %50, %51
%53 = sitofp i32 %52 to double
%54 = fdiv double %49, %53
%55 = call double @atan(double %54)
%56 = fadd double 0x400921FB54442D18, %55
store double %56, double* %5
br label %68
57: ; preds = %39, %33
%58 = load i32, i32* %9
%59 = load i32, i32* %7
%60 = sub nsw i32 %58, %59
%61 = sitofp i32 %60 to double
%62 = load i32, i32* %8
%63 = load i32, i32* %6
%64 = sub nsw i32 %62, %63
%65 = sitofp i32 %64 to double
%66 = fdiv double %61, %65
%67 = call double @atan(double %66)
store double %67, double* %5
br label %68
68: ; preds = %57, %45, %20, %19
%69 = load double, double* %5
ret double %69
}

LLVM IR (-O2):
define double @md_get_pho(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = sub nsw i32 %2, %0
%6 = icmp eq i32 %5, 0
br i1 %6, label %7, label %10
7: ; preds = %4
%8 = icmp sgt i32 %3, %1
%9 = select i1 %8, double 0x3FF921FB54442D18, double 0xBFF921FB54442D18
br label %19
10: ; preds = %4
%11 = sitofp i32 %5 to double
%12 = icmp slt i32 %5, 0
%13 = sub nsw i32 %3, %1
%14 = sitofp i32 %13 to double
%15 = fdiv double %14, %11
%16 = tail call double @atan(double %15)
%17 = fadd double %16, 0x400921FB54442D18
%18 = select i1 %12, double %17, double %16
br label %19
19: ; preds = %10, %7
%20 = phi double [ %9, %7 ], [ %18, %10 ]
ret double %20
}