void frelon_lines_sub(float *__restrict img, float *__restrict drk, 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++) {
      img[p + j] = img[p + j] - drk[p + 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_sub>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movss  %xmm0,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    126b <frelon_lines_sub+0x16b>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11f1 <frelon_lines_sub+0xf1>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x1c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    11de <frelon_lines_sub+0xde>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11e3 <frelon_lines_sub+0xe3>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1157 <frelon_lines_sub+0x57>
cmpl   $0x0,-0x2c(%rbp)
jle    120e <frelon_lines_sub+0x10e>
movss  -0x30(%rbp),%xmm0
cvtsi2ssl -0x2c(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1258 <frelon_lines_sub+0x158>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x34(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1215 <frelon_lines_sub+0x115>
jmp    125d <frelon_lines_sub+0x15d>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112b <frelon_lines_sub+0x2b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<frelon_lines_sub>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,-0x38(%rsp)
test   %edx,%edx
jle    1357 <frelon_lines_sub+0x257>
mov    %rdi,%r10
movss  (%rdi),%xmm1
movslq %ecx,%rsi
mov    %edx,%eax
mov    %rax,-0x8(%rsp)
mov    %esi,%r15d
mov    %r15d,%eax
and    $0xfffffff8,%eax
mov    %rax,-0x18(%rsp)
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x20(%rsp)
mov    %rax,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
mov    %r15d,%r13d
and    $0xfffffffe,%r13d
mov    %rdx,%rax
mov    %rdx,-0x28(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x30(%rsp)
mov    -0x38(%rsp),%rax
lea    0x4(%rax),%rdx
mov    %rsi,-0x10(%rsp)
lea    0x0(,%rsi,4),%r12
lea    0x4(%rdi),%rax
lea    0x30(%rdi),%rbp
xor    %r14d,%r14d
movss  0xe75(%rip),%xmm2
jmp    11ab <frelon_lines_sub+0xab>
nopl   (%rax)
add    $0x1,%r14
add    %r12,%rdx
add    %r12,%rax
add    %r12,%rbp
add    %r12,%rdi
cmp    -0x8(%rsp),%r14
je     1357 <frelon_lines_sub+0x257>
test   %ecx,%ecx
jle    1190 <frelon_lines_sub+0x90>
xor    %r9d,%r9d
xorps  %xmm3,%xmm3
xor    %r11d,%r11d
cmp    $0x1,%ecx
jne    124b <frelon_lines_sub+0x14b>
mov    %r14,%r8
imul   -0x10(%rsp),%r8
test   $0x1,%r15b
je     1201 <frelon_lines_sub+0x101>
add    %r8,%r9
movss  (%r10,%r9,4),%xmm4
mov    -0x38(%rsp),%rsi
subss  (%rsi,%r9,4),%xmm4
movss  %xmm4,(%r10,%r9,4)
xor    %esi,%esi
ucomiss %xmm4,%xmm0
seta   %bl
ja     11f7 <frelon_lines_sub+0xf7>
movaps %xmm2,%xmm4
addss  %xmm4,%xmm3
mov    %bl,%sil
add    %esi,%r11d
test   %r11d,%r11d
je     1215 <frelon_lines_sub+0x115>
xorps  %xmm1,%xmm1
cvtsi2ss %r11d,%xmm1
divss  %xmm1,%xmm3
movaps %xmm3,%xmm1
test   %ecx,%ecx
jle    1190 <frelon_lines_sub+0x90>
cmp    $0x8,%ecx
jae    1294 <frelon_lines_sub+0x194>
xor    %esi,%esi
jmp    1330 <frelon_lines_sub+0x230>
nopl   0x0(%rax)
addss  %xmm5,%xmm3
addss  %xmm4,%xmm3
mov    %bl,%sil
add    %esi,%r11d
add    $0x2,%r9
cmp    %r9,%r13
je     11c1 <frelon_lines_sub+0xc1>
movss  -0x4(%rax,%r9,4),%xmm5
movss  (%rax,%r9,4),%xmm4
subss  -0x4(%rdx,%r9,4),%xmm5
movss  %xmm5,-0x4(%rax,%r9,4)
xor    %esi,%esi
ucomiss %xmm5,%xmm0
seta   %bl
ja     1273 <frelon_lines_sub+0x173>
movaps %xmm2,%xmm5
mov    %bl,%sil
add    %esi,%r11d
subss  (%rdx,%r9,4),%xmm4
movss  %xmm4,(%rax,%r9,4)
xor    %esi,%esi
ucomiss %xmm4,%xmm0
seta   %bl
ja     1230 <frelon_lines_sub+0x130>
movaps %xmm2,%xmm4
jmp    1230 <frelon_lines_sub+0x130>
movaps %xmm1,%xmm3
shufps $0x0,%xmm1,%xmm3
cmpq   $0x0,-0x20(%rsp)
je     134c <frelon_lines_sub+0x24c>
mov    -0x30(%rsp),%rbx
xor    %esi,%esi
xchg   %ax,%ax
movups -0x30(%rbp,%rsi,4),%xmm4
movups -0x20(%rbp,%rsi,4),%xmm5
movups -0x10(%rbp,%rsi,4),%xmm6
movups 0x0(%rbp,%rsi,4),%xmm7
subps  %xmm3,%xmm4
subps  %xmm3,%xmm5
movups %xmm4,-0x30(%rbp,%rsi,4)
movups %xmm5,-0x20(%rbp,%rsi,4)
subps  %xmm3,%xmm6
subps  %xmm3,%xmm7
movups %xmm6,-0x10(%rbp,%rsi,4)
movups %xmm7,0x0(%rbp,%rsi,4)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rbx
jne    12b0 <frelon_lines_sub+0x1b0>
testb  $0x1,-0x28(%rsp)
je     1314 <frelon_lines_sub+0x214>
add    %r8,%rsi
movups (%r10,%rsi,4),%xmm4
movups 0x10(%r10,%rsi,4),%xmm5
subps  %xmm3,%xmm4
subps  %xmm3,%xmm5
movups %xmm4,(%r10,%rsi,4)
movups %xmm5,0x10(%r10,%rsi,4)
mov    -0x18(%rsp),%rbx
mov    %rbx,%rsi
cmp    %r15,%rbx
je     1190 <frelon_lines_sub+0x90>
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rdi,%rsi,4),%xmm3
subss  %xmm1,%xmm3
movss  %xmm3,(%rdi,%rsi,4)
add    $0x1,%rsi
cmp    %rsi,%r15
jne    1330 <frelon_lines_sub+0x230>
jmp    1190 <frelon_lines_sub+0x90>
xor    %esi,%esi
testb  $0x1,-0x28(%rsp)
jne    12f5 <frelon_lines_sub+0x1f5>
jmp    1314 <frelon_lines_sub+0x214>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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