#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 */
scalar_t__ floor(double);
double sqrt(double);
int hex_xy2nm(double p, double x, double y, int *nx, int *my) {
  int m, m2, n;
  double d, dx, dy, xr, yr2;
  double x0 = 0.0, y0 = 0.0;
  d = sqrt(3.0) / 2.0 * p;
  dx = 3.0 / 4.0 * d;
  dy = p;
  m2 = (int)floor((y - y0) * 2.0 / dy);
  n = (int)floor((x - x0) / dx);
  xr = x - x0 - n * dx;
  yr2 = y - y0 - m2 * dy / 2.0;
  if (n % 2 == 0 && m2 % 2 == 0) {
    if ((yr2 + sqrt(3.0) * (xr - 0.5 * d)) < 0) {
      *nx = n;
      m = m2;
    } else {
      *nx = n + 1;
      m = m2;
    }
  }
  if (n % 2 == 0 && m2 % 2 != 0) {
    if ((yr2 - sqrt(3.0) * (xr - 0.25 * d)) > 0) {
      *nx = n;
      m = m2 + 1;
    } else {
      *nx = n + 1;
      m = m2 - 1;
    }
  }
  if (n % 2 != 0 && m2 % 2 == 0) {
    if ((yr2 - sqrt(3.0) * (xr - 0.25 * d)) > 0) {
      *nx = n;
      m = m2;
    } else {
      *nx = n + 1;
      m = m2;
    }
  }
  if (n % 2 != 0 && m2 % 2 != 0) {
    if ((yr2 + sqrt(3.0) * (xr - 0.5 * d)) < 0) {
      *nx = n;
      m = m2 - 1;
    } else {
      *nx = n + 1;
      m = m2 + 1;
    }
  }
  *my = m >> 1;
  return (*nx) * 10000 + (*my);
}


x86 assembly (clang, -O0):
<hex_xy2nm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
mov    %rdi,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0xeb6(%rip),%xmm0
call   1030 <sqrt@plt>
movsd  0xe99(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe8b(%rip),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
movsd  0xe5d(%rip),%xmm1
mulsd  %xmm1,%xmm0
divsd  -0x50(%rbp),%xmm0
call   1040 <floor@plt>
mov    %eax,-0x30(%rbp)
movsd  -0x10(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
call   1040 <floor@plt>
mov    %eax,-0x34(%rbp)
movsd  -0x10(%rbp),%xmm1
subsd  -0x68(%rbp),%xmm1
cvtsi2sdl -0x34(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x18(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
movsd  0xde1(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x34(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    12cd <hex_xy2nm+0x1ad>
mov    -0x30(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    12cd <hex_xy2nm+0x1ad>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xdaa(%rip),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm3
movsd  -0x40(%rbp),%xmm2
movsd  0xd93(%rip),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    12b6 <hex_xy2nm+0x196>
mov    -0x34(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    12c8 <hex_xy2nm+0x1a8>
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    12cd <hex_xy2nm+0x1ad>
mov    -0x34(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1388 <hex_xy2nm+0x268>
mov    -0x30(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
je     1388 <hex_xy2nm+0x268>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  0xd09(%rip),%xmm0
call   1030 <sqrt@plt>
movsd  -0x80(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm3
movsd  -0x40(%rbp),%xmm2
movsd  0xcfd(%rip),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    136e <hex_xy2nm+0x24e>
mov    -0x34(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1383 <hex_xy2nm+0x263>
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1388 <hex_xy2nm+0x268>
mov    -0x34(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
je     1443 <hex_xy2nm+0x323>
mov    -0x30(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1443 <hex_xy2nm+0x323>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xc4b(%rip),%xmm0
call   1030 <sqrt@plt>
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x58(%rbp),%xmm3
movsd  -0x40(%rbp),%xmm2
movsd  0xc3c(%rip),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    142c <hex_xy2nm+0x30c>
mov    -0x34(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    143e <hex_xy2nm+0x31e>
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    1443 <hex_xy2nm+0x323>
mov    -0x34(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
je     14f0 <hex_xy2nm+0x3d0>
mov    -0x30(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
je     14f0 <hex_xy2nm+0x3d0>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0xb90(%rip),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x90(%rbp),%xmm0
ff
movsd  -0x58(%rbp),%xmm3
movsd  -0x40(%rbp),%xmm2
movsd  0xb76(%rip),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    14d6 <hex_xy2nm+0x3b6>
mov    -0x34(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    14eb <hex_xy2nm+0x3cb>
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    14f0 <hex_xy2nm+0x3d0>
mov    -0x2c(%rbp),%ecx
sar    $0x1,%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x20(%rbp),%rax
imul   $0x2710,(%rax),%eax
mov    -0x28(%rbp),%rcx
add    (%rcx),%eax
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<hex_xy2nm>:
push   %r15
push   %r14
push   %rbx
sub    $0x30,%rsp
mov    %rsi,%r14
mov    %rdi,%r15
movsd  %xmm2,0x28(%rsp)
movsd  %xmm1,0x8(%rsp)
movapd %xmm0,%xmm1
movsd  %xmm0,0x20(%rsp)
movsd  0xec3(%rip),%xmm3
mulsd  %xmm0,%xmm3
movsd  %xmm3,0x18(%rsp)
movsd  0xeb9(%rip),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,0x10(%rsp)
movapd %xmm2,%xmm0
addsd  %xmm2,%xmm0
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
mov    %rax,%rbx
movsd  0x8(%rsp),%xmm0
divsd  0x10(%rsp),%xmm0
call   1030 <floor@plt>
movsd  0x8(%rsp),%xmm4
movsd  0x18(%rsp),%xmm3
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
mulsd  0x10(%rsp),%xmm1
mulsd  0x20(%rsp),%xmm0
mulsd  0xe64(%rip),%xmm0
subsd  %xmm1,%xmm4
addsd  0x28(%rsp),%xmm0
mov    %eax,%ecx
or     %ebx,%ecx
test   $0x1,%cl
jne    11ed <hex_xy2nm+0xdd>
movsd  0xe49(%rip),%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm4,%xmm1
mulsd  0xe41(%rip),%xmm1
addsd  %xmm0,%xmm1
xorpd  %xmm2,%xmm2
xor    %ecx,%ecx
ucomisd %xmm1,%xmm2
setbe  %cl
add    %eax,%ecx
mov    %ecx,(%r15)
mov    %ebx,%r8d
and    $0x1,%r8d
mov    %eax,%edx
mov    %ebx,%esi
and    $0x1,%dl
jne    123f <hex_xy2nm+0x12f>
mov    %ebx,%esi
test   %r8d,%r8d
je     123f <hex_xy2nm+0x12f>
movsd  0xe14(%rip),%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm4,%xmm1
mulsd  0xe0c(%rip),%xmm1
addsd  %xmm0,%xmm1
xorpd  %xmm2,%xmm2
xor    %esi,%esi
xor    %edi,%edi
ucomisd %xmm2,%xmm1
seta   %sil
setbe  %dil
add    %eax,%edi
mov    %edi,(%r15)
lea    (%rbx,%rsi,2),%esi
add    $0xffffffff,%esi
test   $0x1,%bl
sete   %cl
and    %cl,%dl
cmp    $0x1,%dl
jne    127c <hex_xy2nm+0x16c>
movsd  0xdcc(%rip),%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm4,%xmm1
mulsd  0xdc4(%rip),%xmm1
addsd  %xmm0,%xmm1
xorpd  %xmm2,%xmm2
xor    %ecx,%ecx
ucomisd %xmm2,%xmm1
setbe  %cl
add    %eax,%ecx
mov    %ecx,(%r15)
mov    %ebx,%esi
test   $0x1,%al
je     12b5 <hex_xy2nm+0x1a5>
test   %r8d,%r8d
je     12b5 <hex_xy2nm+0x1a5>
mulsd  0xd83(%rip),%xmm3
addsd  %xmm4,%xmm3
mulsd  0xd7f(%rip),%xmm3
addsd  %xmm0,%xmm3
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
ucomisd %xmm3,%xmm0
setbe  %cl
add    %ecx,%eax
mov    %eax,(%r15)
lea    (%rbx,%rcx,2),%esi
add    $0xffffffff,%esi
sar    %esi
mov    %esi,(%r14)
imul   $0x2710,(%r15),%eax
add    %esi,%eax
add    $0x30,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define i32 @hex_xy2nm(double %0, double %1, double %2, i32* %3, i32* %4) {
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca i32*
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
store double %0, double* %6
store double %1, double* %7
store double %2, double* %8
store i32* %3, i32** %9
store i32* %4, i32** %10
store double 0.000000e+00, double* %19
store double 0.000000e+00, double* %20
%21 = call double @sqrt(double 3.000000e+00)
%22 = fdiv double %21, 2.000000e+00
%23 = load double, double* %6
%24 = fmul double %22, %23
store double %24, double* %14
%25 = load double, double* %14
%26 = fmul double 7.500000e-01, %25
store double %26, double* %15
%27 = load double, double* %6
store double %27, double* %16
%28 = load double, double* %8
%29 = load double, double* %20
%30 = fsub double %28, %29
%31 = fmul double %30, 2.000000e+00
%32 = load double, double* %16
%33 = fdiv double %31, %32
%34 = call i64 @floor(double %33)
%35 = trunc i64 %34 to i32
store i32 %35, i32* %12
%36 = load double, double* %7
%37 = load double, double* %19
%38 = fsub double %36, %37
%39 = load double, double* %15
%40 = fdiv double %38, %39
%41 = call i64 @floor(double %40)
%42 = trunc i64 %41 to i32
store i32 %42, i32* %13
%43 = load double, double* %7
%44 = load double, double* %19
%45 = fsub double %43, %44
%46 = load i32, i32* %13
%47 = sitofp i32 %46 to double
%48 = load double, double* %15
%49 = fneg double %47
%50 = call double @llvm.fmuladd.f64(double %49, double %48, double %45)
store double %50, double* %17
%51 = load double, double* %8
%52 = load double, double* %20
%53 = fsub double %51, %52
%54 = load i32, i32* %12
%55 = sitofp i32 %54 to double
%56 = load double, double* %16
%57 = fmul double %55, %56
%58 = fdiv double %57, 2.000000e+00
%59 = fsub double %53, %58
store double %59, double* %18
%60 = load i32, i32* %13
%61 = srem i32 %60, 2
%62 = icmp eq i32 %61, 0
br i1 %62, label %63, label %85
63: ; preds = %5
%64 = load i32, i32* %12
%65 = srem i32 %64, 2
%66 = icmp eq i32 %65, 0
br i1 %66, label %67, label %85
67: ; preds = %63
%68 = load double, double* %18
%69 = call double @sqrt(double 3.000000e+00)
%70 = load double, double* %17
%71 = load double, double* %14
%72 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %71, double %70)
%73 = call double @llvm.fmuladd.f64(double %69, double %72, double %68)
%74 = fcmp olt double %73, 0.000000e+00
br i1 %74, label %75, label %79
75: ; preds = %67
%76 = load i32, i32* %13
%77 = load i32*, i32** %9
store i32 %76, i32* %77
%78 = load i32, i32* %12
store i32 %78, i32* %11
br label %84
79: ; preds = %67
%80 = load i32, i32* %13
%81 = add nsw i32 %80, 1
%82 = load i32*, i32** %9
store i32 %81, i32* %82
%83 = load i32, i32* %12
store i32 %83, i32* %11
br label %84
84: ; preds = %79, %75
br label %85
85: ; preds = %84, %63, %5
%86 = load i32, i32* %13
%87 = srem i32 %86, 2
%88 = icmp eq i32 %87, 0
br i1 %88, label %89, label %114
89: ; preds = %85
%90 = load i32, i32* %12
%91 = srem i32 %90, 2
%92 = icmp ne i32 %91, 0
br i1 %92, label %93, label %114
93: ; preds = %89
%94 = load double, double* %18
%95 = call double @sqrt(double 3.000000e+00)
%96 = load double, double* %17
%97 = load double, double* %14
%98 = call double @llvm.fmuladd.f64(double -2.500000e-01, double %97, double %96)
%99 = fneg double %95
%100 = call double @llvm.fmuladd.f64(double %99, double %98, double %94)
%101 = fcmp ogt double %100, 0.000000e+00
br i1 %101, label %102, label %107
102: ; preds = %93
%103 = load i32, i32* %13
%104 = load i32*, i32** %9
store i32 %103, i32* %104
%105 = load i32, i32* %12
%106 = add nsw i32 %105, 1
store i32 %106, i32* %11
br label %113
107: ; preds = %93
%108 = load i32, i32* %13
%109 = add nsw i32 %108, 1
%110 = load i32*, i32** %9
store i32 %109, i32* %110
%111 = load i32, i32* %12
%112 = sub nsw i32 %111, 1
store i32 %112, i32* %11
br label %113
113: ; preds = %107, %102
br label %114
114: ; preds = %113, %89, %85
%115 = load i32, i32* %13
%116 = srem i32 %115, 2
%117 = icmp ne i32 %116, 0
br i1 %117, label %118, label %141
118: ; preds = %114
%119 = load i32, i32* %12
%120 = srem i32 %119, 2
%121 = icmp eq i32 %120, 0
br i1 %121, label %122, label %141
122: ; preds = %118
%123 = load double, double* %18
%124 = call double @sqrt(double 3.000000e+00)
%125 = load double, double* %17
%126 = load double, double* %14
%127 = call double @llvm.fmuladd.f64(double -2.500000e-01, double %126, double %125)
%128 = fneg double %124
%129 = call double @llvm.fmuladd.f64(double %128, double %127, double %123)
%130 = fcmp ogt double %129, 0.000000e+00
br i1 %130, label %131, label %135
131: ; preds = %122
%132 = load i32, i32* %13
%133 = load i32*, i32** %9
store i32 %132, i32* %133
%134 = load i32, i32* %12
store i32 %134, i32* %11
br label %140
135: ; preds = %122
%136 = load i32, i32* %13
%137 = add nsw i32 %136, 1
%138 = load i32*, i32** %9
store i32 %137, i32* %138
%139 = load i32, i32* %12
store i32 %139, i32* %11
br label %140
140: ; preds = %135, %131
br label %141
141: ; preds = %140, %118, %114
%142 = load i32, i32* %13
%143 = srem i32 %142, 2
%144 = icmp ne i32 %143, 0
br i1 %144, label %145, label %169
145: ; preds = %141
%146 = load i32, i32* %12
%147 = srem i32 %146, 2
%148 = icmp ne i32 %147, 0
br i1 %148, label %149, label %169
149: ; preds = %145
%150 = load double, double* %18
%151 = call double @sqrt(double 3.000000e+00)
%152 = load double, double* %17
%153 = load double, double* %14
%154 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %153, double %152)
%155 = call double @llvm.fmuladd.f64(double %151, double %154, double %150)
%156 = fcmp olt double %155, 0.000000e+00
br i1 %156, label %157, label %162
157: ; preds = %149
%158 = load i32, i32* %13
%159 = load i32*, i32** %9
store i32 %158, i32* %159
%160 = load i32, i32* %12
%161 = sub nsw i32 %160, 1
store i32 %161, i32* %11
br label %168
162: ; preds = %149
%163 = load i32, i32* %13
%164 = add nsw i32 %163, 1
%165 = load i32*, i32** %9
store i32 %164, i32* %165
%166 = load i32, i32* %12
%167 = add nsw i32 %166, 1
store i32 %167, i32* %11
br label %168
168: ; preds = %162, %157
br label %169
169: ; preds = %168, %145, %141
%170 = load i32, i32* %11
%171 = ashr i32 %170, 1
%172 = load i32*, i32** %10
store i32 %171, i32* %172
%173 = load i32*, i32** %9
%174 = load i32, i32* %173
%175 = mul nsw i32 %174, 10000
%176 = load i32*, i32** %10
%177 = load i32, i32* %176
%178 = add nsw i32 %175, %177
ret i32 %178
}

LLVM IR (-O2):
define i32 @hex_xy2nm(double %0, double %1, double %2, i32* nocapture %3, i32* nocapture writeonly %4) {
%6 = fmul double %0, 0x3FEBB67AE8584CAA
%7 = fmul double %6, 7.500000e-01
%8 = fmul double %2, 2.000000e+00
%9 = fdiv double %8, %0
%10 = tail call i64 @floor(double %9)
%11 = trunc i64 %10 to i32
%12 = fdiv double %1, %7
%13 = tail call i64 @floor(double %12)
%14 = trunc i64 %13 to i32
%15 = sitofp i32 %14 to double
%16 = fneg double %15
%17 = tail call double @llvm.fmuladd.f64(double %16, double %7, double %1)
%18 = sitofp i32 %11 to double
%19 = fmul double %18, %0
%20 = fmul double %19, 5.000000e-01
%21 = fsub double %2, %20
%22 = and i32 %14, 1
%23 = icmp eq i32 %22, 0
%24 = and i32 %11, 1
%25 = icmp eq i32 %24, 0
%26 = select i1 %23, i1 %25, i1 false
br i1 %26, label %27, label %33
27: ; preds = %5
%28 = tail call double @llvm.fmuladd.f64(double %6, double -5.000000e-01, double %17)
%29 = tail call double @llvm.fmuladd.f64(double %28, double 0x3FFBB67AE8584CAA, double %21)
%30 = fcmp uge double %29, 0.000000e+00
%31 = zext i1 %30 to i32
%32 = add nuw nsw i32 %14, %31
store i32 %32, i32* %3
br label %33
33: ; preds = %27, %5
%34 = xor i1 %23, true
%35 = select i1 %34, i1 true, i1 %25
br i1 %35, label %45, label %36
36: ; preds = %33
%37 = tail call double @llvm.fmuladd.f64(double %6, double -2.500000e-01, double %17)
%38 = tail call double @llvm.fmuladd.f64(double %37, double 0xBFFBB67AE8584CAA, double %21)
%39 = fcmp ogt double %38, 0.000000e+00
%40 = xor i1 %39, true
%41 = zext i1 %40 to i32
%42 = add nsw i32 %14, %41
%43 = select i1 %39, i32 1, i32 -1
store i32 %42, i32* %3
%44 = add nsw i32 %43, %11
br label %45
45: ; preds = %36, %33
%46 = phi i32 [ %11, %33 ], [ %44, %36 ]
%47 = select i1 %34, i1 %25, i1 false
br i1 %47, label %48, label %54
48: ; preds = %45
%49 = tail call double @llvm.fmuladd.f64(double %6, double -2.500000e-01, double %17)
%50 = tail call double @llvm.fmuladd.f64(double %49, double 0xBFFBB67AE8584CAA, double %21)
%51 = fcmp ule double %50, 0.000000e+00
%52 = zext i1 %51 to i32
%53 = add nsw i32 %14, %52
store i32 %53, i32* %3
br label %54
54: ; preds = %48, %45
%55 = phi i32 [ %46, %45 ], [ %11, %48 ]
%56 = select i1 %23, i1 true, i1 %25
br i1 %56, label %66, label %57
57: ; preds = %54
%58 = tail call double @llvm.fmuladd.f64(double %6, double -5.000000e-01, double %17)
%59 = tail call double @llvm.fmuladd.f64(double %58, double 0x3FFBB67AE8584CAA, double %21)
%60 = fcmp olt double %59, 0.000000e+00
%61 = xor i1 %60, true
%62 = zext i1 %61 to i32
%63 = add nsw i32 %14, %62
%64 = select i1 %60, i32 -1, i32 1
store i32 %63, i32* %3
%65 = add nsw i32 %64, %11
br label %66
66: ; preds = %57, %54
%67 = phi i32 [ %55, %54 ], [ %65, %57 ]
%68 = ashr i32 %67, 1
store i32 %68, i32* %4
%69 = load i32, i32* %3
%70 = mul nsw i32 %69, 10000
%71 = add nsw i32 %70, %68
ret i32 %71
}