#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 */
float exp(float);
float fgamma_0f(float);
float fgammaf(int, float);
void fgamma_setf(int max, float x, float *fset) {
  float a, b, c;
  int i, j;
  if (x < 0.0005) {
    for (i = 0; i <= max; i++)
      fset[i] = 1.0 / (i + i + 1.0) - x / (i + i + 3.0) +
                x * x / (i + i + i + i + 10.0);
    return;
  }
  fset[0] = fgamma_0f(x);
  if (max < 1)
    return;
  a = 1.0 / (x + x);
  c = a * exp(-x);
  for (i = 1; i <= max; i++) {
    if ((b = (i + i - 1) * a) > 1.0)
      break;
    fset[i] = b * fset[i - 1] - c;
  }
  if (max < i)
    return;
  fset[max] = fgammaf(max, x);
  a = 1.0 / a;
  c = c * a;
  for (j = max - 1; j >= i; j--)
    fset[j] = (a * fset[j + 1] + c) / (j + j + 1);
}


x86 assembly (clang, -O0):
<fgamma_setf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
movss  %xmm0,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xeab(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1211 <fgamma_setf+0xe1>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     120c <fgamma_setf+0xdc>
mov    -0x20(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xe84(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xe78(%rip),%xmm0
divsd  %xmm1,%xmm0
movss  -0x8(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    -0x20(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm2
movsd  0xe69(%rip),%xmm3
addsd  %xmm3,%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movss  -0x8(%rbp),%xmm1
mulss  -0x8(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    -0x20(%rbp),%eax
add    -0x20(%rbp),%eax
add    -0x20(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm2
movsd  0xe2f(%rip),%xmm3
addsd  %xmm3,%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1166 <fgamma_setf+0x36>
jmp    13d5 <fgamma_setf+0x2a5>
movss  -0x8(%rbp),%xmm0
call   1040 <fgamma_0f@plt>
mov    -0x10(%rbp),%rax
movss  %xmm0,(%rax)
cmpl   $0x1,-0x4(%rbp)
jge    1232 <fgamma_setf+0x102>
jmp    13d5 <fgamma_setf+0x2a5>
movss  -0x8(%rbp),%xmm0
addss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdc0(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x8(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
call   1050 <exp@plt>
movaps %xmm0,%xmm1
movss  -0x28(%rbp),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1318 <fgamma_setf+0x1e8>
mov    -0x20(%rbp),%eax
add    -0x20(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
cvtss2sd %xmm0,%xmm0
movsd  0xd4b(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12cc <fgamma_setf+0x19c>
jmp    1318 <fgamma_setf+0x1e8>
movss  -0x18(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x1c(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    128e <fgamma_setf+0x15e>
mov    -0x4(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1329 <fgamma_setf+0x1f9>
jmp    13d5 <fgamma_setf+0x2a5>
mov    -0x4(%rbp),%edi
movss  -0x8(%rbp),%xmm0
call   1030 <fgammaf@plt>
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xcb4(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x1c(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jl     13d5 <fgamma_setf+0x2a5>
movss  -0x14(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x1c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x24(%rbp),%eax
add    -0x24(%rbp),%eax
add    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1379 <fgamma_setf+0x249>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fgamma_setf>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x28,%rsp
mov    %rsi,%r15
mov    %edi,%r14d
cvtss2sd %xmm0,%xmm9
movsd  0xeb3(%rip),%xmm1
ucomisd %xmm9,%xmm1
jbe    117b <fgamma_setf+0x4b>
test   %r14d,%r14d
js     149b <fgamma_setf+0x36b>
mulss  %xmm0,%xmm0
cvtss2sd %xmm0,%xmm8
lea    0x1(%r14),%eax
cmp    $0x3,%r14d
jae    12f8 <fgamma_setf+0x1c8>
xor    %ecx,%ecx
jmp    1417 <fgamma_setf+0x2e7>
movaps %xmm0,0x10(%rsp)
call   1040 <fgamma_0f@plt>
movss  %xmm0,(%r15)
test   %r14d,%r14d
jle    149b <fgamma_setf+0x36b>
movaps 0x10(%rsp),%xmm1
movaps %xmm1,%xmm0
addss  %xmm1,%xmm0
movss  0xe79(%rip),%xmm2
divss  %xmm0,%xmm2
movss  %xmm2,0xc(%rsp)
movaps %xmm1,%xmm0
xorps  0xe75(%rip),%xmm0
call   1050 <exp@plt>
movaps 0xe69(%rip),%xmm4
movss  0xc(%rsp),%xmm3
movss  0xe4b(%rip),%xmm1
movaps %xmm0,%xmm2
mulss  %xmm3,%xmm2
xorps  %xmm2,%xmm4
mov    %r14d,%eax
mov    $0x1,%ecx
xor    %ebx,%ebx
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
mulss  %xmm3,%xmm0
ucomiss %xmm1,%xmm0
ja     1222 <fgamma_setf+0xf2>
mulss  (%r15,%rbx,4),%xmm0
addss  %xmm4,%xmm0
movss  %xmm0,0x4(%r15,%rbx,4)
add    $0x1,%rbx
add    $0x2,%ecx
cmp    %rbx,%rax
jne    11f0 <fgamma_setf+0xc0>
jmp    149b <fgamma_setf+0x36b>
lea    0x1(%rbx),%ebp
mov    %r14d,%edi
movaps 0x10(%rsp),%xmm0
movaps %xmm2,0x10(%rsp)
call   1030 <fgammaf@plt>
movss  0xc(%rsp),%xmm1
movaps 0x10(%rsp),%xmm3
movss  0xdd6(%rip),%xmm2
movslq %r14d,%rcx
movss  %xmm0,(%r15,%rcx,4)
cmp    %r14d,%ebp
jge    149b <fgamma_setf+0x36b>
divss  %xmm1,%xmm2
mulss  %xmm2,%xmm3
mov    %ebx,%edx
not    %edx
add    %r14d,%edx
mov    %rcx,%rax
test   $0x1,%dl
je     1297 <fgamma_setf+0x167>
lea    -0x1(%rcx),%rax
lea    (%rax,%rax,1),%edx
add    $0x1,%edx
xorps  %xmm1,%xmm1
cvtsi2ss %edx,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm3,%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x4(%r15,%rcx,4)
add    $0xfffffffffffffffe,%rcx
cmp    %rbx,%rcx
je     149b <fgamma_setf+0x36b>
add    $0x1,%rbx
lea    (%rax,%rax,1),%ecx
add    $0xfffffffd,%ecx
xchg   %ax,%ax
mulss  %xmm2,%xmm0
addss  %xmm3,%xmm0
lea    0x2(%rcx),%edx
xorps  %xmm1,%xmm1
cvtsi2ss %edx,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x4(%r15,%rax,4)
mulss  %xmm2,%xmm0
addss  %xmm3,%xmm0
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x8(%r15,%rax,4)
add    $0xfffffffffffffffe,%rax
add    $0xfffffffc,%ecx
cmp    %rbx,%rax
jg     12b0 <fgamma_setf+0x180>
jmp    149b <fgamma_setf+0x36b>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm9,%xmm13
unpcklpd %xmm9,%xmm13
movaps %xmm8,%xmm14
movlhps %xmm8,%xmm14
lea    0x0(,%rax,4),%rdx
and    $0xfffffffffffffff0,%rdx
movdqa 0xd1d(%rip),%xmm6
xor    %esi,%esi
movapd 0xd22(%rip),%xmm15
movapd 0xd29(%rip),%xmm10
movapd 0xd30(%rip),%xmm11
movdqa 0xd37(%rip),%xmm12
movdqa %xmm6,%xmm4
nopl   (%rax)
movdqa %xmm6,%xmm0
paddd  %xmm6,%xmm0
pshufd $0xee,%xmm0,%xmm7
cvtdq2pd %xmm0,%xmm1
cvtdq2pd %xmm7,%xmm2
movapd %xmm1,%xmm0
addpd  %xmm15,%xmm0
movapd %xmm2,%xmm7
addpd  %xmm15,%xmm7
movapd %xmm15,%xmm5
movapd %xmm15,%xmm3
divpd  %xmm0,%xmm3
divpd  %xmm7,%xmm5
addpd  %xmm10,%xmm1
addpd  %xmm10,%xmm2
movapd %xmm13,%xmm0
movapd %xmm13,%xmm7
divpd  %xmm1,%xmm7
movdqa %xmm4,%xmm1
divpd  %xmm2,%xmm0
pslld  $0x2,%xmm1
cvtdq2pd %xmm1,%xmm2
subpd  %xmm7,%xmm3
pshufd $0xee,%xmm1,%xmm1
cvtdq2pd %xmm1,%xmm1
subpd  %xmm0,%xmm5
addpd  %xmm11,%xmm1
addpd  %xmm11,%xmm2
movaps %xmm14,%xmm0
divpd  %xmm2,%xmm0
movaps %xmm14,%xmm2
addpd  %xmm3,%xmm0
divpd  %xmm1,%xmm2
addpd  %xmm5,%xmm2
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm2,%xmm1
unpcklpd %xmm1,%xmm0
movupd %xmm0,(%r15,%rsi,1)
paddd  %xmm12,%xmm6
paddd  %xmm12,%xmm4
add    $0x10,%rsi
cmp    %rsi,%rdx
jne    1350 <fgamma_setf+0x220>
cmp    %rax,%rcx
je     149b <fgamma_setf+0x36b>
lea    (%rcx,%rcx,1),%edx
lea    0x0(,%rcx,4),%esi
movsd  0xbdf(%rip),%xmm0
movsd  0xbdf(%rip),%xmm1
movsd  0xbdf(%rip),%xmm2
nopl   0x0(%rax)
xorps  %xmm3,%xmm3
cvtsi2sd %edx,%xmm3
movapd %xmm3,%xmm4
addsd  %xmm0,%xmm4
movapd %xmm0,%xmm5
divsd  %xmm4,%xmm5
addsd  %xmm1,%xmm3
movapd %xmm9,%xmm4
divsd  %xmm3,%xmm4
subsd  %xmm4,%xmm5
xorps  %xmm3,%xmm3
cvtsi2sd %esi,%xmm3
addsd  %xmm2,%xmm3
movaps %xmm8,%xmm4
divsd  %xmm3,%xmm4
addsd  %xmm5,%xmm4
xorps  %xmm3,%xmm3
cvtsd2ss %xmm4,%xmm3
movss  %xmm3,(%r15,%rcx,4)
add    $0x1,%rcx
add    $0x2,%edx
add    $0x4,%esi
cmp    %rcx,%rax
jne    1440 <fgamma_setf+0x310>
add    $0x28,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @fgamma_setf(i32 %0, float %1, float* nocapture %2) {
%4 = fpext float %1 to double
%5 = fcmp olt double %4, 5.000000e-04
br i1 %5, label %6, label %67
6: ; preds = %3
%7 = icmp slt i32 %0, 0
br i1 %7, label %144, label %8
8: ; preds = %6
%9 = fmul float %1, %1
%10 = fpext float %9 to double
%11 = add nuw i32 %0, 1
%12 = zext i32 %11 to i64
%13 = icmp ult i32 %0, 3
br i1 %13, label %45, label %14
14: ; preds = %8
%15 = and i64 %12, 4294967292
%16 = insertelement <4 x double> poison, double %4, i64 0
%17 = shufflevector <4 x double> %16, <4 x double> poison, <4 x i32> zeroinitializer
%18 = insertelement <4 x double> poison, double %10, i64 0
%19 = shufflevector <4 x double> %18, <4 x double> poison, <4 x i32> zeroinitializer
br label %20
20: ; preds = %20, %14
%21 = phi i64 [ 0, %14 ], [ %39, %20 ]
%22 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %14 ], [ %40, %20 ]
%23 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %14 ], [ %41, %20 ]
%24 = shl <4 x i32> %22, <i32 1, i32 1, i32 1, i32 1>
%25 = sitofp <4 x i32> %24 to <4 x double>
%26 = fadd <4 x double> %25, <double 1.000000e+00, double 1.000000e+00, double 1.000000e+00, double 1.000000e+00>
%27 = fdiv <4 x double> <double 1.000000e+00, double 1.000000e+00, double 1.000000e+00, double 1.000000e+00>, %26
%28 = fadd <4 x double> %25, <double 3.000000e+00, double 3.000000e+00, double 3.000000e+00, double 3.000000e+00>
%29 = fdiv <4 x double> %17, %28
%30 = fsub <4 x double> %27, %29
%31 = shl nsw <4 x i32> %23, <i32 2, i32 2, i32 2, i32 2>
%32 = sitofp <4 x i32> %31 to <4 x double>
%33 = fadd <4 x double> %32, <double 1.000000e+01, double 1.000000e+01, double 1.000000e+01, double 1.000000e+01>
%34 = fdiv <4 x double> %19, %33
%35 = fadd <4 x double> %34, %30
%36 = fptrunc <4 x double> %35 to <4 x float>
%37 = getelementptr inbounds float, float* %2, i64 %21
%38 = bitcast float* %37 to <4 x float>*
store <4 x float> %36, <4 x float>* %38
%39 = add nuw i64 %21, 4
%40 = add <4 x i32> %22, <i32 4, i32 4, i32 4, i32 4>
%41 = add <4 x i32> %23, <i32 4, i32 4, i32 4, i32 4>
%42 = icmp eq i64 %39, %15
br i1 %42, label %43, label %20
43: ; preds = %20
%44 = icmp eq i64 %15, %12
br i1 %44, label %144, label %45
45: ; preds = %8, %43
%46 = phi i64 [ 0, %8 ], [ %15, %43 ]
br label %47
47: ; preds = %45, %47
%48 = phi i64 [ %65, %47 ], [ %46, %45 ]
%49 = trunc i64 %48 to i32
%50 = shl i32 %49, 1
%51 = sitofp i32 %50 to double
%52 = fadd double %51, 1.000000e+00
%53 = fdiv double 1.000000e+00, %52
%54 = fadd double %51, 3.000000e+00
%55 = fdiv double %4, %54
%56 = fsub double %53, %55
%57 = trunc i64 %48 to i32
%58 = shl nsw i32 %57, 2
%59 = sitofp i32 %58 to double
%60 = fadd double %59, 1.000000e+01
%61 = fdiv double %10, %60
%62 = fadd double %61, %56
%63 = fptrunc double %62 to float
%64 = getelementptr inbounds float, float* %2, i64 %48
store float %63, float* %64
%65 = add nuw nsw i64 %48, 1
%66 = icmp eq i64 %65, %12
br i1 %66, label %144, label %47
67: ; preds = %3
%68 = tail call float @fgamma_0f(float %1)
store float %68, float* %2
%69 = icmp slt i32 %0, 1
br i1 %69, label %144, label %70
70: ; preds = %67
%71 = fadd float %1, %1
%72 = fdiv float 1.000000e+00, %71
%73 = fneg float %1
%74 = tail call float @exp(float %73)
%75 = fmul float %72, %74
%76 = fneg float %75
%77 = add nuw i32 %0, 1
%78 = zext i32 %77 to i64
br label %79
79: ; preds = %70, %87
%80 = phi i64 [ 1, %70 ], [ %93, %87 ]
%81 = trunc i64 %80 to i32
%82 = shl i32 %81, 1
%83 = add i32 %82, -1
%84 = sitofp i32 %83 to float
%85 = fmul float %72, %84
%86 = fcmp ogt float %85, 1.000000e+00
br i1 %86, label %95, label %87
87: ; preds = %79
%88 = add nsw i64 %80, -1
%89 = getelementptr inbounds float, float* %2, i64 %88
%90 = load float, float* %89
%91 = tail call float @llvm.fmuladd.f32(float %85, float %90, float %76)
%92 = getelementptr inbounds float, float* %2, i64 %80
store float %91, float* %92
%93 = add nuw nsw i64 %80, 1
%94 = icmp eq i64 %93, %78
br i1 %94, label %144, label %79
95: ; preds = %79
%96 = trunc i64 %80 to i32
%97 = tail call float @fgammaf(i32 %0, float %1)
%98 = sext i32 %0 to i64
%99 = getelementptr inbounds float, float* %2, i64 %98
store float %97, float* %99
%100 = fdiv float 1.000000e+00, %72
%101 = fmul float %100, %75
%102 = icmp slt i32 %96, %0
br i1 %102, label %103, label %144
103: ; preds = %95
%104 = shl i64 %80, 32
%105 = ashr exact i64 %104, 32
%106 = sub i64 %98, %80
%107 = xor i64 %105, -1
%108 = and i64 %106, 1
%109 = icmp eq i64 %108, 0
br i1 %109, label %119, label %110
110: ; preds = %103
%111 = add nsw i64 %98, -1
%112 = tail call float @llvm.fmuladd.f32(float %100, float %97, float %101)
%113 = trunc i64 %111 to i32
%114 = shl nsw i32 %113, 1
%115 = or i32 %114, 1
%116 = sitofp i32 %115 to float
%117 = fdiv float %112, %116
%118 = getelementptr inbounds float, float* %2, i64 %111
store float %117, float* %118
br label %119
119: ; preds = %110, %103
%120 = phi float [ %97, %103 ], [ %117, %110 ]
%121 = phi i64 [ %98, %103 ], [ %111, %110 ]
%122 = sub nsw i64 0, %98
%123 = icmp eq i64 %107, %122
br i1 %123, label %144, label %124
124: ; preds = %119, %124
%125 = phi float [ %141, %124 ], [ %120, %119 ]
%126 = phi i64 [ %135, %124 ], [ %121, %119 ]
%127 = add nsw i64 %126, -1
%128 = tail call float @llvm.fmuladd.f32(float %100, float %125, float %101)
%129 = trunc i64 %127 to i32
%130 = shl nsw i32 %129, 1
%131 = or i32 %130, 1
%132 = sitofp i32 %131 to float
%133 = fdiv float %128, %132
%134 = getelementptr inbounds float, float* %2, i64 %127
store float %133, float* %134
%135 = add nsw i64 %126, -2
%136 = tail call float @llvm.fmuladd.f32(float %100, float %133, float %101)
%137 = trunc i64 %135 to i32
%138 = shl nsw i32 %137, 1
%139 = or i32 %138, 1
%140 = sitofp i32 %139 to float
%141 = fdiv float %136, %140
%142 = getelementptr inbounds float, float* %2, i64 %135
store float %141, float* %142
%143 = icmp sgt i64 %135, %105
br i1 %143, label %124, label %144
144: ; preds = %87, %119, %124, %47, %43, %95, %6, %67
ret void
}