float bessel_Kn_scaled_small_x(int n) {
  int k;
  float k_term, sum1;
  for (k = 1; k <= n - 1; k++) {
    k_term *= -1 / (k * (n - k));
    sum1 += k_term;
  }
  return sum1;
}


x86 assembly (clang, -O0):
<bessel_Kn_scaled_small_x>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movl   $0x1,-0x8(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     115e <bessel_Kn_scaled_small_x+0x5e>
mov    -0x8(%rbp),%ecx
mov    -0x4(%rbp),%eax
sub    -0x8(%rbp),%eax
imul   %eax,%ecx
mov    $0xffffffff,%eax
cltd
idiv   %ecx
cvtsi2ss %eax,%xmm0
mulss  -0xc(%rbp),%xmm0
movss  %xmm0,-0xc(%rbp)
movss  -0xc(%rbp),%xmm0
addss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    110e <bessel_Kn_scaled_small_x+0xe>
movss  -0x10(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<bessel_Kn_scaled_small_x>:
cmp    $0x2,%edi
jl     11a6 <bessel_Kn_scaled_small_x+0xa6>
lea    -0x1(%rdi),%r8d
cmp    $0x2,%edi
jne    1119 <bessel_Kn_scaled_small_x+0x19>
mov    $0x1,%eax
jmp    1184 <bessel_Kn_scaled_small_x+0x84>
lea    -0x2(%rdi),%esi
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
neg    %r9d
mov    $0x1,%r10d
nopl   0x0(%rax)
lea    0x1(%rsi),%ecx
imul   %r10d,%ecx
mov    $0xffffffff,%eax
cltd
idiv   %ecx
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
mulss  %xmm1,%xmm2
lea    0x1(%r10),%ecx
imul   %esi,%ecx
mov    $0xffffffff,%eax
cltd
idiv   %ecx
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
addss  %xmm2,%xmm0
mulss  %xmm2,%xmm1
addss  %xmm1,%xmm0
add    $0xfffffffe,%esi
mov    %r10d,%eax
add    $0x2,%eax
lea    (%r9,%r10,1),%ecx
add    $0x2,%ecx
mov    %eax,%r10d
cmp    $0x1,%ecx
jne    1130 <bessel_Kn_scaled_small_x+0x30>
test   $0x1,%r8b
je     11a6 <bessel_Kn_scaled_small_x+0xa6>
sub    %eax,%edi
imul   %eax,%edi
mov    $0xffffffff,%eax
cltd
idiv   %edi
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
mulss  %xmm2,%xmm1
addss  %xmm1,%xmm0
ret

LLVM IR (-O0):
define float @bessel_Kn_scaled_small_x(i32 %0) {
%2 = alloca i32
%3 = alloca i32
%4 = alloca float
%5 = alloca float
store i32 %0, i32* %2
store i32 1, i32* %3
br label %6
6: ; preds = %24, %1
%7 = load i32, i32* %3
%8 = load i32, i32* %2
%9 = sub nsw i32 %8, 1
%10 = icmp sle i32 %7, %9
br i1 %10, label %11, label %27
11: ; preds = %6
%12 = load i32, i32* %3
%13 = load i32, i32* %2
%14 = load i32, i32* %3
%15 = sub nsw i32 %13, %14
%16 = mul nsw i32 %12, %15
%17 = sdiv i32 -1, %16
%18 = sitofp i32 %17 to float
%19 = load float, float* %4
%20 = fmul float %19, %18
store float %20, float* %4
%21 = load float, float* %4
%22 = load float, float* %5
%23 = fadd float %22, %21
store float %23, float* %5
br label %24
24: ; preds = %11
%25 = load i32, i32* %3
%26 = add nsw i32 %25, 1
store i32 %26, i32* %3
br label %6
27: ; preds = %6
%28 = load float, float* %5
ret float %28
}

LLVM IR (-O2):
define float @bessel_Kn_scaled_small_x(i32 %0) {
%2 = icmp sgt i32 %0, 1
br i1 %2, label %3, label %42
3: ; preds = %1
%4 = add i32 %0, -1
%5 = and i32 %4, 1
%6 = icmp eq i32 %0, 2
br i1 %6, label %30, label %7
7: ; preds = %3
%8 = and i32 %4, -2
br label %9
9: ; preds = %9, %7
%10 = phi float [ undef, %7 ], [ %26, %9 ]
%11 = phi float [ undef, %7 ], [ %25, %9 ]
%12 = phi i32 [ 1, %7 ], [ %27, %9 ]
%13 = phi i32 [ 0, %7 ], [ %28, %9 ]
%14 = sub nsw i32 %0, %12
%15 = mul nsw i32 %14, %12
%16 = sdiv i32 -1, %15
%17 = sitofp i32 %16 to float
%18 = fmul float %11, %17
%19 = fadd float %10, %18
%20 = add nuw nsw i32 %12, 1
%21 = sub nsw i32 %0, %20
%22 = mul nsw i32 %21, %20
%23 = sdiv i32 -1, %22
%24 = sitofp i32 %23 to float
%25 = fmul float %18, %24
%26 = fadd float %19, %25
%27 = add nuw nsw i32 %12, 2
%28 = add nuw i32 %13, 2
%29 = icmp eq i32 %28, %8
br i1 %29, label %30, label %9
30: ; preds = %9, %3
%31 = phi float [ undef, %3 ], [ %26, %9 ]
%32 = phi float [ undef, %3 ], [ %25, %9 ]
%33 = phi i32 [ 1, %3 ], [ %27, %9 ]
%34 = icmp eq i32 %5, 0
br i1 %34, label %42, label %35
35: ; preds = %30
%36 = sub nsw i32 %0, %33
%37 = mul nsw i32 %36, %33
%38 = sdiv i32 -1, %37
%39 = sitofp i32 %38 to float
%40 = fmul float %32, %39
%41 = fadd float %31, %40
br label %42
42: ; preds = %35, %30, %1
%43 = phi float [ undef, %1 ], [ %31, %30 ], [ %41, %35 ]
ret float %43
}