float bessy(int n, float x) {
  float bessy0(float x);
  float bessy1(float x);
  void nrerror(char error_text[]);
  int j;
  float by, bym, byp, tox;
  if (n < 2)
    nrerror("Index n less than 2 in bessy");
  tox = 2.0 / x;
  by = bessy1(x);
  bym = bessy0(x);
  for (j = 1; j < n; j++) {
    byp = j * tox * by - bym;
    bym = by;
    by = byp;
  }
  return by;
}


x86 assembly (clang, -O0):
<bessy>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x4(%rbp)
movss  %xmm0,-0x8(%rbp)
cmpl   $0x2,-0x4(%rbp)
jge    1156 <bessy+0x26>
lea    0xeb7(%rip),%rdi
call   1040 <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   1030 <bessy1@plt>
movss  %xmm0,-0x10(%rbp)
movss  -0x8(%rbp),%xmm0
call   1050 <bessy0@plt>
movss  %xmm0,-0x14(%rbp)
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11f5 <bessy+0xc5>
cvtsi2ssl -0xc(%rbp),%xmm0
mulss  -0x1c(%rbp),%xmm0
movss  -0x10(%rbp),%xmm2
movss  -0x14(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
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 <bessy+0x69>
movss  -0x10(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessy>:
push   %rbx
sub    $0x10,%rsp
mov    %edi,%ebx
cmp    $0x1,%edi
movss  %xmm0,0x8(%rsp)
jg     1154 <bessy+0x24>
lea    0xebb(%rip),%rdi
call   1040 <nrerror@plt>
movss  0x8(%rsp),%xmm0
call   1030 <bessy1@plt>
movss  %xmm0,0xc(%rsp)
movss  0x8(%rsp),%xmm0
call   1050 <bessy0@plt>
cmp    $0x2,%ebx
jl     119a <bessy+0x6a>
movss  0xe89(%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    11a5 <bessy+0x75>
mov    $0x1,%ecx
jmp    120c <bessy+0xdc>
movss  0xc(%rsp),%xmm4
jmp    1231 <bessy+0x101>
and    $0xfffffffc,%edx
mov    $0x1,%ecx
nopl   (%rax)
xorps  %xmm2,%xmm2
cvtsi2ss %ecx,%xmm2
mulss  %xmm1,%xmm2
mulss  %xmm4,%xmm2
lea    0x1(%rcx),%esi
xorps  %xmm3,%xmm3
cvtsi2ss %esi,%xmm3
subss  %xmm0,%xmm2
mulss  %xmm1,%xmm3
mulss  %xmm2,%xmm3
subss  %xmm4,%xmm3
lea    0x2(%rcx),%esi
xorps  %xmm0,%xmm0
cvtsi2ss %esi,%xmm0
mulss  %xmm1,%xmm0
mulss  %xmm3,%xmm0
subss  %xmm2,%xmm0
lea    0x3(%rcx),%esi
xorps  %xmm4,%xmm4
cvtsi2ss %esi,%xmm4
mulss  %xmm1,%xmm4
mulss  %xmm0,%xmm4
subss  %xmm3,%xmm4
add    $0x4,%ecx
cmp    %edx,%esi
jne    11b0 <bessy+0x80>
test   %eax,%eax
je     1231 <bessy+0x101>
movaps %xmm4,%xmm2
xorps  %xmm4,%xmm4
cvtsi2ss %ecx,%xmm4
mulss  %xmm1,%xmm4
mulss  %xmm2,%xmm4
subss  %xmm0,%xmm4
add    $0x1,%ecx
movaps %xmm2,%xmm0
add    $0xffffffff,%eax
jne    1210 <bessy+0xe0>
movaps %xmm4,%xmm0
add    $0x10,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@.str = [29 x i8] c"Index n less than 2 in bessy\00"
define float @bessy(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 @bessy1(float %18)
store float %19, float* %6
%20 = load float, float* %4
%21 = call float @bessy0(float %20)
store float %21, float* %7
store i32 1, i32* %5
br label %22
22: ; preds = %37, %13
%23 = load i32, i32* %5
%24 = load i32, i32* %3
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %40
26: ; preds = %22
%27 = load i32, i32* %5
%28 = sitofp i32 %27 to float
%29 = load float, float* %9
%30 = fmul float %28, %29
%31 = load float, float* %6
%32 = load float, float* %7
%33 = fneg float %32
%34 = call float @llvm.fmuladd.f32(float %30, float %31, float %33)
store float %34, float* %8
%35 = load float, float* %6
store float %35, float* %7
%36 = load float, float* %8
store float %36, float* %6
br label %37
37: ; preds = %26
%38 = load i32, i32* %5
%39 = add nsw i32 %38, 1
store i32 %39, i32* %5
br label %22
40: ; preds = %22
%41 = load float, float* %6
ret float %41
}

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