void frelon_lines(float *img, int ns, int nf, float cut) {
  int i, j, p, npx;
  float rowsum, avg;
  avg = img[0];
#pragma omp parallel for private(i, j, rowsum, npx, p) firstprivate(avg)
  for (i = 0; i < ns; i++) {
    rowsum = 0.;
    npx = 0;
    p = i * nf;
    for (j = 0; j < nf; j++) {
      if (img[p + j] < cut) {
        rowsum += img[p + j];
        npx++;
      }
    }
    if (npx > 0)
      avg = rowsum / npx;
    for (j = 0; j < nf; j++)
      img[p + j] = img[p + j] - avg;
  }
}


x86 assembly (clang, -O0):
<frelon_lines>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1231 <frelon_lines+0x131>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x18(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x20(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11b7 <frelon_lines+0xb7>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x14(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    11a4 <frelon_lines+0xa4>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a9 <frelon_lines+0xa9>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1153 <frelon_lines+0x53>
cmpl   $0x0,-0x24(%rbp)
jle    11d4 <frelon_lines+0xd4>
movss  -0x28(%rbp),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    121e <frelon_lines+0x11e>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11db <frelon_lines+0xdb>
jmp    1223 <frelon_lines+0x123>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1127 <frelon_lines+0x27>
pop    %rbp
ret

x86 assembly (clang, -O2):
<frelon_lines>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdi,-0x28(%rsp)
test   %esi,%esi
jle    1319 <frelon_lines+0x219>
mov    -0x28(%rsp),%rbp
movss  0x0(%rbp),%xmm1
movslq %edx,%rax
mov    %esi,%edi
mov    %rax,%rsi
mov    %esi,%r11d
mov    %r11d,%eax
and    $0xfffffff8,%eax
mov    %rax,-0x8(%rsp)
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x10(%rsp)
mov    %rax,%rcx
shr    $0x3,%rcx
add    $0x1,%rcx
mov    %r11d,%r8d
and    $0xfffffffe,%r8d
mov    %rcx,%rax
mov    %rcx,-0x18(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x20(%rsp)
lea    0x4(%rbp),%r15
lea    0x0(,%rsi,4),%r12
lea    0x30(%rbp),%rcx
xor    %r13d,%r13d
movss  0xe82(%rip),%xmm2
jmp    1196 <frelon_lines+0x96>
add    $0x1,%r13
add    %r12,%r15
add    %r12,%rcx
add    %r12,%rbp
cmp    %rdi,%r13
je     1319 <frelon_lines+0x219>
test   %edx,%edx
jle    1180 <frelon_lines+0x80>
xor    %r10d,%r10d
xorps  %xmm3,%xmm3
xor    %r9d,%r9d
cmp    $0x1,%edx
jne    1216 <frelon_lines+0x116>
mov    %r13,%r14
imul   %rsi,%r14
test   $0x1,%r11b
je     11d9 <frelon_lines+0xd9>
add    %r14,%r10
mov    -0x28(%rsp),%rax
movss  (%rax,%r10,4),%xmm4
xor    %eax,%eax
ucomiss %xmm4,%xmm0
seta   %bl
ja     11d0 <frelon_lines+0xd0>
movaps %xmm2,%xmm4
addss  %xmm4,%xmm3
mov    %bl,%al
add    %eax,%r9d
test   %r9d,%r9d
je     11ed <frelon_lines+0xed>
xorps  %xmm1,%xmm1
cvtsi2ss %r9d,%xmm1
divss  %xmm1,%xmm3
movaps %xmm3,%xmm1
test   %edx,%edx
jle    1180 <frelon_lines+0x80>
cmp    $0x8,%edx
jae    1244 <frelon_lines+0x144>
xor    %eax,%eax
jmp    12f0 <frelon_lines+0x1f0>
nopl   (%rax)
addss  %xmm4,%xmm3
addss  %xmm5,%xmm3
mov    %bl,%al
add    %eax,%r9d
add    $0x2,%r10
cmp    %r10,%r8
je     11a8 <frelon_lines+0xa8>
movss  -0x4(%r15,%r10,4),%xmm4
xor    %eax,%eax
ucomiss %xmm4,%xmm0
seta   %bl
ja     122a <frelon_lines+0x12a>
movaps %xmm2,%xmm4
movss  (%r15,%r10,4),%xmm5
mov    %bl,%al
add    %eax,%r9d
xor    %eax,%eax
ucomiss %xmm5,%xmm0
seta   %bl
ja     1200 <frelon_lines+0x100>
movaps %xmm2,%xmm5
jmp    1200 <frelon_lines+0x100>
mov    %rsi,%rbx
movaps %xmm1,%xmm3
shufps $0x0,%xmm1,%xmm3
cmpq   $0x0,-0x10(%rsp)
je     130e <frelon_lines+0x20e>
mov    -0x20(%rsp),%rsi
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movups -0x30(%rcx,%rax,4),%xmm4
movups -0x20(%rcx,%rax,4),%xmm5
movups -0x10(%rcx,%rax,4),%xmm6
movups (%rcx,%rax,4),%xmm7
subps  %xmm3,%xmm4
subps  %xmm3,%xmm5
movups %xmm4,-0x30(%rcx,%rax,4)
movups %xmm5,-0x20(%rcx,%rax,4)
subps  %xmm3,%xmm6
subps  %xmm3,%xmm7
movups %xmm6,-0x10(%rcx,%rax,4)
movups %xmm7,(%rcx,%rax,4)
add    $0x10,%rax
add    $0xfffffffffffffffe,%rsi
jne    1270 <frelon_lines+0x170>
testb  $0x1,-0x18(%rsp)
je     12d3 <frelon_lines+0x1d3>
add    %r14,%rax
mov    -0x28(%rsp),%rsi
movups (%rsi,%rax,4),%xmm4
movups 0x10(%rsi,%rax,4),%xmm5
subps  %xmm3,%xmm4
subps  %xmm3,%xmm5
movups %xmm4,(%rsi,%rax,4)
movups %xmm5,0x10(%rsi,%rax,4)
mov    -0x8(%rsp),%rsi
mov    %rsi,%rax
cmp    %r11,%rsi
mov    %rbx,%rsi
je     1180 <frelon_lines+0x80>
nopw   0x0(%rax,%rax,1)
movss  0x0(%rbp,%rax,4),%xmm3
subss  %xmm1,%xmm3
movss  %xmm3,0x0(%rbp,%rax,4)
add    $0x1,%rax
cmp    %rax,%r11
jne    12f0 <frelon_lines+0x1f0>
jmp    1180 <frelon_lines+0x80>
xor    %eax,%eax
testb  $0x1,-0x18(%rsp)
jne    12b3 <frelon_lines+0x1b3>
jmp    12d3 <frelon_lines+0x1d3>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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