float bessel(float x) {
  float xh, sum, pow, ds;
  int k;
  xh = (float)0.5 * x;
  sum = 1.0;
  pow = 1.0;
  k = 0;
  ds = 1.0;
  while (ds > sum * 1E-16) {
    ++k;
    pow = pow * (xh / k);
    ds = pow * pow;
    sum = sum + ds;
  }
  return sum;
}


x86 assembly (clang, -O0):
<bessel>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
movss  0xef3(%rip),%xmm0
mulss  -0x4(%rbp),%xmm0
movss  %xmm0,-0x8(%rbp)
movss  0xedd(%rip),%xmm0
movss  %xmm0,-0xc(%rbp)
movss  0xed0(%rip),%xmm0
movss  %xmm0,-0x10(%rbp)
movl   $0x0,-0x18(%rbp)
movss  0xebc(%rip),%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0xc(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xea5(%rip),%xmm2
mulsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    11b9 <bessel+0xb9>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
movss  -0x10(%rbp),%xmm0
movss  -0x8(%rbp),%xmm1
cvtsi2ssl -0x18(%rbp),%xmm2
divss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0xc(%rbp),%xmm0
addss  -0x14(%rbp),%xmm0
movss  %xmm0,-0xc(%rbp)
jmp    1149 <bessel+0x49>
movss  -0xc(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessel>:
mulss  0xef8(%rip),%xmm0
movss  0xef4(%rip),%xmm2
mov    $0x1,%eax
movsd  0xeeb(%rip),%xmm3
movaps %xmm2,%xmm1
xorps  %xmm4,%xmm4
cvtsi2ss %eax,%xmm4
movaps %xmm0,%xmm5
divss  %xmm4,%xmm5
mulss  %xmm5,%xmm2
movaps %xmm2,%xmm4
mulss  %xmm2,%xmm4
addss  %xmm4,%xmm1
cvtss2sd %xmm4,%xmm4
xorps  %xmm5,%xmm5
cvtss2sd %xmm1,%xmm5
mulsd  %xmm3,%xmm5
add    $0x1,%eax
ucomisd %xmm5,%xmm4
ja     1120 <bessel+0x20>
movaps %xmm1,%xmm0
ret

LLVM IR (-O0):
define float @bessel(float %0) {
%2 = alloca float
%3 = alloca float
%4 = alloca float
%5 = alloca float
%6 = alloca float
%7 = alloca i32
store float %0, float* %2
%8 = load float, float* %2
%9 = fmul float 5.000000e-01, %8
store float %9, float* %3
store float 1.000000e+00, float* %4
store float 1.000000e+00, float* %5
store i32 0, i32* %7
store float 1.000000e+00, float* %6
br label %10
10: ; preds = %17, %1
%11 = load float, float* %6
%12 = fpext float %11 to double
%13 = load float, float* %4
%14 = fpext float %13 to double
%15 = fmul double %14, 0x3C9CD2B297D889BC
%16 = fcmp ogt double %12, %15
br i1 %16, label %17, label %32
17: ; preds = %10
%18 = load i32, i32* %7
%19 = add nsw i32 %18, 1
store i32 %19, i32* %7
%20 = load float, float* %5
%21 = load float, float* %3
%22 = load i32, i32* %7
%23 = sitofp i32 %22 to float
%24 = fdiv float %21, %23
%25 = fmul float %20, %24
store float %25, float* %5
%26 = load float, float* %5
%27 = load float, float* %5
%28 = fmul float %26, %27
store float %28, float* %6
%29 = load float, float* %4
%30 = load float, float* %6
%31 = fadd float %29, %30
store float %31, float* %4
br label %10
32: ; preds = %10
%33 = load float, float* %4
ret float %33
}

LLVM IR (-O2):
define float @bessel(float %0) {
%2 = fmul float %0, 5.000000e-01
br label %3
3: ; preds = %1, %3
%4 = phi i32 [ 0, %1 ], [ %7, %3 ]
%5 = phi float [ 1.000000e+00, %1 ], [ %10, %3 ]
%6 = phi float [ 1.000000e+00, %1 ], [ %12, %3 ]
%7 = add nuw nsw i32 %4, 1
%8 = sitofp i32 %7 to float
%9 = fdiv float %2, %8
%10 = fmul float %5, %9
%11 = fmul float %10, %10
%12 = fadd float %6, %11
%13 = fpext float %11 to double
%14 = fpext float %12 to double
%15 = fmul double %14, 0x3C9CD2B297D889BC
%16 = fcmp olt double %15, %13
br i1 %16, label %3, label %17
17: ; preds = %3
ret float %12
}