double qselect(double *v, int len, int k) {
  int i, st, tmp;
  for (st = i = 0; i < len - 1; i++) {
    if (v[i] > v[len - 1])
      continue;
    tmp = v[i];
    v[i] = v[st];
    v[st] = tmp;
    st++;
  }
  tmp = v[len - 1];
  v[len - 1] = v[st];
  v[st] = tmp;
  return k == st  ? v[st]
         : st > k ? qselect(v, st, k)
                  : qselect(v + st, len - st, k - st);
}


x86 assembly (clang, -O0):
<qselect>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11be <qselect+0xae>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    116b <qselect+0x5b>
jmp    11b0 <qselect+0xa0>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1130 <qselect+0x20>
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cvttsd2si (%rax,%rcx,8),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jne    1227 <qselect+0x117>
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    127b <qselect+0x16b>
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
jle    124c <qselect+0x13c>
mov    -0x8(%rbp),%rdi
mov    -0x18(%rbp),%esi
mov    -0x10(%rbp),%edx
call   1030 <qselect@plt>
movsd  %xmm0,-0x30(%rbp)
jmp    1271 <qselect+0x161>
mov    -0x8(%rbp),%rdi
movslq -0x18(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0xc(%rbp),%esi
sub    -0x18(%rbp),%esi
mov    -0x10(%rbp),%edx
sub    -0x18(%rbp),%edx
call   1030 <qselect@plt>
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qselect>:
mov    %esi,%r8d
lea    -0x1(%r8),%eax
cmp    $0x1,%esi
jg     1150 <qselect+0x50>
movslq %eax,%rcx
xor    %esi,%esi
movslq %esi,%rax
movsd  (%rdi,%rax,8),%xmm0
movsd  (%rdi,%rcx,8),%xmm1
movsd  %xmm0,(%rdi,%rcx,8)
cvttpd2dq %xmm1,%xmm0
cvtdq2pd %xmm0,%xmm0
movlpd %xmm0,(%rdi,%rax,8)
mov    %edx,%ecx
sub    %eax,%ecx
je     1194 <qselect+0x94>
jl     1100 <qselect>
lea    (%rdi,%rax,8),%rdi
sub    %eax,%r8d
mov    %r8d,%esi
mov    %ecx,%edx
jmp    1100 <qselect>
cs nopw 0x0(%rax,%rax,1)
mov    %eax,%ecx
xor    %eax,%eax
xor    %esi,%esi
jmp    1169 <qselect+0x69>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
cmp    %rax,%rcx
je     1111 <qselect+0x11>
movsd  (%rdi,%rax,8),%xmm0
ucomisd (%rdi,%rcx,8),%xmm0
ja     1160 <qselect+0x60>
movslq %esi,%rsi
movsd  (%rdi,%rsi,8),%xmm1
movsd  %xmm1,(%rdi,%rax,8)
cvttpd2dq %xmm0,%xmm0
cvtdq2pd %xmm0,%xmm0
movlpd %xmm0,(%rdi,%rsi,8)
add    $0x1,%esi
jmp    1160 <qselect+0x60>
ret

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

LLVM IR (-O2):
define double @qselect(double* %0, i32 %1, i32 %2) {
br label %4
4: ; preds = %49, %3
%5 = phi double* [ %50, %49 ], [ %0, %3 ]
%6 = phi i32 [ %51, %49 ], [ %1, %3 ]
%7 = phi i32 [ %52, %49 ], [ %2, %3 ]
br label %8
8: ; preds = %4, %47
%9 = phi i32 [ %38, %47 ], [ %6, %4 ]
%10 = add i32 %9, -1
%11 = icmp sgt i32 %9, 1
br i1 %11, label %14, label %12
12: ; preds = %8
%13 = sext i32 %10 to i64
br label %36
14: ; preds = %8
%15 = zext i32 %10 to i64
%16 = getelementptr inbounds double, double* %5, i64 %15
%17 = zext i32 %10 to i64
br label %18
18: ; preds = %14, %32
%19 = phi i64 [ 0, %14 ], [ %34, %32 ]
%20 = phi i32 [ 0, %14 ], [ %33, %32 ]
%21 = getelementptr inbounds double, double* %5, i64 %19
%22 = load double, double* %21
%23 = load double, double* %16
%24 = fcmp ogt double %22, %23
br i1 %24, label %32, label %25
25: ; preds = %18
%26 = fptosi double %22 to i32
%27 = sext i32 %20 to i64
%28 = getelementptr inbounds double, double* %5, i64 %27
%29 = load double, double* %28
store double %29, double* %21
%30 = sitofp i32 %26 to double
store double %30, double* %28
%31 = add nsw i32 %20, 1
br label %32
32: ; preds = %18, %25
%33 = phi i32 [ %20, %18 ], [ %31, %25 ]
%34 = add nuw nsw i64 %19, 1
%35 = icmp eq i64 %34, %17
br i1 %35, label %36, label %18
36: ; preds = %32, %12
%37 = phi i64 [ %13, %12 ], [ %15, %32 ]
%38 = phi i32 [ 0, %12 ], [ %33, %32 ]
%39 = getelementptr inbounds double, double* %5, i64 %37
%40 = load double, double* %39
%41 = fptosi double %40 to i32
%42 = sext i32 %38 to i64
%43 = getelementptr inbounds double, double* %5, i64 %42
%44 = load double, double* %43
store double %44, double* %39
%45 = sitofp i32 %41 to double
store double %45, double* %43
%46 = icmp eq i32 %7, %38
br i1 %46, label %53, label %47
47: ; preds = %36
%48 = icmp slt i32 %7, %38
br i1 %48, label %8, label %49
49: ; preds = %47
%50 = getelementptr inbounds double, double* %5, i64 %42
%51 = sub nsw i32 %9, %38
%52 = sub nsw i32 %7, %38
br label %4
53: ; preds = %36
ret double %45
}