#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 */
typedef int uint_fast16_t;

/* Variables and functions */
double circularMvDistance(int, double *, double *);
double k2nn_dist_inSet(int p, int idx, int n, double *set) {
  uint_fast16_t i, j;
  double lower = -1;
  double dist = 0;
  double v;
  for (i = 0, j = 0; j < 2; i++) {
    if (i != idx) {
      j++;
      v = circularMvDistance(p, set + idx * p, set + i * p);
      if (v > dist) {
        lower = dist;
        dist = v;
      } else if (v > lower)
        lower = v;
    }
  }
  for (; i < n; i++) {
    if (i != idx) {
      v = circularMvDistance(p, set + idx * p, set + i * p);
      if (v > lower) {
        if ((v < dist))
          dist = v;
      } else {
        dist = lower;
        lower = v;
      }
    }
  }
  return dist;
}


x86 assembly (clang, -O0):
<k2nn_dist_inSet>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
movsd  0xed3(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
cmpl   $0x2,-0x20(%rbp)
jge    11fc <k2nn_dist_inSet+0xec>
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
je     11e9 <k2nn_dist_inSet+0xd9>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x8(%rbp),%eax
imul   -0x4(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%eax
imul   -0x4(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
call   1030 <circularMvDistance@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    11c5 <k2nn_dist_inSet+0xb5>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    11e4 <k2nn_dist_inSet+0xd4>
movsd  -0x38(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    11df <k2nn_dist_inSet+0xcf>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11e4 <k2nn_dist_inSet+0xd4>
jmp    11e9 <k2nn_dist_inSet+0xd9>
jmp    11ee <k2nn_dist_inSet+0xde>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1148 <k2nn_dist_inSet+0x38>
jmp    1201 <k2nn_dist_inSet+0xf1>
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12ad <k2nn_dist_inSet+0x19d>
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
je     129a <k2nn_dist_inSet+0x18a>
mov    -0x4(%rbp),%edi
mov    -0x18(%rbp),%rsi
mov    -0x8(%rbp),%eax
imul   -0x4(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%eax
imul   -0x4(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
call   1030 <circularMvDistance@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    1281 <k2nn_dist_inSet+0x171>
movsd  -0x38(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    127c <k2nn_dist_inSet+0x16c>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    1295 <k2nn_dist_inSet+0x185>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    129a <k2nn_dist_inSet+0x18a>
jmp    129f <k2nn_dist_inSet+0x18f>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1201 <k2nn_dist_inSet+0xf1>
movsd  -0x30(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<k2nn_dist_inSet>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rcx,%rbx
mov    %rdx,0x30(%rsp)
mov    %edi,%r13d
mov    %esi,%eax
imul   %edi,%eax
cltq
lea    (%rcx,%rax,8),%rax
mov    %rax,0x18(%rsp)
movslq %edi,%rax
mov    %esi,0x14(%rsp)
mov    %esi,%r14d
neg    %r14
mov    %rax,0x20(%rsp)
lea    0x0(,%rax,8),%r15
xor    %r12d,%r12d
movsd  0xea2(%rip),%xmm1
xorpd  %xmm2,%xmm2
mov    %rcx,0x28(%rsp)
xor    %ebp,%ebp
jmp    1184 <k2nn_dist_inSet+0x74>
nopl   0x0(%rax,%rax,1)
movapd %xmm2,%xmm1
movapd %xmm0,%xmm2
add    $0xffffffffffffffff,%r12
add    %r15,%rbx
cmp    $0x2,%ebp
jge    11c4 <k2nn_dist_inSet+0xb4>
cmp    %r12,%r14
je     1178 <k2nn_dist_inSet+0x68>
movsd  %xmm1,(%rsp)
add    $0x1,%ebp
mov    %r13d,%edi
mov    0x18(%rsp),%rsi
mov    %rbx,%rdx
movsd  %xmm2,0x8(%rsp)
call   1030 <circularMvDistance@plt>
movsd  0x8(%rsp),%xmm2
ucomisd %xmm2,%xmm0
ja     1170 <k2nn_dist_inSet+0x60>
movsd  (%rsp),%xmm1
ucomisd %xmm1,%xmm0
jbe    1178 <k2nn_dist_inSet+0x68>
movapd %xmm0,%xmm1
jmp    1178 <k2nn_dist_inSet+0x68>
mov    %r12d,%eax
neg    %eax
mov    0x30(%rsp),%rdx
cmp    %edx,%eax
jge    127a <k2nn_dist_inSet+0x16a>
movslq 0x14(%rsp),%rax
mov    %r12d,%ecx
neg    %ecx
movslq %ecx,%r14
mov    0x20(%rsp),%rcx
imul   %r14,%rcx
sub    %rax,%r14
mov    0x28(%rsp),%rax
lea    (%rax,%rcx,8),%rbx
neg    %edx
sub    %r12,%rdx
xor    %ebp,%ebp
mov    %rdx,%r12
jmp    1224 <k2nn_dist_inSet+0x114>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm2,%xmm0
add    $0xffffffffffffffff,%rbp
add    %r15,%rbx
movapd %xmm0,%xmm2
cmp    %ebp,%r12d
je     127e <k2nn_dist_inSet+0x16e>
cmp    %rbp,%r14
je     1210 <k2nn_dist_inSet+0x100>
movsd  %xmm2,0x8(%rsp)
mov    %r13d,%edi
mov    0x18(%rsp),%rsi
mov    %rbx,%rdx
movsd  %xmm1,(%rsp)
call   1030 <circularMvDistance@plt>
movsd  (%rsp),%xmm3
ucomisd %xmm3,%xmm0
jbe    1270 <k2nn_dist_inSet+0x160>
movsd  0x8(%rsp),%xmm2
movapd %xmm3,%xmm1
ucomisd %xmm0,%xmm2
jbe    1210 <k2nn_dist_inSet+0x100>
jmp    1214 <k2nn_dist_inSet+0x104>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd %xmm0,%xmm1
movapd %xmm3,%xmm0
jmp    1214 <k2nn_dist_inSet+0x104>
movapd %xmm2,%xmm0
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @k2nn_dist_inSet(i32 %0, i32 %1, i32 %2, double* %3) {
%5 = mul nsw i32 %1, %0
%6 = sext i32 %5 to i64
%7 = getelementptr inbounds double, double* %3, i64 %6
%8 = sext i32 %0 to i64
%9 = zext i32 %1 to i64
br label %17
10: ; preds = %32
%11 = trunc i64 %36 to i32
%12 = icmp slt i32 %11, %2
br i1 %12, label %13, label %57
13: ; preds = %10
%14 = shl i64 %36, 32
%15 = ashr exact i64 %14, 32
%16 = sext i32 %1 to i64
br label %38
17: ; preds = %4, %32
%18 = phi i64 [ 0, %4 ], [ %36, %32 ]
%19 = phi double [ 0.000000e+00, %4 ], [ %35, %32 ]
%20 = phi double [ -1.000000e+00, %4 ], [ %34, %32 ]
%21 = phi i32 [ 0, %4 ], [ %33, %32 ]
%22 = icmp eq i64 %18, %9
br i1 %22, label %32, label %23
23: ; preds = %17
%24 = add nsw i32 %21, 1
%25 = mul nsw i64 %18, %8
%26 = getelementptr inbounds double, double* %3, i64 %25
%27 = tail call double @circularMvDistance(i32 %0, double* %7, double* %26)
%28 = fcmp ogt double %27, %19
br i1 %28, label %32, label %29
29: ; preds = %23
%30 = fcmp ogt double %27, %20
br i1 %30, label %31, label %32
31: ; preds = %29
br label %32
32: ; preds = %23, %17, %29, %31
%33 = phi i32 [ %24, %31 ], [ %24, %29 ], [ %21, %17 ], [ %24, %23 ]
%34 = phi double [ %27, %31 ], [ %20, %29 ], [ %20, %17 ], [ %19, %23 ]
%35 = phi double [ %19, %31 ], [ %19, %29 ], [ %19, %17 ], [ %27, %23 ]
%36 = add nuw i64 %18, 1
%37 = icmp slt i32 %33, 2
br i1 %37, label %17, label %10
38: ; preds = %13, %51
%39 = phi i64 [ %15, %13 ], [ %54, %51 ]
%40 = phi double [ %35, %13 ], [ %53, %51 ]
%41 = phi double [ %34, %13 ], [ %52, %51 ]
%42 = icmp eq i64 %39, %16
br i1 %42, label %51, label %43
43: ; preds = %38
%44 = mul nsw i64 %39, %8
%45 = getelementptr inbounds double, double* %3, i64 %44
%46 = tail call double @circularMvDistance(i32 %0, double* %7, double* %45)
%47 = fcmp ogt double %46, %41
br i1 %47, label %48, label %51
48: ; preds = %43
%49 = fcmp olt double %46, %40
br i1 %49, label %50, label %51
50: ; preds = %48
br label %51
51: ; preds = %43, %38, %48, %50
%52 = phi double [ %41, %50 ], [ %41, %48 ], [ %41, %38 ], [ %46, %43 ]
%53 = phi double [ %46, %50 ], [ %40, %48 ], [ %40, %38 ], [ %41, %43 ]
%54 = add nsw i64 %39, 1
%55 = trunc i64 %54 to i32
%56 = icmp eq i32 %55, %2
br i1 %56, label %57, label %38
57: ; preds = %51, %10
%58 = phi double [ %35, %10 ], [ %53, %51 ]
ret double %58
}