void CCutil_rselect(arr, l, r, m, coord) int *arr;
int l;
int r;
int m;
double *coord;
{
  double samplevals[3];
  int i;
  int st, en;
  int n;
  arr += l;
  n = r - l + 1;
  m -= l;
  while (n > 20) {
    for (i = 0; i < 3; i++) {
      samplevals[i] = coord[arr[CCutil_lprand() % n]];
    }
    select_sort_dsample(samplevals, 3);
    select_split(arr, n, samplevals[(3 - 1) / 2], &st, &en, coord);
    if (st > m) {
      n = st;
    } else if (en <= m) {
      arr += en;
      n -= en;
      m -= en;
    } else {
      return;
    }
  }
  select_sort(arr, n, coord);
  return;
}


x86 assembly (clang, -O0):
<CCutil_rselect>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x14(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x14,-0x50(%rbp)
jle    1278 <CCutil_rselect+0x138>
movl   $0x0,-0x44(%rbp)
cmpl   $0x3,-0x44(%rbp)
jge    11e7 <CCutil_rselect+0xa7>
mov    -0x20(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    $0x0,%al
call   1040 <CCutil_lprand@plt>
mov    -0x60(%rbp),%rcx
cltd
idivl  -0x50(%rbp)
mov    -0x58(%rbp),%rax
movslq %edx,%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0x44(%rbp),%rax
movsd  %xmm0,-0x40(%rbp,%rax,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1196 <CCutil_rselect+0x56>
lea    -0x40(%rbp),%rdi
mov    $0x3,%esi
mov    $0x0,%al
call   1030 <select_sort_dsample@plt>
mov    -0x8(%rbp),%rdi
mov    -0x50(%rbp),%esi
movsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%r8
lea    -0x48(%rbp),%rdx
lea    -0x4c(%rbp),%rcx
mov    $0x1,%al
call   1050 <select_split@plt>
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jle    122d <CCutil_rselect+0xed>
mov    -0x48(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    1273 <CCutil_rselect+0x133>
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1269 <CCutil_rselect+0x129>
mov    -0x4c(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x4c(%rbp),%ecx
mov    -0x50(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%ecx
mov    -0x14(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x14(%rbp)
jmp    126e <CCutil_rselect+0x12e>
jmp    128a <CCutil_rselect+0x14a>
jmp    1273 <CCutil_rselect+0x133>
jmp    1185 <CCutil_rselect+0x45>
mov    -0x8(%rbp),%rdi
mov    -0x50(%rbp),%esi
mov    -0x20(%rbp),%rdx
mov    $0x0,%al
call   1060 <select_sort@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<CCutil_rselect>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8,%r14
mov    %edx,%ebp
movslq %esi,%rax
lea    (%rdi,%rax,4),%rbx
sub    %eax,%ebp
add    $0x1,%ebp
cmp    $0x15,%ebp
jl     123b <CCutil_rselect+0xfb>
mov    %ecx,%r15d
sub    %esi,%r15d
lea    0xc(%rsp),%r13
lea    0x8(%rsp),%r12
jmp    118b <CCutil_rselect+0x4b>
nopw   0x0(%rax,%rax,1)
mov    %esi,%ebp
cmp    $0x14,%esi
jle    123d <CCutil_rselect+0xfd>
xor    %eax,%eax
call   1040 <CCutil_lprand@plt>
cltd
idiv   %ebp
movslq %edx,%rax
movslq (%rbx,%rax,4),%rax
movsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,0x10(%rsp)
xor    %eax,%eax
call   1040 <CCutil_lprand@plt>
cltd
idiv   %ebp
movslq %edx,%rax
movslq (%rbx,%rax,4),%rax
movsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,0x18(%rsp)
xor    %eax,%eax
call   1040 <CCutil_lprand@plt>
cltd
idiv   %ebp
movslq %edx,%rax
movslq (%rbx,%rax,4),%rax
movsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,0x20(%rsp)
lea    0x10(%rsp),%rdi
mov    $0x3,%esi
xor    %eax,%eax
call   1030 <select_sort_dsample@plt>
movsd  0x18(%rsp),%xmm0
mov    %rbx,%rdi
mov    %ebp,%esi
mov    %r13,%rdx
mov    %r12,%rcx
mov    %r14,%r8
mov    $0x1,%al
call   1050 <select_split@plt>
mov    0xc(%rsp),%esi
cmp    %r15d,%esi
jg     1180 <CCutil_rselect+0x40>
movslq 0x8(%rsp),%rax
cmp    %eax,%r15d
jl     124a <CCutil_rselect+0x10a>
lea    (%rbx,%rax,4),%rbx
sub    %eax,%ebp
sub    %eax,%r15d
mov    %ebp,%esi
cmp    $0x14,%esi
jg     118b <CCutil_rselect+0x4b>
jmp    123d <CCutil_rselect+0xfd>
mov    %ebp,%esi
mov    %rbx,%rdi
mov    %r14,%rdx
xor    %eax,%eax
call   1060 <select_sort@plt>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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