#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);
float sqrt(float);
void pearsn(float x[], float y[], unsigned long n, float *r, float *prob,
            float *z) {
  float betai(float a, float b, float x);
  float erfcc(float x);
  unsigned long j;
  float yt, xt, t, df;
  float syy = 0.0, sxy = 0.0, sxx = 0.0, ay = 0.0, ax = 0.0;
  for (j = 1; j <= n; j++) {
    ax += x[j];
    ay += y[j];
  }
  ax /= n;
  ay /= n;
  for (j = 1; j <= n; j++) {
    xt = x[j] - ax;
    yt = y[j] - ay;
    sxx += xt * xt;
    syy += yt * yt;
    sxy += xt * yt;
  }
  *r = sxy / sqrt(sxx * syy);
  *z = 0.5 * log((1.0 + (*r) + 1.0e-20) / (1.0 - (*r) + 1.0e-20));
  df = n - 2;
  t = (*r) * sqrt(df / ((1.0 - (*r) + 1.0e-20) * (1.0 + (*r) + 1.0e-20)));
  *prob = betai(0.5 * df, 0.5, df / (df + t * t));
}


x86 assembly (clang, -O0):
<pearsn>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x5c(%rbp)
movq   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     11d0 <pearsn+0xa0>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x5c(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1183 <pearsn+0x53>
mov    -0x18(%rbp),%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x64(%rbp)
test   %rax,%rax
movss  %xmm0,-0x60(%rbp)
js     120d <pearsn+0xdd>
movss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x60(%rbp),%xmm1
movss  -0x5c(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x6c(%rbp)
test   %rax,%rax
movss  %xmm0,-0x68(%rbp)
js     125d <pearsn+0x12d>
movss  -0x6c(%rbp),%xmm0
movss  %xmm0,-0x68(%rbp)
movss  -0x68(%rbp),%xmm1
movss  -0x58(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x58(%rbp)
movq   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     1319 <pearsn+0x1e9>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x5c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x40(%rbp),%xmm0
movss  -0x40(%rbp),%xmm2
movss  -0x54(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x4c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x40(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x50(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1278 <pearsn+0x148>
movss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x78(%rbp)
movss  -0x54(%rbp),%xmm0
movss  -0x4c(%rbp),%xmm1
mulss  %xmm1,%xmm0
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x78(%rbp),%xmm0
divss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm3
movsd  0xcba(%rip),%xmm1
movaps %xmm3,%xmm0
addsd  %xmm1,%xmm0
movsd  0xca3(%rip),%xmm2
addsd  %xmm2,%xmm0
subsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1040 <log@plt>
movsd  0xc7e(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x30(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x74(%rbp)
test   %rax,%rax
movss  %xmm0,-0x70(%rbp)
js     13db <pearsn+0x2ab>
movss  -0x74(%rbp),%xmm0
movss  %xmm0,-0x70(%rbp)
movss  -0x70(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x7c(%rbp)
movss  -0x48(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xc09(%rip),%xmm1
subsd  %xmm2,%xmm1
movsd  0xbf5(%rip),%xmm2
addsd  %xmm2,%xmm1
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm2
cvtss2sd %xmm2,%xmm3
movsd  0xbe5(%rip),%xmm2
addsd  %xmm3,%xmm2
movsd  0xbd1(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x7c(%rbp),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x48(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xb92(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  -0x48(%rbp),%xmm2
movss  -0x48(%rbp),%xmm3
movss  -0x44(%rbp),%xmm1
movss  -0x44(%rbp),%xmm4
mulss  %xmm4,%xmm1
addss  %xmm3,%xmm1
divss  %xmm1,%xmm2
movss  0xb5a(%rip),%xmm1
call   1030 <betai@plt>
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pearsn>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %r9,%r12
mov    %r8,%r14
mov    %rcx,%r15
mov    %rdx,%r13
xorps  %xmm1,%xmm1
xorps  %xmm6,%xmm6
xorps  %xmm0,%xmm0
test   %rdx,%rdx
je     1314 <pearsn+0x1e4>
lea    -0x1(%r13),%r8
mov    %r13d,%ecx
and    $0x3,%ecx
cmp    $0x3,%r8
jae    117d <pearsn+0x4d>
xorps  %xmm2,%xmm2
mov    $0x1,%edx
test   %rcx,%rcx
jne    11ea <pearsn+0xba>
jmp    1219 <pearsn+0xe9>
mov    %r13,%rbx
and    $0xfffffffffffffffc,%rbx
xorps  %xmm2,%xmm2
xor    %edx,%edx
nopl   0x0(%rax)
movss  0x4(%rsi,%rdx,4),%xmm0
movss  0x8(%rsi,%rdx,4),%xmm3
movss  0xc(%rsi,%rdx,4),%xmm4
movss  0x10(%rsi,%rdx,4),%xmm5
movss  0x4(%rdi,%rdx,4),%xmm6
unpcklps %xmm0,%xmm6
addps  %xmm2,%xmm6
movss  0x8(%rdi,%rdx,4),%xmm0
unpcklps %xmm3,%xmm0
addps  %xmm6,%xmm0
movss  0xc(%rdi,%rdx,4),%xmm3
unpcklps %xmm4,%xmm3
addps  %xmm0,%xmm3
movss  0x10(%rdi,%rdx,4),%xmm2
unpcklps %xmm5,%xmm2
addps  %xmm3,%xmm2
add    $0x4,%rdx
cmp    %rdx,%rbx
jne    1190 <pearsn+0x60>
add    $0x1,%rdx
test   %rcx,%rcx
je     1219 <pearsn+0xe9>
lea    (%rsi,%rdx,4),%rbx
lea    (%rdi,%rdx,4),%rdx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rbx,%rax,4),%xmm0
movss  (%rdx,%rax,4),%xmm3
unpcklps %xmm0,%xmm3
addps  %xmm3,%xmm2
add    $0x1,%rax
cmp    %rax,%rcx
jne    1200 <pearsn+0xd0>
mov    %r13,%rcx
shr    %rcx
mov    %r13d,%eax
and    $0x1,%eax
or     %rcx,%rax
test   %r13,%r13
js     1237 <pearsn+0x107>
xorps  %xmm3,%xmm3
cvtsi2ss %r13,%xmm3
jmp    1243 <pearsn+0x113>
xorps  %xmm3,%xmm3
cvtsi2ss %rax,%xmm3
addss  %xmm3,%xmm3
xorps  %xmm6,%xmm6
xorps  %xmm0,%xmm0
je     1314 <pearsn+0x1e4>
shufps $0x0,%xmm3,%xmm3
divps  %xmm3,%xmm2
test   %r8,%r8
je     144c <pearsn+0x31c>
mov    %r13,%rcx
and    $0xfffffffffffffffe,%rcx
neg    %rcx
xorps  %xmm6,%xmm6
mov    $0x1,%eax
xorps  %xmm3,%xmm3
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rsi,%rax,4),%xmm0
movss  0x4(%rsi,%rax,4),%xmm4
movss  (%rdi,%rax,4),%xmm5
unpcklps %xmm0,%xmm5
subps  %xmm2,%xmm5
movaps %xmm5,%xmm0
shufps $0x55,%xmm5,%xmm0
mulps  %xmm5,%xmm0
mulss  %xmm5,%xmm5
addss  %xmm3,%xmm5
movss  0x4(%rdi,%rax,4),%xmm3
unpcklps %xmm4,%xmm3
addps  %xmm6,%xmm0
subps  %xmm2,%xmm3
movaps %xmm3,%xmm6
shufps $0x55,%xmm3,%xmm6
mulps  %xmm3,%xmm6
mulss  %xmm3,%xmm3
addss  %xmm5,%xmm3
addps  %xmm0,%xmm6
lea    (%rcx,%rax,1),%rdx
add    $0x2,%rdx
add    $0x2,%rax
cmp    $0x1,%rdx
jne    1280 <pearsn+0x150>
test   $0x1,%r13b
je     1309 <pearsn+0x1d9>
movss  (%rsi,%rax,4),%xmm0
movss  (%rdi,%rax,4),%xmm4
unpcklps %xmm0,%xmm4
subps  %xmm2,%xmm4
movaps %xmm4,%xmm0
shufps $0x55,%xmm4,%xmm0
mulps  %xmm4,%xmm0
mulss  %xmm4,%xmm4
addss  %xmm4,%xmm3
addps  %xmm0,%xmm6
movaps %xmm6,%xmm0
shufps $0x55,%xmm6,%xmm0
mulss  %xmm3,%xmm0
ucomiss %xmm1,%xmm0
jb     131f <pearsn+0x1ef>
sqrtss %xmm0,%xmm0
jmp    132e <pearsn+0x1fe>
movaps %xmm6,0x10(%rsp)
call   1050 <sqrt@plt>
movaps 0x10(%rsp),%xmm6
divss  %xmm0,%xmm6
movss  %xmm6,(%r15)
xorps  %xmm0,%xmm0
cvtss2sd %xmm6,%xmm0
movsd  0xcba(%rip),%xmm2
movapd %xmm2,%xmm1
subsd  %xmm0,%xmm1
addsd  %xmm2,%xmm0
movsd  0xcae(%rip),%xmm2
addsd  %xmm2,%xmm0
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1040 <log@plt>
mulsd  0xc9d(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r12)
add    $0xfffffffffffffffe,%r13
js     138d <pearsn+0x25d>
xorps  %xmm2,%xmm2
cvtsi2ss %r13,%xmm2
jmp    13a6 <pearsn+0x276>
mov    %r13,%rax
shr    %rax
and    $0x1,%r13d
or     %rax,%r13
xorps  %xmm2,%xmm2
cvtsi2ss %r13,%xmm2
addss  %xmm2,%xmm2
movss  (%r15),%xmm6
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm6,%xmm1
movsd  0xc3f(%rip),%xmm4
movapd %xmm4,%xmm3
subsd  %xmm1,%xmm3
movsd  0xc37(%rip),%xmm5
addsd  %xmm5,%xmm3
addsd  %xmm4,%xmm1
addsd  %xmm5,%xmm1
mulsd  %xmm3,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
xorpd  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     13f8 <pearsn+0x2c8>
sqrtss %xmm0,%xmm0
jmp    1415 <pearsn+0x2e5>
movss  %xmm2,0x10(%rsp)
movss  %xmm6,0xc(%rsp)
call   1050 <sqrt@plt>
movss  0xc(%rsp),%xmm6
movss  0x10(%rsp),%xmm2
mulss  %xmm0,%xmm6
movss  0xbf7(%rip),%xmm1
movaps %xmm2,%xmm0
mulss  %xmm1,%xmm0
mulss  %xmm6,%xmm6
addss  %xmm2,%xmm6
divss  %xmm6,%xmm2
call   1030 <betai@plt>
movss  %xmm0,(%r14)
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
xorps  %xmm6,%xmm6
mov    $0x1,%eax
xorps  %xmm3,%xmm3
test   $0x1,%r13b
jne    12e4 <pearsn+0x1b4>
jmp    1309 <pearsn+0x1d9>

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

LLVM IR (-O2):
define void @pearsn(float* nocapture %0, float* nocapture %1, i64 %2, float* nocapture %3, float* nocapture writeonly %4, float* nocapture writeonly %5) {
%7 = icmp eq i64 %2, 0
br i1 %7, label %138, label %8
8: ; preds = %6
%9 = add i64 %2, -1
%10 = and i64 %2, 3
%11 = icmp ult i64 %9, 3
br i1 %11, label %52, label %12
12: ; preds = %8
%13 = and i64 %2, -4
br label %14
14: ; preds = %14, %12
%15 = phi i64 [ 1, %12 ], [ %49, %14 ]
%16 = phi <2 x float> [ zeroinitializer, %12 ], [ %48, %14 ]
%17 = phi i64 [ 0, %12 ], [ %50, %14 ]
%18 = getelementptr inbounds float, float* %0, i64 %15
%19 = load float, float* %18
%20 = getelementptr inbounds float, float* %1, i64 %15
%21 = load float, float* %20
%22 = insertelement <2 x float> poison, float %19, i64 0
%23 = insertelement <2 x float> %22, float %21, i64 1
%24 = fadd <2 x float> %16, %23
%25 = add nuw nsw i64 %15, 1
%26 = getelementptr inbounds float, float* %0, i64 %25
%27 = load float, float* %26
%28 = getelementptr inbounds float, float* %1, i64 %25
%29 = load float, float* %28
%30 = insertelement <2 x float> poison, float %27, i64 0
%31 = insertelement <2 x float> %30, float %29, i64 1
%32 = fadd <2 x float> %24, %31
%33 = add nuw nsw i64 %15, 2
%34 = getelementptr inbounds float, float* %0, i64 %33
%35 = load float, float* %34
%36 = getelementptr inbounds float, float* %1, i64 %33
%37 = load float, float* %36
%38 = insertelement <2 x float> poison, float %35, i64 0
%39 = insertelement <2 x float> %38, float %37, i64 1
%40 = fadd <2 x float> %32, %39
%41 = add nuw i64 %15, 3
%42 = getelementptr inbounds float, float* %0, i64 %41
%43 = load float, float* %42
%44 = getelementptr inbounds float, float* %1, i64 %41
%45 = load float, float* %44
%46 = insertelement <2 x float> poison, float %43, i64 0
%47 = insertelement <2 x float> %46, float %45, i64 1
%48 = fadd <2 x float> %40, %47
%49 = add nuw i64 %15, 4
%50 = add i64 %17, 4
%51 = icmp eq i64 %50, %13
br i1 %51, label %52, label %14
52: ; preds = %14, %8
%53 = phi <2 x float> [ undef, %8 ], [ %48, %14 ]
%54 = phi i64 [ 1, %8 ], [ %49, %14 ]
%55 = phi <2 x float> [ zeroinitializer, %8 ], [ %48, %14 ]
%56 = icmp eq i64 %10, 0
br i1 %56, label %71, label %57
57: ; preds = %52, %57
%58 = phi i64 [ %68, %57 ], [ %54, %52 ]
%59 = phi <2 x float> [ %67, %57 ], [ %55, %52 ]
%60 = phi i64 [ %69, %57 ], [ 0, %52 ]
%61 = getelementptr inbounds float, float* %0, i64 %58
%62 = load float, float* %61
%63 = getelementptr inbounds float, float* %1, i64 %58
%64 = load float, float* %63
%65 = insertelement <2 x float> poison, float %62, i64 0
%66 = insertelement <2 x float> %65, float %64, i64 1
%67 = fadd <2 x float> %59, %66
%68 = add nuw i64 %58, 1
%69 = add i64 %60, 1
%70 = icmp eq i64 %69, %10
br i1 %70, label %71, label %57
71: ; preds = %57, %52
%72 = phi <2 x float> [ %53, %52 ], [ %67, %57 ]
%73 = uitofp i64 %2 to float
%74 = insertelement <2 x float> poison, float %73, i64 0
%75 = shufflevector <2 x float> %74, <2 x float> poison, <2 x i32> zeroinitializer
%76 = fdiv <2 x float> %72, %75
br i1 %7, label %138, label %77
77: ; preds = %71
%78 = and i64 %2, 1
%79 = icmp eq i64 %9, 0
br i1 %79, label %113, label %80
80: ; preds = %77
%81 = and i64 %2, -2
br label %82
82: ; preds = %82, %80
%83 = phi float [ 0.000000e+00, %80 ], [ %107, %82 ]
%84 = phi i64 [ 1, %80 ], [ %110, %82 ]
%85 = phi <2 x float> [ zeroinitializer, %80 ], [ %109, %82 ]
%86 = phi i64 [ 0, %80 ], [ %111, %82 ]
%87 = getelementptr inbounds float, float* %0, i64 %84
%88 = load float, float* %87
%89 = getelementptr inbounds float, float* %1, i64 %84
%90 = load float, float* %89
%91 = insertelement <2 x float> poison, float %88, i64 0
%92 = insertelement <2 x float> %91, float %90, i64 1
%93 = fsub <2 x float> %92, %76
%94 = extractelement <2 x float> %93, i64 0
%95 = tail call float @llvm.fmuladd.f32(float %94, float %94, float %83)
%96 = shufflevector <2 x float> %93, <2 x float> undef, <2 x i32> <i32 1, i32 1>
%97 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %93, <2 x float> %96, <2 x float> %85)
%98 = add nuw i64 %84, 1
%99 = getelementptr inbounds float, float* %0, i64 %98
%100 = load float, float* %99
%101 = getelementptr inbounds float, float* %1, i64 %98
%102 = load float, float* %101
%103 = insertelement <2 x float> poison, float %100, i64 0
%104 = insertelement <2 x float> %103, float %102, i64 1
%105 = fsub <2 x float> %104, %76
%106 = extractelement <2 x float> %105, i64 0
%107 = tail call float @llvm.fmuladd.f32(float %106, float %106, float %95)
%108 = shufflevector <2 x float> %105, <2 x float> undef, <2 x i32> <i32 1, i32 1>
%109 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %105, <2 x float> %108, <2 x float> %97)
%110 = add nuw i64 %84, 2
%111 = add i64 %86, 2
%112 = icmp eq i64 %111, %81
br i1 %112, label %113, label %82
113: ; preds = %82, %77
%114 = phi float [ undef, %77 ], [ %107, %82 ]
%115 = phi <2 x float> [ undef, %77 ], [ %109, %82 ]
%116 = phi float [ 0.000000e+00, %77 ], [ %107, %82 ]
%117 = phi i64 [ 1, %77 ], [ %110, %82 ]
%118 = phi <2 x float> [ zeroinitializer, %77 ], [ %109, %82 ]
%119 = icmp eq i64 %78, 0
br i1 %119, label %132, label %120
120: ; preds = %113
%121 = getelementptr inbounds float, float* %0, i64 %117
%122 = load float, float* %121
%123 = getelementptr inbounds float, float* %1, i64 %117
%124 = load float, float* %123
%125 = insertelement <2 x float> poison, float %122, i64 0
%126 = insertelement <2 x float> %125, float %124, i64 1
%127 = fsub <2 x float> %126, %76
%128 = extractelement <2 x float> %127, i64 0
%129 = tail call float @llvm.fmuladd.f32(float %128, float %128, float %116)
%130 = shufflevector <2 x float> %127, <2 x float> undef, <2 x i32> <i32 1, i32 1>
%131 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %127, <2 x float> %130, <2 x float> %118)
br label %132
132: ; preds = %113, %120
%133 = phi float [ %114, %113 ], [ %129, %120 ]
%134 = phi <2 x float> [ %115, %113 ], [ %131, %120 ]
%135 = extractelement <2 x float> %134, i64 1
%136 = fmul float %135, %133
%137 = extractelement <2 x float> %134, i64 0
br label %138
138: ; preds = %6, %132, %71
%139 = phi float [ 0.000000e+00, %71 ], [ %137, %132 ], [ 0.000000e+00, %6 ]
%140 = phi float [ 0.000000e+00, %71 ], [ %136, %132 ], [ 0.000000e+00, %6 ]
%141 = tail call float @sqrt(float %140)
%142 = fdiv float %139, %141
store float %142, float* %3
%143 = fpext float %142 to double
%144 = fadd double %143, 1.000000e+00
%145 = fadd double %144, 0x3BC79CA10C924223
%146 = fsub double 1.000000e+00, %143
%147 = fadd double %146, 0x3BC79CA10C924223
%148 = fdiv double %145, %147
%149 = tail call double @log(double %148)
%150 = fmul double %149, 5.000000e-01
%151 = fptrunc double %150 to float
store float %151, float* %5
%152 = add i64 %2, -2
%153 = uitofp i64 %152 to float
%154 = load float, float* %3
%155 = fpext float %153 to double
%156 = fpext float %154 to double
%157 = fsub double 1.000000e+00, %156
%158 = fadd double %157, 0x3BC79CA10C924223
%159 = fadd double %156, 1.000000e+00
%160 = fadd double %159, 0x3BC79CA10C924223
%161 = fmul double %158, %160
%162 = fdiv double %155, %161
%163 = fptrunc double %162 to float
%164 = tail call float @sqrt(float %163)
%165 = fmul float %154, %164
%166 = fmul float %153, 5.000000e-01
%167 = tail call float @llvm.fmuladd.f32(float %165, float %165, float %153)
%168 = fdiv float %153, %167
%169 = tail call float @betai(float %166, float 5.000000e-01, float %168)
store float %169, float* %4
ret void
}