#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double sqrt(float);
void arrayf_mean(float *a, int n, double *mean, double *stddev) {
  int i;
  double sum, avg, sum_sq;
  if (n < 1)
    return;
  for (sum = sum_sq = i = 0; i < n; i++)
    sum += a[i];
  avg = sum / n;
  if (mean != NULL)
    (*mean) = avg;
  if (stddev != NULL) {
    float sum_sq;
    for (sum_sq = i = 0; i < n; i++)
      sum_sq += (a[i] - avg) * (a[i] - avg);
    (*stddev) = sqrt(sum_sq / n);
  }
}


x86 assembly (clang, -O0):
<arrayf_mean>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
cmpl   $0x1,-0xc(%rbp)
jge    1136 <arrayf_mean+0x26>
jmp    1242 <arrayf_mean+0x132>
movl   $0x0,-0x24(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1182 <arrayf_mean+0x72>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    114d <arrayf_mean+0x3d>
movsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
cmpq   $0x0,-0x18(%rbp)
je     11ad <arrayf_mean+0x9d>
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
cmpq   $0x0,-0x20(%rbp)
je     1242 <arrayf_mean+0x132>
movl   $0x0,-0x24(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1227 <arrayf_mean+0x117>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
subsd  -0x38(%rbp),%xmm2
movss  -0x44(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11c7 <arrayf_mean+0xb7>
movss  -0x44(%rbp),%xmm0
cvtsi2ssl -0xc(%rbp),%xmm1
divss  %xmm1,%xmm0
call   1030 <sqrt@plt>
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<arrayf_mean>:
push   %r14
push   %rbx
push   %rax
test   %esi,%esi
jle    128b <arrayf_mean+0x17b>
mov    %rcx,%r14
mov    %esi,%r9d
lea    -0x1(%r9),%r8
mov    %r9d,%eax
and    $0x3,%eax
cmp    $0x3,%r8
jae    113a <arrayf_mean+0x2a>
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
jmp    1193 <arrayf_mean+0x83>
mov    %r9d,%ebx
and    $0xfffffffc,%ebx
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movss  (%rdi,%rcx,4),%xmm0
movss  0x4(%rdi,%rcx,4),%xmm2
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
addsd  %xmm0,%xmm1
movss  0x8(%rdi,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
movss  0xc(%rdi,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm0,%xmm1
add    $0x4,%rcx
cmp    %rcx,%rbx
jne    1150 <arrayf_mean+0x40>
test   %rax,%rax
je     11b6 <arrayf_mean+0xa6>
lea    (%rdi,%rcx,4),%rcx
xor    %ebx,%ebx
xchg   %ax,%ax
movss  (%rcx,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm0,%xmm1
add    $0x1,%rbx
cmp    %rbx,%rax
jne    11a0 <arrayf_mean+0x90>
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
divsd  %xmm0,%xmm1
test   %rdx,%rdx
je     11ca <arrayf_mean+0xba>
movsd  %xmm1,(%rdx)
test   %r14,%r14
je     128b <arrayf_mean+0x17b>
test   %esi,%esi
jle    124c <arrayf_mean+0x13c>
test   %r8,%r8
je     1252 <arrayf_mean+0x142>
mov    %r9d,%ecx
and    $0xfffffffe,%ecx
xorpd  %xmm0,%xmm0
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movss  (%rdi,%rax,4),%xmm2
movss  0x4(%rdi,%rax,4),%xmm3
cvtss2sd %xmm2,%xmm2
subsd  %xmm1,%xmm2
cvtss2sd %xmm0,%xmm0
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsd2ss %xmm2,%xmm0
xorps  %xmm2,%xmm2
cvtss2sd %xmm3,%xmm2
subsd  %xmm1,%xmm2
cvtss2sd %xmm0,%xmm0
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsd2ss %xmm2,%xmm0
add    $0x2,%rax
cmp    %rax,%rcx
jne    11f0 <arrayf_mean+0xe0>
xorps  %xmm2,%xmm2
cvtss2sd %xmm0,%xmm2
test   $0x1,%r9b
jne    125d <arrayf_mean+0x14d>
jmp    1276 <arrayf_mean+0x166>
xorpd  %xmm0,%xmm0
jmp    1276 <arrayf_mean+0x166>
xorps  %xmm2,%xmm2
xor    %eax,%eax
test   $0x1,%r9b
je     1276 <arrayf_mean+0x166>
movss  (%rdi,%rax,4),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
cvtsd2ss %xmm0,%xmm0
xorps  %xmm1,%xmm1
cvtsi2ss %esi,%xmm1
divss  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,(%r14)
add    $0x8,%rsp
pop    %rbx
pop    %r14
ret

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

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