void remove_dc_bias(float speech[], int len, float sigout[]) {
  float bias;
  int i;
  bias = 0.0f;
  for (i = 0; i < len; i++)
    bias += speech[i];
  bias /= len;
  for (i = 0; i < len; i++)
    sigout[i] = speech[i] - bias;
}


x86 assembly (clang, -O0):
<remove_dc_bias>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    114f <remove_dc_bias+0x4f>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    111e <remove_dc_bias+0x1e>
cvtsi2ssl -0xc(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11a2 <remove_dc_bias+0xa2>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x1c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1169 <remove_dc_bias+0x69>
pop    %rbp
ret

x86 assembly (clang, -O2):
<remove_dc_bias>:
test   %esi,%esi
jle    12f3 <remove_dc_bias+0x1f3>
mov    %esi,%r9d
lea    -0x1(%r9),%rax
mov    %r9d,%r8d
and    $0x7,%r8d
cmp    $0x7,%rax
jae    1123 <remove_dc_bias+0x23>
xorps  %xmm0,%xmm0
xor    %eax,%eax
jmp    1168 <remove_dc_bias+0x68>
mov    %r9d,%ecx
and    $0xfffffff8,%ecx
xorps  %xmm0,%xmm0
xor    %eax,%eax
xchg   %ax,%ax
addss  (%rdi,%rax,4),%xmm0
addss  0x4(%rdi,%rax,4),%xmm0
addss  0x8(%rdi,%rax,4),%xmm0
addss  0xc(%rdi,%rax,4),%xmm0
addss  0x10(%rdi,%rax,4),%xmm0
addss  0x14(%rdi,%rax,4),%xmm0
addss  0x18(%rdi,%rax,4),%xmm0
addss  0x1c(%rdi,%rax,4),%xmm0
add    $0x8,%rax
cmp    %rax,%rcx
jne    1130 <remove_dc_bias+0x30>
test   %r8,%r8
je     118e <remove_dc_bias+0x8e>
lea    (%rdi,%rax,4),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addss  (%rax,%rcx,4),%xmm0
add    $0x1,%rcx
cmp    %rcx,%r8
jne    1180 <remove_dc_bias+0x80>
test   %esi,%esi
jle    12f3 <remove_dc_bias+0x1f3>
cvtsi2ss %esi,%xmm1
divss  %xmm1,%xmm0
cmp    $0x8,%esi
jb     11bd <remove_dc_bias+0xbd>
lea    (%rdi,%r9,4),%rax
cmp    %rdx,%rax
jbe    125c <remove_dc_bias+0x15c>
lea    (%rdx,%r9,4),%rax
cmp    %rdi,%rax
jbe    125c <remove_dc_bias+0x15c>
xor    %ecx,%ecx
mov    %rcx,%rax
not    %rax
add    %r9,%rax
mov    %r9,%rsi
and    $0x3,%rsi
je     11f8 <remove_dc_bias+0xf8>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  (%rdi,%rcx,4),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,(%rdx,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    11e0 <remove_dc_bias+0xe0>
cmp    $0x3,%rax
jb     12f3 <remove_dc_bias+0x1f3>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rdi,%rcx,4),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,(%rdx,%rcx,4)
movss  0x4(%rdi,%rcx,4),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,0x4(%rdx,%rcx,4)
movss  0x8(%rdi,%rcx,4),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,0x8(%rdx,%rcx,4)
movss  0xc(%rdi,%rcx,4),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,0xc(%rdx,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%r9
jne    1210 <remove_dc_bias+0x110>
jmp    12f3 <remove_dc_bias+0x1f3>
mov    %r9d,%ecx
and    $0xfffffff8,%ecx
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
lea    -0x8(%rcx),%rax
mov    %rax,%r8
shr    $0x3,%r8
add    $0x1,%r8
test   %rax,%rax
je     12f4 <remove_dc_bias+0x1f4>
mov    %r8,%rax
and    $0xfffffffffffffffe,%rax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
movups (%rdi,%rsi,4),%xmm2
movups 0x10(%rdi,%rsi,4),%xmm3
subps  %xmm1,%xmm2
subps  %xmm1,%xmm3
movups %xmm2,(%rdx,%rsi,4)
movups %xmm3,0x10(%rdx,%rsi,4)
movups 0x20(%rdi,%rsi,4),%xmm2
movups 0x30(%rdi,%rsi,4),%xmm3
subps  %xmm1,%xmm2
subps  %xmm1,%xmm3
movups %xmm2,0x20(%rdx,%rsi,4)
movups %xmm3,0x30(%rdx,%rsi,4)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rax
jne    1290 <remove_dc_bias+0x190>
test   $0x1,%r8b
je     12ea <remove_dc_bias+0x1ea>
movups (%rdi,%rsi,4),%xmm2
movups 0x10(%rdi,%rsi,4),%xmm3
subps  %xmm1,%xmm2
subps  %xmm1,%xmm3
movups %xmm2,(%rdx,%rsi,4)
movups %xmm3,0x10(%rdx,%rsi,4)
cmp    %r9,%rcx
jne    11bf <remove_dc_bias+0xbf>
ret
xor    %esi,%esi
test   $0x1,%r8b
jne    12d2 <remove_dc_bias+0x1d2>
jmp    12ea <remove_dc_bias+0x1ea>

LLVM IR (-O0):
define void @remove_dc_bias(float* %0, i32 %1, float* %2) {
%4 = alloca float*
%5 = alloca i32
%6 = alloca float*
%7 = alloca float
%8 = alloca i32
store float* %0, float** %4
store i32 %1, i32* %5
store float* %2, float** %6
store float 0.000000e+00, float* %7
store i32 0, i32* %8
br label %9
9: ; preds = %21, %3
%10 = load i32, i32* %8
%11 = load i32, i32* %5
%12 = icmp slt i32 %10, %11
br i1 %12, label %13, label %24
13: ; preds = %9
%14 = load float*, float** %4
%15 = load i32, i32* %8
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds float, float* %14, i64 %16
%18 = load float, float* %17
%19 = load float, float* %7
%20 = fadd float %19, %18
store float %20, float* %7
br label %21
21: ; preds = %13
%22 = load i32, i32* %8
%23 = add nsw i32 %22, 1
store i32 %23, i32* %8
br label %9
24: ; preds = %9
%25 = load i32, i32* %5
%26 = sitofp i32 %25 to float
%27 = load float, float* %7
%28 = fdiv float %27, %26
store float %28, float* %7
store i32 0, i32* %8
br label %29
29: ; preds = %45, %24
%30 = load i32, i32* %8
%31 = load i32, i32* %5
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %48
33: ; preds = %29
%34 = load float*, float** %4
%35 = load i32, i32* %8
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds float, float* %34, i64 %36
%38 = load float, float* %37
%39 = load float, float* %7
%40 = fsub float %38, %39
%41 = load float*, float** %6
%42 = load i32, i32* %8
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
store float %40, float* %44
br label %45
45: ; preds = %33
%46 = load i32, i32* %8
%47 = add nsw i32 %46, 1
store i32 %47, i32* %8
br label %29
48: ; preds = %29
ret void
}

LLVM IR (-O2):
define void @remove_dc_bias(float* nocapture %0, i32 %1, float* nocapture writeonly %2) {
%4 = icmp sgt i32 %1, 0
br i1 %4, label %5, label %182
5: ; preds = %3
%6 = zext i32 %1 to i64
%7 = add nsw i64 %6, -1
%8 = and i64 %6, 7
%9 = icmp ult i64 %7, 7
br i1 %9, label %50, label %10
10: ; preds = %5
%11 = and i64 %6, 4294967288
br label %12
12: ; preds = %12, %10
%13 = phi i64 [ 0, %10 ], [ %47, %12 ]
%14 = phi float [ 0.000000e+00, %10 ], [ %46, %12 ]
%15 = phi i64 [ 0, %10 ], [ %48, %12 ]
%16 = getelementptr inbounds float, float* %0, i64 %13
%17 = load float, float* %16
%18 = fadd float %14, %17
%19 = or i64 %13, 1
%20 = getelementptr inbounds float, float* %0, i64 %19
%21 = load float, float* %20
%22 = fadd float %18, %21
%23 = or i64 %13, 2
%24 = getelementptr inbounds float, float* %0, i64 %23
%25 = load float, float* %24
%26 = fadd float %22, %25
%27 = or i64 %13, 3
%28 = getelementptr inbounds float, float* %0, i64 %27
%29 = load float, float* %28
%30 = fadd float %26, %29
%31 = or i64 %13, 4
%32 = getelementptr inbounds float, float* %0, i64 %31
%33 = load float, float* %32
%34 = fadd float %30, %33
%35 = or i64 %13, 5
%36 = getelementptr inbounds float, float* %0, i64 %35
%37 = load float, float* %36
%38 = fadd float %34, %37
%39 = or i64 %13, 6
%40 = getelementptr inbounds float, float* %0, i64 %39
%41 = load float, float* %40
%42 = fadd float %38, %41
%43 = or i64 %13, 7
%44 = getelementptr inbounds float, float* %0, i64 %43
%45 = load float, float* %44
%46 = fadd float %42, %45
%47 = add nuw nsw i64 %13, 8
%48 = add i64 %15, 8
%49 = icmp eq i64 %48, %11
br i1 %49, label %50, label %12
50: ; preds = %12, %5
%51 = phi float [ undef, %5 ], [ %46, %12 ]
%52 = phi i64 [ 0, %5 ], [ %47, %12 ]
%53 = phi float [ 0.000000e+00, %5 ], [ %46, %12 ]
%54 = icmp eq i64 %8, 0
br i1 %54, label %65, label %55
55: ; preds = %50, %55
%56 = phi i64 [ %62, %55 ], [ %52, %50 ]
%57 = phi float [ %61, %55 ], [ %53, %50 ]
%58 = phi i64 [ %63, %55 ], [ 0, %50 ]
%59 = getelementptr inbounds float, float* %0, i64 %56
%60 = load float, float* %59
%61 = fadd float %57, %60
%62 = add nuw nsw i64 %56, 1
%63 = add i64 %58, 1
%64 = icmp eq i64 %63, %8
br i1 %64, label %65, label %55
65: ; preds = %55, %50
%66 = phi float [ %51, %50 ], [ %61, %55 ]
%67 = sitofp i32 %1 to float
%68 = fdiv float %66, %67
br i1 %4, label %69, label %182
69: ; preds = %65
%70 = zext i32 %1 to i64
%71 = icmp ult i32 %1, 8
br i1 %71, label %140, label %72
72: ; preds = %69
%73 = getelementptr float, float* %2, i64 %6
%74 = getelementptr float, float* %0, i64 %6
%75 = icmp ugt float* %74, %2
%76 = icmp ugt float* %73, %0
%77 = and i1 %75, %76
br i1 %77, label %140, label %78
78: ; preds = %72
%79 = and i64 %6, 4294967288
%80 = insertelement <4 x float> poison, float %68, i64 0
%81 = shufflevector <4 x float> %80, <4 x float> poison, <4 x i32> zeroinitializer
%82 = insertelement <4 x float> poison, float %68, i64 0
%83 = shufflevector <4 x float> %82, <4 x float> poison, <4 x i32> zeroinitializer
%84 = add nsw i64 %79, -8
%85 = lshr exact i64 %84, 3
%86 = add nuw nsw i64 %85, 1
%87 = and i64 %86, 1
%88 = icmp eq i64 %84, 0
br i1 %88, label %122, label %89
89: ; preds = %78
%90 = and i64 %86, 4611686018427387902
br label %91
91: ; preds = %91, %89
%92 = phi i64 [ 0, %89 ], [ %119, %91 ]
%93 = phi i64 [ 0, %89 ], [ %120, %91 ]
%94 = getelementptr inbounds float, float* %0, i64 %92
%95 = bitcast float* %94 to <4 x float>*
%96 = load <4 x float>, <4 x float>* %95
%97 = getelementptr inbounds float, float* %94, i64 4
%98 = bitcast float* %97 to <4 x float>*
%99 = load <4 x float>, <4 x float>* %98
%100 = fsub <4 x float> %96, %81
%101 = fsub <4 x float> %99, %83
%102 = getelementptr inbounds float, float* %2, i64 %92
%103 = bitcast float* %102 to <4 x float>*
store <4 x float> %100, <4 x float>* %103
%104 = getelementptr inbounds float, float* %102, i64 4
%105 = bitcast float* %104 to <4 x float>*
store <4 x float> %101, <4 x float>* %105
%106 = or i64 %92, 8
%107 = getelementptr inbounds float, float* %0, i64 %106
%108 = bitcast float* %107 to <4 x float>*
%109 = load <4 x float>, <4 x float>* %108
%110 = getelementptr inbounds float, float* %107, i64 4
%111 = bitcast float* %110 to <4 x float>*
%112 = load <4 x float>, <4 x float>* %111
%113 = fsub <4 x float> %109, %81
%114 = fsub <4 x float> %112, %83
%115 = getelementptr inbounds float, float* %2, i64 %106
%116 = bitcast float* %115 to <4 x float>*
store <4 x float> %113, <4 x float>* %116
%117 = getelementptr inbounds float, float* %115, i64 4
%118 = bitcast float* %117 to <4 x float>*
store <4 x float> %114, <4 x float>* %118
%119 = add nuw i64 %92, 16
%120 = add i64 %93, 2
%121 = icmp eq i64 %120, %90
br i1 %121, label %122, label %91
122: ; preds = %91, %78
%123 = phi i64 [ 0, %78 ], [ %119, %91 ]
%124 = icmp eq i64 %87, 0
br i1 %124, label %138, label %125
125: ; preds = %122
%126 = getelementptr inbounds float, float* %0, i64 %123
%127 = bitcast float* %126 to <4 x float>*
%128 = load <4 x float>, <4 x float>* %127
%129 = getelementptr inbounds float, float* %126, i64 4
%130 = bitcast float* %129 to <4 x float>*
%131 = load <4 x float>, <4 x float>* %130
%132 = fsub <4 x float> %128, %81
%133 = fsub <4 x float> %131, %83
%134 = getelementptr inbounds float, float* %2, i64 %123
%135 = bitcast float* %134 to <4 x float>*
store <4 x float> %132, <4 x float>* %135
%136 = getelementptr inbounds float, float* %134, i64 4
%137 = bitcast float* %136 to <4 x float>*
store <4 x float> %133, <4 x float>* %137
br label %138
138: ; preds = %122, %125
%139 = icmp eq i64 %79, %6
br i1 %139, label %182, label %140
140: ; preds = %72, %69, %138
%141 = phi i64 [ 0, %72 ], [ 0, %69 ], [ %79, %138 ]
%142 = xor i64 %141, -1
%143 = add nsw i64 %142, %6
%144 = and i64 %6, 3
%145 = icmp eq i64 %144, 0
br i1 %145, label %156, label %146
146: ; preds = %140, %146
%147 = phi i64 [ %153, %146 ], [ %141, %140 ]
%148 = phi i64 [ %154, %146 ], [ 0, %140 ]
%149 = getelementptr inbounds float, float* %0, i64 %147
%150 = load float, float* %149
%151 = fsub float %150, %68
%152 = getelementptr inbounds float, float* %2, i64 %147
store float %151, float* %152
%153 = add nuw nsw i64 %147, 1
%154 = add i64 %148, 1
%155 = icmp eq i64 %154, %144
br i1 %155, label %156, label %146
156: ; preds = %146, %140
%157 = phi i64 [ %141, %140 ], [ %153, %146 ]
%158 = icmp ult i64 %143, 3
br i1 %158, label %182, label %159
159: ; preds = %156, %159
%160 = phi i64 [ %180, %159 ], [ %157, %156 ]
%161 = getelementptr inbounds float, float* %0, i64 %160
%162 = load float, float* %161
%163 = fsub float %162, %68
%164 = getelementptr inbounds float, float* %2, i64 %160
store float %163, float* %164
%165 = add nuw nsw i64 %160, 1
%166 = getelementptr inbounds float, float* %0, i64 %165
%167 = load float, float* %166
%168 = fsub float %167, %68
%169 = getelementptr inbounds float, float* %2, i64 %165
store float %168, float* %169
%170 = add nuw nsw i64 %160, 2
%171 = getelementptr inbounds float, float* %0, i64 %170
%172 = load float, float* %171
%173 = fsub float %172, %68
%174 = getelementptr inbounds float, float* %2, i64 %170
store float %173, float* %174
%175 = add nuw nsw i64 %160, 3
%176 = getelementptr inbounds float, float* %0, i64 %175
%177 = load float, float* %176
%178 = fsub float %177, %68
%179 = getelementptr inbounds float, float* %2, i64 %175
store float %178, float* %179
%180 = add nuw nsw i64 %160, 4
%181 = icmp eq i64 %180, %70
br i1 %181, label %182, label %159
182: ; preds = %156, %159, %138, %3, %65
ret void
}