void closest_vec(double x[], int dim, int nv, int closest[]) {
  int i, j, k, ib;
  double scor, best, t;
#pragma omp parallel for private(i, j, k, ib, scor, best, t)
  for (i = 0; i < nv; i++) {
    j = (i + 1) % nv;
    best = 0.;
    for (k = 0; k < dim; k++) {
      t = x[i * dim + k] - x[j * dim + k];
      best += t * t;
    }
    ib = j;
    for (j = 0; j < nv; j++) {
      if (i == j)
        continue;
      scor = 0.;
      for (k = 0; k < dim; k++) {
        t = x[i * dim + k] - x[j * dim + k];
        scor += t * t;
      }
      if (scor < best) {
        ib = j;
        best = scor;
      }
    }
    closest[i] = ib;
  }
}


x86 assembly (clang, -O0):
<closest_vec>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    129b <closest_vec+0x19b>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
cltd
idivl  -0x10(%rbp)
mov    %edx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11a8 <closest_vec+0xa8>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1141 <closest_vec+0x41>
mov    -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    127f <closest_vec+0x17f>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jne    11d2 <closest_vec+0xd2>
jmp    1271 <closest_vec+0x171>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1248 <closest_vec+0x148>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11e1 <closest_vec+0xe1>
movsd  -0x30(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    126c <closest_vec+0x16c>
mov    -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    1271 <closest_vec+0x171>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11b5 <closest_vec+0xb5>
mov    -0x28(%rbp),%edx
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <closest_vec+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<closest_vec>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x18(%rsp)
mov    %rdi,-0x30(%rsp)
test   %edx,%edx
jle    1321 <closest_vec+0x221>
movslq %esi,%rdi
mov    %edx,%r15d
mov    %edi,%r10d
lea    -0x1(%r10),%r9
mov    %r10d,%eax
and    $0xfffffffe,%eax
mov    -0x30(%rsp),%rcx
lea    0x8(%rcx),%rdx
mov    %rdi,-0x28(%rsp)
lea    0x0(,%rdi,8),%r12
mov    %rdx,-0x20(%rsp)
xor    %r14d,%r14d
jmp    116d <closest_vec+0x6d>
nop
mov    -0x18(%rsp),%rcx
mov    %r8d,(%rcx,%r14,4)
add    %r12,%rdx
mov    -0x10(%rsp),%rcx
mov    %rcx,%r14
cmp    %r15,%rcx
je     1321 <closest_vec+0x221>
lea    0x1(%r14),%rdi
cmp    %r15,%rdi
mov    %edi,%r8d
mov    $0x0,%ecx
cmove  %ecx,%r8d
test   %esi,%esi
mov    %rdi,-0x10(%rsp)
jle    11f0 <closest_vec+0xf0>
mov    %r8d,%ecx
imul   %esi,%ecx
movslq %ecx,%r11
test   %r9,%r9
je     1200 <closest_vec+0x100>
mov    -0x20(%rsp),%rcx
lea    (%rcx,%r11,8),%rbp
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%rdx,%rcx,8),%xmm2
movsd  (%rdx,%rcx,8),%xmm0
subsd  -0x8(%rbp,%rcx,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
subsd  0x0(%rbp,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
add    $0x2,%rcx
movapd %xmm0,%xmm1
cmp    %rcx,%rax
jne    11b0 <closest_vec+0xb0>
test   $0x1,%r10b
jne    120c <closest_vec+0x10c>
jmp    1235 <closest_vec+0x135>
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
jmp    1235 <closest_vec+0x135>
cs nopw 0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
test   $0x1,%r10b
je     1235 <closest_vec+0x135>
mov    %r14,%rbx
imul   -0x28(%rsp),%rbx
add    %rcx,%rbx
mov    -0x30(%rsp),%rbp
movsd  0x0(%rbp,%rbx,8),%xmm1
add    %rcx,%r11
subsd  0x0(%rbp,%r11,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
mov    %r14,%rcx
imul   -0x28(%rsp),%rcx
mov    %rcx,-0x8(%rsp)
mov    -0x20(%rsp),%r11
xor    %ecx,%ecx
jmp    1260 <closest_vec+0x160>
nopl   0x0(%rax)
add    $0x1,%rcx
add    %r12,%r11
cmp    %r15,%rcx
je     1150 <closest_vec+0x50>
cmp    %rcx,%r14
je     1250 <closest_vec+0x150>
test   %esi,%esi
jle    12bf <closest_vec+0x1bf>
test   %r9,%r9
je     12c5 <closest_vec+0x1c5>
xorpd  %xmm2,%xmm2
xor    %r13d,%r13d
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x8(%rdx,%r13,8),%xmm3
movsd  (%rdx,%r13,8),%xmm1
subsd  -0x8(%r11,%r13,8),%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm2,%xmm3
subsd  (%r11,%r13,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm3,%xmm1
add    $0x2,%r13
movapd %xmm1,%xmm2
cmp    %r13,%rax
jne    1280 <closest_vec+0x180>
test   $0x1,%r10b
jne    12d2 <closest_vec+0x1d2>
jmp    130b <closest_vec+0x20b>
xorpd  %xmm1,%xmm1
jmp    130b <closest_vec+0x20b>
xorpd  %xmm1,%xmm1
xor    %r13d,%r13d
test   $0x1,%r10b
je     130b <closest_vec+0x20b>
mov    %rcx,%rbp
imul   -0x28(%rsp),%rbp
mov    %r10,%rbx
mov    %r9,%r10
mov    -0x8(%rsp),%rdi
lea    (%rdi,%r13,1),%r9
mov    -0x30(%rsp),%rdi
movsd  (%rdi,%r9,8),%xmm2
mov    %r10,%r9
mov    %rbx,%r10
add    %r13,%rbp
subsd  (%rdi,%rbp,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    1250 <closest_vec+0x150>
mov    %ecx,%r8d
movapd %xmm1,%xmm0
jmp    1250 <closest_vec+0x150>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @closest_vec(double* nocapture %0, i32 %1, i32 %2, i32* nocapture writeonly %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %128
6: ; preds = %4
%7 = icmp sgt i32 %1, 0
%8 = sext i32 %1 to i64
%9 = zext i32 %2 to i64
%10 = zext i32 %1 to i64
%11 = add nsw i64 %10, -1
%12 = and i64 %10, 1
%13 = icmp eq i64 %11, 0
%14 = and i64 %10, 4294967294
%15 = icmp eq i64 %12, 0
%16 = and i64 %10, 1
%17 = icmp eq i64 %11, 0
%18 = and i64 %10, 4294967294
%19 = icmp eq i64 %16, 0
br label %20
20: ; preds = %6, %125
%21 = phi i64 [ 0, %6 ], [ %22, %125 ]
%22 = add nuw nsw i64 %21, 1
%23 = icmp eq i64 %22, %9
%24 = trunc i64 %22 to i32
%25 = select i1 %23, i32 0, i32 %24
br i1 %7, label %26, label %43
26: ; preds = %20
%27 = mul nsw i64 %21, %8
%28 = mul nsw i32 %25, %1
%29 = sext i32 %28 to i64
br i1 %13, label %30, label %46
30: ; preds = %46, %26
%31 = phi double [ undef, %26 ], [ %66, %46 ]
%32 = phi i64 [ 0, %26 ], [ %67, %46 ]
%33 = phi double [ 0.000000e+00, %26 ], [ %66, %46 ]
br i1 %15, label %43, label %34
34: ; preds = %30
%35 = add nsw i64 %32, %27
%36 = getelementptr inbounds double, double* %0, i64 %35
%37 = load double, double* %36
%38 = add nsw i64 %32, %29
%39 = getelementptr inbounds double, double* %0, i64 %38
%40 = load double, double* %39
%41 = fsub double %37, %40
%42 = tail call double @llvm.fmuladd.f64(double %41, double %41, double %33)
br label %43
43: ; preds = %34, %30, %20
%44 = phi double [ 0.000000e+00, %20 ], [ %31, %30 ], [ %42, %34 ]
%45 = mul nsw i64 %21, %8
br label %70
46: ; preds = %26, %46
%47 = phi i64 [ %67, %46 ], [ 0, %26 ]
%48 = phi double [ %66, %46 ], [ 0.000000e+00, %26 ]
%49 = phi i64 [ %68, %46 ], [ 0, %26 ]
%50 = add nsw i64 %47, %27
%51 = getelementptr inbounds double, double* %0, i64 %50
%52 = load double, double* %51
%53 = add nsw i64 %47, %29
%54 = getelementptr inbounds double, double* %0, i64 %53
%55 = load double, double* %54
%56 = fsub double %52, %55
%57 = tail call double @llvm.fmuladd.f64(double %56, double %56, double %48)
%58 = or i64 %47, 1
%59 = add nsw i64 %58, %27
%60 = getelementptr inbounds double, double* %0, i64 %59
%61 = load double, double* %60
%62 = add nsw i64 %58, %29
%63 = getelementptr inbounds double, double* %0, i64 %62
%64 = load double, double* %63
%65 = fsub double %61, %64
%66 = tail call double @llvm.fmuladd.f64(double %65, double %65, double %57)
%67 = add nuw nsw i64 %47, 2
%68 = add i64 %49, 2
%69 = icmp eq i64 %68, %14
br i1 %69, label %30, label %46
70: ; preds = %43, %120
%71 = phi i64 [ 0, %43 ], [ %123, %120 ]
%72 = phi double [ %44, %43 ], [ %122, %120 ]
%73 = phi i32 [ %25, %43 ], [ %121, %120 ]
%74 = icmp eq i64 %21, %71
br i1 %74, label %120, label %75
75: ; preds = %70
br i1 %7, label %76, label %115
76: ; preds = %75
%77 = mul nsw i64 %71, %8
br i1 %17, label %102, label %78
78: ; preds = %76, %78
%79 = phi i64 [ %99, %78 ], [ 0, %76 ]
%80 = phi double [ %98, %78 ], [ 0.000000e+00, %76 ]
%81 = phi i64 [ %100, %78 ], [ 0, %76 ]
%82 = add nsw i64 %79, %45
%83 = getelementptr inbounds double, double* %0, i64 %82
%84 = load double, double* %83
%85 = add nsw i64 %79, %77
%86 = getelementptr inbounds double, double* %0, i64 %85
%87 = load double, double* %86
%88 = fsub double %84, %87
%89 = tail call double @llvm.fmuladd.f64(double %88, double %88, double %80)
%90 = or i64 %79, 1
%91 = add nsw i64 %90, %45
%92 = getelementptr inbounds double, double* %0, i64 %91
%93 = load double, double* %92
%94 = add nsw i64 %90, %77
%95 = getelementptr inbounds double, double* %0, i64 %94
%96 = load double, double* %95
%97 = fsub double %93, %96
%98 = tail call double @llvm.fmuladd.f64(double %97, double %97, double %89)
%99 = add nuw nsw i64 %79, 2
%100 = add i64 %81, 2
%101 = icmp eq i64 %100, %18
br i1 %101, label %102, label %78
102: ; preds = %78, %76
%103 = phi double [ undef, %76 ], [ %98, %78 ]
%104 = phi i64 [ 0, %76 ], [ %99, %78 ]
%105 = phi double [ 0.000000e+00, %76 ], [ %98, %78 ]
br i1 %19, label %115, label %106
106: ; preds = %102
%107 = add nsw i64 %104, %45
%108 = getelementptr inbounds double, double* %0, i64 %107
%109 = load double, double* %108
%110 = add nsw i64 %104, %77
%111 = getelementptr inbounds double, double* %0, i64 %110
%112 = load double, double* %111
%113 = fsub double %109, %112
%114 = tail call double @llvm.fmuladd.f64(double %113, double %113, double %105)
br label %115
115: ; preds = %106, %102, %75
%116 = phi double [ 0.000000e+00, %75 ], [ %103, %102 ], [ %114, %106 ]
%117 = fcmp olt double %116, %72
br i1 %117, label %118, label %120
118: ; preds = %115
%119 = trunc i64 %71 to i32
br label %120
120: ; preds = %115, %118, %70
%121 = phi i32 [ %73, %70 ], [ %119, %118 ], [ %73, %115 ]
%122 = phi double [ %72, %70 ], [ %116, %118 ], [ %72, %115 ]
%123 = add nuw nsw i64 %71, 1
%124 = icmp eq i64 %123, %9
br i1 %124, label %125, label %70
125: ; preds = %120
%126 = getelementptr inbounds i32, i32* %3, i64 %21
store i32 %121, i32* %126
%127 = icmp eq i64 %22, %9
br i1 %127, label %128, label %20
128: ; preds = %125, %4
ret void
}