int MeanVar(double x[], int n, double *m, double *v) {
  int i;
  for (i = 0, *m = 0; i < n; i++)
    *m += x[i];
  *m /= n;
  for (i = 0, *v = 0; i < n; i++)
    *v += ((x[i] - *m) * (x[i] - *m));
  if (n > 1)
    *v /= (n - 1.);
  return (0);
}


x86 assembly (clang, -O0):
<MeanVar>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1158 <MeanVar+0x58>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1125 <MeanVar+0x25>
cvtsi2sdl -0xc(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x24(%rbp)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11d7 <MeanVar+0xd7>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
subsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
subsd  (%rax),%xmm2
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    117f <MeanVar+0x7f>
cmpl   $0x1,-0xc(%rbp)
jle    1202 <MeanVar+0x102>
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  0xe12(%rip),%xmm0
subsd  %xmm0,%xmm1
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<MeanVar>:
movq   $0x0,(%rdx)
test   %esi,%esi
jle    1127 <MeanVar+0x27>
mov    %esi,%r9d
lea    -0x1(%r9),%r8
mov    %r9d,%r10d
and    $0x3,%r10d
cmp    $0x3,%r8
jae    1140 <MeanVar+0x40>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
jmp    1180 <MeanVar+0x80>
cvtsi2sd %esi,%xmm0
xorpd  %xmm1,%xmm1
xorpd  %xmm2,%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,(%rdx)
jmp    1230 <MeanVar+0x130>
mov    %r9d,%r11d
and    $0xfffffffc,%r11d
xorpd  %xmm1,%xmm1
xor    %eax,%eax
nopl   (%rax)
addsd  (%rdi,%rax,8),%xmm1
movsd  %xmm1,(%rdx)
addsd  0x8(%rdi,%rax,8),%xmm1
movsd  %xmm1,(%rdx)
addsd  0x10(%rdi,%rax,8),%xmm1
movsd  %xmm1,(%rdx)
addsd  0x18(%rdi,%rax,8),%xmm1
movsd  %xmm1,(%rdx)
add    $0x4,%rax
cmp    %rax,%r11
jne    1150 <MeanVar+0x50>
test   %r10,%r10
je     11a3 <MeanVar+0xa3>
lea    (%rdi,%rax,8),%r11
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
addsd  (%r11,%rax,8),%xmm1
movsd  %xmm1,(%rdx)
add    $0x1,%rax
cmp    %rax,%r10
jne    1190 <MeanVar+0x90>
cvtsi2sd %esi,%xmm0
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx)
movq   $0x0,(%rcx)
test   %esi,%esi
jle    1234 <MeanVar+0x134>
test   %r8,%r8
je     1237 <MeanVar+0x137>
mov    %r9d,%r8d
and    $0xfffffffe,%r8d
xorpd  %xmm1,%xmm1
xor    %eax,%eax
nopl   0x0(%rax)
movsd  (%rdi,%rax,8),%xmm2
subsd  (%rdx),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rcx)
movsd  0x8(%rdi,%rax,8),%xmm1
subsd  (%rdx),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rcx)
add    $0x2,%rax
cmp    %rax,%r8
jne    11d0 <MeanVar+0xd0>
test   $0x1,%r9b
je     121f <MeanVar+0x11f>
movsd  (%rdi,%rax,8),%xmm2
subsd  (%rdx),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rcx)
cmp    $0x2,%esi
jl     1234 <MeanVar+0x134>
addsd  0xdd4(%rip),%xmm0
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
xor    %eax,%eax
ret
xorpd  %xmm1,%xmm1
xor    %eax,%eax
test   $0x1,%r9b
jne    120a <MeanVar+0x10a>
jmp    121f <MeanVar+0x11f>

LLVM IR (-O0):
define i32 @MeanVar(double* %0, i32 %1, double* %2, double* %3) {
%5 = alloca double*
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
store double* %0, double** %5
store i32 %1, i32* %6
store double* %2, double** %7
store double* %3, double** %8
store i32 0, i32* %9
%10 = load double*, double** %7
store double 0.000000e+00, double* %10
br label %11
11: ; preds = %24, %4
%12 = load i32, i32* %9
%13 = load i32, i32* %6
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %27
15: ; preds = %11
%16 = load double*, double** %5
%17 = load i32, i32* %9
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds double, double* %16, i64 %18
%20 = load double, double* %19
%21 = load double*, double** %7
%22 = load double, double* %21
%23 = fadd double %22, %20
store double %23, double* %21
br label %24
24: ; preds = %15
%25 = load i32, i32* %9
%26 = add nsw i32 %25, 1
store i32 %26, i32* %9
br label %11
27: ; preds = %11
%28 = load i32, i32* %6
%29 = sitofp i32 %28 to double
%30 = load double*, double** %7
%31 = load double, double* %30
%32 = fdiv double %31, %29
store double %32, double* %30
store i32 0, i32* %9
%33 = load double*, double** %8
store double 0.000000e+00, double* %33
br label %34
34: ; preds = %58, %27
%35 = load i32, i32* %9
%36 = load i32, i32* %6
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %61
38: ; preds = %34
%39 = load double*, double** %5
%40 = load i32, i32* %9
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %39, i64 %41
%43 = load double, double* %42
%44 = load double*, double** %7
%45 = load double, double* %44
%46 = fsub double %43, %45
%47 = load double*, double** %5
%48 = load i32, i32* %9
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
%52 = load double*, double** %7
%53 = load double, double* %52
%54 = fsub double %51, %53
%55 = load double*, double** %8
%56 = load double, double* %55
%57 = call double @llvm.fmuladd.f64(double %46, double %54, double %56)
store double %57, double* %55
br label %58
58: ; preds = %38
%59 = load i32, i32* %9
%60 = add nsw i32 %59, 1
store i32 %60, i32* %9
br label %34
61: ; preds = %34
%62 = load i32, i32* %6
%63 = icmp sgt i32 %62, 1
br i1 %63, label %64, label %71
64: ; preds = %61
%65 = load i32, i32* %6
%66 = sitofp i32 %65 to double
%67 = fsub double %66, 1.000000e+00
%68 = load double*, double** %8
%69 = load double, double* %68
%70 = fdiv double %69, %67
store double %70, double* %68
br label %71
71: ; preds = %64, %61
ret i32 0
}

LLVM IR (-O2):
define i32 @MeanVar(double* nocapture %0, i32 %1, double* nocapture %2, double* nocapture %3) {
store double 0.000000e+00, double* %2
%5 = icmp sgt i32 %1, 0
br i1 %5, label %9, label %6
6: ; preds = %4
%7 = sitofp i32 %1 to double
%8 = fdiv double 0.000000e+00, %7
store double %8, double* %2
br label %98
9: ; preds = %4
%10 = zext i32 %1 to i64
%11 = add nsw i64 %10, -1
%12 = and i64 %10, 3
%13 = icmp ult i64 %11, 3
br i1 %13, label %38, label %14
14: ; preds = %9
%15 = and i64 %10, 4294967292
br label %16
16: ; preds = %16, %14
%17 = phi double [ 0.000000e+00, %14 ], [ %34, %16 ]
%18 = phi i64 [ 0, %14 ], [ %35, %16 ]
%19 = phi i64 [ 0, %14 ], [ %36, %16 ]
%20 = getelementptr inbounds double, double* %0, i64 %18
%21 = load double, double* %20
%22 = fadd double %21, %17
store double %22, double* %2
%23 = or i64 %18, 1
%24 = getelementptr inbounds double, double* %0, i64 %23
%25 = load double, double* %24
%26 = fadd double %25, %22
store double %26, double* %2
%27 = or i64 %18, 2
%28 = getelementptr inbounds double, double* %0, i64 %27
%29 = load double, double* %28
%30 = fadd double %29, %26
store double %30, double* %2
%31 = or i64 %18, 3
%32 = getelementptr inbounds double, double* %0, i64 %31
%33 = load double, double* %32
%34 = fadd double %33, %30
store double %34, double* %2
%35 = add nuw nsw i64 %18, 4
%36 = add i64 %19, 4
%37 = icmp eq i64 %36, %15
br i1 %37, label %38, label %16
38: ; preds = %16, %9
%39 = phi double [ undef, %9 ], [ %34, %16 ]
%40 = phi double [ 0.000000e+00, %9 ], [ %34, %16 ]
%41 = phi i64 [ 0, %9 ], [ %35, %16 ]
%42 = icmp eq i64 %12, 0
br i1 %42, label %53, label %43
43: ; preds = %38, %43
%44 = phi double [ %49, %43 ], [ %40, %38 ]
%45 = phi i64 [ %50, %43 ], [ %41, %38 ]
%46 = phi i64 [ %51, %43 ], [ 0, %38 ]
%47 = getelementptr inbounds double, double* %0, i64 %45
%48 = load double, double* %47
%49 = fadd double %48, %44
store double %49, double* %2
%50 = add nuw nsw i64 %45, 1
%51 = add i64 %46, 1
%52 = icmp eq i64 %51, %12
br i1 %52, label %53, label %43
53: ; preds = %43, %38
%54 = phi double [ %39, %38 ], [ %49, %43 ]
%55 = sitofp i32 %1 to double
%56 = fdiv double %54, %55
store double %56, double* %2
store double 0.000000e+00, double* %3
br i1 %5, label %57, label %100
57: ; preds = %53
%58 = zext i32 %1 to i64
%59 = and i64 %58, 1
%60 = icmp eq i64 %11, 0
br i1 %60, label %81, label %61
61: ; preds = %57
%62 = and i64 %58, 4294967294
br label %63
63: ; preds = %63, %61
%64 = phi double [ 0.000000e+00, %61 ], [ %77, %63 ]
%65 = phi i64 [ 0, %61 ], [ %78, %63 ]
%66 = phi i64 [ 0, %61 ], [ %79, %63 ]
%67 = getelementptr inbounds double, double* %0, i64 %65
%68 = load double, double* %67
%69 = load double, double* %2
%70 = fsub double %68, %69
%71 = tail call double @llvm.fmuladd.f64(double %70, double %70, double %64)
store double %71, double* %3
%72 = or i64 %65, 1
%73 = getelementptr inbounds double, double* %0, i64 %72
%74 = load double, double* %73
%75 = load double, double* %2
%76 = fsub double %74, %75
%77 = tail call double @llvm.fmuladd.f64(double %76, double %76, double %71)
store double %77, double* %3
%78 = add nuw nsw i64 %65, 2
%79 = add i64 %66, 2
%80 = icmp eq i64 %79, %62
br i1 %80, label %81, label %63
81: ; preds = %63, %57
%82 = phi double [ undef, %57 ], [ %77, %63 ]
%83 = phi double [ 0.000000e+00, %57 ], [ %77, %63 ]
%84 = phi i64 [ 0, %57 ], [ %78, %63 ]
%85 = icmp eq i64 %59, 0
br i1 %85, label %92, label %86
86: ; preds = %81
%87 = getelementptr inbounds double, double* %0, i64 %84
%88 = load double, double* %87
%89 = load double, double* %2
%90 = fsub double %88, %89
%91 = tail call double @llvm.fmuladd.f64(double %90, double %90, double %83)
store double %91, double* %3
br label %92
92: ; preds = %81, %86
%93 = phi double [ %82, %81 ], [ %91, %86 ]
%94 = icmp sgt i32 %1, 1
br i1 %94, label %95, label %100
95: ; preds = %92
%96 = fadd double %55, -1.000000e+00
%97 = fdiv double %93, %96
br label %98
98: ; preds = %95, %6
%99 = phi double [ 0.000000e+00, %6 ], [ %97, %95 ]
store double %99, double* %3
br label %100
100: ; preds = %98, %53, %92
ret i32 0
}