float bessk(int n, float x) {
  float bessk0(float x);
  float bessk1(float x);
  void nrerror(char error_text[]);
  int j;
  float bk, bkm, bkp, tox;
  if (n < 2)
    nrerror("Index n less than 2 in bessk");
  tox = 2.0 / x;
  bkm = bessk0(x);
  bk = bessk1(x);
  for (j = 1; j < n; j++) {
    bkp = bkm + j * tox * bk;
    bkm = bk;
    bk = bkp;
  }
  return bk;
}


x86 assembly (clang, -O0):
<bessk>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x4(%rbp)
movss  %xmm0,-0x8(%rbp)
cmpl   $0x2,-0x4(%rbp)
jge    1156 <bessk+0x26>
lea    0xeb7(%rip),%rdi
call   1030 <nrerror@plt>
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe99(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x8(%rbp),%xmm0
call   1040 <bessk0@plt>
movss  %xmm0,-0x14(%rbp)
movss  -0x8(%rbp),%xmm0
call   1050 <bessk1@plt>
movss  %xmm0,-0x10(%rbp)
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11e8 <bessk+0xb8>
movss  -0x14(%rbp),%xmm1
cvtsi2ssl -0xc(%rbp),%xmm0
mulss  -0x1c(%rbp),%xmm0
movss  -0x10(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    1199 <bessk+0x69>
movss  -0x10(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessk>:
push   %rbx
sub    $0x10,%rsp
mov    %edi,%ebx
cmp    $0x1,%edi
movss  %xmm0,0x8(%rsp)
jg     1154 <bessk+0x24>
lea    0xebb(%rip),%rdi
call   1030 <nrerror@plt>
movss  0x8(%rsp),%xmm0
call   1040 <bessk0@plt>
movss  %xmm0,0xc(%rsp)
movss  0x8(%rsp),%xmm0
call   1050 <bessk1@plt>
cmp    $0x2,%ebx
jl     1231 <bessk+0x101>
movss  0xe85(%rip),%xmm1
divss  0x8(%rsp),%xmm1
lea    -0x1(%rbx),%edx
add    $0xfffffffe,%ebx
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%ebx
movss  0xc(%rsp),%xmm4
jae    119e <bessk+0x6e>
mov    $0x1,%ecx
jmp    120c <bessk+0xdc>
and    $0xfffffffc,%edx
mov    $0x1,%ecx
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm2,%xmm2
cvtsi2ss %ecx,%xmm2
mulss  %xmm1,%xmm2
mulss  %xmm0,%xmm2
lea    0x1(%rcx),%esi
xorps  %xmm3,%xmm3
cvtsi2ss %esi,%xmm3
addss  %xmm4,%xmm2
mulss  %xmm1,%xmm3
mulss  %xmm2,%xmm3
addss  %xmm0,%xmm3
lea    0x2(%rcx),%esi
xorps  %xmm4,%xmm4
cvtsi2ss %esi,%xmm4
mulss  %xmm1,%xmm4
mulss  %xmm3,%xmm4
addss  %xmm2,%xmm4
lea    0x3(%rcx),%esi
xorps  %xmm0,%xmm0
cvtsi2ss %esi,%xmm0
mulss  %xmm1,%xmm0
mulss  %xmm4,%xmm0
addss  %xmm3,%xmm0
add    $0x4,%ecx
cmp    %edx,%esi
jne    11b0 <bessk+0x80>
test   %eax,%eax
je     1231 <bessk+0x101>
movaps %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
mulss  %xmm1,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm4,%xmm0
add    $0x1,%ecx
movaps %xmm2,%xmm4
add    $0xffffffff,%eax
jne    1210 <bessk+0xe0>
add    $0x10,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@.str = [29 x i8] c"Index n less than 2 in bessk\00"
define float @bessk(i32 %0, float %1) {
%3 = alloca i32
%4 = alloca float
%5 = alloca i32
%6 = alloca float
%7 = alloca float
%8 = alloca float
%9 = alloca float
store i32 %0, i32* %3
store float %1, float* %4
%10 = load i32, i32* %3
%11 = icmp slt i32 %10, 2
br i1 %11, label %12, label %13
12: ; preds = %2
br label %13
13: ; preds = %12, %2
%14 = load float, float* %4
%15 = fpext float %14 to double
%16 = fdiv double 2.000000e+00, %15
%17 = fptrunc double %16 to float
store float %17, float* %9
%18 = load float, float* %4
%19 = call float @bessk0(float %18)
store float %19, float* %7
%20 = load float, float* %4
%21 = call float @bessk1(float %20)
store float %21, float* %6
store i32 1, i32* %5
br label %22
22: ; preds = %36, %13
%23 = load i32, i32* %5
%24 = load i32, i32* %3
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %39
26: ; preds = %22
%27 = load float, float* %7
%28 = load i32, i32* %5
%29 = sitofp i32 %28 to float
%30 = load float, float* %9
%31 = fmul float %29, %30
%32 = load float, float* %6
%33 = call float @llvm.fmuladd.f32(float %31, float %32, float %27)
store float %33, float* %8
%34 = load float, float* %6
store float %34, float* %7
%35 = load float, float* %8
store float %35, float* %6
br label %36
36: ; preds = %26
%37 = load i32, i32* %5
%38 = add nsw i32 %37, 1
store i32 %38, i32* %5
br label %22
39: ; preds = %22
%40 = load float, float* %6
ret float %40
}

LLVM IR (-O2):
@.str = [29 x i8] c"Index n less than 2 in bessk\00"
define float @bessk(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 = fdiv float 2.000000e+00, %1
%7 = tail call float @bessk0(float %1)
%8 = tail call float @bessk1(float %1)
%9 = icmp sgt i32 %0, 1
br i1 %9, label %10, label %57
10: ; preds = %5
%11 = add i32 %0, -1
%12 = add i32 %0, -2
%13 = and i32 %11, 3
%14 = icmp ult i32 %12, 3
br i1 %14, label %40, label %15
15: ; preds = %10
%16 = and i32 %11, -4
br label %17
17: ; preds = %17, %15
%18 = phi i32 [ 1, %15 ], [ %37, %17 ]
%19 = phi float [ %8, %15 ], [ %36, %17 ]
%20 = phi float [ %7, %15 ], [ %32, %17 ]
%21 = phi i32 [ 0, %15 ], [ %38, %17 ]
%22 = sitofp i32 %18 to float
%23 = fmul float %6, %22
%24 = tail call float @llvm.fmuladd.f32(float %23, float %19, float %20)
%25 = add nuw nsw i32 %18, 1
%26 = sitofp i32 %25 to float
%27 = fmul float %6, %26
%28 = tail call float @llvm.fmuladd.f32(float %27, float %24, float %19)
%29 = add nuw nsw i32 %18, 2
%30 = sitofp i32 %29 to float
%31 = fmul float %6, %30
%32 = tail call float @llvm.fmuladd.f32(float %31, float %28, float %24)
%33 = add nuw nsw i32 %18, 3
%34 = sitofp i32 %33 to float
%35 = fmul float %6, %34
%36 = tail call float @llvm.fmuladd.f32(float %35, float %32, float %28)
%37 = add nuw nsw i32 %18, 4
%38 = add i32 %21, 4
%39 = icmp eq i32 %38, %16
br i1 %39, label %40, label %17
40: ; preds = %17, %10
%41 = phi float [ undef, %10 ], [ %36, %17 ]
%42 = phi i32 [ 1, %10 ], [ %37, %17 ]
%43 = phi float [ %8, %10 ], [ %36, %17 ]
%44 = phi float [ %7, %10 ], [ %32, %17 ]
%45 = icmp eq i32 %13, 0
br i1 %45, label %57, label %46
46: ; preds = %40, %46
%47 = phi i32 [ %54, %46 ], [ %42, %40 ]
%48 = phi float [ %53, %46 ], [ %43, %40 ]
%49 = phi float [ %48, %46 ], [ %44, %40 ]
%50 = phi i32 [ %55, %46 ], [ 0, %40 ]
%51 = sitofp i32 %47 to float
%52 = fmul float %6, %51
%53 = tail call float @llvm.fmuladd.f32(float %52, float %48, float %49)
%54 = add nuw nsw i32 %47, 1
%55 = add i32 %50, 1
%56 = icmp eq i32 %55, %13
br i1 %56, label %57, label %46
57: ; preds = %40, %46, %5
%58 = phi float [ %8, %5 ], [ %41, %40 ], [ %53, %46 ]
ret float %58
}