double bessel0(double x) {
  int i;
  double d, y, d2, sum;
  y = x / 2.0;
  d = 1.0;
  sum = 1.0;
  for (i = 1; i <= 25; i++) {
    d = d * y / i;
    d2 = d * d;
    sum = sum + d2;
    if (d2 < sum * (1.0e-8))
      break;
  }
  return (sum);
}


x86 assembly (clang, -O0):
<bessel0>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  0xef2(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  0xed9(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  0xecc(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0xc(%rbp)
cmpl   $0x19,-0xc(%rbp)
jg     11b4 <bessel0+0xb4>
movsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  0xe83(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11a1 <bessel0+0xa1>
jmp    11b4 <bessel0+0xb4>
jmp    11a6 <bessel0+0xa6>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    1140 <bessel0+0x40>
movsd  -0x30(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessel0>:
mulsd  0xef8(%rip),%xmm0
mov    $0x1,%eax
movsd  0xef3(%rip),%xmm1
movsd  0xef3(%rip),%xmm2
movapd %xmm1,%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mulsd  %xmm0,%xmm3
xorps  %xmm4,%xmm4
cvtsi2sd %eax,%xmm4
divsd  %xmm4,%xmm3
movapd %xmm3,%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm4,%xmm1
movapd %xmm1,%xmm5
mulsd  %xmm2,%xmm5
ucomisd %xmm4,%xmm5
ja     1163 <bessel0+0x63>
lea    0x1(%rax),%ecx
cmp    $0x19,%eax
mov    %ecx,%eax
jne    1130 <bessel0+0x30>
movapd %xmm1,%xmm0
ret

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

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