#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 floor(double);
double bessel(double *v, int n, double t) {
  int i, iy, k;
  double ans, p, B, d[6];
  if (t <= 0) {
    ans = v[0];
    goto done;
  }
  if (t >= n - 1) {
    ans = v[n - 1];
    goto done;
  }
  p = floor(t);
  iy = (int)t;
  ans = v[iy];
  k = iy + 1;
  if (k >= n)
    goto done;
  p = t - p;
  ans += p * (v[k] - v[iy]);
  if ((iy - 1 < 0) || (iy + 2 >= n))
    goto done;
  k = iy - 2;
  for (i = 0; i < 5; i++) {
    if ((k < 0) || (k + 1 >= n))
      d[i] = 0;
    else
      d[i] = v[k + 1] - v[k];
    k += 1;
  }
  for (i = 0; i < 4; i++)
    d[i] = d[i + 1] - d[i];
  B = 0.25 * p * (p - 1.0);
  ans += B * (d[1] + d[2]);
  if (iy + 2 >= n)
    goto done;
  for (i = 0; i < 3; i++)
    d[i] = d[i + 1] - d[i];
  B = 2.0 * B / 3.0;
  ans += (p - 0.5) * B * d[1];
  if ((iy - 2 < 0) || (iy + 3 > n))
    goto done;
  for (i = 0; i < 2; i++)
    d[i] = d[i + 1] - d[i];
  B = 0.125 * B * (p + 1.0) * (p - 2.0);
  ans += B * (d[0] + d[1]);
done:
  return ans;
}


x86 assembly (clang, -O0):
<bessel>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
jb     1144 <bessel+0x34>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    148a <bessel+0x37a>
movsd  -0x18(%rbp),%xmm0
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
ucomisd %xmm1,%xmm0
jb     1179 <bessel+0x69>
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    148a <bessel+0x37a>
movsd  -0x18(%rbp),%xmm0
call   1030 <floor@plt>
movsd  %xmm0,-0x38(%rbp)
cvttsd2si -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     11bc <bessel+0xac>
jmp    148a <bessel+0x37a>
movsd  -0x18(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%eax
sub    $0x1,%eax
cmp    $0x0,%eax
jl     121a <bessel+0x10a>
mov    -0x20(%rbp),%eax
add    $0x2,%eax
cmp    -0xc(%rbp),%eax
jl     121f <bessel+0x10f>
jmp    148a <bessel+0x37a>
mov    -0x20(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x5,-0x1c(%rbp)
jge    12a4 <bessel+0x194>
cmpl   $0x0,-0x24(%rbp)
jl     1252 <bessel+0x142>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
cmp    -0xc(%rbp),%eax
jl     1264 <bessel+0x154>
movslq -0x1c(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp,%rax,8)
jmp    128d <bessel+0x17d>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movslq -0x1c(%rbp),%rax
movsd  %xmm0,-0x70(%rbp,%rax,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    122f <bessel+0x11f>
movl   $0x0,-0x1c(%rbp)
cmpl   $0x4,-0x1c(%rbp)
jge    12e5 <bessel+0x1d5>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
cltq
movsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x1c(%rbp),%rax
subsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x1c(%rbp),%rax
movsd  %xmm0,-0x70(%rbp,%rax,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    12ab <bessel+0x19b>
movsd  0xd1b(%rip),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
movsd  0xd01(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
addsd  -0x60(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%eax
add    $0x2,%eax
cmp    -0xc(%rbp),%eax
jl     1341 <bessel+0x231>
jmp    148a <bessel+0x37a>
movl   $0x0,-0x1c(%rbp)
cmpl   $0x3,-0x1c(%rbp)
jge    1382 <bessel+0x272>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
cltq
movsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x1c(%rbp),%rax
subsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x1c(%rbp),%rax
movsd  %xmm0,-0x70(%rbp,%rax,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1348 <bessel+0x238>
movsd  0xc96(%rip),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  0xc81(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xc63(%rip),%xmm1
subsd  %xmm1,%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%eax
sub    $0x2,%eax
cmp    $0x0,%eax
jl     13eb <bessel+0x2db>
mov    -0x20(%rbp),%eax
add    $0x3,%eax
cmp    -0xc(%rbp),%eax
jle    13f0 <bessel+0x2e0>
jmp    148a <bessel+0x37a>
movl   $0x0,-0x1c(%rbp)
cmpl   $0x2,-0x1c(%rbp)
jge    1431 <bessel+0x321>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
cltq
movsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x1c(%rbp),%rax
subsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0x1c(%rbp),%rax
movsd  %xmm0,-0x70(%rbp,%rax,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    13f7 <bessel+0x2e7>
movsd  0xbef(%rip),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  0xbba(%rip),%xmm1
addsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x38(%rbp),%xmm1
movsd  0xbc4(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
addsd  -0x68(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessel>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rdi,%rbx
xorpd  %xmm1,%xmm1
ucomisd %xmm0,%xmm1
jae    11c2 <bessel+0xb2>
mov    %esi,%r14d
lea    -0x1(%r14),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
ucomisd %xmm1,%xmm0
jae    11cc <bessel+0xbc>
cvttsd2si %xmm0,%r15d
movslq %r15d,%r12
movsd  (%rbx,%r12,8),%xmm3
add    $0x1,%r12d
cmp    %r14d,%r12d
jge    11c6 <bessel+0xb6>
movsd  %xmm0,0x10(%rsp)
movsd  %xmm3,0x8(%rsp)
call   1030 <floor@plt>
movsd  0x8(%rsp),%xmm3
movsd  0x10(%rsp),%xmm4
subsd  %xmm0,%xmm4
movslq %r12d,%rax
movsd  (%rbx,%rax,8),%xmm0
subsd  %xmm3,%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm0,%xmm3
test   %r15d,%r15d
jle    11c6 <bessel+0xb6>
lea    0x2(%r15),%eax
cmp    %r14d,%eax
jge    11c6 <bessel+0xb6>
xor    %ecx,%ecx
mov    %r15d,%edx
xorpd  %xmm9,%xmm9
xorpd  %xmm8,%xmm8
sub    $0x2,%edx
jb     11ef <bessel+0xdf>
lea    -0x1(%r15),%ecx
cmp    %r14d,%r15d
jle    11df <bessel+0xcf>
xorpd  %xmm8,%xmm8
jmp    11ef <bessel+0xdf>
movsd  (%rbx),%xmm3
movapd %xmm3,%xmm0
jmp    11d3 <bessel+0xc3>
cltq
movsd  (%rbx,%rax,8),%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
mov    %ecx,%esi
movsd  (%rbx,%rsi,8),%xmm8
mov    %edx,%edx
subsd  (%rbx,%rdx,8),%xmm8
mov    %r15d,%edx
cmp    %r14d,%r15d
jge    1205 <bessel+0xf5>
movsd  (%rbx,%rdx,8),%xmm9
mov    %ecx,%ecx
subsd  (%rbx,%rcx,8),%xmm9
mov    %r12d,%ecx
movsd  (%rbx,%rcx,8),%xmm0
xorpd  %xmm11,%xmm11
xorpd  %xmm10,%xmm10
cmp    %r14d,%eax
jge    122a <bessel+0x11a>
mov    %eax,%esi
movsd  (%rbx,%rsi,8),%xmm10
subsd  (%rbx,%rcx,8),%xmm10
subsd  (%rbx,%rdx,8),%xmm0
lea    0x3(%r15),%ecx
cmp    %r14d,%ecx
jge    1248 <bessel+0x138>
mov    %ecx,%edx
movsd  (%rbx,%rdx,8),%xmm11
mov    %eax,%eax
subsd  (%rbx,%rax,8),%xmm11
movapd %xmm10,%xmm7
subsd  %xmm0,%xmm7
subsd  %xmm9,%xmm0
movsd  0xda2(%rip),%xmm1
addsd  %xmm4,%xmm1
movapd %xmm0,%xmm2
addsd  %xmm7,%xmm2
movapd %xmm7,%xmm5
subsd  %xmm0,%xmm5
movsd  0xd8e(%rip),%xmm6
mulsd  %xmm4,%xmm6
mulsd  %xmm1,%xmm6
mulsd  %xmm6,%xmm2
addsd  %xmm3,%xmm2
addsd  %xmm6,%xmm6
divsd  0xd7a(%rip),%xmm6
movsd  0xd7a(%rip),%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm6,%xmm3
mulsd  %xmm5,%xmm3
addsd  %xmm2,%xmm3
cmp    $0x2,%r15d
jb     11c6 <bessel+0xb6>
cmp    %r14d,%ecx
jg     11c6 <bessel+0xb6>
subsd  %xmm10,%xmm11
subsd  %xmm7,%xmm11
subsd  %xmm8,%xmm9
subsd  %xmm9,%xmm0
subsd  %xmm5,%xmm11
subsd  %xmm0,%xmm5
addsd  %xmm11,%xmm5
mulsd  0xd35(%rip),%xmm6
movsd  0xd35(%rip),%xmm0
addsd  %xmm4,%xmm0
addsd  0xd31(%rip),%xmm4
mulsd  %xmm6,%xmm0
mulsd  %xmm0,%xmm4
mulsd  %xmm5,%xmm4
addsd  %xmm3,%xmm4
movapd %xmm4,%xmm0
jmp    11d3 <bessel+0xc3>

LLVM IR (-O0):
define double @bessel(double* %0, i32 %1, double %2) {
%4 = alloca double*
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca [6 x double]
store double* %0, double** %4
store i32 %1, i32* %5
store double %2, double* %6
%14 = load double, double* %6
%15 = fcmp ole double %14, 0.000000e+00
br i1 %15, label %16, label %20
16: ; preds = %3
%17 = load double*, double** %4
%18 = getelementptr inbounds double, double* %17, i64 0
%19 = load double, double* %18
store double %19, double* %10
br label %236
20: ; preds = %3
%21 = load double, double* %6
%22 = load i32, i32* %5
%23 = sub nsw i32 %22, 1
%24 = sitofp i32 %23 to double
%25 = fcmp oge double %21, %24
br i1 %25, label %26, label %33
26: ; preds = %20
%27 = load double*, double** %4
%28 = load i32, i32* %5
%29 = sub nsw i32 %28, 1
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %27, i64 %30
%32 = load double, double* %31
store double %32, double* %10
br label %236
33: ; preds = %20
%34 = load double, double* %6
%35 = call double @llvm.floor.f64(double %34)
store double %35, double* %11
%36 = load double, double* %6
%37 = fptosi double %36 to i32
store i32 %37, i32* %8
%38 = load double*, double** %4
%39 = load i32, i32* %8
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
store double %42, double* %10
%43 = load i32, i32* %8
%44 = add nsw i32 %43, 1
store i32 %44, i32* %9
%45 = load i32, i32* %9
%46 = load i32, i32* %5
%47 = icmp sge i32 %45, %46
br i1 %47, label %48, label %49
48: ; preds = %33
br label %236
49: ; preds = %33
%50 = load double, double* %6
%51 = load double, double* %11
%52 = fsub double %50, %51
store double %52, double* %11
%53 = load double, double* %11
%54 = load double*, double** %4
%55 = load i32, i32* %9
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %54, i64 %56
%58 = load double, double* %57
%59 = load double*, double** %4
%60 = load i32, i32* %8
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %59, i64 %61
%63 = load double, double* %62
%64 = fsub double %58, %63
%65 = load double, double* %10
%66 = call double @llvm.fmuladd.f64(double %53, double %64, double %65)
store double %66, double* %10
%67 = load i32, i32* %8
%68 = sub nsw i32 %67, 1
%69 = icmp slt i32 %68, 0
br i1 %69, label %75, label %70
70: ; preds = %49
%71 = load i32, i32* %8
%72 = add nsw i32 %71, 2
%73 = load i32, i32* %5
%74 = icmp sge i32 %72, %73
br i1 %74, label %75, label %76
75: ; preds = %70, %49
br label %236
76: ; preds = %70
%77 = load i32, i32* %8
%78 = sub nsw i32 %77, 2
store i32 %78, i32* %9
store i32 0, i32* %7
br label %79
79: ; preds = %113, %76
%80 = load i32, i32* %7
%81 = icmp slt i32 %80, 5
br i1 %81, label %82, label %116
82: ; preds = %79
%83 = load i32, i32* %9
%84 = icmp slt i32 %83, 0
br i1 %84, label %90, label %85
85: ; preds = %82
%86 = load i32, i32* %9
%87 = add nsw i32 %86, 1
%88 = load i32, i32* %5
%89 = icmp sge i32 %87, %88
br i1 %89, label %90, label %94
90: ; preds = %85, %82
%91 = load i32, i32* %7
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %92
store double 0.000000e+00, double* %93
br label %110
94: ; preds = %85
%95 = load double*, double** %4
%96 = load i32, i32* %9
%97 = add nsw i32 %96, 1
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %95, i64 %98
%100 = load double, double* %99
%101 = load double*, double** %4
%102 = load i32, i32* %9
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %101, i64 %103
%105 = load double, double* %104
%106 = fsub double %100, %105
%107 = load i32, i32* %7
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %108
store double %106, double* %109
br label %110
110: ; preds = %94, %90
%111 = load i32, i32* %9
%112 = add nsw i32 %111, 1
store i32 %112, i32* %9
br label %113
113: ; preds = %110
%114 = load i32, i32* %7
%115 = add nsw i32 %114, 1
store i32 %115, i32* %7
br label %79
116: ; preds = %79
store i32 0, i32* %7
br label %117
117: ; preds = %134, %116
%118 = load i32, i32* %7
%119 = icmp slt i32 %118, 4
br i1 %119, label %120, label %137
120: ; preds = %117
%121 = load i32, i32* %7
%122 = add nsw i32 %121, 1
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %123
%125 = load double, double* %124
%126 = load i32, i32* %7
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %127
%129 = load double, double* %128
%130 = fsub double %125, %129
%131 = load i32, i32* %7
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %132
store double %130, double* %133
br label %134
134: ; preds = %120
%135 = load i32, i32* %7
%136 = add nsw i32 %135, 1
store i32 %136, i32* %7
br label %117
137: ; preds = %117
%138 = load double, double* %11
%139 = fmul double 2.500000e-01, %138
%140 = load double, double* %11
%141 = fsub double %140, 1.000000e+00
%142 = fmul double %139, %141
store double %142, double* %12
%143 = load double, double* %12
%144 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 1
%145 = load double, double* %144
%146 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 2
%147 = load double, double* %146
%148 = fadd double %145, %147
%149 = load double, double* %10
%150 = call double @llvm.fmuladd.f64(double %143, double %148, double %149)
store double %150, double* %10
%151 = load i32, i32* %8
%152 = add nsw i32 %151, 2
%153 = load i32, i32* %5
%154 = icmp sge i32 %152, %153
br i1 %154, label %155, label %156
155: ; preds = %137
br label %236
156: ; preds = %137
store i32 0, i32* %7
br label %157
157: ; preds = %174, %156
%158 = load i32, i32* %7
%159 = icmp slt i32 %158, 3
br i1 %159, label %160, label %177
160: ; preds = %157
%161 = load i32, i32* %7
%162 = add nsw i32 %161, 1
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %163
%165 = load double, double* %164
%166 = load i32, i32* %7
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %167
%169 = load double, double* %168
%170 = fsub double %165, %169
%171 = load i32, i32* %7
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %172
store double %170, double* %173
br label %174
174: ; preds = %160
%175 = load i32, i32* %7
%176 = add nsw i32 %175, 1
store i32 %176, i32* %7
br label %157
177: ; preds = %157
%178 = load double, double* %12
%179 = fmul double 2.000000e+00, %178
%180 = fdiv double %179, 3.000000e+00
store double %180, double* %12
%181 = load double, double* %11
%182 = fsub double %181, 5.000000e-01
%183 = load double, double* %12
%184 = fmul double %182, %183
%185 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 1
%186 = load double, double* %185
%187 = load double, double* %10
%188 = call double @llvm.fmuladd.f64(double %184, double %186, double %187)
store double %188, double* %10
%189 = load i32, i32* %8
%190 = sub nsw i32 %189, 2
%191 = icmp slt i32 %190, 0
br i1 %191, label %197, label %192
192: ; preds = %177
%193 = load i32, i32* %8
%194 = add nsw i32 %193, 3
%195 = load i32, i32* %5
%196 = icmp sgt i32 %194, %195
br i1 %196, label %197, label %198
197: ; preds = %192, %177
br label %236
198: ; preds = %192
store i32 0, i32* %7
br label %199
199: ; preds = %216, %198
%200 = load i32, i32* %7
%201 = icmp slt i32 %200, 2
br i1 %201, label %202, label %219
202: ; preds = %199
%203 = load i32, i32* %7
%204 = add nsw i32 %203, 1
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %205
%207 = load double, double* %206
%208 = load i32, i32* %7
%209 = sext i32 %208 to i64
%210 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %209
%211 = load double, double* %210
%212 = fsub double %207, %211
%213 = load i32, i32* %7
%214 = sext i32 %213 to i64
%215 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 %214
store double %212, double* %215
br label %216
216: ; preds = %202
%217 = load i32, i32* %7
%218 = add nsw i32 %217, 1
store i32 %218, i32* %7
br label %199
219: ; preds = %199
%220 = load double, double* %12
%221 = fmul double 1.250000e-01, %220
%222 = load double, double* %11
%223 = fadd double %222, 1.000000e+00
%224 = fmul double %221, %223
%225 = load double, double* %11
%226 = fsub double %225, 2.000000e+00
%227 = fmul double %224, %226
store double %227, double* %12
%228 = load double, double* %12
%229 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 0
%230 = load double, double* %229
%231 = getelementptr inbounds [6 x double], [6 x double]* %13, i64 0, i64 1
%232 = load double, double* %231
%233 = fadd double %230, %232
%234 = load double, double* %10
%235 = call double @llvm.fmuladd.f64(double %228, double %233, double %234)
store double %235, double* %10
br label %236
236: ; preds = %219, %197, %155, %75, %48, %26, %16
%237 = load double, double* %10
ret double %237
}

LLVM IR (-O2):
define double @bessel(double* nocapture %0, i32 %1, double %2) {
%4 = fcmp ugt double %2, 0.000000e+00
br i1 %4, label %7, label %5
5: ; preds = %3
%6 = load double, double* %0
br label %124
7: ; preds = %3
%8 = add nsw i32 %1, -1
%9 = sitofp i32 %8 to double
%10 = fcmp ugt double %9, %2
br i1 %10, label %15, label %11
11: ; preds = %7
%12 = sext i32 %8 to i64
%13 = getelementptr inbounds double, double* %0, i64 %12
%14 = load double, double* %13
br label %124
15: ; preds = %7
%16 = fptosi double %2 to i32
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds double, double* %0, i64 %17
%19 = load double, double* %18
%20 = add nsw i32 %16, 1
%21 = icmp slt i32 %20, %1
br i1 %21, label %22, label %124
22: ; preds = %15
%23 = tail call double @llvm.floor.f64(double %2)
%24 = fsub double %2, %23
%25 = sext i32 %20 to i64
%26 = getelementptr inbounds double, double* %0, i64 %25
%27 = load double, double* %26
%28 = fsub double %27, %19
%29 = tail call double @llvm.fmuladd.f64(double %24, double %28, double %19)
%30 = icmp sgt i32 %16, 0
%31 = add nsw i32 %16, 2
%32 = icmp slt i32 %31, %1
%33 = select i1 %30, i1 %32, i1 false
br i1 %33, label %34, label %124
34: ; preds = %22
%35 = add nsw i32 %16, -2
%36 = icmp ult i32 %16, 2
br i1 %36, label %48, label %37
37: ; preds = %34
%38 = add nsw i32 %16, -1
%39 = icmp sgt i32 %16, %1
br i1 %39, label %48, label %40
40: ; preds = %37
%41 = zext i32 %38 to i64
%42 = getelementptr inbounds double, double* %0, i64 %41
%43 = load double, double* %42
%44 = zext i32 %35 to i64
%45 = getelementptr inbounds double, double* %0, i64 %44
%46 = load double, double* %45
%47 = fsub double %43, %46
br label %48
48: ; preds = %34, %40, %37
%49 = phi i32 [ %38, %40 ], [ %38, %37 ], [ 0, %34 ]
%50 = phi double [ %47, %40 ], [ 0.000000e+00, %37 ], [ 0.000000e+00, %34 ]
%51 = icmp slt i32 %16, %1
br i1 %51, label %52, label %60
52: ; preds = %48
%53 = zext i32 %16 to i64
%54 = getelementptr inbounds double, double* %0, i64 %53
%55 = load double, double* %54
%56 = zext i32 %49 to i64
%57 = getelementptr inbounds double, double* %0, i64 %56
%58 = load double, double* %57
%59 = fsub double %55, %58
br label %60
60: ; preds = %48, %52
%61 = phi double [ %59, %52 ], [ 0.000000e+00, %48 ]
%62 = zext i32 %20 to i64
%63 = getelementptr inbounds double, double* %0, i64 %62
%64 = load double, double* %63
%65 = zext i32 %16 to i64
%66 = getelementptr inbounds double, double* %0, i64 %65
%67 = load double, double* %66
%68 = fsub double %64, %67
%69 = add nuw nsw i32 %16, 2
%70 = icmp slt i32 %69, %1
br i1 %70, label %71, label %79
71: ; preds = %60
%72 = zext i32 %69 to i64
%73 = getelementptr inbounds double, double* %0, i64 %72
%74 = load double, double* %73
%75 = zext i32 %20 to i64
%76 = getelementptr inbounds double, double* %0, i64 %75
%77 = load double, double* %76
%78 = fsub double %74, %77
br label %79
79: ; preds = %71, %60
%80 = phi double [ %78, %71 ], [ 0.000000e+00, %60 ]
%81 = add nuw nsw i32 %16, 3
%82 = icmp slt i32 %81, %1
br i1 %82, label %83, label %91
83: ; preds = %79
%84 = zext i32 %81 to i64
%85 = getelementptr inbounds double, double* %0, i64 %84
%86 = load double, double* %85
%87 = zext i32 %69 to i64
%88 = getelementptr inbounds double, double* %0, i64 %87
%89 = load double, double* %88
%90 = fsub double %86, %89
br label %91
91: ; preds = %79, %83
%92 = phi double [ %90, %83 ], [ 0.000000e+00, %79 ]
%93 = fsub double %68, %61
%94 = fsub double %80, %68
%95 = fadd double %24, -1.000000e+00
%96 = fadd double %93, %94
%97 = fsub double %94, %93
%98 = fmul double %24, 2.500000e-01
%99 = fmul double %98, %95
%100 = tail call double @llvm.fmuladd.f64(double %99, double %96, double %29)
%101 = fmul double %99, 2.000000e+00
%102 = fdiv double %101, 3.000000e+00
%103 = fadd double %24, -5.000000e-01
%104 = fmul double %103, %102
%105 = tail call double @llvm.fmuladd.f64(double %104, double %97, double %100)
%106 = icmp ult i32 %16, 2
%107 = add nuw nsw i32 %16, 3
%108 = icmp sgt i32 %107, %1
%109 = select i1 %106, i1 true, i1 %108
br i1 %109, label %124, label %110
110: ; preds = %91
%111 = fsub double %92, %80
%112 = fsub double %111, %94
%113 = fsub double %61, %50
%114 = fsub double %93, %113
%115 = fsub double %97, %114
%116 = fsub double %112, %97
%117 = fmul double %102, 1.250000e-01
%118 = fadd double %24, 1.000000e+00
%119 = fmul double %118, %117
%120 = fadd double %24, -2.000000e+00
%121 = fmul double %120, %119
%122 = fadd double %115, %116
%123 = tail call double @llvm.fmuladd.f64(double %121, double %122, double %105)
br label %124
124: ; preds = %91, %22, %15, %110, %11, %5
%125 = phi double [ %6, %5 ], [ %14, %11 ], [ %19, %15 ], [ %29, %22 ], [ %105, %91 ], [ %123, %110 ]
ret double %125
}