#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 fabs(double);
scalar_t__ log(double);
int BinomialK(double alpha, int n, double C[], double S[]) {
  int i, nround = n, alphaint = (int)alpha;
  double c = 1, large = 1E200;
  if (alpha > 0 && fabs(alpha - alphaint) < 1e-100) {
    nround = ((n) < (alphaint / 2) ? (n) : (alphaint / 2));
  }
  C[0] = 1;
  S[0] = 0;
  for (i = 1; i <= nround; i++) {
    c *= (alpha - i + 1) / i;
    S[i] = S[i - 1];
    if (c > large) {
      S[i] += log(c);
      c = 1;
    }
    C[i] = c;
  }
  for (; i <= ((n) < (alphaint) ? (n) : (alphaint)); i++) {
    C[i] = C[alphaint - i];
    S[i] = S[alphaint - i];
  }
  for (; i <= n; i++) {
    C[i] = 0;
    S[i] = 0;
  }
  return (0);
}


x86 assembly (clang, -O0):
<BinomialK>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x28(%rbp)
cvttsd2si -0x8(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movsd  0xeba(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xea5(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    11d8 <BinomialK+0xb8>
movsd  -0x8(%rbp),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm1
subsd  %xmm1,%xmm0
call   1030 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xe7f(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    11d8 <BinomialK+0xb8>
mov    -0xc(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x2c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x44(%rbp),%eax
cmp    %ecx,%eax
jge    11c4 <BinomialK+0xa4>
mov    -0xc(%rbp),%eax
mov    %eax,-0x48(%rbp)
jmp    11d2 <BinomialK+0xb2>
mov    -0x2c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movsd  0xe24(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     12b0 <BinomialK+0x190>
movsd  -0x8(%rbp),%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  0xdec(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1290 <BinomialK+0x170>
movsd  -0x38(%rbp),%xmm0
call   1040 <log@plt>
cvtsi2sd %rax,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  0xd7d(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11fa <BinomialK+0xda>
jmp    12b5 <BinomialK+0x195>
mov    -0x24(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    12d2 <BinomialK+0x1b2>
mov    -0xc(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    12d8 <BinomialK+0x1b8>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x50(%rbp),%ecx
cmp    %ecx,%eax
jg     1332 <BinomialK+0x212>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12b5 <BinomialK+0x195>
jmp    1337 <BinomialK+0x217>
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1371 <BinomialK+0x251>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1337 <BinomialK+0x217>
xor    %eax,%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<BinomialK>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
movapd %xmm0,%xmm4
cvttsd2si %xmm0,%ebx
mov    %rdx,%rbp
mov    %rsi,%r13
mov    %edi,%r14d
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm4
mov    %edi,%r15d
movsd  %xmm4,0x10(%rsp)
jbe    1184 <BinomialK+0x64>
cvtsi2sd %ebx,%xmm1
subsd  %xmm1,%xmm4
movapd %xmm4,%xmm0
call   1030 <fabs@plt>
movsd  0x10(%rsp),%xmm4
mov    %r14d,%r15d
test   %eax,%eax
jg     1184 <BinomialK+0x64>
mov    %ebx,%r15d
shr    $0x1f,%r15d
add    %ebx,%r15d
sar    %r15d
cmp    %r14d,%r15d
cmovg  %r14d,%r15d
movabs $0x3ff0000000000000,%rax
mov    %rax,0x0(%r13)
movq   $0x0,0x0(%rbp)
test   %r15d,%r15d
jle    1266 <BinomialK+0x146>
mov    %ebx,0xc(%rsp)
add    $0x1,%r15d
mov    $0x2,%ebx
sub    %r15,%rbx
movsd  0xe45(%rip),%xmm5
mov    $0x1,%r12d
xorpd  %xmm1,%xmm1
movsd  0xe3b(%rip),%xmm6
movapd %xmm5,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm2,%xmm2
cvtsi2sd %r12d,%xmm2
movapd %xmm0,%xmm3
movapd %xmm4,%xmm0
subsd  %xmm2,%xmm0
addsd  %xmm5,%xmm0
divsd  %xmm2,%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm1,0x0(%rbp,%r12,8)
ucomisd %xmm6,%xmm0
jbe    1242 <BinomialK+0x122>
call   1040 <log@plt>
movsd  0xdee(%rip),%xmm6
movsd  0xdde(%rip),%xmm5
movsd  0x10(%rsp),%xmm4
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x0(%rbp,%r12,8),%xmm0
movsd  %xmm0,0x0(%rbp,%r12,8)
movapd %xmm5,%xmm0
movsd  %xmm0,0x0(%r13,%r12,8)
lea    (%rbx,%r12,1),%rax
cmp    $0x1,%rax
je     1260 <BinomialK+0x140>
movsd  0x0(%rbp,%r12,8),%xmm1
add    $0x1,%r12
jmp    11e0 <BinomialK+0xc0>
mov    0xc(%rsp),%ebx
jmp    126c <BinomialK+0x14c>
mov    $0x1,%r15d
cmp    %r14d,%ebx
mov    %ebx,%esi
cmovg  %r14d,%esi
cmp    %r15d,%esi
jl     131c <BinomialK+0x1fc>
movslq %r15d,%rax
movslq %ebx,%rdx
lea    0x1(%rsi),%ecx
mov    %ecx,%edi
sub    %r15d,%edi
test   $0x1,%dil
je     12b9 <BinomialK+0x199>
mov    %rdx,%rdi
sub    %rax,%rdi
movsd  0x0(%r13,%rdi,8),%xmm0
mov    %r15d,%ebx
movsd  %xmm0,0x0(%r13,%rbx,8)
movsd  0x0(%rbp,%rdi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rbx,8)
add    $0x1,%rax
cmp    %r15d,%esi
je     1319 <BinomialK+0x1f9>
mov    %ecx,%esi
sub    %rax,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x0(%r13,%rdx,8),%xmm0
mov    %eax,%edi
shl    $0x3,%rdi
movsd  %xmm0,0x0(%r13,%rdi,1)
movsd  0x0(%rbp,%rdx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdi,1)
lea    0x1(%rax),%edi
movsd  -0x8(%r13,%rdx,8),%xmm0
movsd  %xmm0,0x0(%r13,%rdi,8)
movsd  -0x8(%rbp,%rdx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdi,8)
add    $0x2,%rax
add    $0xfffffffffffffffe,%rdx
cmp    %eax,%esi
jne    12d0 <BinomialK+0x1b0>
mov    %ecx,%r15d
cmp    %r14d,%r15d
jg     14cc <BinomialK+0x3ac>
mov    %r15d,%r9d
mov    %r14d,%ecx
sub    %r15d,%ecx
cmp    $0x3,%ecx
jb     1487 <BinomialK+0x367>
lea    0x0(,%r9,8),%rdx
add    %r13,%rdx
lea    (%rcx,%r9,1),%rax
lea    0x8(,%rax,8),%rsi
add    %rbp,%rsi
cmp    %rsi,%rdx
jae    1375 <BinomialK+0x255>
lea    0x8(,%rax,8),%rax
add    %r13,%rax
lea    0x0(,%r9,8),%rdx
add    %rbp,%rdx
cmp    %rax,%rdx
jb     1487 <BinomialK+0x367>
add    $0x1,%rcx
mov    %rcx,%r8
and    $0xfffffffffffffffc,%r8
lea    -0x4(%r8),%rax
mov    %rax,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%esi
and    $0x3,%esi
cmp    $0xc,%rax
jae    13a1 <BinomialK+0x281>
xor    %ebx,%ebx
jmp    1438 <BinomialK+0x318>
and    $0xfffffffffffffffc,%rdi
lea    0x70(,%r9,8),%rdx
add    %rbp,%rdx
lea    0x70(,%r9,8),%rax
add    %r13,%rax
xor    %ebx,%ebx
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movupd %xmm0,-0x70(%rax,%rbx,8)
movupd %xmm0,-0x60(%rax,%rbx,8)
movupd %xmm0,-0x70(%rdx,%rbx,8)
movupd %xmm0,-0x60(%rdx,%rbx,8)
movupd %xmm0,-0x50(%rax,%rbx,8)
movupd %xmm0,-0x40(%rax,%rbx,8)
movupd %xmm0,-0x50(%rdx,%rbx,8)
movupd %xmm0,-0x40(%rdx,%rbx,8)
movupd %xmm0,-0x30(%rax,%rbx,8)
movupd %xmm0,-0x20(%rax,%rbx,8)
movupd %xmm0,-0x30(%rdx,%rbx,8)
movupd %xmm0,-0x20(%rdx,%rbx,8)
movupd %xmm0,-0x10(%rax,%rbx,8)
movupd %xmm0,(%rax,%rbx,8)
movupd %xmm0,-0x10(%rdx,%rbx,8)
movupd %xmm0,(%rdx,%rbx,8)
add    $0x10,%rbx
add    $0xfffffffffffffffc,%rdi
jne    13d0 <BinomialK+0x2b0>
test   %rsi,%rsi
je     147f <BinomialK+0x35f>
add    %r9,%rbx
lea    0x10(,%rbx,8),%rax
add    %rbp,%rax
lea    0x10(,%rbx,8),%rdx
add    %r13,%rdx
shl    $0x5,%rsi
xor    %edi,%edi
xorpd  %xmm0,%xmm0
movupd %xmm0,-0x10(%rdx,%rdi,1)
movupd %xmm0,(%rdx,%rdi,1)
movupd %xmm0,-0x10(%rax,%rdi,1)
movupd %xmm0,(%rax,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rsi
jne    1460 <BinomialK+0x340>
cmp    %r8,%rcx
je     14cc <BinomialK+0x3ac>
add    %r8,%r9
lea    -0x1(%r9),%eax
lea    0x0(,%r9,8),%rcx
add    %rbp,%rcx
lea    0x0(,%r9,8),%rdx
add    %r13,%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movq   $0x0,(%rdx,%rsi,8)
movq   $0x0,(%rcx,%rsi,8)
add    $0x1,%rsi
lea    (%rax,%rsi,1),%edi
cmp    %r14d,%edi
jl     14b0 <BinomialK+0x390>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @BinomialK(double %0, i32 %1, double* %2, double* %3) {
%5 = alloca double
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca double
store double %0, double* %5
store i32 %1, i32* %6
store double* %2, double** %7
store double* %3, double** %8
%14 = load i32, i32* %6
store i32 %14, i32* %10
%15 = load double, double* %5
%16 = fptosi double %15 to i32
store i32 %16, i32* %11
store double 1.000000e+00, double* %12
store double 0x6974E718D7D7625A, double* %13
%17 = load double, double* %5
%18 = fcmp ogt double %17, 0.000000e+00
br i1 %18, label %19, label %39
19: ; preds = %4
%20 = load double, double* %5
%21 = load i32, i32* %11
%22 = sitofp i32 %21 to double
%23 = fsub double %20, %22
%24 = call i32 @fabs(double %23)
%25 = sitofp i32 %24 to double
%26 = fcmp olt double %25, 1.000000e-100
br i1 %26, label %27, label %39
27: ; preds = %19
%28 = load i32, i32* %6
%29 = load i32, i32* %11
%30 = sdiv i32 %29, 2
%31 = icmp slt i32 %28, %30
br i1 %31, label %32, label %34
32: ; preds = %27
%33 = load i32, i32* %6
br label %37
34: ; preds = %27
%35 = load i32, i32* %11
%36 = sdiv i32 %35, 2
br label %37
37: ; preds = %34, %32
%38 = phi i32 [ %33, %32 ], [ %36, %34 ]
store i32 %38, i32* %10
br label %39
39: ; preds = %37, %19, %4
%40 = load double*, double** %7
%41 = getelementptr inbounds double, double* %40, i64 0
store double 1.000000e+00, double* %41
%42 = load double*, double** %8
%43 = getelementptr inbounds double, double* %42, i64 0
store double 0.000000e+00, double* %43
store i32 1, i32* %9
br label %44
44: ; preds = %88, %39
%45 = load i32, i32* %9
%46 = load i32, i32* %10
%47 = icmp sle i32 %45, %46
br i1 %47, label %48, label %91
48: ; preds = %44
%49 = load double, double* %5
%50 = load i32, i32* %9
%51 = sitofp i32 %50 to double
%52 = fsub double %49, %51
%53 = fadd double %52, 1.000000e+00
%54 = load i32, i32* %9
%55 = sitofp i32 %54 to double
%56 = fdiv double %53, %55
%57 = load double, double* %12
%58 = fmul double %57, %56
store double %58, double* %12
%59 = load double*, double** %8
%60 = load i32, i32* %9
%61 = sub nsw i32 %60, 1
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %59, i64 %62
%64 = load double, double* %63
%65 = load double*, double** %8
%66 = load i32, i32* %9
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %65, i64 %67
store double %64, double* %68
%69 = load double, double* %12
%70 = load double, double* %13
%71 = fcmp ogt double %69, %70
br i1 %71, label %72, label %82
72: ; preds = %48
%73 = load double, double* %12
%74 = call i64 @log(double %73)
%75 = sitofp i64 %74 to double
%76 = load double*, double** %8
%77 = load i32, i32* %9
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %76, i64 %78
%80 = load double, double* %79
%81 = fadd double %80, %75
store double %81, double* %79
store double 1.000000e+00, double* %12
br label %82
82: ; preds = %72, %48
%83 = load double, double* %12
%84 = load double*, double** %7
%85 = load i32, i32* %9
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
store double %83, double* %87
br label %88
88: ; preds = %82
%89 = load i32, i32* %9
%90 = add nsw i32 %89, 1
store i32 %90, i32* %9
br label %44
91: ; preds = %44
br label %92
92: ; preds = %127, %91
%93 = load i32, i32* %9
%94 = load i32, i32* %6
%95 = load i32, i32* %11
%96 = icmp slt i32 %94, %95
br i1 %96, label %97, label %99
97: ; preds = %92
%98 = load i32, i32* %6
br label %101
99: ; preds = %92
%100 = load i32, i32* %11
br label %101
101: ; preds = %99, %97
%102 = phi i32 [ %98, %97 ], [ %100, %99 ]
%103 = icmp sle i32 %93, %102
br i1 %103, label %104, label %130
104: ; preds = %101
%105 = load double*, double** %7
%106 = load i32, i32* %11
%107 = load i32, i32* %9
%108 = sub nsw i32 %106, %107
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %105, i64 %109
%111 = load double, double* %110
%112 = load double*, double** %7
%113 = load i32, i32* %9
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %112, i64 %114
store double %111, double* %115
%116 = load double*, double** %8
%117 = load i32, i32* %11
%118 = load i32, i32* %9
%119 = sub nsw i32 %117, %118
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %116, i64 %120
%122 = load double, double* %121
%123 = load double*, double** %8
%124 = load i32, i32* %9
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double, double* %123, i64 %125
store double %122, double* %126
br label %127
127: ; preds = %104
%128 = load i32, i32* %9
%129 = add nsw i32 %128, 1
store i32 %129, i32* %9
br label %92
130: ; preds = %101
br label %131
131: ; preds = %144, %130
%132 = load i32, i32* %9
%133 = load i32, i32* %6
%134 = icmp sle i32 %132, %133
br i1 %134, label %135, label %147
135: ; preds = %131
%136 = load double*, double** %7
%137 = load i32, i32* %9
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds double, double* %136, i64 %138
store double 0.000000e+00, double* %139
%140 = load double*, double** %8
%141 = load i32, i32* %9
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %140, i64 %142
store double 0.000000e+00, double* %143
br label %144
144: ; preds = %135
%145 = load i32, i32* %9
%146 = add nsw i32 %145, 1
store i32 %146, i32* %9
br label %131
147: ; preds = %131
ret i32 0
}

LLVM IR (-O2):
define i32 @BinomialK(double %0, i32 %1, double* nocapture %2, double* nocapture %3) {
%5 = fptosi double %0 to i32
%6 = fcmp ogt double %0, 0.000000e+00
br i1 %6, label %7, label %16
7: ; preds = %4
%8 = sitofp i32 %5 to double
%9 = fsub double %0, %8
%10 = tail call i32 @fabs(double %9)
%11 = icmp slt i32 %10, 1
br i1 %11, label %12, label %16
12: ; preds = %7
%13 = sdiv i32 %5, 2
%14 = icmp sgt i32 %13, %1
%15 = select i1 %14, i32 %1, i32 %13
br label %16
16: ; preds = %12, %7, %4
%17 = phi i32 [ %15, %12 ], [ %1, %7 ], [ %1, %4 ]
store double 1.000000e+00, double* %2
store double 0.000000e+00, double* %3
%18 = icmp slt i32 %17, 1
br i1 %18, label %22, label %19
19: ; preds = %16
%20 = add nuw i32 %17, 1
%21 = zext i32 %20 to i64
br label %47
22: ; preds = %64, %16
%23 = phi i32 [ 1, %16 ], [ %20, %64 ]
%24 = icmp sgt i32 %5, %1
%25 = select i1 %24, i32 %1, i32 %5
%26 = icmp slt i32 %25, %23
br i1 %26, label %71, label %27
27: ; preds = %22
%28 = sext i32 %23 to i64
%29 = sext i32 %5 to i64
%30 = add i32 %25, 1
%31 = sub i32 %30, %23
%32 = and i32 %31, 1
%33 = icmp eq i32 %32, 0
br i1 %33, label %44, label %34
34: ; preds = %27
%35 = sub nsw i64 %29, %28
%36 = getelementptr inbounds double, double* %2, i64 %35
%37 = load double, double* %36
%38 = and i64 %28, 4294967295
%39 = getelementptr inbounds double, double* %2, i64 %38
store double %37, double* %39
%40 = getelementptr inbounds double, double* %3, i64 %35
%41 = load double, double* %40
%42 = getelementptr inbounds double, double* %3, i64 %38
store double %41, double* %42
%43 = add nsw i64 %28, 1
br label %44
44: ; preds = %34, %27
%45 = phi i64 [ %28, %27 ], [ %43, %34 ]
%46 = icmp eq i32 %25, %23
br i1 %46, label %71, label %169
47: ; preds = %69, %19
%48 = phi double [ 0.000000e+00, %19 ], [ %70, %69 ]
%49 = phi i64 [ 1, %19 ], [ %67, %69 ]
%50 = phi double [ 1.000000e+00, %19 ], [ %65, %69 ]
%51 = trunc i64 %49 to i32
%52 = sitofp i32 %51 to double
%53 = fsub double %0, %52
%54 = fadd double %53, 1.000000e+00
%55 = fdiv double %54, %52
%56 = fmul double %50, %55
%57 = getelementptr inbounds double, double* %3, i64 %49
store double %48, double* %57
%58 = fcmp ogt double %56, 0x6974E718D7D7625A
br i1 %58, label %59, label %64
59: ; preds = %47
%60 = tail call i64 @log(double %56)
%61 = sitofp i64 %60 to double
%62 = load double, double* %57
%63 = fadd double %62, %61
store double %63, double* %57
br label %64
64: ; preds = %59, %47
%65 = phi double [ 1.000000e+00, %59 ], [ %56, %47 ]
%66 = getelementptr inbounds double, double* %2, i64 %49
store double %65, double* %66
%67 = add nuw nsw i64 %49, 1
%68 = icmp eq i64 %67, %21
br i1 %68, label %22, label %69
69: ; preds = %64
%70 = load double, double* %57
br label %47
71: ; preds = %44, %169, %22
%72 = phi i32 [ %23, %22 ], [ %30, %169 ], [ %30, %44 ]
%73 = icmp sgt i32 %72, %1
br i1 %73, label %198, label %74
74: ; preds = %71
%75 = zext i32 %72 to i64
%76 = sub i32 %1, %72
%77 = zext i32 %76 to i64
%78 = add nuw nsw i64 %77, 1
%79 = icmp ult i32 %76, 3
br i1 %79, label %167, label %80
80: ; preds = %74
%81 = getelementptr double, double* %2, i64 %75
%82 = sub i32 %1, %72
%83 = zext i32 %82 to i64
%84 = add nuw nsw i64 %75, %83
%85 = add nuw nsw i64 %84, 1
%86 = getelementptr double, double* %2, i64 %85
%87 = getelementptr double, double* %3, i64 %75
%88 = getelementptr double, double* %3, i64 %85
%89 = icmp ult double* %81, %88
%90 = icmp ult double* %87, %86
%91 = and i1 %89, %90
br i1 %91, label %167, label %92
92: ; preds = %80
%93 = and i64 %78, 8589934588
%94 = add nuw nsw i64 %93, %75
%95 = add nsw i64 %93, -4
%96 = lshr exact i64 %95, 2
%97 = add nuw nsw i64 %96, 1
%98 = and i64 %97, 3
%99 = icmp ult i64 %95, 12
br i1 %99, label %147, label %100
100: ; preds = %92
%101 = and i64 %97, 9223372036854775804
br label %102
102: ; preds = %102, %100
%103 = phi i64 [ 0, %100 ], [ %144, %102 ]
%104 = phi i64 [ 0, %100 ], [ %145, %102 ]
%105 = add i64 %103, %75
%106 = getelementptr inbounds double, double* %2, i64 %105
%107 = bitcast double* %106 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %107
%108 = getelementptr inbounds double, double* %106, i64 2
%109 = bitcast double* %108 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %109
%110 = getelementptr inbounds double, double* %3, i64 %105
%111 = bitcast double* %110 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %111
%112 = getelementptr inbounds double, double* %110, i64 2
%113 = bitcast double* %112 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %113
%114 = or i64 %103, 4
%115 = add i64 %114, %75
%116 = getelementptr inbounds double, double* %2, i64 %115
%117 = bitcast double* %116 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %117
%118 = getelementptr inbounds double, double* %116, i64 2
%119 = bitcast double* %118 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %119
%120 = getelementptr inbounds double, double* %3, i64 %115
%121 = bitcast double* %120 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %121
%122 = getelementptr inbounds double, double* %120, i64 2
%123 = bitcast double* %122 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %123
%124 = or i64 %103, 8
%125 = add i64 %124, %75
%126 = getelementptr inbounds double, double* %2, i64 %125
%127 = bitcast double* %126 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %127
%128 = getelementptr inbounds double, double* %126, i64 2
%129 = bitcast double* %128 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %129
%130 = getelementptr inbounds double, double* %3, i64 %125
%131 = bitcast double* %130 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %131
%132 = getelementptr inbounds double, double* %130, i64 2
%133 = bitcast double* %132 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %133
%134 = or i64 %103, 12
%135 = add i64 %134, %75
%136 = getelementptr inbounds double, double* %2, i64 %135
%137 = bitcast double* %136 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %137
%138 = getelementptr inbounds double, double* %136, i64 2
%139 = bitcast double* %138 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %139
%140 = getelementptr inbounds double, double* %3, i64 %135
%141 = bitcast double* %140 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %141
%142 = getelementptr inbounds double, double* %140, i64 2
%143 = bitcast double* %142 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %143
%144 = add nuw i64 %103, 16
%145 = add nuw i64 %104, 4
%146 = icmp eq i64 %145, %101
br i1 %146, label %147, label %102
147: ; preds = %102, %92
%148 = phi i64 [ 0, %92 ], [ %144, %102 ]
%149 = icmp eq i64 %98, 0
br i1 %149, label %165, label %150
150: ; preds = %147, %150
%151 = phi i64 [ %162, %150 ], [ %148, %147 ]
%152 = phi i64 [ %163, %150 ], [ 0, %147 ]
%153 = add i64 %151, %75
%154 = getelementptr inbounds double, double* %2, i64 %153
%155 = bitcast double* %154 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %155
%156 = getelementptr inbounds double, double* %154, i64 2
%157 = bitcast double* %156 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %157
%158 = getelementptr inbounds double, double* %3, i64 %153
%159 = bitcast double* %158 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %159
%160 = getelementptr inbounds double, double* %158, i64 2
%161 = bitcast double* %160 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %161
%162 = add nuw i64 %151, 4
%163 = add i64 %152, 1
%164 = icmp eq i64 %163, %98
br i1 %164, label %165, label %150
165: ; preds = %150, %147
%166 = icmp eq i64 %78, %93
br i1 %166, label %198, label %167
167: ; preds = %80, %74, %165
%168 = phi i64 [ %75, %80 ], [ %75, %74 ], [ %94, %165 ]
br label %191
169: ; preds = %44, %169
%170 = phi i64 [ %188, %169 ], [ %45, %44 ]
%171 = sub nsw i64 %29, %170
%172 = getelementptr inbounds double, double* %2, i64 %171
%173 = load double, double* %172
%174 = and i64 %170, 4294967295
%175 = getelementptr inbounds double, double* %2, i64 %174
store double %173, double* %175
%176 = getelementptr inbounds double, double* %3, i64 %171
%177 = load double, double* %176
%178 = getelementptr inbounds double, double* %3, i64 %174
store double %177, double* %178
%179 = add nsw i64 %170, 1
%180 = sub nsw i64 %29, %179
%181 = getelementptr inbounds double, double* %2, i64 %180
%182 = load double, double* %181
%183 = and i64 %179, 4294967295
%184 = getelementptr inbounds double, double* %2, i64 %183
store double %182, double* %184
%185 = getelementptr inbounds double, double* %3, i64 %180
%186 = load double, double* %185
%187 = getelementptr inbounds double, double* %3, i64 %183
store double %186, double* %187
%188 = add nsw i64 %170, 2
%189 = trunc i64 %188 to i32
%190 = icmp eq i32 %30, %189
br i1 %190, label %71, label %169
191: ; preds = %167, %191
%192 = phi i64 [ %195, %191 ], [ %168, %167 ]
%193 = getelementptr inbounds double, double* %2, i64 %192
store double 0.000000e+00, double* %193
%194 = getelementptr inbounds double, double* %3, i64 %192
store double 0.000000e+00, double* %194
%195 = add nuw nsw i64 %192, 1
%196 = trunc i64 %192 to i32
%197 = icmp slt i32 %196, %1
br i1 %197, label %191, label %198
198: ; preds = %191, %165, %71
ret i32 0
}