#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 log(double);
double pow(double, double);
void pmfBartzero(double A, double k, int n, double *prob) {
  double eta[n + 1], xi[n + 1];
  int i, j;
  for (i = 0; i <= n; i++)
    eta[i] = 0;
  for (i = 0; i <= n; i++)
    xi[i] = 0;
  for (i = 0; i <= n; i++)
    prob[i] = 0;
  eta[0] = 1.0 + A;
  for (i = 1; i <= n; i++)
    eta[i] = -A / i / (i + 1.0);
  xi[0] = log(1.0 + A);
  for (i = 1; i <= n; i++) {
    for (j = 1; j <= i - 1; j++) {
      xi[i] += j * xi[j] * eta[i - j];
    }
    xi[i] = (i * eta[i] - xi[i]) / (i * (1.0 + A));
  }
  prob[0] = 1.0 / pow(1.0 + A, k);
  for (i = 1; i <= n; i++) {
    for (j = 1; j <= i; j++) {
      prob[i] += j * xi[j] * prob[i - j];
    }
    prob[i] *= (-k / i);
  }
}


x86 assembly (clang, -O0):
<pmfBartzero>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
mov    %rsi,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %rsp,%rcx
mov    %rcx,-0x28(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x50(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x30(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x48(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11b7 <pmfBartzero+0x97>
mov    -0x50(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    118d <pmfBartzero+0x6d>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11e8 <pmfBartzero+0xc8>
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11be <pmfBartzero+0x9e>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1219 <pmfBartzero+0xf9>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11ef <pmfBartzero+0xcf>
mov    -0x50(%rbp),%rax
movsd  0xddb(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1296 <pmfBartzero+0x176>
mov    -0x50(%rbp),%rax
movsd  -0x8(%rbp),%xmm0
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
divsd  %xmm1,%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
movsd  0xd89(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1235 <pmfBartzero+0x115>
movsd  0xd62(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
call   1030 <log@plt>
mov    -0x48(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1392 <pmfBartzero+0x272>
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x3c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1323 <pmfBartzero+0x203>
mov    -0x48(%rbp),%rax
mov    -0x50(%rbp),%rcx
cvtsi2sdl -0x40(%rbp),%xmm0
movslq -0x40(%rbp),%rdx
mulsd  (%rax,%rdx,8),%xmm0
mov    -0x3c(%rbp),%edx
sub    -0x40(%rbp),%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm2
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    12ca <pmfBartzero+0x1aa>
mov    -0x48(%rbp),%rax
mov    -0x50(%rbp),%rcx
cvtsi2sdl -0x3c(%rbp),%xmm0
movslq -0x3c(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm2
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
movsd  0xc92(%rip),%xmm2
addsd  -0x8(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12b7 <pmfBartzero+0x197>
movsd  0xc66(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xc4c(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1477 <pmfBartzero+0x357>
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     1432 <pmfBartzero+0x312>
mov    -0x48(%rbp),%rax
cvtsi2sdl -0x40(%rbp),%xmm0
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
sub    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    13da <pmfBartzero+0x2ba>
movsd  -0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    13c7 <pmfBartzero+0x2a7>
mov    -0x28(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pmfBartzero>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rsi,%r14
movdqa %xmm1,-0x40(%rbp)
lea    0x1(%rdi),%r12d
mov    %rsp,%r15
lea    0xf(,%r12,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%r15
mov    %r15,%rsp
mov    %rsp,%rbx
sub    %rax,%rbx
mov    %rbx,%rsp
test   %edi,%edi
js     1385 <pmfBartzero+0x255>
mov    %edi,%eax
mov    %rdi,-0x58(%rbp)
lea    0x8(,%rax,8),%r13
mov    %r15,%rdi
xor    %esi,%esi
mov    %r13,%rdx
movaps %xmm0,-0x50(%rbp)
call   1030 <memset@plt>
mov    %rbx,%rdi
xor    %esi,%esi
mov    %r13,%rdx
call   1030 <memset@plt>
mov    %r14,%rdi
xor    %esi,%esi
mov    %r13,%rdx
mov    -0x58(%rbp),%r13
call   1030 <memset@plt>
movapd -0x50(%rbp),%xmm7
movsd  0xe43(%rip),%xmm0
addsd  %xmm7,%xmm0
movsd  %xmm0,(%r15)
test   %r13d,%r13d
movsd  %xmm0,-0x30(%rbp)
je     14c6 <pmfBartzero+0x396>
xorpd  0xe34(%rip),%xmm7
lea    -0x1(%r12),%rcx
mov    $0x1,%eax
cmp    $0x2,%rcx
jb     124d <pmfBartzero+0x11d>
mov    %rcx,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rcx,%rax
or     $0x1,%rax
movapd %xmm7,%xmm6
unpcklpd %xmm7,%xmm6
movdqa 0xe16(%rip),%xmm1
xor    %esi,%esi
movapd 0xe1c(%rip),%xmm2
movdqa 0xe24(%rip),%xmm3
nopl   0x0(%rax)
cvtdq2pd %xmm1,%xmm4
movapd %xmm6,%xmm5
divpd  %xmm4,%xmm5
addpd  %xmm2,%xmm4
divpd  %xmm4,%xmm5
movupd %xmm5,0x8(%r15,%rsi,8)
add    $0x2,%rsi
paddd  %xmm3,%xmm1
cmp    %rsi,%rdx
jne    1220 <pmfBartzero+0xf0>
cmp    %rdx,%rcx
je     1286 <pmfBartzero+0x156>
movsd  0xdab(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
nop
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movapd %xmm7,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,(%r15,%rax,8)
add    $0x1,%rax
cmp    %rax,%r12
jne    1260 <pmfBartzero+0x130>
call   1040 <log@plt>
movsd  %xmm0,(%rbx)
test   %r13d,%r13d
je     14da <pmfBartzero+0x3aa>
mov    $0x1,%eax
xor    %r9d,%r9d
mov    %r15,%r8
movaps -0x40(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
jmp    12e7 <pmfBartzero+0x1b7>
xchg   %ax,%ax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
movapd %xmm0,%xmm2
mulsd  %xmm3,%xmm2
mulsd  (%r15,%rax,8),%xmm3
subsd  %xmm4,%xmm3
divsd  %xmm2,%xmm3
movsd  %xmm3,(%rbx,%rax,8)
add    $0x1,%rax
add    $0x1,%r9
add    $0x8,%r8
cmp    %r12,%rax
je     13aa <pmfBartzero+0x27a>
movsd  (%rbx,%rax,8),%xmm4
cmp    $0x2,%rax
jb     12b0 <pmfBartzero+0x180>
mov    $0x1,%esi
cmp    $0x1,%r9
je     1355 <pmfBartzero+0x225>
mov    %r9,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0x1,%esi
mov    %r8,%rdx
nopl   0x0(%rax)
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  (%rbx,%rsi,8),%xmm2
mulsd  (%rdx),%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,(%rbx,%rax,8)
lea    0x1(%rsi),%rcx
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
mulsd  0x8(%rbx,%rsi,8),%xmm4
mulsd  -0x8(%rdx),%xmm4
addsd  %xmm2,%xmm4
movsd  %xmm4,(%rbx,%rax,8)
add    $0x2,%rsi
add    $0xfffffffffffffff0,%rdx
cmp    %rdi,%rcx
jne    1310 <pmfBartzero+0x1e0>
test   $0x1,%r9b
je     12b0 <pmfBartzero+0x180>
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  (%rbx,%rsi,8),%xmm2
mov    %rax,%rcx
sub    %rsi,%rcx
mulsd  (%r15,%rcx,8),%xmm2
addsd  %xmm2,%xmm4
movsd  %xmm4,(%rbx,%rax,8)
jmp    12b0 <pmfBartzero+0x180>
addsd  0xc73(%rip),%xmm0
movapd %xmm0,-0x50(%rbp)
movsd  %xmm0,(%r15)
call   1040 <log@plt>
movsd  %xmm0,(%rbx)
movapd -0x50(%rbp),%xmm2
jmp    14d4 <pmfBartzero+0x3a4>
call   1050 <pow@plt>
movapd -0x40(%rbp),%xmm3
movsd  0xc44(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r14)
xorpd  0xc43(%rip),%xmm3
mov    $0x1,%eax
xor    %r8d,%r8d
mov    %r14,%r9
jmp    140a <pmfBartzero+0x2da>
nopw   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movapd %xmm3,%xmm2
divsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%r14,%rax,8)
add    $0x1,%rax
add    $0x8,%r9
cmp    %r12,%rax
je     14fd <pmfBartzero+0x3cd>
mov    %r8,%rcx
add    $0x1,%r8
movsd  (%r14,%rax,8),%xmm0
test   %rcx,%rcx
je     1490 <pmfBartzero+0x360>
mov    %r8,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0x1,%esi
mov    %r9,%rcx
nopl   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mulsd  (%rbx,%rsi,8),%xmm1
mulsd  (%rcx),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rax,8)
lea    0x1(%rsi),%rdx
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
mulsd  0x8(%rbx,%rsi,8),%xmm0
mulsd  -0x8(%rcx),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r14,%rax,8)
add    $0x2,%rsi
add    $0xfffffffffffffff0,%rcx
cmp    %rdi,%rdx
jne    1430 <pmfBartzero+0x300>
test   $0x1,%r8b
je     13e0 <pmfBartzero+0x2b0>
jmp    149f <pmfBartzero+0x36f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    $0x1,%esi
test   $0x1,%r8b
je     13e0 <pmfBartzero+0x2b0>
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mulsd  (%rbx,%rsi,8),%xmm1
mov    %rax,%rcx
sub    %rsi,%rcx
mulsd  (%r14,%rcx,8),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r14,%rax,8)
jmp    13e0 <pmfBartzero+0x2b0>
call   1040 <log@plt>
movsd  -0x30(%rbp),%xmm2
movsd  %xmm0,(%rbx)
movaps -0x40(%rbp),%xmm1
jmp    14e3 <pmfBartzero+0x3b3>
movaps -0x40(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
movapd %xmm2,%xmm0
call   1050 <pow@plt>
movsd  0xb0c(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r14)
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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