#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 */
int nBinDigits(long);
void getPos(long i, double p[3]) {
  int n = nBinDigits(i);
  double h = 0.25;
  p[0] = p[1] = p[2] = 0.5;
  for (int j = n - 2 * 3; j >= 0; j -= 3) {
    for (int k = 0; k < 3; k++) {
      int l = j + k;
      p[k] += h * (double)(2 * ((i & (1 << l)) >> l) - 1);
    }
    h *= 0.5;
  }
  return;
}


x86 assembly (clang, -O0):
<getPos>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x8(%rbp),%rdi
call   1030 <nBinDigits@plt>
mov    %eax,-0x14(%rbp)
movsd  0xed4(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movsd  0xebb(%rip),%xmm0
movsd  %xmm0,0x10(%rax)
mov    -0x10(%rbp),%rax
movsd  0xeaa(%rip),%xmm0
movsd  %xmm0,0x8(%rax)
mov    -0x10(%rbp),%rax
movsd  0xe99(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x14(%rbp),%eax
sub    $0x6,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     120e <getPos+0xfe>
movl   $0x0,-0x28(%rbp)
cmpl   $0x3,-0x28(%rbp)
jge    11ee <getPos+0xde>
mov    -0x24(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %edx,%ecx
movslq %ecx,%rcx
and    %rcx,%rax
mov    -0x2c(%rbp),%ecx
sar    %cl,%rax
shl    $0x1,%rax
sub    $0x1,%rax
cvtsi2sd %rax,%xmm2
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1185 <getPos+0x75>
movsd  0xe0a(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x3,%eax
mov    %eax,-0x24(%rbp)
jmp    1174 <getPos+0x64>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<getPos>:
push   %r15
push   %r14
push   %rbx
mov    %rsi,%r14
mov    %rdi,%r15
call   1030 <nBinDigits@plt>
movabs $0x3fe0000000000000,%rcx
mov    %rcx,0x10(%r14)
movapd 0xeca(%rip),%xmm0
movupd %xmm0,(%r14)
cmp    $0x6,%eax
jl     1222 <getPos+0x112>
lea    -0x6(%rax),%edx
lea    -0x4(%rax),%esi
add    $0xfffffffb,%eax
movsd  0xebb(%rip),%xmm0
movsd  0xebb(%rip),%xmm1
movapd %xmm1,%xmm2
movapd %xmm1,%xmm3
movapd %xmm1,%xmm4
nopl   0x0(%rax)
mov    $0x1,%edi
mov    %edx,%ecx
shl    %cl,%edi
movslq %edi,%rdi
and    %r15,%rdi
sar    %cl,%rdi
lea    (%rdi,%rdi,1),%rcx
add    $0xffffffffffffffff,%rcx
xorps  %xmm5,%xmm5
cvtsi2sd %rcx,%xmm5
lea    0x1(%rdx),%ecx
mov    $0x1,%edi
shl    %cl,%edi
movslq %edi,%rdi
and    %r15,%rdi
mov    %eax,%ecx
sar    %cl,%rdi
mulsd  %xmm0,%xmm5
lea    (%rdi,%rdi,1),%rcx
add    $0xffffffffffffffff,%rcx
xorps  %xmm6,%xmm6
cvtsi2sd %rcx,%xmm6
addsd  %xmm5,%xmm4
mulsd  %xmm0,%xmm6
lea    0x2(%rdx),%rdi
mov    $0x1,%ebx
mov    %edi,%ecx
shl    %cl,%ebx
movslq %ebx,%rbx
and    %r15,%rbx
mov    %esi,%ecx
sar    %cl,%rbx
addsd  %xmm6,%xmm3
lea    (%rbx,%rbx,1),%rcx
add    $0xffffffffffffffff,%rcx
xorps  %xmm5,%xmm5
cvtsi2sd %rcx,%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm5,%xmm2
mulsd  %xmm1,%xmm0
add    $0xfffffffd,%esi
add    $0xfffffffd,%eax
add    $0xfffffffe,%edi
add    $0xfffffffffffffffd,%rdx
cmp    $0x2,%edi
ja     1170 <getPos+0x60>
movsd  %xmm4,(%r14)
movsd  %xmm3,0x8(%r14)
movsd  %xmm2,0x10(%r14)
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @getPos(i64 %0, double* %1) {
%3 = alloca i64
%4 = alloca double*
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store i64 %0, i64* %3
store double* %1, double** %4
%10 = load i64, i64* %3
%11 = call i32 @nBinDigits(i64 %10)
store i32 %11, i32* %5
store double 2.500000e-01, double* %6
%12 = load double*, double** %4
%13 = getelementptr inbounds double, double* %12, i64 2
store double 5.000000e-01, double* %13
%14 = load double*, double** %4
%15 = getelementptr inbounds double, double* %14, i64 1
store double 5.000000e-01, double* %15
%16 = load double*, double** %4
%17 = getelementptr inbounds double, double* %16, i64 0
store double 5.000000e-01, double* %17
%18 = load i32, i32* %5
%19 = sub nsw i32 %18, 6
store i32 %19, i32* %7
br label %20
20: ; preds = %55, %2
%21 = load i32, i32* %7
%22 = icmp sge i32 %21, 0
br i1 %22, label %23, label %58
23: ; preds = %20
store i32 0, i32* %8
br label %24
24: ; preds = %49, %23
%25 = load i32, i32* %8
%26 = icmp slt i32 %25, 3
br i1 %26, label %27, label %52
27: ; preds = %24
%28 = load i32, i32* %7
%29 = load i32, i32* %8
%30 = add nsw i32 %28, %29
store i32 %30, i32* %9
%31 = load double, double* %6
%32 = load i64, i64* %3
%33 = load i32, i32* %9
%34 = shl i32 1, %33
%35 = sext i32 %34 to i64
%36 = and i64 %32, %35
%37 = load i32, i32* %9
%38 = zext i32 %37 to i64
%39 = ashr i64 %36, %38
%40 = mul nsw i64 2, %39
%41 = sub nsw i64 %40, 1
%42 = sitofp i64 %41 to double
%43 = load double*, double** %4
%44 = load i32, i32* %8
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %43, i64 %45
%47 = load double, double* %46
%48 = call double @llvm.fmuladd.f64(double %31, double %42, double %47)
store double %48, double* %46
br label %49
49: ; preds = %27
%50 = load i32, i32* %8
%51 = add nsw i32 %50, 1
store i32 %51, i32* %8
br label %24
52: ; preds = %24
%53 = load double, double* %6
%54 = fmul double %53, 5.000000e-01
store double %54, double* %6
br label %55
55: ; preds = %52
%56 = load i32, i32* %7
%57 = sub nsw i32 %56, 3
store i32 %57, i32* %7
br label %20
58: ; preds = %20
ret void
}

LLVM IR (-O2):
define void @getPos(i64 %0, double* nocapture %1) {
%3 = tail call i32 @nBinDigits(i64 %0)
%4 = getelementptr inbounds double, double* %1, i64 2
store double 5.000000e-01, double* %4
%5 = getelementptr inbounds double, double* %1, i64 1
%6 = bitcast double* %1 to <2 x double>*
store <2 x double> <double 5.000000e-01, double 5.000000e-01>, <2 x double>* %6
%7 = icmp sgt i32 %3, 5
br i1 %7, label %8, label %49
8: ; preds = %2
%9 = add nsw i32 %3, -6
br label %10
10: ; preds = %8, %10
%11 = phi double [ %44, %10 ], [ 5.000000e-01, %8 ]
%12 = phi double [ %34, %10 ], [ 5.000000e-01, %8 ]
%13 = phi double [ %24, %10 ], [ 5.000000e-01, %8 ]
%14 = phi double [ %45, %10 ], [ 2.500000e-01, %8 ]
%15 = phi i32 [ %46, %10 ], [ %9, %8 ]
%16 = shl nuw i32 1, %15
%17 = sext i32 %16 to i64
%18 = and i64 %17, %0
%19 = zext i32 %15 to i64
%20 = ashr i64 %18, %19
%21 = shl nsw i64 %20, 1
%22 = add nsw i64 %21, -1
%23 = sitofp i64 %22 to double
%24 = tail call double @llvm.fmuladd.f64(double %14, double %23, double %13)
%25 = add nuw nsw i32 %15, 1
%26 = shl nuw i32 1, %25
%27 = sext i32 %26 to i64
%28 = and i64 %27, %0
%29 = zext i32 %25 to i64
%30 = ashr i64 %28, %29
%31 = shl nsw i64 %30, 1
%32 = add nsw i64 %31, -1
%33 = sitofp i64 %32 to double
%34 = tail call double @llvm.fmuladd.f64(double %14, double %33, double %12)
%35 = add nuw nsw i32 %15, 2
%36 = shl nuw i32 1, %35
%37 = sext i32 %36 to i64
%38 = and i64 %37, %0
%39 = zext i32 %35 to i64
%40 = ashr i64 %38, %39
%41 = shl nsw i64 %40, 1
%42 = add nsw i64 %41, -1
%43 = sitofp i64 %42 to double
%44 = tail call double @llvm.fmuladd.f64(double %14, double %43, double %11)
%45 = fmul double %14, 5.000000e-01
%46 = add nsw i32 %15, -3
%47 = icmp ugt i32 %15, 2
br i1 %47, label %10, label %48
48: ; preds = %10
store double %24, double* %1
store double %34, double* %5
store double %44, double* %4
br label %49
49: ; preds = %48, %2
ret void
}