int restore_viirs_bowtie(float **img, float **mask, int nx, int ny,
                         float fillval) {
  int ix, ntot = 0;
  if (nx != 3200) {
    return -1;
  }
  for (ix = 0; ix < nx * ny; ix++) {
    if (fabs((mask[0][ix]) - fillval) > 0.001) {
      img[0][ix] = mask[0][ix];
      ntot++;
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<restore_viirs_bowtie>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
movss  %xmm0,-0x24(%rbp)
movl   $0x0,-0x2c(%rbp)
cmpl   $0xc80,-0x1c(%rbp)
je     1137 <restore_viirs_bowtie+0x37>
movl   $0xffffffff,-0x4(%rbp)
jmp    11c9 <restore_viirs_bowtie+0xc9>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x1c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
cmp    %ecx,%eax
jge    11c2 <restore_viirs_bowtie+0xc2>
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xea0(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0xe84(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11af <restore_viirs_bowtie+0xaf>
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b4 <restore_viirs_bowtie+0xb4>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113e <restore_viirs_bowtie+0x3e>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<restore_viirs_bowtie>:
mov    $0xffffffff,%eax
cmp    $0xc80,%edx
jne    115f <restore_viirs_bowtie+0x5f>
test   %ecx,%ecx
jle    1160 <restore_viirs_bowtie+0x60>
imul   $0xc80,%ecx,%edx
mov    (%rsi),%rcx
xor    %eax,%eax
mov    $0x1,%r8d
cmp    $0x2,%edx
cmovge %edx,%r8d
jge    1163 <restore_viirs_bowtie+0x63>
xor    %edx,%edx
test   $0x1,%r8b
je     115f <restore_viirs_bowtie+0x5f>
movss  (%rcx,%rdx,4),%xmm1
movaps %xmm1,%xmm2
subss  %xmm0,%xmm2
andps  0xeba(%rip),%xmm2
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
ucomisd 0xebb(%rip),%xmm0
jbe    115f <restore_viirs_bowtie+0x5f>
mov    (%rdi),%rcx
movss  %xmm1,(%rcx,%rdx,4)
ret
xor    %eax,%eax
ret
mov    %r8d,%r9d
and    $0x7fffff80,%r9d
xor    %edx,%edx
movaps 0xe8a(%rip),%xmm1
movsd  0xe92(%rip),%xmm2
jmp    1189 <restore_viirs_bowtie+0x89>
add    $0x2,%rdx
cmp    %rdx,%r9
je     112d <restore_viirs_bowtie+0x2d>
movss  (%rcx,%rdx,4),%xmm3
movaps %xmm3,%xmm4
subss  %xmm0,%xmm4
andps  %xmm1,%xmm4
cvtss2sd %xmm4,%xmm4
ucomisd %xmm2,%xmm4
jbe    11aa <restore_viirs_bowtie+0xaa>
mov    (%rdi),%rsi
movss  %xmm3,(%rsi,%rdx,4)
movss  0x4(%rcx,%rdx,4),%xmm3
movaps %xmm3,%xmm4
subss  %xmm0,%xmm4
andps  %xmm1,%xmm4
cvtss2sd %xmm4,%xmm4
ucomisd %xmm2,%xmm4
jbe    1180 <restore_viirs_bowtie+0x80>
mov    (%rdi),%rsi
movss  %xmm3,0x4(%rsi,%rdx,4)
jmp    1180 <restore_viirs_bowtie+0x80>

LLVM IR (-O0):
define i32 @restore_viirs_bowtie(float** %0, float** %1, i32 %2, i32 %3, float %4) {
%6 = alloca i32
%7 = alloca float**
%8 = alloca float**
%9 = alloca i32
%10 = alloca i32
%11 = alloca float
%12 = alloca i32
%13 = alloca i32
store float** %0, float*** %7
store float** %1, float*** %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store float %4, float* %11
store i32 0, i32* %13
%14 = load i32, i32* %9
%15 = icmp ne i32 %14, 3200
br i1 %15, label %16, label %17
16: ; preds = %5
store i32 -1, i32* %6
br label %58
17: ; preds = %5
store i32 0, i32* %12
br label %18
18: ; preds = %54, %17
%19 = load i32, i32* %12
%20 = load i32, i32* %9
%21 = load i32, i32* %10
%22 = mul nsw i32 %20, %21
%23 = icmp slt i32 %19, %22
br i1 %23, label %24, label %57
24: ; preds = %18
%25 = load float**, float*** %8
%26 = getelementptr inbounds float*, float** %25, i64 0
%27 = load float*, float** %26
%28 = load i32, i32* %12
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %27, i64 %29
%31 = load float, float* %30
%32 = load float, float* %11
%33 = fsub float %31, %32
%34 = fpext float %33 to double
%35 = call double @llvm.fabs.f64(double %34)
%36 = fcmp ogt double %35, 1.000000e-03
br i1 %36, label %37, label %53
37: ; preds = %24
%38 = load float**, float*** %8
%39 = getelementptr inbounds float*, float** %38, i64 0
%40 = load float*, float** %39
%41 = load i32, i32* %12
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %40, i64 %42
%44 = load float, float* %43
%45 = load float**, float*** %7
%46 = getelementptr inbounds float*, float** %45, i64 0
%47 = load float*, float** %46
%48 = load i32, i32* %12
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds float, float* %47, i64 %49
store float %44, float* %50
%51 = load i32, i32* %13
%52 = add nsw i32 %51, 1
store i32 %52, i32* %13
br label %53
53: ; preds = %37, %24
br label %54
54: ; preds = %53
%55 = load i32, i32* %12
%56 = add nsw i32 %55, 1
store i32 %56, i32* %12
br label %18
57: ; preds = %18
store i32 0, i32* %6
br label %58
58: ; preds = %57, %16
%59 = load i32, i32* %6
ret i32 %59
}

LLVM IR (-O2):
define i32 @restore_viirs_bowtie(float** nocapture %0, float** nocapture %1, i32 %2, i32 %3, float %4) {
%6 = icmp eq i32 %2, 3200
br i1 %6, label %7, label %58
7: ; preds = %5
%8 = icmp sgt i32 %3, 0
br i1 %8, label %9, label %58
9: ; preds = %7
%10 = mul i32 %3, 3200
%11 = load float*, float** %1
%12 = call i32 @llvm.smax.i32(i32 %10, i32 1)
%13 = zext i32 %12 to i64
%14 = and i64 %13, 1
%15 = icmp slt i32 %10, 2
br i1 %15, label %45, label %16
16: ; preds = %9
%17 = and i64 %13, 2147483520
br label %18
18: ; preds = %41, %16
%19 = phi i64 [ 0, %16 ], [ %42, %41 ]
%20 = phi i64 [ 0, %16 ], [ %43, %41 ]
%21 = getelementptr inbounds float, float* %11, i64 %19
%22 = load float, float* %21
%23 = fsub float %22, %4
%24 = tail call float @llvm.fabs.f32(float %23)
%25 = fpext float %24 to double
%26 = fcmp ogt double %25, 1.000000e-03
br i1 %26, label %27, label %30
27: ; preds = %18
%28 = load float*, float** %0
%29 = getelementptr inbounds float, float* %28, i64 %19
store float %22, float* %29
br label %30
30: ; preds = %18, %27
%31 = or i64 %19, 1
%32 = getelementptr inbounds float, float* %11, i64 %31
%33 = load float, float* %32
%34 = fsub float %33, %4
%35 = tail call float @llvm.fabs.f32(float %34)
%36 = fpext float %35 to double
%37 = fcmp ogt double %36, 1.000000e-03
br i1 %37, label %38, label %41
38: ; preds = %30
%39 = load float*, float** %0
%40 = getelementptr inbounds float, float* %39, i64 %31
store float %33, float* %40
br label %41
41: ; preds = %38, %30
%42 = add nuw nsw i64 %19, 2
%43 = add i64 %20, 2
%44 = icmp eq i64 %43, %17
br i1 %44, label %45, label %18
45: ; preds = %41, %9
%46 = phi i64 [ 0, %9 ], [ %42, %41 ]
%47 = icmp eq i64 %14, 0
br i1 %47, label %58, label %48
48: ; preds = %45
%49 = getelementptr inbounds float, float* %11, i64 %46
%50 = load float, float* %49
%51 = fsub float %50, %4
%52 = tail call float @llvm.fabs.f32(float %51)
%53 = fpext float %52 to double
%54 = fcmp ogt double %53, 1.000000e-03
br i1 %54, label %55, label %58
55: ; preds = %48
%56 = load float*, float** %0
%57 = getelementptr inbounds float, float* %56, i64 %46
store float %50, float* %57
br label %58
58: ; preds = %45, %55, %48, %7, %5
%59 = phi i32 [ -1, %5 ], [ 0, %7 ], [ 0, %48 ], [ 0, %55 ], [ 0, %45 ]
ret i32 %59
}