void sigs2weights(float *sigs, int nchan, float thresh) {
  float tot = 0;
  for (int i = 0; i < nchan; i++) {
    if (sigs[i] > 0) {
      sigs[i] = 1. / sigs[i] / sigs[i];
      tot += sigs[i];
    }
  }
  tot /= nchan;
  if (thresh > 0)
    for (int i = 0; i < nchan; i++)
      if (sigs[i] < tot * thresh)
        sigs[i] = 0;
}


x86 assembly (clang, -O0):
<sigs2weights>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11b1 <sigs2weights+0xb1>
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    119e <sigs2weights+0x9e>
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xea3(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
jmp    11a3 <sigs2weights+0xa3>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    111f <sigs2weights+0x1f>
cvtsi2ssl -0xc(%rbp),%xmm1
movss  -0x14(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    1230 <sigs2weights+0x130>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    122b <sigs2weights+0x12b>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x14(%rbp),%xmm0
mulss  -0x10(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1218 <sigs2weights+0x118>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    121d <sigs2weights+0x11d>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11dc <sigs2weights+0xdc>
jmp    1230 <sigs2weights+0x130>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sigs2weights>:
test   %esi,%esi
jle    1251 <sigs2weights+0x151>
mov    %esi,%eax
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
movsd  0xee9(%rip),%xmm3
xorps  %xmm1,%xmm1
jmp    1129 <sigs2weights+0x29>
nopl   0x0(%rax)
add    $0x1,%rcx
cmp    %rcx,%rax
je     1155 <sigs2weights+0x55>
movss  (%rdi,%rcx,4),%xmm4
ucomiss %xmm2,%xmm4
jbe    1120 <sigs2weights+0x20>
cvtss2sd %xmm4,%xmm4
movapd %xmm3,%xmm5
divsd  %xmm4,%xmm5
divsd  %xmm4,%xmm5
xorps  %xmm4,%xmm4
cvtsd2ss %xmm5,%xmm4
movss  %xmm4,(%rdi,%rcx,4)
addss  %xmm4,%xmm1
jmp    1120 <sigs2weights+0x20>
xorps  %xmm2,%xmm2
ucomiss %xmm2,%xmm0
jbe    1251 <sigs2weights+0x151>
test   %esi,%esi
jle    1251 <sigs2weights+0x151>
xorps  %xmm2,%xmm2
cvtsi2ss %esi,%xmm2
divss  %xmm2,%xmm1
mulss  %xmm0,%xmm1
cmp    $0x8,%esi
jae    1184 <sigs2weights+0x84>
xor    %ecx,%ecx
jmp    1269 <sigs2weights+0x169>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
movaps %xmm1,%xmm0
shufps $0x0,%xmm1,%xmm0
lea    0x0(,%rax,4),%r8
and    $0xffffffffffffffe0,%r8
xor    %esi,%esi
jmp    11ad <sigs2weights+0xad>
add    $0x20,%rsi
cmp    %rsi,%r8
je     124c <sigs2weights+0x14c>
movups (%rdi,%rsi,1),%xmm2
ucomiss %xmm2,%xmm0
jbe    11bd <sigs2weights+0xbd>
movl   $0x0,(%rdi,%rsi,1)
cmpltps %xmm0,%xmm2
pextrw $0x2,%xmm2,%edx
test   $0x1,%dl
je     11d3 <sigs2weights+0xd3>
movl   $0x0,0x4(%rdi,%rsi,1)
pextrw $0x4,%xmm2,%edx
test   $0x1,%dl
je     11e5 <sigs2weights+0xe5>
movl   $0x0,0x8(%rdi,%rsi,1)
pextrw $0x6,%xmm2,%edx
test   $0x1,%dl
je     11f7 <sigs2weights+0xf7>
movl   $0x0,0xc(%rdi,%rsi,1)
movups 0x10(%rdi,%rsi,1),%xmm2
ucomiss %xmm2,%xmm0
jbe    1209 <sigs2weights+0x109>
movl   $0x0,0x10(%rdi,%rsi,1)
cmpltps %xmm0,%xmm2
pextrw $0x2,%xmm2,%edx
test   $0x1,%dl
je     121f <sigs2weights+0x11f>
movl   $0x0,0x14(%rdi,%rsi,1)
pextrw $0x4,%xmm2,%edx
test   $0x1,%dl
je     1231 <sigs2weights+0x131>
movl   $0x0,0x18(%rdi,%rsi,1)
pextrw $0x6,%xmm2,%edx
test   $0x1,%dl
je     11a0 <sigs2weights+0xa0>
movl   $0x0,0x1c(%rdi,%rsi,1)
jmp    11a0 <sigs2weights+0xa0>
cmp    %rax,%rcx
jne    1269 <sigs2weights+0x169>
ret
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rcx
cmp    %rcx,%rax
je     1251 <sigs2weights+0x151>
ucomiss (%rdi,%rcx,4),%xmm1
jbe    1260 <sigs2weights+0x160>
movl   $0x0,(%rdi,%rcx,4)
jmp    1260 <sigs2weights+0x160>

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

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