void addQual(double *qual, double *qmean, double *qm2, int total, int n) {
  int i;
  double d;
  for (i = 0; i < n; i++) {
    d = qual[i] - qmean[i];
    qmean[i] = qmean[i] + (d / total);
    qm2[i] = qm2[i] + d * (qual[i] - qmean[i]);
  }
}


x86 assembly (clang, -O0):
<addQual>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11c4 <addQual+0xc4>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x30(%rbp),%xmm1
cvtsi2sdl -0x1c(%rbp),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111e <addQual+0x1e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<addQual>:
push   %rbp
push   %r15
push   %r14
push   %rbx
test   %r8d,%r8d
jle    120f <addQual+0x10f>
cvtsi2sd %ecx,%xmm0
mov    %r8d,%r9d
cmp    $0x1,%r8d
jne    1123 <addQual+0x23>
xor    %ecx,%ecx
jmp    11d0 <addQual+0xd0>
lea    (%rdx,%r9,8),%rcx
cmp    %rsi,%rcx
seta   %r15b
lea    (%rsi,%r9,8),%rbx
cmp    %rdx,%rbx
seta   %r14b
lea    (%rdi,%r9,8),%rbp
cmp    %rsi,%rbp
seta   %al
cmp    %rdi,%rbx
seta   %r10b
cmp    %rdx,%rbp
seta   %r8b
cmp    %rdi,%rcx
seta   %r11b
xor    %ecx,%ecx
test   %r14b,%r15b
jne    11d0 <addQual+0xd0>
and    %r10b,%al
jne    11d0 <addQual+0xd0>
and    %r11b,%r8b
jne    11d0 <addQual+0xd0>
mov    %r9d,%ecx
and    $0xfffffffe,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %eax,%eax
nopl   0x0(%rax)
movupd (%rdi,%rax,8),%xmm2
movupd (%rsi,%rax,8),%xmm3
subpd  %xmm3,%xmm2
movapd %xmm2,%xmm4
divpd  %xmm1,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%rsi,%rax,8)
movupd (%rdx,%rax,8),%xmm3
movupd (%rdi,%rax,8),%xmm5
subpd  %xmm4,%xmm5
mulpd  %xmm2,%xmm5
addpd  %xmm3,%xmm5
movupd %xmm5,(%rdx,%rax,8)
add    $0x2,%rax
cmp    %rax,%rcx
jne    1180 <addQual+0x80>
cmp    %r9,%rcx
je     120f <addQual+0x10f>
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rcx,8),%xmm1
movsd  (%rsi,%rcx,8),%xmm2
subsd  %xmm2,%xmm1
movapd %xmm1,%xmm3
divsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
movsd  (%rdi,%rcx,8),%xmm2
subsd  %xmm3,%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rdx,%rcx,8),%xmm2
movsd  %xmm2,(%rdx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r9
jne    11d0 <addQual+0xd0>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @addQual(double* nocapture %0, double* nocapture %1, double* nocapture %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %72
7: ; preds = %5
%8 = sitofp i32 %3 to double
%9 = zext i32 %4 to i64
%10 = icmp eq i32 %4, 1
br i1 %10, label %54, label %11
11: ; preds = %7
%12 = getelementptr double, double* %1, i64 %9
%13 = getelementptr double, double* %2, i64 %9
%14 = getelementptr double, double* %0, i64 %9
%15 = icmp ugt double* %13, %1
%16 = icmp ugt double* %12, %2
%17 = and i1 %15, %16
%18 = icmp ugt double* %14, %1
%19 = icmp ugt double* %12, %0
%20 = and i1 %18, %19
%21 = or i1 %17, %20
%22 = icmp ugt double* %14, %2
%23 = icmp ugt double* %13, %0
%24 = and i1 %22, %23
%25 = or i1 %21, %24
br i1 %25, label %54, label %26
26: ; preds = %11
%27 = and i64 %9, 4294967294
%28 = insertelement <2 x double> poison, double %8, i64 0
%29 = shufflevector <2 x double> %28, <2 x double> poison, <2 x i32> zeroinitializer
br label %30
30: ; preds = %30, %26
%31 = phi i64 [ 0, %26 ], [ %50, %30 ]
%32 = getelementptr inbounds double, double* %0, i64 %31
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = getelementptr inbounds double, double* %1, i64 %31
%36 = bitcast double* %35 to <2 x double>*
%37 = load <2 x double>, <2 x double>* %36
%38 = fsub <2 x double> %34, %37
%39 = fdiv <2 x double> %38, %29
%40 = fadd <2 x double> %37, %39
%41 = bitcast double* %35 to <2 x double>*
store <2 x double> %40, <2 x double>* %41
%42 = getelementptr inbounds double, double* %2, i64 %31
%43 = bitcast double* %42 to <2 x double>*
%44 = load <2 x double>, <2 x double>* %43
%45 = bitcast double* %32 to <2 x double>*
%46 = load <2 x double>, <2 x double>* %45
%47 = fsub <2 x double> %46, %40
%48 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %38, <2 x double> %47, <2 x double> %44)
%49 = bitcast double* %42 to <2 x double>*
store <2 x double> %48, <2 x double>* %49
%50 = add nuw i64 %31, 2
%51 = icmp eq i64 %50, %27
br i1 %51, label %52, label %30
52: ; preds = %30
%53 = icmp eq i64 %27, %9
br i1 %53, label %72, label %54
54: ; preds = %11, %7, %52
%55 = phi i64 [ 0, %11 ], [ 0, %7 ], [ %27, %52 ]
br label %56
56: ; preds = %54, %56
%57 = phi i64 [ %70, %56 ], [ %55, %54 ]
%58 = getelementptr inbounds double, double* %0, i64 %57
%59 = load double, double* %58
%60 = getelementptr inbounds double, double* %1, i64 %57
%61 = load double, double* %60
%62 = fsub double %59, %61
%63 = fdiv double %62, %8
%64 = fadd double %61, %63
store double %64, double* %60
%65 = getelementptr inbounds double, double* %2, i64 %57
%66 = load double, double* %65
%67 = load double, double* %58
%68 = fsub double %67, %64
%69 = tail call double @llvm.fmuladd.f64(double %62, double %68, double %66)
store double %69, double* %65
%70 = add nuw nsw i64 %57, 1
%71 = icmp eq i64 %70, %9
br i1 %71, label %72, label %56
72: ; preds = %56, %52, %5
ret void
}