#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 fabs(float);
scalar_t__ sqrt(double);
float bessj(int n, float x) {
  float bessj0(float x);
  float bessj1(float x);
  void nrerror(char error_text[]);
  int j, jsum, m;
  float ax, bj, bjm, bjp, sum, tox, ans;
  if (n < 2)
    nrerror("Index n less than 2 in bessj");
  ax = fabs(x);
  if (ax == 0.0)
    return 0.0;
  else if (ax > (float)n) {
    tox = 2.0 / ax;
    bjm = bessj0(ax);
    bj = bessj1(ax);
    for (j = 1; j < n; j++) {
      bjp = j * tox * bj - bjm;
      bjm = bj;
      bj = bjp;
    }
    ans = bj;
  } else {
    tox = 2.0 / ax;
    m = 2 * ((n + (int)sqrt(40.0 * n)) / 2);
    jsum = 0;
    bjp = ans = sum = 0.0;
    bj = 1.0;
    for (j = m; j > 0; j--) {
      bjm = j * tox * bj - bjp;
      bjp = bj;
      bj = bjm;
      if (fabs(bj) > 1.0e10) {
        bj *= 1.0e-10;
        bjp *= 1.0e-10;
        ans *= 1.0e-10;
        sum *= 1.0e-10;
      }
      if (jsum)
        sum += bj;
      jsum = !jsum;
      if (j == n)
        ans = bjp;
    }
    sum = 2.0 * sum - bj;
    ans /= sum;
  }
  return x < 0.0 && (n & 1) ? -ans : ans;
}


x86 assembly (clang, -O0):
<bessj>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
cmpl   $0x2,-0x8(%rbp)
jge    1176 <bessj+0x26>
lea    0xeb7(%rip),%rdi
call   1040 <nrerror@plt>
movss  -0xc(%rbp),%xmm0
call   1050 <fabs@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11b2 <bessj+0x62>
jp     11b2 <bessj+0x62>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    14cf <bessj+0x37f>
movss  -0x1c(%rbp),%xmm0
cvtsi2ssl -0x8(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jbe    1273 <bessj+0x123>
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe3a(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x1c(%rbp),%xmm0
call   1060 <bessj0@plt>
movss  %xmm0,-0x24(%rbp)
movss  -0x1c(%rbp),%xmm0
call   1030 <bessj1@plt>
movss  %xmm0,-0x20(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1264 <bessj+0x114>
cvtsi2ssl -0x10(%rbp),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  -0x20(%rbp),%xmm2
movss  -0x24(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1208 <bessj+0xb8>
movss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
jmp    1475 <bessj+0x325>
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd8c(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x38(%rbp)
cvtsi2sdl -0x8(%rbp),%xmm1
movsd  0xd64(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1070 <sqrt@plt>
mov    %rax,%rcx
mov    -0x38(%rbp),%eax
add    %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
shl    $0x1,%eax
mov    %eax,-0x18(%rbp)
movl   $0x0,-0x14(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  0xd16(%rip),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x10(%rbp)
cmpl   $0x0,-0x10(%rbp)
jle    1428 <bessj+0x2d8>
cvtsi2ssl -0x10(%rbp),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  -0x20(%rbp),%xmm2
movss  -0x28(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
call   1050 <fabs@plt>
movsd  0xcc5(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    13d5 <bessj+0x285>
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xcb2(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc94(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc76(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc58(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
cmpl   $0x0,-0x14(%rbp)
je     13ee <bessj+0x29e>
movss  -0x20(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
cmpl   $0x0,-0x14(%rbp)
setne  %al
xor    $0xff,%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x14(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x8(%rbp),%eax
jne    1415 <bessj+0x2c5>
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
jmp    141a <bessj+0x2ca>
mov    -0x10(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x10(%rbp)
jmp    12f5 <bessj+0x1a5>
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x20(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm1
movss  -0x34(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
jmp    147a <bessj+0x32a>
movss  -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    14bb <bessj+0x36b>
mov    -0x8(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     14bb <bessj+0x36b>
movss  -0x34(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x3c(%rbp)
jmp    14c5 <bessj+0x375>
movss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessj>:
push   %rbp
push   %r14
push   %rbx
sub    $0x40,%rsp
mov    %edi,%r14d
cmp    $0x1,%edi
movaps %xmm0,0x10(%rsp)
jg     1176 <bessj+0x26>
lea    0xef4(%rip),%rdi
call   1040 <nrerror@plt>
movaps 0x10(%rsp),%xmm0
call   1050 <fabs@plt>
cvtsd2ss %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm1
jne    118d <bessj+0x3d>
jnp    144b <bessj+0x2fb>
xorps  %xmm0,%xmm0
cvtsi2ss %r14d,%xmm0
ucomiss %xmm0,%xmm1
movss  0xe60(%rip),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,0x8(%rsp)
jbe    1206 <bessj+0xb6>
movaps %xmm1,%xmm0
movss  %xmm1,0xc(%rsp)
call   1060 <bessj0@plt>
movss  %xmm0,0x20(%rsp)
movss  0xc(%rsp),%xmm0
call   1030 <bessj1@plt>
cmp    $0x2,%r14d
jl     1332 <bessj+0x1e2>
movss  0x20(%rsp),%xmm5
lea    -0x1(%r14),%edx
lea    -0x2(%r14),%ecx
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%ecx
jae    1378 <bessj+0x228>
mov    $0x1,%ecx
movaps 0x10(%rsp),%xmm3
movss  0x8(%rsp),%xmm4
jmp    13f1 <bessj+0x2a1>
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
mulsd  0xdfa(%rip),%xmm0
call   1070 <sqrt@plt>
mov    %rax,%rbx
add    %r14d,%ebx
cmp    $0x2,%ebx
jl     133f <bessj+0x1ef>
and    $0xfffffffe,%ebx
movaps 0xdfc(%rip),%xmm0
xor    %ebp,%ebp
xorps  %xmm6,%xmm6
xorps  %xmm7,%xmm7
movss  0x8(%rsp),%xmm3
jmp    126b <bessj+0x11b>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addss  %xmm2,%xmm6
xor    $0x1,%ebp
lea    -0x1(%rbx),%eax
movaps %xmm1,%xmm0
unpcklps %xmm4,%xmm0
cmp    $0x1,%ebx
mov    %eax,%ebx
jle    1350 <bessj+0x200>
movaps %xmm0,0x30(%rsp)
movss  %xmm6,0xc(%rsp)
xorps  %xmm1,%xmm1
cvtsi2ss %ebx,%xmm1
mulss  %xmm3,%xmm1
shufps $0x55,%xmm0,%xmm0
mulss  %xmm1,%xmm0
subss  %xmm7,%xmm0
movaps %xmm0,0x20(%rsp)
call   1050 <fabs@plt>
ucomisd 0xd79(%rip),%xmm0
jbe    12f0 <bessj+0x1a0>
movaps 0x20(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd6e(%rip),%xmm1
mulsd  %xmm1,%xmm0
xorps  %xmm4,%xmm4
cvtsd2ss %xmm0,%xmm4
cvtps2pd 0x30(%rsp),%xmm0
mulpd  0xd76(%rip),%xmm0
cvtpd2ps %xmm0,%xmm5
movss  0xc(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
xorps  %xmm2,%xmm2
cvtsd2ss %xmm0,%xmm2
movss  0x8(%rsp),%xmm3
jmp    1307 <bessj+0x1b7>
nopl   0x0(%rax,%rax,1)
movss  0xc(%rsp),%xmm2
movss  0x8(%rsp),%xmm3
movaps 0x20(%rsp),%xmm4
movapd 0x30(%rsp),%xmm5
movss  0xcf9(%rip),%xmm6
test   %ebp,%ebp
je     1316 <bessj+0x1c6>
movaps %xmm4,%xmm6
movapd %xmm5,%xmm7
shufps $0x55,%xmm5,%xmm7
movaps %xmm7,%xmm1
cmp    %ebx,%r14d
je     1250 <bessj+0x100>
movaps %xmm5,%xmm1
jmp    1250 <bessj+0x100>
movaps %xmm0,%xmm1
movaps 0x10(%rsp),%xmm3
jmp    1421 <bessj+0x2d1>
xorps  %xmm1,%xmm1
xorpd  %xmm0,%xmm0
movss  0xcb6(%rip),%xmm4
jmp    1357 <bessj+0x207>
xorps  %xmm0,%xmm0
cvtss2sd %xmm6,%xmm0
movaps 0x10(%rsp),%xmm3
xorps  %xmm2,%xmm2
cvtss2sd %xmm4,%xmm2
addsd  %xmm0,%xmm0
subsd  %xmm2,%xmm0
cvtsd2ss %xmm0,%xmm0
divss  %xmm0,%xmm1
jmp    1421 <bessj+0x2d1>
and    $0xfffffffc,%edx
mov    $0x1,%ecx
movss  0x8(%rsp),%xmm4
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
mulss  %xmm4,%xmm1
mulss  %xmm0,%xmm1
lea    0x1(%rcx),%esi
xorps  %xmm2,%xmm2
cvtsi2ss %esi,%xmm2
subss  %xmm5,%xmm1
mulss  %xmm4,%xmm2
mulss  %xmm1,%xmm2
subss  %xmm0,%xmm2
lea    0x2(%rcx),%esi
xorps  %xmm5,%xmm5
cvtsi2ss %esi,%xmm5
mulss  %xmm4,%xmm5
mulss  %xmm2,%xmm5
subss  %xmm1,%xmm5
lea    0x3(%rcx),%esi
xorps  %xmm0,%xmm0
cvtsi2ss %esi,%xmm0
mulss  %xmm4,%xmm0
mulss  %xmm5,%xmm0
subss  %xmm2,%xmm0
add    $0x4,%ecx
cmp    %edx,%esi
jne    1390 <bessj+0x240>
movaps 0x10(%rsp),%xmm3
movaps %xmm0,%xmm1
test   %eax,%eax
je     1421 <bessj+0x2d1>
nopl   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
mulss  %xmm4,%xmm1
mulss  %xmm0,%xmm1
subss  %xmm5,%xmm1
add    $0x1,%ecx
movaps %xmm0,%xmm5
movaps %xmm1,%xmm0
add    $0xffffffff,%eax
jne    1400 <bessj+0x2b0>
movaps %xmm1,%xmm0
test   $0x1,%r14b
je     1434 <bessj+0x2e4>
movaps 0xc1f(%rip),%xmm0
xorps  %xmm1,%xmm0
xorps  %xmm2,%xmm2
cmpnltss %xmm2,%xmm3
movaps %xmm3,%xmm2
andnps %xmm0,%xmm2
andps  %xmm1,%xmm3
orps   %xmm2,%xmm3
movaps %xmm3,%xmm0
add    $0x40,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
@.str = [29 x i8] c"Index n less than 2 in bessj\00"
define float @bessj(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
store i32 %0, i32* %4
store float %1, float* %5
%16 = load i32, i32* %4
%17 = icmp slt i32 %16, 2
br i1 %17, label %18, label %19
18: ; preds = %2
br label %19
19: ; preds = %18, %2
%20 = load float, float* %5
%21 = call double @fabs(float %20)
%22 = fptrunc double %21 to float
store float %22, float* %9
%23 = load float, float* %9
%24 = fpext float %23 to double
%25 = fcmp oeq double %24, 0.000000e+00
br i1 %25, label %26, label %27
26: ; preds = %19
store float 0.000000e+00, float* %3
br label %158
27: ; preds = %19
%28 = load float, float* %9
%29 = load i32, i32* %4
%30 = sitofp i32 %29 to float
%31 = fcmp ogt float %28, %30
br i1 %31, label %32, label %61
32: ; preds = %27
%33 = load float, float* %9
%34 = fpext float %33 to double
%35 = fdiv double 2.000000e+00, %34
%36 = fptrunc double %35 to float
store float %36, float* %14
%37 = load float, float* %9
%38 = call float @bessj0(float %37)
store float %38, float* %11
%39 = load float, float* %9
%40 = call float @bessj1(float %39)
store float %40, float* %10
store i32 1, i32* %6
br label %41
41: ; preds = %56, %32
%42 = load i32, i32* %6
%43 = load i32, i32* %4
%44 = icmp slt i32 %42, %43
br i1 %44, label %45, label %59
45: ; preds = %41
%46 = load i32, i32* %6
%47 = sitofp i32 %46 to float
%48 = load float, float* %14
%49 = fmul float %47, %48
%50 = load float, float* %10
%51 = load float, float* %11
%52 = fneg float %51
%53 = call float @llvm.fmuladd.f32(float %49, float %50, float %52)
store float %53, float* %12
%54 = load float, float* %10
store float %54, float* %11
%55 = load float, float* %12
store float %55, float* %10
br label %56
56: ; preds = %45
%57 = load i32, i32* %6
%58 = add nsw i32 %57, 1
store i32 %58, i32* %6
br label %41
59: ; preds = %41
%60 = load float, float* %10
store float %60, float* %15
br label %142
61: ; preds = %27
%62 = load float, float* %9
%63 = fpext float %62 to double
%64 = fdiv double 2.000000e+00, %63
%65 = fptrunc double %64 to float
store float %65, float* %14
%66 = load i32, i32* %4
%67 = load i32, i32* %4
%68 = sitofp i32 %67 to double
%69 = fmul double 4.000000e+01, %68
%70 = call i64 @sqrt(double %69)
%71 = trunc i64 %70 to i32
%72 = add nsw i32 %66, %71
%73 = sdiv i32 %72, 2
%74 = mul nsw i32 2, %73
store i32 %74, i32* %8
store i32 0, i32* %7
store float 0.000000e+00, float* %13
store float 0.000000e+00, float* %15
store float 0.000000e+00, float* %12
store float 1.000000e+00, float* %10
%75 = load i32, i32* %8
store i32 %75, i32* %6
br label %76
76: ; preds = %128, %61
%77 = load i32, i32* %6
%78 = icmp sgt i32 %77, 0
br i1 %78, label %79, label %131
79: ; preds = %76
%80 = load i32, i32* %6
%81 = sitofp i32 %80 to float
%82 = load float, float* %14
%83 = fmul float %81, %82
%84 = load float, float* %10
%85 = load float, float* %12
%86 = fneg float %85
%87 = call float @llvm.fmuladd.f32(float %83, float %84, float %86)
store float %87, float* %11
%88 = load float, float* %10
store float %88, float* %12
%89 = load float, float* %11
store float %89, float* %10
%90 = load float, float* %10
%91 = call double @fabs(float %90)
%92 = fcmp ogt double %91, 1.000000e+10
br i1 %92, label %93, label %110
93: ; preds = %79
%94 = load float, float* %10
%95 = fpext float %94 to double
%96 = fmul double %95, 1.000000e-10
%97 = fptrunc double %96 to float
store float %97, float* %10
%98 = load float, float* %12
%99 = fpext float %98 to double
%100 = fmul double %99, 1.000000e-10
%101 = fptrunc double %100 to float
store float %101, float* %12
%102 = load float, float* %15
%103 = fpext float %102 to double
%104 = fmul double %103, 1.000000e-10
%105 = fptrunc double %104 to float
store float %105, float* %15
%106 = load float, float* %13
%107 = fpext float %106 to double
%108 = fmul double %107, 1.000000e-10
%109 = fptrunc double %108 to float
store float %109, float* %13
br label %110
110: ; preds = %93, %79
%111 = load i32, i32* %7
%112 = icmp ne i32 %111, 0
br i1 %112, label %113, label %117
113: ; preds = %110
%114 = load float, float* %10
%115 = load float, float* %13
%116 = fadd float %115, %114
store float %116, float* %13
br label %117
117: ; preds = %113, %110
%118 = load i32, i32* %7
%119 = icmp ne i32 %118, 0
%120 = xor i1 %119, true
%121 = zext i1 %120 to i32
store i32 %121, i32* %7
%122 = load i32, i32* %6
%123 = load i32, i32* %4
%124 = icmp eq i32 %122, %123
br i1 %124, label %125, label %127
125: ; preds = %117
%126 = load float, float* %12
store float %126, float* %15
br label %127
127: ; preds = %125, %117
br label %128
128: ; preds = %127
%129 = load i32, i32* %6
%130 = add nsw i32 %129, -1
store i32 %130, i32* %6
br label %76
131: ; preds = %76
%132 = load float, float* %13
%133 = fpext float %132 to double
%134 = load float, float* %10
%135 = fpext float %134 to double
%136 = fneg double %135
%137 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %133, double %136)
%138 = fptrunc double %137 to float
store float %138, float* %13
%139 = load float, float* %13
%140 = load float, float* %15
%141 = fdiv float %140, %139
store float %141, float* %15
br label %142
142: ; preds = %131, %59
br label %143
143: ; preds = %142
%144 = load float, float* %5
%145 = fpext float %144 to double
%146 = fcmp olt double %145, 0.000000e+00
br i1 %146, label %147, label %154
147: ; preds = %143
%148 = load i32, i32* %4
%149 = and i32 %148, 1
%150 = icmp ne i32 %149, 0
br i1 %150, label %151, label %154
151: ; preds = %147
%152 = load float, float* %15
%153 = fneg float %152
br label %156
154: ; preds = %147, %143
%155 = load float, float* %15
br label %156
156: ; preds = %154, %151
%157 = phi float [ %153, %151 ], [ %155, %154 ]
store float %157, float* %3
br label %158
158: ; preds = %156, %26
%159 = load float, float* %3
ret float %159
}

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