void pitch_yin_getcum(float *diff, int N) {
  int tau;
  float tmp;
  tmp = 0.;
  diff[0] = 1.;
  for (tau = 1; tau < N; tau++) {
    tmp += diff[tau];
    diff[tau] *= tau / tmp;
  }
}


x86 assembly (clang, -O0):
<pitch_yin_getcum>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
mov    -0x8(%rbp),%rax
movss  0xee1(%rip),%xmm0
movss  %xmm0,(%rax)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1177 <pitch_yin_getcum+0x77>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
cvtsi2ssl -0x10(%rbp),%xmm0
divss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    112a <pitch_yin_getcum+0x2a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<pitch_yin_getcum>:
movl   $0x3f800000,(%rdi)
cmp    $0x2,%esi
jl     11ab <pitch_yin_getcum+0xab>
mov    %esi,%ecx
add    $0xffffffffffffffff,%rcx
cmp    $0x2,%esi
jne    1124 <pitch_yin_getcum+0x24>
xorps  %xmm0,%xmm0
mov    $0x1,%eax
jmp    1189 <pitch_yin_getcum+0x89>
mov    %rcx,%rdx
and    $0xfffffffffffffffe,%rdx
xorps  %xmm0,%xmm0
mov    $0x1,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%rdi,%rax,4),%xmm1
movss  0x4(%rdi,%rax,4),%xmm2
addss  %xmm1,%xmm0
xorps  %xmm3,%xmm3
cvtsi2ss %eax,%xmm3
divss  %xmm0,%xmm3
mulss  %xmm1,%xmm3
movss  %xmm3,(%rdi,%rax,4)
addss  %xmm2,%xmm0
lea    0x1(%rax),%rsi
xorps  %xmm1,%xmm1
cvtsi2ss %esi,%xmm1
divss  %xmm0,%xmm1
mulss  %xmm2,%xmm1
movss  %xmm1,0x4(%rdi,%rax,4)
add    $0x2,%rax
cmp    %rdx,%rsi
jne    1140 <pitch_yin_getcum+0x40>
test   $0x1,%cl
je     11ab <pitch_yin_getcum+0xab>
movss  (%rdi,%rax,4),%xmm1
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
addss  %xmm1,%xmm0
divss  %xmm0,%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,(%rdi,%rax,4)
ret

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

LLVM IR (-O2):
define void @pitch_yin_getcum(float* nocapture %0, i32 %1) {
store float 1.000000e+00, float* %0
%3 = icmp sgt i32 %1, 1
br i1 %3, label %4, label %45
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = add nsw i64 %5, -1
%7 = and i64 %6, 1
%8 = icmp eq i32 %1, 2
br i1 %8, label %33, label %9
9: ; preds = %4
%10 = and i64 %6, -2
br label %11
11: ; preds = %11, %9
%12 = phi i64 [ 1, %9 ], [ %30, %11 ]
%13 = phi float [ 0.000000e+00, %9 ], [ %25, %11 ]
%14 = phi i64 [ 0, %9 ], [ %31, %11 ]
%15 = getelementptr inbounds float, float* %0, i64 %12
%16 = load float, float* %15
%17 = fadd float %13, %16
%18 = trunc i64 %12 to i32
%19 = sitofp i32 %18 to float
%20 = fdiv float %19, %17
%21 = fmul float %16, %20
store float %21, float* %15
%22 = add nuw nsw i64 %12, 1
%23 = getelementptr inbounds float, float* %0, i64 %22
%24 = load float, float* %23
%25 = fadd float %17, %24
%26 = trunc i64 %22 to i32
%27 = sitofp i32 %26 to float
%28 = fdiv float %27, %25
%29 = fmul float %24, %28
store float %29, float* %23
%30 = add nuw nsw i64 %12, 2
%31 = add i64 %14, 2
%32 = icmp eq i64 %31, %10
br i1 %32, label %33, label %11
33: ; preds = %11, %4
%34 = phi i64 [ 1, %4 ], [ %30, %11 ]
%35 = phi float [ 0.000000e+00, %4 ], [ %25, %11 ]
%36 = icmp eq i64 %7, 0
br i1 %36, label %45, label %37
37: ; preds = %33
%38 = getelementptr inbounds float, float* %0, i64 %34
%39 = load float, float* %38
%40 = fadd float %35, %39
%41 = trunc i64 %34 to i32
%42 = sitofp i32 %41 to float
%43 = fdiv float %42, %40
%44 = fmul float %39, %43
store float %44, float* %38
br label %45
45: ; preds = %37, %33, %2
ret void
}