void add_knn(int k, const void **nn, double *nndist, const void *ndp,
             double d) {
  int i, j;
  if (d >= nndist[k - 1]) {
    return;
  }
  for (i = 0; i < k && nndist[i] < d; i++)
    ;
  for (j = k - 1; j > i; j--) {
    nn[j] = nn[j - 1];
    nndist[j] = nndist[j - 1];
  }
  nn[i] = ndp;
  nndist[i] = d;
}


x86 assembly (clang, -O0):
<add_knn>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     113a <add_knn+0x3a>
jmp    1212 <add_knn+0x112>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%ecx
xor    %eax,%eax
cmp    -0x4(%rbp),%ecx
mov    %al,-0x31(%rbp)
jge    116e <add_knn+0x6e>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x31(%rbp)
mov    -0x31(%rbp),%al
test   $0x1,%al
jne    117e <add_knn+0x7e>
jmp    1191 <add_knn+0x91>
jmp    1183 <add_knn+0x83>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1141 <add_knn+0x41>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jle    11f0 <add_knn+0xf0>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    119a <add_knn+0x9a>
mov    -0x20(%rbp),%rdx
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<add_knn>:
movslq %edi,%r8
ucomisd -0x8(%rdx,%r8,8),%xmm0
jae    11c5 <add_knn+0xc5>
test   %edi,%edi
jle    1138 <add_knn+0x38>
lea    -0x1(%rdi),%r9d
mov    %edi,%r10d
xor    %eax,%eax
nopl   (%rax)
ucomisd (%rdx,%rax,8),%xmm0
jbe    113f <add_knn+0x3f>
add    $0x1,%rax
cmp    %rax,%r10
jne    1120 <add_knn+0x20>
mov    %r10d,%eax
jmp    11ba <add_knn+0xba>
xor    %eax,%eax
jmp    11ba <add_knn+0xba>
cmp    %eax,%r9d
jle    11ba <add_knn+0xba>
mov    %eax,%r9d
not    %r9d
add    %edi,%r9d
lea    -0x2(%r8),%r10
test   $0x1,%r9b
jne    115d <add_knn+0x5d>
add    $0xffffffffffffffff,%r8
jmp    1178 <add_knn+0x78>
mov    -0x10(%rsi,%r8,8),%rdi
mov    %rdi,-0x8(%rsi,%r8,8)
movsd  -0x10(%rdx,%r8,8),%xmm1
movsd  %xmm1,-0x8(%rdx,%r8,8)
mov    %r10,%r8
cmp    %rax,%r10
je     11ba <add_knn+0xba>
nopl   (%rax)
mov    -0x8(%rsi,%r8,8),%rdi
mov    %rdi,(%rsi,%r8,8)
movsd  -0x8(%rdx,%r8,8),%xmm1
movsd  %xmm1,(%rdx,%r8,8)
mov    -0x10(%rsi,%r8,8),%rdi
mov    %rdi,-0x8(%rsi,%r8,8)
movsd  -0x10(%rdx,%r8,8),%xmm1
movsd  %xmm1,-0x8(%rdx,%r8,8)
lea    -0x2(%r8),%rdi
mov    %rdi,%r8
cmp    %rax,%rdi
jg     1180 <add_knn+0x80>
mov    %eax,%eax
mov    %rcx,(%rsi,%rax,8)
movsd  %xmm0,(%rdx,%rax,8)
ret

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

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