#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
float dist_calc(float, float, float, float, float, float);
int remove_unclustered(float *latitudes, float *longitudes, float *depths,
                       long n_locs, unsigned char *mask, float distance_cutoff,
                       int n_threads) {
  int out = 0;
  long i;
#pragma omp parallel for num_threads(n_threads)
  for (i = 0; i < n_locs; ++i) {
    long j;
    float dist;
    if (mask[i] != 0) {
      continue;
    }
    for (j = 0; j < n_locs; ++j) {
      dist = dist_calc(latitudes[j], longitudes[j], depths[j], latitudes[i],
                       longitudes[i], depths[i]);
      if (j != i && dist < distance_cutoff) {
        mask[i] = 1;
        mask[j] = 1;
        break;
      }
    }
  }
  return out;
}


x86 assembly (clang, -O0):
<remove_unclustered>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movss  %xmm0,-0x2c(%rbp)
mov    %r9d,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    1244 <remove_unclustered+0x134>
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movzbl (%rax,%rcx,1),%eax
cmp    $0x0,%eax
je     116c <remove_unclustered+0x5c>
jmp    1233 <remove_unclustered+0x123>
movq   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    122e <remove_unclustered+0x11e>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm3
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm5
call   1030 <dist_calc@plt>
movss  %xmm0,-0x4c(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x40(%rbp),%rax
je     1218 <remove_unclustered+0x108>
movss  -0x4c(%rbp),%xmm1
movss  -0x2c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1218 <remove_unclustered+0x108>
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movb   $0x1,(%rax,%rcx,1)
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
movb   $0x1,(%rax,%rcx,1)
jmp    122e <remove_unclustered+0x11e>
jmp    121d <remove_unclustered+0x10d>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1174 <remove_unclustered+0x64>
jmp    1233 <remove_unclustered+0x123>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1144 <remove_unclustered+0x34>
mov    -0x34(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<remove_unclustered>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
movss  %xmm0,0xc(%rsp)
mov    %r8,0x10(%rsp)
test   %rcx,%rcx
jle    11c3 <remove_unclustered+0xb3>
mov    %rcx,%r15
mov    %rdx,%r12
mov    %rsi,%r13
mov    %rdi,%rbx
xor    %ebp,%ebp
jmp    1159 <remove_unclustered+0x49>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rbp
cmp    %r15,%rbp
je     11c3 <remove_unclustered+0xb3>
mov    0x10(%rsp),%rax
cmpb   $0x0,(%rax,%rbp,1)
jne    1150 <remove_unclustered+0x40>
xor    %r14d,%r14d
jmp    1179 <remove_unclustered+0x69>
nopl   0x0(%rax)
add    $0x1,%r14
cmp    %r14,%r15
je     1150 <remove_unclustered+0x40>
movss  (%rbx,%r14,4),%xmm0
movss  0x0(%r13,%r14,4),%xmm1
movss  (%r12,%r14,4),%xmm2
movss  (%rbx,%rbp,4),%xmm3
movss  0x0(%r13,%rbp,4),%xmm4
movss  (%r12,%rbp,4),%xmm5
call   1030 <dist_calc@plt>
cmp    %r14,%rbp
je     1170 <remove_unclustered+0x60>
movss  0xc(%rsp),%xmm1
ucomiss %xmm0,%xmm1
jbe    1170 <remove_unclustered+0x60>
mov    0x10(%rsp),%rax
movb   $0x1,(%rax,%rbp,1)
movb   $0x1,(%rax,%r14,1)
jmp    1150 <remove_unclustered+0x40>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @remove_unclustered(float* nocapture %0, float* nocapture %1, float* nocapture %2, i64 %3, i8* nocapture %4, float %5, i32 %6) {
%8 = icmp sgt i64 %3, 0
br i1 %8, label %9, label %41
9: ; preds = %7, %38
%10 = phi i64 [ %39, %38 ], [ 0, %7 ]
%11 = getelementptr inbounds i8, i8* %4, i64 %10
%12 = load i8, i8* %11
%13 = icmp eq i8 %12, 0
br i1 %13, label %14, label %38
14: ; preds = %9
%15 = getelementptr inbounds float, float* %0, i64 %10
%16 = getelementptr inbounds float, float* %1, i64 %10
%17 = getelementptr inbounds float, float* %2, i64 %10
br label %18
18: ; preds = %14, %35
%19 = phi i64 [ 0, %14 ], [ %36, %35 ]
%20 = getelementptr inbounds float, float* %0, i64 %19
%21 = load float, float* %20
%22 = getelementptr inbounds float, float* %1, i64 %19
%23 = load float, float* %22
%24 = getelementptr inbounds float, float* %2, i64 %19
%25 = load float, float* %24
%26 = load float, float* %15
%27 = load float, float* %16
%28 = load float, float* %17
%29 = tail call float @dist_calc(float %21, float %23, float %25, float %26, float %27, float %28)
%30 = icmp ne i64 %19, %10
%31 = fcmp olt float %29, %5
%32 = select i1 %30, i1 %31, i1 false
br i1 %32, label %33, label %35
33: ; preds = %18
store i8 1, i8* %11
%34 = getelementptr inbounds i8, i8* %4, i64 %19
store i8 1, i8* %34
br label %38
35: ; preds = %18
%36 = add nuw nsw i64 %19, 1
%37 = icmp eq i64 %36, %3
br i1 %37, label %38, label %18
38: ; preds = %35, %9, %33
%39 = add nuw nsw i64 %10, 1
%40 = icmp eq i64 %39, %3
br i1 %40, label %41, label %9
41: ; preds = %38, %7
ret i32 0
}