double owl_stats_mean_stable(double *x, int n) {
  double t = 0;
  for (int i = 0; i < n; i++)
    t += (*(x + i) - t) / (i + 1);
  return t;
}


x86 assembly (clang, -O0):
<owl_stats_mean_stable>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    115e <owl_stats_mean_stable+0x5e>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x18(%rbp),%xmm0
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    111a <owl_stats_mean_stable+0x1a>
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<owl_stats_mean_stable>:
test   %esi,%esi
jle    1113 <owl_stats_mean_stable+0x13>
mov    %esi,%eax
cmp    $0x1,%esi
jne    1117 <owl_stats_mean_stable+0x17>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    1171 <owl_stats_mean_stable+0x71>
xorps  %xmm0,%xmm0
ret
mov    %eax,%edx
and    $0xfffffffe,%edx
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rdi,%rcx,8),%xmm2
movsd  0x8(%rdi,%rcx,8),%xmm0
subsd  %xmm1,%xmm2
lea    0x1(%rcx),%esi
xorps  %xmm3,%xmm3
cvtsi2sd %esi,%xmm3
divsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
add    $0x2,%rcx
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
subsd  %xmm2,%xmm0
divsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
movapd %xmm0,%xmm1
cmp    %rcx,%rdx
jne    1130 <owl_stats_mean_stable+0x30>
test   $0x1,%al
je     1190 <owl_stats_mean_stable+0x90>
movsd  (%rdi,%rcx,8),%xmm1
subsd  %xmm0,%xmm1
add    $0x1,%ecx
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
ret

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

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