#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 exp(double);
double log(double);
scalar_t__ modf(double, double *);
double expint(double n, double z) {
  double y;
  {
    double junk;
    if (n < 0 || z < 0 || modf(n, &junk))
      return -1.0;
    if (z == 0 && n < 2)
      return -1.0;
  }
  if (n == 0)
    return exp(-z) / z;
  if (z == 0)
    return 1 / (n - 1);
  if (z > 3) {
    double n0, n1, n2, n3, d0, d1, d2, d3, y_prev = 1;
    int i;
    n0 = 0;
    n1 = 1;
    d0 = 1;
    d1 = z;
    for (i = 0; i < 333; i++) {
      n2 = n0 * (n + i) + n1;
      d2 = d0 * (n + i) + d1;
      n3 = n1 * (1 + i) + n2 * z;
      d3 = d1 * (1 + i) + d2 * z;
      y = n3 / d3;
      if (y == y_prev)
        break;
      y_prev = y;
      n0 = n2;
      n1 = n3;
      d0 = d2;
      d1 = d3;
      if (d3 >= 1.0E+37) {
        n0 /= 1.0E+37;
        n1 /= 1.0E+37;
        d0 /= 1.0E+37;
        d1 /= 1.0E+37;
      }
    }
    y = exp(-z) * y;
  } else {
    double euler_constant = .577215664901532860606512;
    double y_prev = 0;
    double t, m;
    y = -euler_constant - log(z);
    t = 1;
    for (m = 1; m < 333; m++) {
      t = -t * z / m;
      y = y - t / m;
      if (y == y_prev)
        break;
      y_prev = y;
    }
    for (m = 1; m < n; m++)
      y = (exp(-z) - z * y) / m;
  }
  return y;
}


x86 assembly (clang, -O0):
<expint>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
ja     1179 <expint+0x49>
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
ja     1179 <expint+0x49>
movsd  -0x10(%rbp),%xmm0
lea    -0x28(%rbp),%rdi
call   1030 <modf@plt>
cmp    $0x0,%rax
je     118b <expint+0x5b>
movsd  0xeaf(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    162e <expint+0x4fe>
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11c8 <expint+0x98>
jp     11c8 <expint+0x98>
movsd  0xe55(%rip),%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    11c8 <expint+0x98>
movsd  0xe72(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    162e <expint+0x4fe>
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1210 <expint+0xe0>
jp     1210 <expint+0xe0>
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1040 <exp@plt>
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    162e <expint+0x4fe>
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    124f <expint+0x11f>
jp     124f <expint+0x11f>
movsd  -0x10(%rbp),%xmm1
movsd  0xddb(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  0xdcf(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    162e <expint+0x4fe>
movsd  -0x18(%rbp),%xmm0
movsd  0xdac(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1457 <expint+0x327>
movsd  0xda2(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xd8d(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xd80(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x74(%rbp)
cmpl   $0x14d,-0x74(%rbp)
jge    1427 <expint+0x2f7>
movsd  -0x30(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm2
cvtsi2sdl -0x74(%rbp),%xmm1
addsd  %xmm1,%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm2
cvtsi2sdl -0x74(%rbp),%xmm1
addsd  %xmm1,%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x74(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x74(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x48(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
jne    1373 <expint+0x243>
jp     1373 <expint+0x243>
jmp    1427 <expint+0x2f7>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  0xc76(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1414 <expint+0x2e4>
movsd  -0x30(%rbp),%xmm0
movsd  0xc5f(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xc49(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  0xc33(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  0xc1d(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    1419 <expint+0x2e9>
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    12a6 <expint+0x176>
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1040 <exp@plt>
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1624 <expint+0x4f4>
movsd  0xbb9(%rip),%xmm0
movsd  %xmm0,-0x80(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
call   1050 <log@plt>
movaps %xmm0,%xmm1
movsd  -0xa0(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  0xb57(%rip),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0xb47(%rip),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xb47(%rip),%xmm0
ucomisd -0x98(%rbp),%xmm0
ff
jbe    1581 <expint+0x451>
movsd  -0x90(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x18(%rbp),%xmm0
divsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
movsd  -0x90(%rbp),%xmm1
ff
divsd  -0x98(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
ucomisd -0x88(%rbp),%xmm0
ff
jne    1557 <expint+0x427>
jp     1557 <expint+0x427>
jmp    1581 <expint+0x451>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xaa4(%rip),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
jmp    14d1 <expint+0x3a1>
movsd  0xa87(%rip),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm1
ff
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    161f <expint+0x4ef>
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1040 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x20(%rbp)
movsd  0xa06(%rip),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
jmp    1591 <expint+0x461>
jmp    1624 <expint+0x4f4>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<expint>:
sub    $0x48,%rsp
xorpd  %xmm2,%xmm2
ucomisd %xmm0,%xmm2
ja     1189 <expint+0x59>
ucomisd %xmm1,%xmm2
ja     1189 <expint+0x59>
movapd %xmm1,0x10(%rsp)
lea    0x40(%rsp),%rdi
movsd  %xmm0,0x8(%rsp)
call   1030 <modf@plt>
test   %rax,%rax
jne    1189 <expint+0x59>
movsd  0xea1(%rip),%xmm0
movsd  0x8(%rsp),%xmm1
ucomisd %xmm1,%xmm0
movapd 0x10(%rsp),%xmm2
movapd %xmm1,%xmm3
jbe    119a <expint+0x6a>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm2
jne    119a <expint+0x6a>
jp     119a <expint+0x6a>
movsd  0xe6f(%rip),%xmm5
movapd %xmm5,%xmm0
add    $0x48,%rsp
ret
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm3
jne    11c3 <expint+0x93>
jp     11c3 <expint+0x93>
movapd 0xe92(%rip),%xmm0
xorpd  %xmm2,%xmm0
call   1040 <exp@plt>
movapd %xmm0,%xmm5
divsd  0x10(%rsp),%xmm5
jmp    1191 <expint+0x61>
ucomisd %xmm0,%xmm2
jne    11e1 <expint+0xb1>
jp     11e1 <expint+0xb1>
addsd  0xe2d(%rip),%xmm3
movsd  0xe45(%rip),%xmm5
divsd  %xmm3,%xmm5
jmp    1191 <expint+0x61>
ucomisd 0xe27(%rip),%xmm2
jbe    12c2 <expint+0x192>
movapd 0xe59(%rip),%xmm4
unpcklpd %xmm2,%xmm4
movsd  0xe1d(%rip),%xmm7
xorpd  %xmm1,%xmm1
unpcklpd %xmm7,%xmm1
movapd %xmm2,%xmm0
unpcklpd %xmm2,%xmm0
mov    $0x1,%eax
movsd  0xe0f(%rip),%xmm8
movapd 0xe36(%rip),%xmm9
jmp    1247 <expint+0x117>
movapd %xmm6,%xmm4
movsd  %xmm1,%xmm4
shufpd $0x1,%xmm6,%xmm1
add    $0x1,%eax
cmp    $0x14e,%eax
je     13be <expint+0x28e>
movapd %xmm7,%xmm6
lea    -0x1(%rax),%ecx
xorps  %xmm5,%xmm5
cvtsi2sd %ecx,%xmm5
addsd  %xmm3,%xmm5
unpcklpd %xmm5,%xmm5
mulpd  %xmm1,%xmm5
addpd  %xmm4,%xmm5
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movapd %xmm5,%xmm7
mulpd  %xmm0,%xmm7
unpcklpd %xmm1,%xmm1
mulpd  %xmm4,%xmm1
addpd  %xmm7,%xmm1
movapd %xmm1,%xmm4
unpckhpd %xmm1,%xmm4
movapd %xmm1,%xmm7
divsd  %xmm4,%xmm7
ucomisd %xmm6,%xmm7
jne    129c <expint+0x16c>
jnp    13be <expint+0x28e>
ucomisd %xmm8,%xmm4
movapd %xmm5,%xmm6
unpckhpd %xmm1,%xmm6
unpcklpd %xmm5,%xmm1
jb     122c <expint+0xfc>
divpd  %xmm9,%xmm1
divpd  %xmm9,%xmm6
jmp    122c <expint+0xfc>
movapd %xmm2,%xmm0
call   1050 <log@plt>
movapd 0x10(%rsp),%xmm1
movsd  0xd3f(%rip),%xmm2
subsd  %xmm0,%xmm2
xorpd  %xmm5,%xmm5
movsd  0xd37(%rip),%xmm0
movapd 0xd4e(%rip),%xmm8
movsd  0xd2d(%rip),%xmm9
movapd %xmm0,%xmm3
movapd %xmm0,%xmm4
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd %xmm5,%xmm6
xorpd  %xmm8,%xmm4
mulsd  %xmm1,%xmm4
divsd  %xmm3,%xmm4
movapd %xmm4,%xmm7
divsd  %xmm3,%xmm7
movapd %xmm2,%xmm5
subsd  %xmm7,%xmm5
ucomisd %xmm6,%xmm5
addsd  %xmm0,%xmm3
jne    133d <expint+0x20d>
jnp    1348 <expint+0x218>
ucomisd %xmm3,%xmm9
movapd %xmm5,%xmm2
ja     1310 <expint+0x1e0>
movsd  0x8(%rsp),%xmm0
ucomisd 0xcca(%rip),%xmm0
jbe    1191 <expint+0x61>
xorpd  0xcdc(%rip),%xmm1
movsd  0xcb4(%rip),%xmm2
movapd %xmm1,0x10(%rsp)
movsd  %xmm2,0x38(%rsp)
movapd %xmm5,0x20(%rsp)
movapd %xmm1,%xmm0
call   1040 <exp@plt>
movsd  0x38(%rsp),%xmm2
movapd 0x20(%rsp),%xmm5
movapd 0x10(%rsp),%xmm1
mulsd  %xmm1,%xmm5
addsd  %xmm0,%xmm5
divsd  %xmm2,%xmm5
addsd  0xc73(%rip),%xmm2
movsd  0x8(%rsp),%xmm0
ucomisd %xmm2,%xmm0
ja     1372 <expint+0x242>
jmp    1191 <expint+0x61>
xorpd  0xc7a(%rip),%xmm2
movapd %xmm2,%xmm0
movapd %xmm7,0x20(%rsp)
call   1040 <exp@plt>
movapd 0x20(%rsp),%xmm5
mulsd  %xmm0,%xmm5
jmp    1191 <expint+0x61>

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

LLVM IR (-O2):
define double @expint(double %0, double %1) {
%3 = alloca double
%4 = bitcast double* %3 to i8*
%5 = fcmp olt double %0, 0.000000e+00
%6 = fcmp olt double %1, 0.000000e+00
%7 = or i1 %5, %6
br i1 %7, label %93, label %8
8: ; preds = %2
%9 = call i64 @modf(double %0, double* %3)
%10 = icmp eq i64 %9, 0
br i1 %10, label %11, label %92
11: ; preds = %8
%12 = fcmp oeq double %1, 0.000000e+00
%13 = fcmp olt double %0, 2.000000e+00
%14 = and i1 %13, %12
br i1 %14, label %94, label %15
15: ; preds = %11
%16 = fcmp oeq double %0, 0.000000e+00
br i1 %16, label %17, label %21
17: ; preds = %15
%18 = fneg double %1
%19 = tail call double @exp(double %18)
%20 = fdiv double %19, %1
br label %94
21: ; preds = %15
br i1 %12, label %22, label %25
22: ; preds = %21
%23 = fadd double %0, -1.000000e+00
%24 = fdiv double 1.000000e+00, %23
br label %94
25: ; preds = %21
%26 = fcmp ogt double %1, 3.000000e+00
br i1 %26, label %27, label %63
27: ; preds = %25
%28 = insertelement <2 x double> <double 1.000000e+00, double poison>, double %1, i64 1
%29 = insertelement <2 x double> poison, double %1, i64 0
%30 = shufflevector <2 x double> %29, <2 x double> poison, <2 x i32> zeroinitializer
br label %31
31: ; preds = %27, %51
%32 = phi i32 [ %41, %51 ], [ 0, %27 ]
%33 = phi double [ %49, %51 ], [ 1.000000e+00, %27 ]
%34 = phi <2 x double> [ %57, %51 ], [ %28, %27 ]
%35 = phi <2 x double> [ %58, %51 ], [ <double 0.000000e+00, double 1.000000e+00>, %27 ]
%36 = sitofp i32 %32 to double
%37 = fadd double %36, %0
%38 = insertelement <2 x double> poison, double %37, i64 0
%39 = shufflevector <2 x double> %38, <2 x double> poison, <2 x i32> zeroinitializer
%40 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %35, <2 x double> %39, <2 x double> %34)
%41 = add nuw nsw i32 %32, 1
%42 = sitofp i32 %41 to double
%43 = fmul <2 x double> %40, %30
%44 = insertelement <2 x double> poison, double %42, i64 0
%45 = shufflevector <2 x double> %44, <2 x double> poison, <2 x i32> zeroinitializer
%46 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %34, <2 x double> %45, <2 x double> %43)
%47 = extractelement <2 x double> %46, i64 0
%48 = extractelement <2 x double> %46, i64 1
%49 = fdiv double %47, %48
%50 = fcmp oeq double %49, %33
br i1 %50, label %59, label %51
51: ; preds = %31
%52 = fcmp ult double %48, 0x479E17B84357691B
%53 = shufflevector <2 x double> %40, <2 x double> %46, <4 x i32> <i32 2, i32 0, i32 1, i32 3>
%54 = fdiv <4 x double> %53, <double 0x479E17B84357691B, double 0x479E17B84357691B, double 0x479E17B84357691B, double 0x479E17B84357691B>
%55 = select i1 %52, <4 x double> %53, <4 x double> %54
%56 = icmp eq i32 %41, 333
%57 = shufflevector <4 x double> %55, <4 x double> undef, <2 x i32> <i32 0, i32 3>
%58 = shufflevector <4 x double> %55, <4 x double> undef, <2 x i32> <i32 1, i32 2>
br i1 %56, label %59, label %31
59: ; preds = %31, %51
%60 = fneg double %1
%61 = tail call double @exp(double %60)
%62 = fmul double %49, %61
br label %94
63: ; preds = %25
%64 = tail call double @log(double %1)
%65 = fsub double 0xBFE2788CFC6FB619, %64
br label %66
66: ; preds = %66, %63
%67 = phi double [ 1.000000e+00, %63 ], [ %77, %66 ]
%68 = phi double [ 1.000000e+00, %63 ], [ %73, %66 ]
%69 = phi double [ 0.000000e+00, %63 ], [ %75, %66 ]
%70 = phi double [ %65, %63 ], [ %75, %66 ]
%71 = fneg double %68
%72 = fmul double %71, %1
%73 = fdiv double %72, %67
%74 = fdiv double %73, %67
%75 = fsub double %70, %74
%76 = fcmp une double %75, %69
%77 = fadd double %67, 1.000000e+00
%78 = fcmp olt double %77, 3.330000e+02
%79 = select i1 %76, i1 %78, i1 false
br i1 %79, label %66, label %80
80: ; preds = %66
%81 = fcmp ogt double %0, 1.000000e+00
br i1 %81, label %82, label %94
82: ; preds = %80
%83 = fneg double %1
br label %84
84: ; preds = %82, %84
%85 = phi double [ 1.000000e+00, %82 ], [ %90, %84 ]
%86 = phi double [ %75, %82 ], [ %89, %84 ]
%87 = tail call double @exp(double %83)
%88 = tail call double @llvm.fmuladd.f64(double %83, double %86, double %87)
%89 = fdiv double %88, %85
%90 = fadd double %85, 1.000000e+00
%91 = fcmp olt double %90, %0
br i1 %91, label %84, label %94
92: ; preds = %8
br label %94
93: ; preds = %2
br label %94
94: ; preds = %84, %80, %59, %93, %92, %11, %22, %17
%95 = phi double [ %20, %17 ], [ %24, %22 ], [ -1.000000e+00, %11 ], [ -1.000000e+00, %92 ], [ -1.000000e+00, %93 ], [ %62, %59 ], [ %75, %80 ], [ %89, %84 ]
ret double %95
}