int prjbchk(double tol, int nphi, int ntheta, int spt, double phi[],
            double theta[], int stat[]) {
  int status = 0;
  register int iphi, itheta, *statp;
  register double *phip, *thetap;
  phip = phi;
  thetap = theta;
  statp = stat;
  for (itheta = 0; itheta < ntheta; itheta++) {
    for (iphi = 0; iphi < nphi; iphi++, phip += spt, thetap += spt, statp++) {
      if (*statp == 0) {
        if (*phip < -180.0) {
          if (*phip < -180.0 - tol) {
            *statp = 1;
            status = 1;
          } else {
            *phip = -180.0;
          }
        } else if (180.0 < *phip) {
          if (180.0 + tol < *phip) {
            *statp = 1;
            status = 1;
          } else {
            *phip = 180.0;
          }
        }
        if (*thetap < -90.0) {
          if (*thetap < -90.0 - tol) {
            *statp = 1;
            status = 1;
          } else {
            *thetap = -90.0;
          }
        } else if (90.0 < *thetap) {
          if (90.0 + tol < *thetap) {
            *statp = 1;
            status = 1;
          } else {
            *thetap = 90.0;
          }
        }
      }
    }
  }
  return status;
}


x86 assembly (clang, -O0):
<prjbchk>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x20(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x30(%rbp),%rax
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1363 <prjbchk+0x263>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1350 <prjbchk+0x250>
mov    -0x48(%rbp),%rax
cmpl   $0x0,(%rax)
jne    1307 <prjbchk+0x207>
mov    -0x50(%rbp),%rax
movsd  0xe84(%rip),%xmm0
ucomisd (%rax),%xmm0
jbe    11d0 <prjbchk+0xd0>
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm1
movsd  0xe6a(%rip),%xmm0
subsd  -0x8(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11bb <prjbchk+0xbb>
mov    -0x48(%rbp),%rax
movl   $0x1,(%rax)
movl   $0x1,-0x34(%rbp)
jmp    11cb <prjbchk+0xcb>
mov    -0x50(%rbp),%rax
movsd  0xe39(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    1239 <prjbchk+0x139>
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xe28(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1234 <prjbchk+0x134>
movsd  0xe16(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    121f <prjbchk+0x11f>
mov    -0x48(%rbp),%rax
movl   $0x1,(%rax)
movl   $0x1,-0x34(%rbp)
jmp    122f <prjbchk+0x12f>
mov    -0x50(%rbp),%rax
movsd  0xddd(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    1234 <prjbchk+0x134>
jmp    1239 <prjbchk+0x139>
mov    -0x58(%rbp),%rax
movsd  0xdcb(%rip),%xmm0
ucomisd (%rax),%xmm0
jbe    1299 <prjbchk+0x199>
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm1
movsd  0xdb1(%rip),%xmm0
subsd  -0x8(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1284 <prjbchk+0x184>
mov    -0x48(%rbp),%rax
movl   $0x1,(%rax)
movl   $0x1,-0x34(%rbp)
jmp    1294 <prjbchk+0x194>
mov    -0x58(%rbp),%rax
movsd  0xd80(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    1302 <prjbchk+0x202>
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xd6f(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12fd <prjbchk+0x1fd>
movsd  0xd5d(%rip),%xmm1
addsd  -0x8(%rbp),%xmm1
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    12e8 <prjbchk+0x1e8>
mov    -0x48(%rbp),%rax
movl   $0x1,(%rax)
movl   $0x1,-0x34(%rbp)
jmp    12f8 <prjbchk+0x1f8>
mov    -0x58(%rbp),%rax
movsd  0xd24(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    12fd <prjbchk+0x1fd>
jmp    1302 <prjbchk+0x202>
jmp    1307 <prjbchk+0x207>
jmp    130c <prjbchk+0x20c>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0x50(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0x58(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x48(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x48(%rbp)
jmp    1157 <prjbchk+0x57>
jmp    1355 <prjbchk+0x255>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1144 <prjbchk+0x44>
mov    -0x34(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<prjbchk>:
push   %rbp
push   %r15
push   %r14
push   %rbx
test   %esi,%esi
jle    123b <prjbchk+0x13b>
movsd  0xee9(%rip),%xmm8
addsd  %xmm0,%xmm8
movsd  0xee4(%rip),%xmm2
movapd %xmm2,%xmm3
subsd  %xmm0,%xmm3
movsd  0xedc(%rip),%xmm4
addsd  %xmm0,%xmm4
movsd  0xed8(%rip),%xmm5
movapd %xmm5,%xmm6
subsd  %xmm0,%xmm6
movslq %edx,%rdx
shl    $0x3,%rdx
xor    %r10d,%r10d
movabs $0xc066800000000000,%r11
80 66 c0
movabs $0xc056800000000000,%rbx
80 56 c0
movsd  0xea2(%rip),%xmm0
movabs $0x4056800000000000,%r14
80 56 40
movsd  0xe80(%rip),%xmm7
movabs $0x4066800000000000,%r15
80 66 40
xor    %eax,%eax
jmp    119d <prjbchk+0x9d>
xchg   %ax,%ax
add    $0x1,%r10d
cmp    %esi,%r10d
je     123d <prjbchk+0x13d>
test   %edi,%edi
jle    1190 <prjbchk+0x90>
mov    %edi,%ebp
jmp    11cf <prjbchk+0xcf>
movl   $0x1,(%r9)
mov    $0x1,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x4,%r9
add    %rdx,%r8
add    %rdx,%rcx
add    $0xffffffff,%ebp
je     1190 <prjbchk+0x90>
cmpl   $0x0,(%r9)
jne    11c0 <prjbchk+0xc0>
movsd  (%rcx),%xmm1
ucomisd %xmm1,%xmm2
jbe    11ea <prjbchk+0xea>
ucomisd %xmm1,%xmm3
ja     11f7 <prjbchk+0xf7>
mov    %r11,(%rcx)
jmp    1210 <prjbchk+0x110>
ucomisd %xmm7,%xmm1
jbe    1210 <prjbchk+0x110>
ucomisd %xmm8,%xmm1
jbe    1205 <prjbchk+0x105>
movl   $0x1,(%r9)
mov    $0x1,%eax
jmp    1210 <prjbchk+0x110>
mov    %r15,(%rcx)
nopl   0x0(%rax,%rax,1)
movsd  (%r8),%xmm1
ucomisd %xmm1,%xmm5
jbe    1226 <prjbchk+0x126>
ucomisd %xmm1,%xmm6
ja     11a5 <prjbchk+0xa5>
mov    %rbx,(%r8)
jmp    11c0 <prjbchk+0xc0>
ucomisd %xmm0,%xmm1
jbe    11c0 <prjbchk+0xc0>
ucomisd %xmm4,%xmm1
ja     11a5 <prjbchk+0xa5>
mov    %r14,(%r8)
jmp    11c0 <prjbchk+0xc0>
xor    %eax,%eax
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @prjbchk(double %0, i32 %1, i32 %2, i32 %3, double* nocapture %4, double* nocapture %5, i32* nocapture %6) {
%8 = icmp sgt i32 %2, 0
br i1 %8, label %9, label %71
9: ; preds = %7
%10 = icmp sgt i32 %1, 0
%11 = fadd double %0, 1.800000e+02
%12 = fsub double -1.800000e+02, %0
%13 = fadd double %0, 9.000000e+01
%14 = fsub double -9.000000e+01, %0
%15 = sext i32 %3 to i64
br label %16
16: ; preds = %9, %64
%17 = phi double* [ %5, %9 ], [ %68, %64 ]
%18 = phi double* [ %4, %9 ], [ %67, %64 ]
%19 = phi i32* [ %6, %9 ], [ %66, %64 ]
%20 = phi i32 [ 0, %9 ], [ %69, %64 ]
%21 = phi i32 [ 0, %9 ], [ %65, %64 ]
br i1 %10, label %22, label %64
22: ; preds = %16, %57
%23 = phi double* [ %61, %57 ], [ %17, %16 ]
%24 = phi double* [ %60, %57 ], [ %18, %16 ]
%25 = phi i32* [ %62, %57 ], [ %19, %16 ]
%26 = phi i32 [ %59, %57 ], [ 0, %16 ]
%27 = phi i32 [ %58, %57 ], [ %21, %16 ]
%28 = load i32, i32* %25
%29 = icmp eq i32 %28, 0
br i1 %29, label %30, label %57
30: ; preds = %22
%31 = load double, double* %24
%32 = fcmp olt double %31, -1.800000e+02
br i1 %32, label %33, label %37
33: ; preds = %30
%34 = fcmp olt double %31, %12
br i1 %34, label %35, label %36
35: ; preds = %33
store i32 1, i32* %25
br label %43
36: ; preds = %33
store double -1.800000e+02, double* %24
br label %43
37: ; preds = %30
%38 = fcmp ogt double %31, 1.800000e+02
br i1 %38, label %39, label %43
39: ; preds = %37
%40 = fcmp olt double %11, %31
br i1 %40, label %41, label %42
41: ; preds = %39
store i32 1, i32* %25
br label %43
42: ; preds = %39
store double 1.800000e+02, double* %24
br label %43
43: ; preds = %37, %42, %41, %35, %36
%44 = phi i32 [ 1, %35 ], [ %27, %36 ], [ 1, %41 ], [ %27, %42 ], [ %27, %37 ]
%45 = load double, double* %23
%46 = fcmp olt double %45, -9.000000e+01
br i1 %46, label %47, label %51
47: ; preds = %43
%48 = fcmp olt double %45, %14
br i1 %48, label %49, label %50
49: ; preds = %47
store i32 1, i32* %25
br label %57
50: ; preds = %47
store double -9.000000e+01, double* %23
br label %57
51: ; preds = %43
%52 = fcmp ogt double %45, 9.000000e+01
br i1 %52, label %53, label %57
53: ; preds = %51
%54 = fcmp olt double %13, %45
br i1 %54, label %55, label %56
55: ; preds = %53
store i32 1, i32* %25
br label %57
56: ; preds = %53
store double 9.000000e+01, double* %23
br label %57
57: ; preds = %22, %51, %56, %55, %49, %50
%58 = phi i32 [ 1, %49 ], [ %44, %50 ], [ 1, %55 ], [ %44, %56 ], [ %44, %51 ], [ %27, %22 ]
%59 = add nuw nsw i32 %26, 1
%60 = getelementptr inbounds double, double* %24, i64 %15
%61 = getelementptr inbounds double, double* %23, i64 %15
%62 = getelementptr inbounds i32, i32* %25, i64 1
%63 = icmp eq i32 %59, %1
br i1 %63, label %64, label %22
64: ; preds = %57, %16
%65 = phi i32 [ %21, %16 ], [ %58, %57 ]
%66 = phi i32* [ %19, %16 ], [ %62, %57 ]
%67 = phi double* [ %18, %16 ], [ %60, %57 ]
%68 = phi double* [ %17, %16 ], [ %61, %57 ]
%69 = add nuw nsw i32 %20, 1
%70 = icmp eq i32 %69, %2
br i1 %70, label %71, label %16
71: ; preds = %64, %7
%72 = phi i32 [ 0, %7 ], [ %65, %64 ]
ret i32 %72
}