#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 */
long double expl(long double);
long double const sqrtl(long double);
void fgammal(long double *fm, const int m, long double t) {
  register int i;
  long double mphalf, t2, expt, term, sum, tcrit;
  const long double sqrpi = 8.8622692545275801365546944e-01L;
  const long double eps = 1.e-14;
  tcrit = (m + m) + 33.;
  mphalf = (long double)m + 0.5;
  if (t > tcrit) {
    fm[0] = sqrpi / sqrtl(t);
    for (i = 1; i <= m; i++)
      fm[i] = fm[i - 1] * (i - 0.5) / t;
    return;
  }
  term = 0.5 / mphalf;
  sum = term;
  for (i = 1; i <= 1000; i++) {
    term *= (t / (mphalf + i));
    sum += term;
    if (term < eps)
      break;
  }
  t2 = t + t;
  expt = expl(-t);
  fm[m] = sum * expt;
  for (i = m - 1; i >= 0; i--)
    fm[i] = (fm[i + 1] * t2 + expt) / (2 * i + 1);
  return;
}


x86 assembly (clang, -O0):
<fgammal>:
push   %rbp
mov    %rsp,%rbp
sub    $0xe0,%rsp
fldt   0x10(%rbp)
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
fstpt  -0x20(%rbp)
fldt   0xec2(%rip)
fstpt  -0xa0(%rbp)
fldl   0xec6(%rip)
fstpt  -0xb0(%rbp)
mov    -0xc(%rbp),%eax
add    %eax,%eax
cvtsi2sd %eax,%xmm0
movsd  0xeb7(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
fldl   -0xc0(%rbp)
fstpt  -0x90(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0xb4(%rbp)
fildl  -0xb4(%rbp)
flds   0xe9a(%rip)
faddp  %st,%st(1)
fstpt  -0x40(%rbp)
fldt   -0x20(%rbp)
fldt   -0x90(%rbp)
fxch   %st(1)
fucomip %st(1),%st
fstp   %st(0)
jbe    1233 <fgammal+0x113>
jmp    11ad <fgammal+0x8d>
fldt   -0x20(%rbp)
mov    %rsp,%rax
fstpt  (%rax)
call   1040 <sqrtl@plt>
fldt   0xe40(%rip)
fdivp  %st,%st(1)
mov    -0x8(%rbp),%rax
fstpt  (%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     122e <fgammal+0x10e>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %ecx,%edx
mov    %edx,%esi
add    $0xffffffff,%esi
movslq %esi,%rsi
shl    $0x4,%rsi
fldt   (%rax,%rsi,1)
cvtsi2sd %edx,%xmm0
movsd  0xe20(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
fldl   -0xd0(%rbp)
fmulp  %st,%st(1)
fldt   -0x20(%rbp)
fdivrp %st,%st(1)
shl    $0x4,%rcx
fstpt  (%rax,%rcx,1)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11cf <fgammal+0xaf>
jmp    1343 <fgammal+0x223>
fldt   -0x40(%rbp)
flds   0xdec(%rip)
fdivp  %st,%st(1)
fstpt  -0x70(%rbp)
fldt   -0x70(%rbp)
fstpt  -0x80(%rbp)
movl   $0x1,-0x24(%rbp)
cmpl   $0x3e8,-0x24(%rbp)
jg     12b7 <fgammal+0x197>
fldt   -0x20(%rbp)
fldt   -0x40(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0xc4(%rbp)
fildl  -0xc4(%rbp)
faddp  %st,%st(1)
fdivrp %st,%st(1)
fldt   -0x70(%rbp)
fmulp  %st,%st(1)
fstpt  -0x70(%rbp)
fldt   -0x70(%rbp)
fldt   -0x80(%rbp)
faddp  %st,%st(1)
fstpt  -0x80(%rbp)
fldt   -0x70(%rbp)
fldl   0xd80(%rip)
fucomip %st(1),%st
fstp   %st(0)
jbe    12a4 <fgammal+0x184>
jmp    129f <fgammal+0x17f>
jmp    12b7 <fgammal+0x197>
jmp    12a9 <fgammal+0x189>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    124e <fgammal+0x12e>
fldt   -0x20(%rbp)
fadd   %st,%st(0)
fstpt  -0x50(%rbp)
fldt   -0x20(%rbp)
fchs
mov    %rsp,%rax
fstpt  (%rax)
call   1030 <expl@plt>
fstpt  -0x60(%rbp)
fldt   -0x80(%rbp)
fldt   -0x60(%rbp)
fmulp  %st,%st(1)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x4,%rcx
fstpt  (%rax,%rcx,1)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     133e <fgammal+0x21e>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %ecx,%esi
shl    $0x4,%rcx
fldt   0x10(%rax,%rcx,1)
fldt   -0x50(%rbp)
fldt   -0x60(%rbp)
fxch   %st(2)
fmulp  %st,%st(1)
faddp  %st,%st(1)
mov    %esi,%edx
lea    0x1(%rdx,%rdx,1),%edx
mov    %edx,-0xc8(%rbp)
fildl  -0xc8(%rbp)
fdivrp %st,%st(1)
fstpt  (%rax,%rcx,1)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    12f1 <fgammal+0x1d1>
jmp    1343 <fgammal+0x223>
add    $0xe0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fgammal>:
push   %r14
push   %rbx
sub    $0x58,%rsp
mov    %esi,%r14d
mov    %rdi,%rbx
lea    (%r14,%r14,1),%eax
cvtsi2sd %eax,%xmm0
addsd  0xec3(%rip),%xmm0
fldt   0x70(%rsp)
movsd  %xmm0,0x40(%rsp)
fldl   0x40(%rsp)
mov    %esi,0x18(%rsp)
fildl  0x18(%rsp)
fadds  0xebf(%rip)
fxch   %st(2)
fucomi %st(1),%st
fstp   %st(1)
jbe    1242 <fgammal+0x122>
fstp   %st(1)
fld    %st(0)
fstpt  (%rsp)
fstpt  0x1c(%rsp)
call   1040 <sqrtl@plt>
fldt   0x1c(%rsp)
fldt   0xea1(%rip)
fdivp  %st,%st(2)
fld    %st(1)
fstpt  (%rbx)
test   %r14d,%r14d
jle    12f9 <fgammal+0x1d9>
mov    %r14d,%ecx
mov    $0x1,%eax
cmp    $0x1,%r14d
je     120c <fgammal+0xec>
mov    %ecx,%edx
and    $0xfffffffe,%edx
lea    0x20(%rbx),%rsi
mov    $0x1,%eax
movsd  0xe5e(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x38(%rsp)
fxch   %st(1)
fmull  0x38(%rsp)
fdiv   %st(1),%st
fld    %st(0)
fstpt  -0x10(%rsi)
lea    0x1(%rax),%rdi
xorps  %xmm1,%xmm1
cvtsi2sd %edi,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x30(%rsp)
fmull  0x30(%rsp)
fdiv   %st(1),%st
fld    %st(0)
fstpt  (%rsi)
add    $0x2,%rax
add    $0x20,%rsi
cmp    %rdx,%rdi
fxch   %st(1)
jne    11c0 <fgammal+0xa0>
test   $0x1,%cl
je     12f9 <fgammal+0x1d9>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
addsd  0xdec(%rip),%xmm0
movsd  %xmm0,0x28(%rsp)
fxch   %st(1)
fmull  0x28(%rsp)
fdivp  %st,%st(1)
shl    $0x4,%rax
fstpt  (%rbx,%rax,1)
fldz
fldz
jmp    12f9 <fgammal+0x1d9>
fld    %st(1)
fdivrs 0xdce(%rip)
mov    $0x1,%eax
fldl   0xdb3(%rip)
fld    %st(1)
nopw   0x0(%rax,%rax,1)
mov    %eax,0x14(%rsp)
fld    %st(4)
fiaddl 0x14(%rsp)
fdivr  %st(4),%st
fmulp  %st,%st(3)
fadd   %st(2),%st
fxch   %st(1)
fucomi %st(2),%st
fxch   %st(1)
ja     1284 <fgammal+0x164>
lea    0x1(%rax),%ecx
cmp    $0x3e8,%eax
mov    %ecx,%eax
jne    1260 <fgammal+0x140>
fstp   %st(4)
fstp   %st(1)
fstp   %st(0)
fld    %st(0)
fadd   %st(1),%st
fstpt  0x1c(%rsp)
fchs
fstpt  (%rsp)
fstpt  0x4c(%rsp)
call   1030 <expl@plt>
fldt   0x1c(%rsp)
fldt   0x4c(%rsp)
fmul   %st(2),%st
movslq %r14d,%rax
shl    $0x4,%rax
fstpt  (%rbx,%rax,1)
test   %r14d,%r14d
jle    12f9 <fgammal+0x1d9>
mov    %r14d,%eax
lea    0x1(%rax),%rcx
shl    $0x4,%rax
add    %rbx,%rax
lea    (%r14,%r14,1),%edx
add    $0xffffffff,%edx
xchg   %ax,%ax
lea    -0x2(%rcx),%esi
fldt   (%rax)
fmul   %st(1),%st
fadd   %st(2),%st
mov    %edx,0x10(%rsp)
fidivl 0x10(%rsp)
shl    $0x4,%rsi
fstpt  (%rbx,%rsi,1)
add    $0xffffffffffffffff,%rcx
add    $0xfffffffffffffff0,%rax
add    $0xfffffffe,%edx
cmp    $0x1,%rcx
ja     12d0 <fgammal+0x1b0>
fstp   %st(1)
fstp   %st(0)
add    $0x58,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @fgammal(x86_fp80* %0, i32 %1, x86_fp80 %2) {
%4 = alloca x86_fp80*
%5 = alloca i32
%6 = alloca x86_fp80
%7 = alloca i32
%8 = alloca x86_fp80
%9 = alloca x86_fp80
%10 = alloca x86_fp80
%11 = alloca x86_fp80
%12 = alloca x86_fp80
%13 = alloca x86_fp80
%14 = alloca x86_fp80
%15 = alloca x86_fp80
store x86_fp80* %0, x86_fp80** %4
store i32 %1, i32* %5
store x86_fp80 %2, x86_fp80* %6
store x86_fp80 0xK3FFEE2DFC48DA77B553D, x86_fp80* %14
store x86_fp80 0xK3FD0B424DC35095CD800, x86_fp80* %15
%16 = load i32, i32* %5
%17 = load i32, i32* %5
%18 = add nsw i32 %16, %17
%19 = sitofp i32 %18 to double
%20 = fadd double %19, 3.300000e+01
%21 = fpext double %20 to x86_fp80
store x86_fp80 %21, x86_fp80* %13
%22 = load i32, i32* %5
%23 = sitofp i32 %22 to x86_fp80
%24 = fadd x86_fp80 %23, 0xK3FFE8000000000000000
store x86_fp80 %24, x86_fp80* %8
%25 = load x86_fp80, x86_fp80* %6
%26 = load x86_fp80, x86_fp80* %13
%27 = fcmp ogt x86_fp80 %25, %26
br i1 %27, label %28, label %60
28: ; preds = %3
%29 = load x86_fp80, x86_fp80* %6
%30 = call x86_fp80 @sqrtl(x86_fp80 %29)
%31 = fdiv x86_fp80 0xK3FFEE2DFC48DA77B553D, %30
%32 = load x86_fp80*, x86_fp80** %4
%33 = getelementptr inbounds x86_fp80, x86_fp80* %32, i64 0
store x86_fp80 %31, x86_fp80* %33
store i32 1, i32* %7
br label %34
34: ; preds = %56, %28
%35 = load i32, i32* %7
%36 = load i32, i32* %5
%37 = icmp sle i32 %35, %36
br i1 %37, label %38, label %59
38: ; preds = %34
%39 = load x86_fp80*, x86_fp80** %4
%40 = load i32, i32* %7
%41 = sub nsw i32 %40, 1
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds x86_fp80, x86_fp80* %39, i64 %42
%44 = load x86_fp80, x86_fp80* %43
%45 = load i32, i32* %7
%46 = sitofp i32 %45 to double
%47 = fsub double %46, 5.000000e-01
%48 = fpext double %47 to x86_fp80
%49 = fmul x86_fp80 %44, %48
%50 = load x86_fp80, x86_fp80* %6
%51 = fdiv x86_fp80 %49, %50
%52 = load x86_fp80*, x86_fp80** %4
%53 = load i32, i32* %7
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds x86_fp80, x86_fp80* %52, i64 %54
store x86_fp80 %51, x86_fp80* %55
br label %56
56: ; preds = %38
%57 = load i32, i32* %7
%58 = add nsw i32 %57, 1
store i32 %58, i32* %7
br label %34
59: ; preds = %34
br label %128
60: ; preds = %3
%61 = load x86_fp80, x86_fp80* %8
%62 = fdiv x86_fp80 0xK3FFE8000000000000000, %61
store x86_fp80 %62, x86_fp80* %11
%63 = load x86_fp80, x86_fp80* %11
store x86_fp80 %63, x86_fp80* %12
store i32 1, i32* %7
br label %64
64: ; preds = %83, %60
%65 = load i32, i32* %7
%66 = icmp sle i32 %65, 1000
br i1 %66, label %67, label %86
67: ; preds = %64
%68 = load x86_fp80, x86_fp80* %6
%69 = load x86_fp80, x86_fp80* %8
%70 = load i32, i32* %7
%71 = sitofp i32 %70 to x86_fp80
%72 = fadd x86_fp80 %69, %71
%73 = fdiv x86_fp80 %68, %72
%74 = load x86_fp80, x86_fp80* %11
%75 = fmul x86_fp80 %74, %73
store x86_fp80 %75, x86_fp80* %11
%76 = load x86_fp80, x86_fp80* %11
%77 = load x86_fp80, x86_fp80* %12
%78 = fadd x86_fp80 %77, %76
store x86_fp80 %78, x86_fp80* %12
%79 = load x86_fp80, x86_fp80* %11
%80 = fcmp olt x86_fp80 %79, 0xK3FD0B424DC35095CD800
br i1 %80, label %81, label %82
81: ; preds = %67
br label %86
82: ; preds = %67
br label %83
83: ; preds = %82
%84 = load i32, i32* %7
%85 = add nsw i32 %84, 1
store i32 %85, i32* %7
br label %64
86: ; preds = %81, %64
%87 = load x86_fp80, x86_fp80* %6
%88 = load x86_fp80, x86_fp80* %6
%89 = fadd x86_fp80 %87, %88
store x86_fp80 %89, x86_fp80* %9
%90 = load x86_fp80, x86_fp80* %6
%91 = fneg x86_fp80 %90
%92 = call x86_fp80 @expl(x86_fp80 %91)
store x86_fp80 %92, x86_fp80* %10
%93 = load x86_fp80, x86_fp80* %12
%94 = load x86_fp80, x86_fp80* %10
%95 = fmul x86_fp80 %93, %94
%96 = load x86_fp80*, x86_fp80** %4
%97 = load i32, i32* %5
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds x86_fp80, x86_fp80* %96, i64 %98
store x86_fp80 %95, x86_fp80* %99
%100 = load i32, i32* %5
%101 = sub nsw i32 %100, 1
store i32 %101, i32* %7
br label %102
102: ; preds = %124, %86
%103 = load i32, i32* %7
%104 = icmp sge i32 %103, 0
br i1 %104, label %105, label %127
105: ; preds = %102
%106 = load x86_fp80*, x86_fp80** %4
%107 = load i32, i32* %7
%108 = add nsw i32 %107, 1
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds x86_fp80, x86_fp80* %106, i64 %109
%111 = load x86_fp80, x86_fp80* %110
%112 = load x86_fp80, x86_fp80* %9
%113 = load x86_fp80, x86_fp80* %10
%114 = call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %111, x86_fp80 %112, x86_fp80 %113)
%115 = load i32, i32* %7
%116 = mul nsw i32 2, %115
%117 = add nsw i32 %116, 1
%118 = sitofp i32 %117 to x86_fp80
%119 = fdiv x86_fp80 %114, %118
%120 = load x86_fp80*, x86_fp80** %4
%121 = load i32, i32* %7
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds x86_fp80, x86_fp80* %120, i64 %122
store x86_fp80 %119, x86_fp80* %123
br label %124
124: ; preds = %105
%125 = load i32, i32* %7
%126 = add nsw i32 %125, -1
store i32 %126, i32* %7
br label %102
127: ; preds = %102
br label %128
128: ; preds = %127, %59
ret void
}

LLVM IR (-O2):
define void @fgammal(x86_fp80* nocapture %0, i32 %1, x86_fp80 %2) {
%4 = shl nsw i32 %1, 1
%5 = sitofp i32 %4 to double
%6 = fadd double %5, 3.300000e+01
%7 = fpext double %6 to x86_fp80
%8 = sitofp i32 %1 to x86_fp80
%9 = fadd x86_fp80 %8, 0xK3FFE8000000000000000
%10 = fcmp olt x86_fp80 %7, %2
br i1 %10, label %11, label %43
11: ; preds = %3
%12 = tail call x86_fp80 @sqrtl(x86_fp80 %2)
%13 = fdiv x86_fp80 0xK3FFEE2DFC48DA77B553D, %12
store x86_fp80 %13, x86_fp80* %0
%14 = icmp slt i32 %1, 1
br i1 %14, label %94, label %15
15: ; preds = %11
%16 = zext i32 %1 to i64
%17 = and i64 %16, 1
%18 = icmp eq i32 %1, 1
br i1 %18, label %82, label %19
19: ; preds = %15
%20 = and i64 %16, 4294967294
br label %21
21: ; preds = %21, %19
%22 = phi x86_fp80 [ %13, %19 ], [ %38, %21 ]
%23 = phi i64 [ 1, %19 ], [ %40, %21 ]
%24 = phi i64 [ 0, %19 ], [ %41, %21 ]
%25 = trunc i64 %23 to i32
%26 = sitofp i32 %25 to double
%27 = fadd double %26, -5.000000e-01
%28 = fpext double %27 to x86_fp80
%29 = fmul x86_fp80 %22, %28
%30 = fdiv x86_fp80 %29, %2
%31 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %23
store x86_fp80 %30, x86_fp80* %31
%32 = add nuw nsw i64 %23, 1
%33 = trunc i64 %32 to i32
%34 = sitofp i32 %33 to double
%35 = fadd double %34, -5.000000e-01
%36 = fpext double %35 to x86_fp80
%37 = fmul x86_fp80 %30, %36
%38 = fdiv x86_fp80 %37, %2
%39 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %32
store x86_fp80 %38, x86_fp80* %39
%40 = add nuw nsw i64 %23, 2
%41 = add i64 %24, 2
%42 = icmp eq i64 %41, %20
br i1 %42, label %82, label %21
43: ; preds = %3
%44 = fdiv x86_fp80 0xK3FFE8000000000000000, %9
br label %45
45: ; preds = %45, %43
%46 = phi i32 [ 1, %43 ], [ %55, %45 ]
%47 = phi x86_fp80 [ %44, %43 ], [ %52, %45 ]
%48 = phi x86_fp80 [ %44, %43 ], [ %53, %45 ]
%49 = sitofp i32 %46 to x86_fp80
%50 = fadd x86_fp80 %9, %49
%51 = fdiv x86_fp80 %2, %50
%52 = fmul x86_fp80 %47, %51
%53 = fadd x86_fp80 %48, %52
%54 = fcmp olt x86_fp80 %52, 0xK3FD0B424DC35095CD800
%55 = add nuw nsw i32 %46, 1
%56 = icmp eq i32 %55, 1001
%57 = select i1 %54, i1 true, i1 %56
br i1 %57, label %58, label %45
58: ; preds = %45
%59 = fadd x86_fp80 %2, %2
%60 = fneg x86_fp80 %2
%61 = tail call x86_fp80 @expl(x86_fp80 %60)
%62 = fmul x86_fp80 %53, %61
%63 = sext i32 %1 to i64
%64 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %63
store x86_fp80 %62, x86_fp80* %64
%65 = icmp sgt i32 %1, 0
br i1 %65, label %66, label %94
66: ; preds = %58
%67 = zext i32 %1 to i64
br label %68
68: ; preds = %66, %68
%69 = phi i64 [ %67, %66 ], [ %70, %68 ]
%70 = add nsw i64 %69, -1
%71 = trunc i64 %70 to i32
%72 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %69
%73 = load x86_fp80, x86_fp80* %72
%74 = tail call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %73, x86_fp80 %59, x86_fp80 %61)
%75 = shl nuw nsw i32 %71, 1
%76 = or i32 %75, 1
%77 = sitofp i32 %76 to x86_fp80
%78 = fdiv x86_fp80 %74, %77
%79 = and i64 %70, 4294967295
%80 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %79
store x86_fp80 %78, x86_fp80* %80
%81 = icmp ugt i64 %69, 1
br i1 %81, label %68, label %94
82: ; preds = %21, %15
%83 = phi x86_fp80 [ %13, %15 ], [ %38, %21 ]
%84 = phi i64 [ 1, %15 ], [ %40, %21 ]
%85 = icmp eq i64 %17, 0
br i1 %85, label %94, label %86
86: ; preds = %82
%87 = trunc i64 %84 to i32
%88 = sitofp i32 %87 to double
%89 = fadd double %88, -5.000000e-01
%90 = fpext double %89 to x86_fp80
%91 = fmul x86_fp80 %83, %90
%92 = fdiv x86_fp80 %91, %2
%93 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %84
store x86_fp80 %92, x86_fp80* %93
br label %94
94: ; preds = %68, %86, %82, %58, %11
ret void
}