float zero_bessel(float x) {
  float num, fact, i;
  float result = 1;
  float x_2 = x / 2;
  num = 1;
  fact = 1;
  for (i = 1; i < 20; i++) {
    num *= x_2 * x_2;
    fact *= i;
    result += num / (fact * fact);
  }
  return result;
}


x86 assembly (clang, -O0):
<zero_bessel>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
movss  0xeef(%rip),%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x4(%rbp),%xmm0
movss  0xee1(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  0xecc(%rip),%xmm0
movss  %xmm0,-0x8(%rbp)
movss  0xebf(%rip),%xmm0
movss  %xmm0,-0xc(%rbp)
movss  0xeb2(%rip),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  0xead(%rip),%xmm0
ucomiss -0x10(%rbp),%xmm0
jbe    11bc <zero_bessel+0xbc>
movss  -0x18(%rbp),%xmm0
mulss  -0x18(%rbp),%xmm0
mulss  -0x8(%rbp),%xmm0
movss  %xmm0,-0x8(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0xc(%rbp),%xmm0
movss  %xmm0,-0xc(%rbp)
movss  -0x8(%rbp),%xmm0
movss  -0xc(%rbp),%xmm1
mulss  -0xc(%rbp),%xmm1
divss  %xmm1,%xmm0
addss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
movss  0xe53(%rip),%xmm0
addss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
jmp    1153 <zero_bessel+0x53>
movss  -0x14(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<zero_bessel>:
mulss  0xef8(%rip),%xmm0
mulss  %xmm0,%xmm0
mov    $0x1,%eax
movss  0xeeb(%rip),%xmm2
movaps %xmm2,%xmm4
movaps %xmm2,%xmm1
nop
xorps  %xmm3,%xmm3
cvtsi2ss %eax,%xmm3
mulss  %xmm1,%xmm3
mulss  %xmm0,%xmm2
movaps %xmm3,%xmm5
mulss  %xmm3,%xmm5
movaps %xmm2,%xmm1
divss  %xmm5,%xmm1
addss  %xmm4,%xmm1
add    $0x1,%eax
cmp    $0x14,%eax
je     1175 <zero_bessel+0x75>
xorps  %xmm4,%xmm4
cvtsi2ss %eax,%xmm4
mulss  %xmm0,%xmm2
mulss  %xmm4,%xmm3
movaps %xmm3,%xmm4
mulss  %xmm3,%xmm4
movaps %xmm2,%xmm5
divss  %xmm4,%xmm5
addss  %xmm5,%xmm1
add    $0x1,%eax
movaps %xmm1,%xmm4
movaps %xmm3,%xmm1
jmp    1120 <zero_bessel+0x20>
movaps %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define float @zero_bessel(float %0) {
%2 = fmul float %0, 5.000000e-01
%3 = fmul float %2, %2
br label %4
4: ; preds = %17, %1
%5 = phi float [ 1.000000e+00, %1 ], [ %19, %17 ]
%6 = phi float [ 1.000000e+00, %1 ], [ %23, %17 ]
%7 = phi i32 [ 1, %1 ], [ %24, %17 ]
%8 = phi float [ 1.000000e+00, %1 ], [ %20, %17 ]
%9 = sitofp i32 %7 to float
%10 = fmul float %3, %5
%11 = fmul float %8, %9
%12 = fmul float %11, %11
%13 = fdiv float %10, %12
%14 = fadd float %6, %13
%15 = add nuw nsw i32 %7, 1
%16 = icmp eq i32 %15, 20
br i1 %16, label %25, label %17
17: ; preds = %4
%18 = sitofp i32 %15 to float
%19 = fmul float %3, %10
%20 = fmul float %11, %18
%21 = fmul float %20, %20
%22 = fdiv float %19, %21
%23 = fadd float %14, %22
%24 = add nuw nsw i32 %7, 2
br label %4
25: ; preds = %4
ret float %14
}