int rejectPixels(data, flat, normx, badpix, npix, sumxsqr, sumxz, sumx, sumz,
                 threshold, ngrow)
float *data;
float *flat;
float *normx;
char *badpix;
int npix;
double *sumxsqr, *sumxz, *sumx, *sumz;
double threshold;
int ngrow;
{
  int ngoodpix, i, j;
  float residual, lcut, hcut;
  double x, z;
  ngoodpix = npix;
  lcut = -threshold;
  hcut = threshold;
  for (i = 0; i < npix; i++) {
    if (badpix[i] == 1)
      ngoodpix = ngoodpix - 1;
    else {
      residual = flat[i];
      if (residual < lcut || residual > hcut) {
        for (j = ((0) > (i - ngrow) ? (0) : (i - ngrow));
             j < ((npix) < (i + ngrow) ? (npix) : (i + ngrow)); j++) {
          if (badpix[j] != 1) {
            if (j <= i) {
              x = (double)normx[j];
              z = (double)data[j];
              *sumxsqr = *sumxsqr - (x * x);
              *sumxz = *sumxz - z * x;
              *sumx = *sumx - x;
              *sumz = *sumz - z;
              badpix[j] = 1;
              ngoodpix = ngoodpix - 1;
            } else
              badpix[j] = 2;
          }
        }
      }
    }
  }
  return (ngoodpix);
}


x86 assembly (clang, -O0):
<rejectPixels>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movsd  -0x38(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
movsd  -0x38(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1379 <rejectPixels+0x279>
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x1,%eax
jne    119f <rejectPixels+0x9f>
mov    -0x3c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1366 <rejectPixels+0x266>
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm1
movss  -0x4c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
ja     11d3 <rejectPixels+0xd3>
movss  -0x48(%rbp),%xmm0
ucomiss -0x50(%rbp),%xmm0
jbe    1361 <rejectPixels+0x261>
mov    -0x40(%rbp),%ecx
sub    0x28(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    11ed <rejectPixels+0xed>
xor    %eax,%eax
mov    %eax,-0x64(%rbp)
jmp    11f6 <rejectPixels+0xf6>
mov    -0x40(%rbp),%eax
sub    0x28(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x40(%rbp),%ecx
add    0x28(%rbp),%ecx
cmp    %ecx,%eax
jge    121e <rejectPixels+0x11e>
mov    -0x24(%rbp),%eax
mov    %eax,-0x6c(%rbp)
jmp    1227 <rejectPixels+0x127>
mov    -0x40(%rbp),%eax
add    0x28(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x68(%rbp),%eax
mov    -0x6c(%rbp),%ecx
cmp    %ecx,%eax
jge    135c <rejectPixels+0x25c>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x1,%eax
je     1349 <rejectPixels+0x249>
mov    -0x44(%rbp),%eax
cmp    -0x40(%rbp),%eax
jg     1338 <rejectPixels+0x238>
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    0x10(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x60(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm0
subsd  -0x58(%rbp),%xmm0
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    0x20(%rbp),%rax
movsd  (%rax),%xmm0
subsd  -0x60(%rbp),%xmm0
mov    0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movb   $0x1,(%rax,%rcx,1)
mov    -0x3c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1344 <rejectPixels+0x244>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movb   $0x2,(%rax,%rcx,1)
jmp    1349 <rejectPixels+0x249>
jmp    134e <rejectPixels+0x24e>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11fc <rejectPixels+0xfc>
jmp    1361 <rejectPixels+0x261>
jmp    1366 <rejectPixels+0x266>
jmp    136b <rejectPixels+0x26b>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1170 <rejectPixels+0x70>
mov    -0x3c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<rejectPixels>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8d,%eax
test   %r8d,%r8d
jle    125e <rejectPixels+0x15e>
mov    0x50(%rsp),%ebx
mov    0x48(%rsp),%r15
mov    0x40(%rsp),%r13
cvtsd2ss %xmm0,%xmm0
movaps 0xed1(%rip),%xmm1
xorps  %xmm0,%xmm1
mov    %ebx,%r14d
neg    %r14d
mov    %eax,%r8d
xor    %ebx,%ebx
mov    %eax,-0xc(%rsp)
mov    %r8,-0x8(%rsp)
jmp    1164 <rejectPixels+0x64>
nopl   0x0(%rax,%rax,1)
add    $0xffffffff,%eax
add    $0x1,%rbx
add    $0x1,%r14d
cmp    %r8,%rbx
je     125e <rejectPixels+0x15e>
test   %r14d,%r14d
mov    $0x0,%r12d
cmovg  %r14d,%r12d
cmpb   $0x1,(%rcx,%rbx,1)
je     1150 <rejectPixels+0x50>
movss  (%rsi,%rbx,4),%xmm2
ucomiss %xmm2,%xmm1
ja     1186 <rejectPixels+0x86>
ucomiss %xmm0,%xmm2
jbe    1153 <rejectPixels+0x53>
mov    %ebx,%r11d
mov    0x50(%rsp),%ebp
mov    %rbp,%r8
sub    %r8d,%r11d
mov    $0x0,%ebp
cmovs  %ebp,%r11d
lea    (%rbx,%r8,1),%r10d
mov    -0xc(%rsp),%r8d
cmp    %r8d,%r10d
cmovg  %r8d,%r10d
cmp    %r10d,%r11d
jge    1254 <rejectPixels+0x154>
movslq %r10d,%r10
mov    0x38(%rsp),%rbp
mov    -0x8(%rsp),%r8
jmp    11e2 <rejectPixels+0xe2>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movb   $0x2,(%rcx,%r12,1)
add    $0x1,%r12
cmp    %r10,%r12
jge    1153 <rejectPixels+0x53>
cmpb   $0x1,(%rcx,%r12,1)
je     11d5 <rejectPixels+0xd5>
cmp    %rbx,%r12
ja     11d0 <rejectPixels+0xd0>
movss  (%rdx,%r12,4),%xmm2
cvtss2sd %xmm2,%xmm2
movss  (%rdi,%r12,4),%xmm3
cvtss2sd %xmm3,%xmm3
movsd  (%r9),%xmm4
movaps %xmm2,%xmm5
mulsd  %xmm2,%xmm5
subsd  %xmm5,%xmm4
movsd  %xmm4,(%r9)
movsd  0x0(%rbp),%xmm4
movaps %xmm3,%xmm5
mulsd  %xmm2,%xmm5
subsd  %xmm5,%xmm4
movsd  %xmm4,0x0(%rbp)
movsd  0x0(%r13),%xmm4
subsd  %xmm2,%xmm4
movsd  %xmm4,0x0(%r13)
movsd  (%r15),%xmm2
subsd  %xmm3,%xmm2
movsd  %xmm2,(%r15)
movb   $0x1,(%rcx,%r12,1)
add    $0xffffffff,%eax
jmp    11d5 <rejectPixels+0xd5>
mov    -0x8(%rsp),%r8
jmp    1153 <rejectPixels+0x53>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @rejectPixels(float* %0, float* %1, float* %2, i8* %3, i32 %4, double* %5, double* %6, double* %7, double* %8, double %9, i32 %10) {
%12 = alloca float*
%13 = alloca float*
%14 = alloca float*
%15 = alloca i8*
%16 = alloca i32
%17 = alloca double*
%18 = alloca double*
%19 = alloca double*
%20 = alloca double*
%21 = alloca double
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca float
%27 = alloca float
%28 = alloca float
%29 = alloca double
%30 = alloca double
store float* %0, float** %12
store float* %1, float** %13
store float* %2, float** %14
store i8* %3, i8** %15
store i32 %4, i32* %16
store double* %5, double** %17
store double* %6, double** %18
store double* %7, double** %19
store double* %8, double** %20
store double %9, double* %21
store i32 %10, i32* %22
%31 = load i32, i32* %16
store i32 %31, i32* %23
%32 = load double, double* %21
%33 = fneg double %32
%34 = fptrunc double %33 to float
store float %34, float* %27
%35 = load double, double* %21
%36 = fptrunc double %35 to float
store float %36, float* %28
store i32 0, i32* %24
br label %37
37: ; preds = %161, %11
%38 = load i32, i32* %24
%39 = load i32, i32* %16
%40 = icmp slt i32 %38, %39
br i1 %40, label %41, label %164
41: ; preds = %37
%42 = load i8*, i8** %15
%43 = load i32, i32* %24
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds i8, i8* %42, i64 %44
%46 = load i8, i8* %45
%47 = sext i8 %46 to i32
%48 = icmp eq i32 %47, 1
br i1 %48, label %49, label %52
49: ; preds = %41
%50 = load i32, i32* %23
%51 = sub nsw i32 %50, 1
store i32 %51, i32* %23
br label %160
52: ; preds = %41
%53 = load float*, float** %13
%54 = load i32, i32* %24
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds float, float* %53, i64 %55
%57 = load float, float* %56
store float %57, float* %26
%58 = load float, float* %26
%59 = load float, float* %27
%60 = fcmp olt float %58, %59
br i1 %60, label %65, label %61
61: ; preds = %52
%62 = load float, float* %26
%63 = load float, float* %28
%64 = fcmp ogt float %62, %63
br i1 %64, label %65, label %159
65: ; preds = %61, %52
%66 = load i32, i32* %24
%67 = load i32, i32* %22
%68 = sub nsw i32 %66, %67
%69 = icmp sgt i32 0, %68
br i1 %69, label %70, label %71
70: ; preds = %65
br label %75
71: ; preds = %65
%72 = load i32, i32* %24
%73 = load i32, i32* %22
%74 = sub nsw i32 %72, %73
br label %75
75: ; preds = %71, %70
%76 = phi i32 [ 0, %70 ], [ %74, %71 ]
store i32 %76, i32* %25
br label %77
77: ; preds = %155, %75
%78 = load i32, i32* %25
%79 = load i32, i32* %16
%80 = load i32, i32* %24
%81 = load i32, i32* %22
%82 = add nsw i32 %80, %81
%83 = icmp slt i32 %79, %82
br i1 %83, label %84, label %86
84: ; preds = %77
%85 = load i32, i32* %16
br label %90
86: ; preds = %77
%87 = load i32, i32* %24
%88 = load i32, i32* %22
%89 = add nsw i32 %87, %88
br label %90
90: ; preds = %86, %84
%91 = phi i32 [ %85, %84 ], [ %89, %86 ]
%92 = icmp slt i32 %78, %91
br i1 %92, label %93, label %158
93: ; preds = %90
%94 = load i8*, i8** %15
%95 = load i32, i32* %25
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds i8, i8* %94, i64 %96
%98 = load i8, i8* %97
%99 = sext i8 %98 to i32
%100 = icmp ne i32 %99, 1
br i1 %100, label %101, label %154
101: ; preds = %93
%102 = load i32, i32* %25
%103 = load i32, i32* %24
%104 = icmp sle i32 %102, %103
br i1 %104, label %105, label %148
105: ; preds = %101
%106 = load float*, float** %14
%107 = load i32, i32* %25
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds float, float* %106, i64 %108
%110 = load float, float* %109
%111 = fpext float %110 to double
store double %111, double* %29
%112 = load float*, float** %12
%113 = load i32, i32* %25
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds float, float* %112, i64 %114
%116 = load float, float* %115
%117 = fpext float %116 to double
store double %117, double* %30
%118 = load double*, double** %17
%119 = load double, double* %118
%120 = load double, double* %29
%121 = load double, double* %29
%122 = fneg double %120
%123 = call double @llvm.fmuladd.f64(double %122, double %121, double %119)
%124 = load double*, double** %17
store double %123, double* %124
%125 = load double*, double** %18
%126 = load double, double* %125
%127 = load double, double* %30
%128 = load double, double* %29
%129 = fneg double %127
%130 = call double @llvm.fmuladd.f64(double %129, double %128, double %126)
%131 = load double*, double** %18
store double %130, double* %131
%132 = load double*, double** %19
%133 = load double, double* %132
%134 = load double, double* %29
%135 = fsub double %133, %134
%136 = load double*, double** %19
store double %135, double* %136
%137 = load double*, double** %20
%138 = load double, double* %137
%139 = load double, double* %30
%140 = fsub double %138, %139
%141 = load double*, double** %20
store double %140, double* %141
%142 = load i8*, i8** %15
%143 = load i32, i32* %25
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds i8, i8* %142, i64 %144
store i8 1, i8* %145
%146 = load i32, i32* %23
%147 = sub nsw i32 %146, 1
store i32 %147, i32* %23
br label %153
148: ; preds = %101
%149 = load i8*, i8** %15
%150 = load i32, i32* %25
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds i8, i8* %149, i64 %151
store i8 2, i8* %152
br label %153
153: ; preds = %148, %105
br label %154
154: ; preds = %153, %93
br label %155
155: ; preds = %154
%156 = load i32, i32* %25
%157 = add nsw i32 %156, 1
store i32 %157, i32* %25
br label %77
158: ; preds = %90
br label %159
159: ; preds = %158, %61
br label %160
160: ; preds = %159, %49
br label %161
161: ; preds = %160
%162 = load i32, i32* %24
%163 = add nsw i32 %162, 1
store i32 %163, i32* %24
br label %37
164: ; preds = %37
%165 = load i32, i32* %23
ret i32 %165
}

LLVM IR (-O2):
define i32 @rejectPixels(float* nocapture %0, float* nocapture %1, float* nocapture %2, i8* nocapture %3, i32 %4, double* nocapture %5, double* nocapture %6, double* nocapture %7, double* nocapture %8, double %9, i32 %10) {
%12 = fptrunc double %9 to float
%13 = fneg float %12
%14 = icmp sgt i32 %4, 0
br i1 %14, label %15, label %82
15: ; preds = %11
%16 = sub i32 0, %10
%17 = zext i32 %4 to i64
br label %18
18: ; preds = %15, %77
%19 = phi i64 [ 0, %15 ], [ %79, %77 ]
%20 = phi i32 [ %16, %15 ], [ %80, %77 ]
%21 = phi i32 [ %4, %15 ], [ %78, %77 ]
%22 = call i32 @llvm.smax.i32(i32 %20, i32 0)
%23 = zext i32 %22 to i64
%24 = getelementptr inbounds i8, i8* %3, i64 %19
%25 = load i8, i8* %24
%26 = icmp eq i8 %25, 1
br i1 %26, label %27, label %29
27: ; preds = %18
%28 = add nsw i32 %21, -1
br label %77
29: ; preds = %18
%30 = getelementptr inbounds float, float* %1, i64 %19
%31 = load float, float* %30
%32 = fcmp olt float %31, %13
%33 = fcmp ogt float %31, %12
%34 = or i1 %32, %33
br i1 %34, label %35, label %77
35: ; preds = %29
%36 = trunc i64 %19 to i32
%37 = sub nsw i32 %36, %10
%38 = icmp sgt i32 %37, 0
%39 = select i1 %38, i32 %37, i32 0
%40 = add nsw i32 %36, %10
%41 = icmp sgt i32 %40, %4
%42 = select i1 %41, i32 %4, i32 %40
%43 = icmp slt i32 %39, %42
br i1 %43, label %44, label %77
44: ; preds = %35
%45 = sext i32 %42 to i64
br label %46
46: ; preds = %44, %73
%47 = phi i64 [ %23, %44 ], [ %75, %73 ]
%48 = phi i32 [ %21, %44 ], [ %74, %73 ]
%49 = getelementptr inbounds i8, i8* %3, i64 %47
%50 = load i8, i8* %49
%51 = icmp eq i8 %50, 1
br i1 %51, label %73, label %52
52: ; preds = %46
%53 = icmp ugt i64 %47, %19
br i1 %53, label %72, label %54
54: ; preds = %52
%55 = getelementptr inbounds float, float* %2, i64 %47
%56 = load float, float* %55
%57 = fpext float %56 to double
%58 = getelementptr inbounds float, float* %0, i64 %47
%59 = load float, float* %58
%60 = fpext float %59 to double
%61 = load double, double* %5
%62 = fneg double %57
%63 = tail call double @llvm.fmuladd.f64(double %62, double %57, double %61)
store double %63, double* %5
%64 = load double, double* %6
%65 = fneg double %60
%66 = tail call double @llvm.fmuladd.f64(double %65, double %57, double %64)
store double %66, double* %6
%67 = load double, double* %7
%68 = fsub double %67, %57
store double %68, double* %7
%69 = load double, double* %8
%70 = fsub double %69, %60
store double %70, double* %8
store i8 1, i8* %49
%71 = add nsw i32 %48, -1
br label %73
72: ; preds = %52
store i8 2, i8* %49
br label %73
73: ; preds = %46, %72, %54
%74 = phi i32 [ %71, %54 ], [ %48, %72 ], [ %48, %46 ]
%75 = add nuw nsw i64 %47, 1
%76 = icmp slt i64 %75, %45
br i1 %76, label %46, label %77
77: ; preds = %73, %35, %29, %27
%78 = phi i32 [ %28, %27 ], [ %21, %29 ], [ %21, %35 ], [ %74, %73 ]
%79 = add nuw nsw i64 %19, 1
%80 = add i32 %20, 1
%81 = icmp eq i64 %79, %17
br i1 %81, label %82, label %18
82: ; preds = %77, %11
%83 = phi i32 [ %4, %11 ], [ %78, %77 ]
ret i32 %83
}