#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);
double fabs(float);
float log(float);
float expint(int n, float x) {
  void nrerror(char error_text[]);
  int i, ii, nm1;
  float a, b, c, d, del, fact, h, psi, ans;
  nm1 = n - 1;
  if (n < 0 || x < 0.0 || (x == 0.0 && (n == 0 || n == 1)))
    nrerror("bad arguments in expint");
  else {
    if (n == 0)
      ans = exp(-x) / x;
    else {
      if (x == 0.0)
        ans = 1.0 / nm1;
      else {
        if (x > 1.0) {
          b = x + n;
          c = 1.0 / 1.0e-30;
          d = 1.0 / b;
          h = d;
          for (i = 1; i <= 100; i++) {
            a = -i * (nm1 + i);
            b += 2.0;
            d = 1.0 / (a * d + b);
            c = b + a / c;
            del = c * d;
            h *= del;
            if (fabs(del - 1.0) < 1.0e-7) {
              ans = h * exp(-x);
              return ans;
            }
          }
          nrerror("continued fraction failed in expint");
        } else {
          ans = (nm1 != 0 ? 1.0 / nm1 : -log(x) - 0.5772156649);
          fact = 1.0;
          for (i = 1; i <= 100; i++) {
            fact *= -x / i;
            if (i != nm1)
              del = -fact / (i - nm1);
            else {
              psi = -0.5772156649;
              for (ii = 1; ii <= nm1; ii++)
                psi += 1.0 / ii;
              del = fact * (-log(x) + psi);
            }
            ans += del;
            if (fabs(del) < fabs(ans) * 1.0e-7)
              return ans;
          }
          nrerror("series failed in expint");
        }
      }
    }
  }
  return ans;
}


x86 assembly (clang, -O0):
<expint>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x0,-0x8(%rbp)
jl     11a9 <expint+0x69>
movss  -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
ja     11a9 <expint+0x69>
movss  -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11ba <expint+0x7a>
jp     11ba <expint+0x7a>
cmpl   $0x0,-0x8(%rbp)
je     11a9 <expint+0x69>
cmpl   $0x1,-0x8(%rbp)
jne    11ba <expint+0x7a>
lea    0xe7c(%rip),%rdi
call   1030 <nrerror@plt>
jmp    15b0 <expint+0x470>
cmpl   $0x0,-0x8(%rbp)
jne    11ea <expint+0xaa>
movss  -0xc(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
call   1050 <exp@plt>
divss  -0xc(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
jmp    15ab <expint+0x46b>
movss  -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1225 <expint+0xe5>
jp     1225 <expint+0xe5>
cvtsi2sdl -0x18(%rbp),%xmm1
movsd  0xded(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
jmp    15a6 <expint+0x466>
movss  -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xdca(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    13cc <expint+0x28c>
movss  -0xc(%rbp),%xmm0
cvtsi2ssl -0x8(%rbp),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  0xdcd(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd8f(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x1,-0x10(%rbp)
cmpl   $0x64,-0x10(%rbp)
jg     13bb <expint+0x27b>
xor    %eax,%eax
sub    -0x10(%rbp),%eax
mov    -0x18(%rbp),%ecx
add    -0x10(%rbp),%ecx
imul   %ecx,%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd57(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x1c(%rbp),%xmm0
movss  -0x28(%rbp),%xmm2
movss  -0x20(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd0f(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x20(%rbp),%xmm0
movss  -0x1c(%rbp),%xmm1
divss  -0x24(%rbp),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
mulss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xcbb(%rip),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xcb3(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    13a8 <expint+0x268>
movss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0xc(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
call   1050 <exp@plt>
movaps %xmm0,%xmm1
movss  -0x40(%rbp),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    15ba <expint+0x47a>
jmp    13ad <expint+0x26d>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    128f <expint+0x14f>
lea    0xc82(%rip),%rdi
call   1030 <nrerror@plt>
jmp    15a1 <expint+0x461>
cmpl   $0x0,-0x18(%rbp)
je     13f1 <expint+0x2b1>
cvtsi2sdl -0x18(%rbp),%xmm1
movsd  0xc1d(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    141d <expint+0x2dd>
movss  -0xc(%rbp),%xmm0
call   1060 <log@plt>
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xbf4(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  0xbed(%rip),%xmm0
movss  %xmm0,-0x30(%rbp)
movl   $0x1,-0x10(%rbp)
cmpl   $0x64,-0x10(%rbp)
jg     1595 <expint+0x455>
movss  -0xc(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
cvtsi2ssl -0x10(%rbp),%xmm1
divss  %xmm1,%xmm0
mulss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
je     14a4 <expint+0x364>
movss  -0x30(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x10(%rbp),%eax
sub    -0x18(%rbp),%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
jmp    1530 <expint+0x3f0>
movss  0xb78(%rip),%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     14f9 <expint+0x3b9>
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  0xb2f(%rip),%xmm1
divsd  %xmm0,%xmm1
movss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    14b8 <expint+0x378>
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0xc(%rbp),%xmm0
call   1060 <log@plt>
movaps %xmm0,%xmm1
movss  -0x4c(%rbp),%xmm0
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
addss  -0x38(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x2c(%rbp),%xmm0
call   1040 <fabs@plt>
movsd  %xmm0,-0x58(%rbp)
movss  -0x3c(%rbp),%xmm0
call   1040 <fabs@plt>
movsd  -0x58(%rbp),%xmm1
movsd  0xaab(%rip),%xmm2
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    1582 <expint+0x442>
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    15ba <expint+0x47a>
jmp    1587 <expint+0x447>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    143f <expint+0x2ff>
lea    0xacc(%rip),%rdi
call   1030 <nrerror@plt>
jmp    15a6 <expint+0x466>
jmp    15ab <expint+0x46b>
jmp    15b0 <expint+0x470>
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<expint>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x58,%rsp
test   %edi,%edi
js     116a <expint+0x2a>
movaps %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
ja     116a <expint+0x2a>
mov    %edi,%ebx
cmp    $0x1,%edi
ja     1181 <expint+0x41>
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm1
jne    1181 <expint+0x41>
jp     1181 <expint+0x41>
lea    0xedf(%rip),%rdi
call   1030 <nrerror@plt>
add    $0x58,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
test   %ebx,%ebx
je     11b0 <expint+0x70>
lea    -0x1(%rbx),%r14d
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm1
jne    11cc <expint+0x8c>
jp     11cc <expint+0x8c>
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
movsd  0xe75(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
jmp    1176 <expint+0x36>
movaps 0xe89(%rip),%xmm0
xorps  %xmm1,%xmm0
movaps %xmm1,0x10(%rsp)
call   1050 <exp@plt>
divss  0x10(%rsp),%xmm0
jmp    1176 <expint+0x36>
ucomiss 0xe2d(%rip),%xmm1
movaps %xmm1,0x10(%rsp)
jbe    12b1 <expint+0x171>
cvtsi2ss %ebx,%xmm2
addss  %xmm1,%xmm2
movss  0xe12(%rip),%xmm0
divss  %xmm2,%xmm0
mov    $0x1,%ebp
movss  0xe09(%rip),%xmm4
movss  %xmm0,0xc(%rsp)
movaps %xmm0,%xmm1
nopl   0x0(%rax,%rax,1)
lea    (%rbx,%rbp,1),%eax
add    $0xffffffff,%eax
imul   %ebp,%eax
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
addss  0xde2(%rip),%xmm2
mulss  %xmm0,%xmm1
addss  %xmm2,%xmm1
movss  0xdc6(%rip),%xmm3
divss  %xmm1,%xmm3
divss  %xmm4,%xmm0
movss  %xmm2,0x30(%rsp)
addss  %xmm2,%xmm0
movss  %xmm0,0x20(%rsp)
movss  %xmm3,0x40(%rsp)
mulss  %xmm3,%xmm0
movss  0xc(%rsp),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0xc(%rsp)
addss  0xd9c(%rip),%xmm0
call   1040 <fabs@plt>
movsd  0xdaf(%rip),%xmm1
ucomisd %xmm0,%xmm1
ja     12cf <expint+0x18f>
add    $0x1,%ebp
cmp    $0x65,%ebp
movss  0x30(%rsp),%xmm2
movss  0x20(%rsp),%xmm4
movss  0x40(%rsp),%xmm1
jne    1210 <expint+0xd0>
lea    0xdbc(%rip),%rdi
jmp    1171 <expint+0x31>
test   %r14d,%r14d
je     12f2 <expint+0x1b2>
movaps %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsi2sd %r14d,%xmm1
movsd  0xd4f(%rip),%xmm0
divsd  %xmm1,%xmm0
jmp    1312 <expint+0x1d2>
movaps 0x10(%rsp),%xmm0
xorps  0xd65(%rip),%xmm0
call   1050 <exp@plt>
movss  0xc(%rsp),%xmm1
mulss  %xmm0,%xmm1
movaps %xmm1,%xmm0
jmp    1176 <expint+0x36>
movaps %xmm1,%xmm0
call   1060 <log@plt>
movaps 0x10(%rsp),%xmm2
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
movsd  0xd12(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0xc(%rsp)
movaps 0xd1d(%rip),%xmm0
xorps  %xmm2,%xmm0
movaps %xmm0,0x40(%rsp)
mov    %r14d,%r15d
and    $0xfffffffe,%r15d
neg    %r15d
mov    $0x1,%ebp
movss  0xcbe(%rip),%xmm1
movsd  0xcce(%rip),%xmm4
movsd  0xcd6(%rip),%xmm5
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
movaps 0x40(%rsp),%xmm2
divss  %xmm0,%xmm2
mulss  %xmm1,%xmm2
mov    %ebp,%eax
sub    %r14d,%eax
movaps %xmm2,0x30(%rsp)
jne    1460 <expint+0x320>
movss  0xc78(%rip),%xmm2
cmp    $0x2,%ebx
jb     1436 <expint+0x2f6>
mov    $0x1,%ecx
movss  0xc62(%rip),%xmm2
je     140f <expint+0x2cf>
mov    $0x1,%eax
movss  0xc53(%rip),%xmm2
movapd %xmm4,%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm4,%xmm1
divsd  %xmm0,%xmm1
cvtss2sd %xmm2,%xmm2
addsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
lea    0x1(%rax),%ecx
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
movapd %xmm4,%xmm3
divsd  %xmm2,%xmm3
cvtss2sd %xmm1,%xmm1
addsd  %xmm3,%xmm1
xorps  %xmm2,%xmm2
cvtsd2ss %xmm1,%xmm2
addsd  %xmm5,%xmm0
mov    %eax,%ecx
add    $0x2,%ecx
lea    (%r15,%rax,1),%edx
add    $0x2,%edx
mov    %ecx,%eax
cmp    $0x1,%edx
jne    13c0 <expint+0x280>
test   $0x1,%r14b
je     1436 <expint+0x2f6>
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
movapd %xmm4,%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
addsd  %xmm1,%xmm0
xorps  %xmm2,%xmm2
cvtsd2ss %xmm0,%xmm2
movaps %xmm2,0x20(%rsp)
movaps 0x10(%rsp),%xmm0
call   1060 <log@plt>
movaps 0x20(%rsp),%xmm2
subss  %xmm0,%xmm2
mulss  0x30(%rsp),%xmm2
jmp    1472 <expint+0x332>
cs nopw 0x0(%rax,%rax,1)
xorps  0xbd9(%rip),%xmm2
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
divss  %xmm0,%xmm2
movss  0xc(%rsp),%xmm0
addss  %xmm2,%xmm0
movss  %xmm0,0xc(%rsp)
movaps %xmm2,%xmm0
call   1040 <fabs@plt>
movsd  %xmm0,0x20(%rsp)
movss  0xc(%rsp),%xmm0
call   1040 <fabs@plt>
mulsd  0xb8d(%rip),%xmm0
ucomisd 0x20(%rsp),%xmm0
ja     14d8 <expint+0x398>
add    $0x1,%ebp
cmp    $0x65,%ebp
movaps 0x30(%rsp),%xmm1
movsd  0xb5a(%rip),%xmm4
movsd  0xb62(%rip),%xmm5
jne    1360 <expint+0x220>
lea    0xbb9(%rip),%rdi
call   1030 <nrerror@plt>
movss  0xc(%rsp),%xmm0
jmp    1176 <expint+0x36>

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

LLVM IR (-O2):
@.str = [24 x i8] c"bad arguments in expint\00"
@.str.1 = [36 x i8] c"continued fraction failed in expint\00"
@.str.2 = [24 x i8] c"series failed in expint\00"
define float @expint(i32 %0, float %1) {
%3 = add i32 %0, -1
%4 = icmp slt i32 %0, 0
%5 = fcmp olt float %1, 0.000000e+00
%6 = or i1 %4, %5
br i1 %6, label %11, label %7
7: ; preds = %2
%8 = fcmp oeq float %1, 0.000000e+00
%9 = icmp ult i32 %0, 2
%10 = and i1 %9, %8
br i1 %10, label %11, label %12
11: ; preds = %7, %2
tail call void @nrerror(i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0))
br label %135
12: ; preds = %7
%13 = icmp eq i32 %0, 0
br i1 %13, label %14, label %18
14: ; preds = %12
%15 = fneg float %1
%16 = tail call float @exp(float %15)
%17 = fdiv float %16, %1
br label %135
18: ; preds = %12
br i1 %8, label %19, label %23
19: ; preds = %18
%20 = sitofp i32 %3 to double
%21 = fdiv double 1.000000e+00, %20
%22 = fptrunc double %21 to float
br label %135
23: ; preds = %18
%24 = fcmp ogt float %1, 1.000000e+00
br i1 %24, label %25, label %57
25: ; preds = %23
%26 = sitofp i32 %0 to float
%27 = fadd float %26, %1
%28 = fdiv float 1.000000e+00, %27
br label %32
29: ; preds = %32
%30 = add nuw nsw i32 %37, 1
%31 = icmp eq i32 %30, 101
br i1 %31, label %56, label %32
32: ; preds = %25, %29
%33 = phi float [ %28, %25 ], [ %48, %29 ]
%34 = phi float [ %28, %25 ], [ %44, %29 ]
%35 = phi float [ 0x46293E5940000000, %25 ], [ %46, %29 ]
%36 = phi float [ %27, %25 ], [ %42, %29 ]
%37 = phi i32 [ 1, %25 ], [ %30, %29 ]
%38 = add nsw i32 %37, %3
%39 = mul i32 %37, %38
%40 = sub i32 0, %39
%41 = sitofp i32 %40 to float
%42 = fadd float %36, 2.000000e+00
%43 = tail call float @llvm.fmuladd.f32(float %41, float %34, float %42)
%44 = fdiv float 1.000000e+00, %43
%45 = fdiv float %41, %35
%46 = fadd float %42, %45
%47 = fmul float %46, %44
%48 = fmul float %33, %47
%49 = fadd float %47, -1.000000e+00
%50 = tail call double @fabs(float %49)
%51 = fcmp olt double %50, 0x3E7AD7F29ABCAF48
br i1 %51, label %52, label %29
52: ; preds = %32
%53 = fneg float %1
%54 = tail call float @exp(float %53)
%55 = fmul float %48, %54
br label %135
56: ; preds = %29
tail call void @nrerror(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.1, i64 0, i64 0))
br label %135
57: ; preds = %23
%58 = icmp eq i32 %3, 0
br i1 %58, label %62, label %59
59: ; preds = %57
%60 = sitofp i32 %3 to double
%61 = fdiv double 1.000000e+00, %60
br label %67
62: ; preds = %57
%63 = tail call float @log(float %1)
%64 = fneg float %63
%65 = fpext float %64 to double
%66 = fadd double %65, 0xBFE2788CFC6F802A
br label %67
67: ; preds = %62, %59
%68 = phi double [ %61, %59 ], [ %66, %62 ]
%69 = fptrunc double %68 to float
%70 = fneg float %1
%71 = icmp ugt i32 %0, 1
%72 = and i32 %3, 1
%73 = icmp eq i32 %0, 2
%74 = and i32 %3, -2
%75 = icmp eq i32 %72, 0
br label %79
76: ; preds = %127
%77 = add nuw nsw i32 %82, 1
%78 = icmp eq i32 %77, 101
br i1 %78, label %134, label %79
79: ; preds = %67, %76
%80 = phi float [ %69, %67 ], [ %129, %76 ]
%81 = phi float [ 1.000000e+00, %67 ], [ %85, %76 ]
%82 = phi i32 [ 1, %67 ], [ %77, %76 ]
%83 = sitofp i32 %82 to float
%84 = fdiv float %70, %83
%85 = fmul float %81, %84
%86 = icmp eq i32 %82, %3
br i1 %86, label %87, label %89
87: ; preds = %79
br i1 %71, label %88, label %122
88: ; preds = %87
br i1 %73, label %112, label %94
89: ; preds = %79
%90 = fneg float %85
%91 = sub nsw i32 %82, %3
%92 = sitofp i32 %91 to float
%93 = fdiv float %90, %92
br label %127
94: ; preds = %88, %94
%95 = phi float [ %108, %94 ], [ 0xBFE2788D00000000, %88 ]
%96 = phi i32 [ %109, %94 ], [ 1, %88 ]
%97 = phi i32 [ %110, %94 ], [ 0, %88 ]
%98 = sitofp i32 %96 to double
%99 = fdiv double 1.000000e+00, %98
%100 = fpext float %95 to double
%101 = fadd double %99, %100
%102 = fptrunc double %101 to float
%103 = add nuw nsw i32 %96, 1
%104 = sitofp i32 %103 to double
%105 = fdiv double 1.000000e+00, %104
%106 = fpext float %102 to double
%107 = fadd double %105, %106
%108 = fptrunc double %107 to float
%109 = add nuw nsw i32 %96, 2
%110 = add i32 %97, 2
%111 = icmp eq i32 %110, %74
br i1 %111, label %112, label %94
112: ; preds = %94, %88
%113 = phi float [ undef, %88 ], [ %108, %94 ]
%114 = phi float [ 0xBFE2788D00000000, %88 ], [ %108, %94 ]
%115 = phi i32 [ 1, %88 ], [ %109, %94 ]
br i1 %75, label %122, label %116
116: ; preds = %112
%117 = sitofp i32 %115 to double
%118 = fdiv double 1.000000e+00, %117
%119 = fpext float %114 to double
%120 = fadd double %118, %119
%121 = fptrunc double %120 to float
br label %122
122: ; preds = %116, %112, %87
%123 = phi float [ 0xBFE2788D00000000, %87 ], [ %113, %112 ], [ %121, %116 ]
%124 = tail call float @log(float %1)
%125 = fsub float %123, %124
%126 = fmul float %85, %125
br label %127
127: ; preds = %122, %89
%128 = phi float [ %93, %89 ], [ %126, %122 ]
%129 = fadd float %80, %128
%130 = tail call double @fabs(float %128)
%131 = tail call double @fabs(float %129)
%132 = fmul double %131, 0x3E7AD7F29ABCAF48
%133 = fcmp olt double %130, %132
br i1 %133, label %135, label %76
134: ; preds = %76
tail call void @nrerror(i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str.2, i64 0, i64 0))
br label %135
135: ; preds = %127, %11, %19, %134, %56, %14, %52
%136 = phi float [ %55, %52 ], [ undef, %11 ], [ %17, %14 ], [ %22, %19 ], [ undef, %56 ], [ %129, %134 ], [ %129, %127 ]
ret float %136
}