void sort_index_dvector_selection(double *key, long nl, long nh, int *index) {
  int i, j, pos, tmp;
  double x, min;
  for (i = nl; i <= nh; i++) {
    index[i] = i - nl;
  }
  for (i = nl; i < nh; i++) {
    min = key[i];
    pos = i;
    for (j = i + 1; j <= nh; j++) {
      x = key[j];
      if (x < min) {
        min = x;
        pos = j;
      }
    }
    key[pos] = key[i];
    key[i] = min;
    tmp = index[i];
    index[i] = index[pos];
    index[pos] = tmp;
  }
}


x86 assembly (clang, -O0):
<sort_index_dvector_selection>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    %eax,-0x24(%rbp)
movslq -0x24(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     114c <sort_index_dvector_selection+0x4c>
movslq -0x24(%rbp),%rax
sub    -0x10(%rbp),%rax
mov    %eax,%edx
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111b <sort_index_dvector_selection+0x1b>
mov    -0x10(%rbp),%rax
mov    %eax,-0x24(%rbp)
movslq -0x24(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    1245 <sort_index_dvector_selection+0x145>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
movslq -0x28(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11d9 <sort_index_dvector_selection+0xd9>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11c6 <sort_index_dvector_selection+0xc6>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    11cb <sort_index_dvector_selection+0xcb>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1182 <sort_index_dvector_selection+0x82>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x30(%rbp),%edx
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1153 <sort_index_dvector_selection+0x53>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sort_index_dvector_selection>:
push   %r14
push   %rbx
movslq %esi,%r8
cmp    %rdx,%r8
jg     1230 <sort_index_dvector_selection+0x130>
lea    0x1(%rdx),%rax
lea    0x1(%r8),%r9
cmp    %r9,%rax
cmovg  %rax,%r9
sub    %r8,%r9
mov    %r8,%rax
cmp    $0x8,%r9
jb     1213 <sort_index_dvector_selection+0x113>
mov    %r9,%rax
and    $0xfffffffffffffff8,%rax
movd   %esi,%xmm0
pshufd $0x0,%xmm0,%xmm0
paddd  0xeba(%rip),%xmm0
movd   %esi,%xmm1
pshufd $0x0,%xmm1,%xmm1
lea    -0x8(%rax),%rbx
mov    %rbx,%r11
shr    $0x3,%r11
add    $0x1,%r11
test   %rbx,%rbx
je     12b3 <sort_index_dvector_selection+0x1b3>
mov    %r11,%rbx
and    $0xfffffffffffffffe,%rbx
lea    (%rcx,%r8,4),%r14
add    $0x30,%r14
xor    %r10d,%r10d
movdqa 0xe8f(%rip),%xmm2
movdqa 0xe97(%rip),%xmm3
movdqa 0xe9f(%rip),%xmm4
movdqa 0xea7(%rip),%xmm5
nopl   0x0(%rax)
movdqa %xmm0,%xmm6
psubd  %xmm1,%xmm6
movdqa %xmm6,%xmm7
paddd  %xmm2,%xmm7
movdqu %xmm6,-0x30(%r14,%r10,4)
movdqu %xmm7,-0x20(%r14,%r10,4)
movdqa %xmm6,%xmm7
paddd  %xmm3,%xmm7
paddd  %xmm4,%xmm6
movdqu %xmm7,-0x10(%r14,%r10,4)
movdqu %xmm6,(%r14,%r10,4)
add    $0x10,%r10
paddd  %xmm5,%xmm0
add    $0xfffffffffffffffe,%rbx
jne    11a0 <sort_index_dvector_selection+0xa0>
test   $0x1,%r11b
je     120b <sort_index_dvector_selection+0x10b>
add    %r8,%r10
psubd  %xmm1,%xmm0
movdqa 0xe16(%rip),%xmm1
paddd  %xmm0,%xmm1
movdqu %xmm0,(%rcx,%r10,4)
movdqu %xmm1,0x10(%rcx,%r10,4)
cmp    %rax,%r9
je     1230 <sort_index_dvector_selection+0x130>
add    %r8,%rax
mov    %esi,%r9d
neg    %r9d
nopl   0x0(%rax)
lea    (%r9,%rax,1),%ebx
mov    %ebx,(%rcx,%rax,4)
add    $0x1,%rax
cmp    %rdx,%rax
jle    1220 <sort_index_dvector_selection+0x120>
cmp    %rdx,%r8
jge    12af <sort_index_dvector_selection+0x1af>
add    $0x1,%esi
movslq %esi,%rsi
jmp    1270 <sort_index_dvector_selection+0x170>
nopl   (%rax)
mov    %r8d,%ebx
movapd %xmm0,%xmm1
movslq %ebx,%rax
movsd  %xmm0,(%rdi,%rax,8)
movsd  %xmm1,(%rdi,%r8,8)
mov    (%rcx,%r8,4),%r10d
mov    (%rcx,%rax,4),%ebx
mov    %ebx,(%rcx,%r8,4)
mov    %r10d,(%rcx,%rax,4)
add    $0x1,%rsi
mov    %r9,%r8
cmp    %rdx,%r9
jge    12af <sort_index_dvector_selection+0x1af>
movsd  (%rdi,%r8,8),%xmm0
lea    0x1(%r8),%r9
cmp    %rdx,%r8
jge    1240 <sort_index_dvector_selection+0x140>
mov    %rsi,%rax
movapd %xmm0,%xmm2
mov    %r8d,%ebx
nopl   0x0(%rax)
movsd  (%rdi,%rax,8),%xmm1
ucomisd %xmm1,%xmm2
cmova  %eax,%ebx
minsd  %xmm2,%xmm1
add    $0x1,%rax
movapd %xmm1,%xmm2
cmp    %rdx,%rax
jle    1290 <sort_index_dvector_selection+0x190>
jmp    1247 <sort_index_dvector_selection+0x147>
pop    %rbx
pop    %r14
ret
xor    %r10d,%r10d
test   $0x1,%r11b
jne    11eb <sort_index_dvector_selection+0xeb>
jmp    120b <sort_index_dvector_selection+0x10b>

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

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