void avevar(float data[], unsigned long n, float *ave, float *var) {
  unsigned long j;
  float s, ep;
  for (*ave = 0.0, j = 1; j <= n; j++)
    *ave += data[j];
  *ave /= n;
  *var = ep = 0.0;
  for (j = 1; j <= n; j++) {
    s = data[j] - (*ave);
    ep += s;
    *var += s * s;
  }
  *var = (*var - ep * ep / n) / (n - 1);
}


x86 assembly (clang, -O0):
<avevar>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movq   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
ja     115f <avevar+0x5f>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1127 <avevar+0x27>
mov    -0x10(%rbp),%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x38(%rbp)
test   %rax,%rax
movss  %xmm0,-0x34(%rbp)
js     119c <avevar+0x9c>
movss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movq   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
ja     1232 <avevar+0x132>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
subss  (%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
addss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x2c(%rbp),%xmm0
movss  -0x2c(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    11cc <avevar+0xcc>
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x30(%rbp),%xmm0
mulss  %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x40(%rbp)
test   %rax,%rax
movss  %xmm0,-0x3c(%rbp)
js     128e <avevar+0x18e>
movss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x48(%rbp),%rax
movss  -0x50(%rbp),%xmm0
movss  -0x4c(%rbp),%xmm1
movss  -0x3c(%rbp),%xmm2
divss  %xmm2,%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x5c(%rbp)
add    $0xffffffffffffffff,%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x58(%rbp)
test   %rax,%rax
movss  %xmm0,-0x54(%rbp)
js     12eb <avevar+0x1eb>
movss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x5c(%rbp),%xmm0
movss  -0x54(%rbp),%xmm1
divss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<avevar>:
movl   $0x0,(%rdx)
test   %rsi,%rsi
je     112b <avevar+0x2b>
lea    -0x1(%rsi),%r8
mov    %esi,%r9d
and    $0x3,%r9d
cmp    $0x3,%r8
jae    113b <avevar+0x3b>
mov    $0x1,%eax
xorps  %xmm1,%xmm1
test   %r9,%r9
jne    118a <avevar+0x8a>
jmp    11a3 <avevar+0xa3>
js     1255 <avevar+0x155>
cvtsi2ss %rsi,%xmm0
jmp    126c <avevar+0x16c>
mov    %rsi,%r10
and    $0xfffffffffffffffc,%r10
xorps  %xmm1,%xmm1
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
addss  0x4(%rdi,%rax,4),%xmm1
movss  %xmm1,(%rdx)
addss  0x8(%rdi,%rax,4),%xmm1
movss  %xmm1,(%rdx)
addss  0xc(%rdi,%rax,4),%xmm1
movss  %xmm1,(%rdx)
addss  0x10(%rdi,%rax,4),%xmm1
movss  %xmm1,(%rdx)
add    $0x4,%rax
cmp    %rax,%r10
jne    1150 <avevar+0x50>
add    $0x1,%rax
test   %r9,%r9
je     11a3 <avevar+0xa3>
lea    (%rdi,%rax,4),%r10
xor    %eax,%eax
addss  (%r10,%rax,4),%xmm1
movss  %xmm1,(%rdx)
add    $0x1,%rax
cmp    %rax,%r9
jne    1190 <avevar+0x90>
mov    %rsi,%r9
shr    %r9
mov    %esi,%eax
and    $0x1,%eax
or     %r9,%rax
test   %rsi,%rsi
js     11bd <avevar+0xbd>
cvtsi2ss %rsi,%xmm0
jmp    11c6 <avevar+0xc6>
cvtsi2ss %rax,%xmm0
addss  %xmm0,%xmm0
divss  %xmm0,%xmm1
movss  %xmm1,(%rdx)
movl   $0x0,(%rcx)
je     124d <avevar+0x14d>
test   %r8,%r8
je     127f <avevar+0x17f>
mov    %rsi,%r8
and    $0xfffffffffffffffe,%r8
neg    %r8
mov    $0x1,%eax
xorps  %xmm1,%xmm1
xorps  %xmm2,%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rdi,%rax,4),%xmm3
subss  (%rdx),%xmm3
addss  %xmm3,%xmm2
mulss  %xmm3,%xmm3
addss  %xmm1,%xmm3
movss  %xmm3,(%rcx)
movss  0x4(%rdi,%rax,4),%xmm1
subss  (%rdx),%xmm1
addss  %xmm1,%xmm2
mulss  %xmm1,%xmm1
addss  %xmm3,%xmm1
movss  %xmm1,(%rcx)
lea    (%r8,%rax,1),%r9
add    $0x2,%r9
add    $0x2,%rax
cmp    $0x1,%r9
jne    1200 <avevar+0x100>
test   $0x1,%sil
jne    1290 <avevar+0x190>
jmp    12a9 <avevar+0x1a9>
xorps  %xmm1,%xmm1
xorps  %xmm2,%xmm2
jmp    12a9 <avevar+0x1a9>
mov    %rsi,%rax
shr    %rax
mov    %esi,%edi
and    $0x1,%edi
or     %rax,%rdi
cvtsi2ss %rdi,%xmm0
addss  %xmm0,%xmm0
xorps  %xmm1,%xmm1
xorps  %xmm2,%xmm2
divss  %xmm0,%xmm2
movss  %xmm2,(%rdx)
xorps  %xmm2,%xmm2
jmp    12a9 <avevar+0x1a9>
mov    $0x1,%eax
xorps  %xmm1,%xmm1
xorps  %xmm2,%xmm2
test   $0x1,%sil
je     12a9 <avevar+0x1a9>
movss  (%rdi,%rax,4),%xmm3
subss  (%rdx),%xmm3
addss  %xmm3,%xmm2
mulss  %xmm3,%xmm3
addss  %xmm3,%xmm1
movss  %xmm1,(%rcx)
mulss  %xmm2,%xmm2
divss  %xmm0,%xmm2
subss  %xmm2,%xmm1
add    $0xffffffffffffffff,%rsi
js     12c5 <avevar+0x1c5>
xorps  %xmm0,%xmm0
cvtsi2ss %rsi,%xmm0
jmp    12dd <avevar+0x1dd>
mov    %rsi,%rax
shr    %rax
and    $0x1,%esi
or     %rax,%rsi
xorps  %xmm0,%xmm0
cvtsi2ss %rsi,%xmm0
addss  %xmm0,%xmm0
divss  %xmm0,%xmm1
movss  %xmm1,(%rcx)
ret

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

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