void hpsort_doubledouble(unsigned long n, double ra[], double ra2[]) {
  unsigned long i, ir, j, l;
  double rra;
  double rra2;
  if (n < 1)
    return;
  l = (n >> 1) + 1;
  ir = n;
  for (;;) {
    if (l > 1) {
      rra = ra[--l - 1];
      rra2 = ra2[l - 1];
    } else {
      rra = ra[ir - 1];
      rra2 = ra2[ir - 1];
      ra[ir - 1] = ra[0];
      ra2[ir - 1] = ra2[0];
      if (--ir == 1) {
        ra[0] = rra;
        ra2[0] = rra2;
        break;
      }
    }
    i = l;
    j = l + l;
    while (j <= ir) {
      if (j < ir && ra[j - 1] < ra[j])
        j++;
      if (rra < ra[j - 1]) {
        ra[i - 1] = ra[j - 1];
        ra2[i - 1] = ra2[j - 1];
        i = j;
        j <<= 1;
      } else
        j = ir + 1;
    }
    ra[i - 1] = rra;
    ra2[i - 1] = rra2;
  }
}


x86 assembly (clang, -O0):
<hpsort_doubledouble>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
cmpq   $0x1,-0x8(%rbp)
jae    1120 <hpsort_doubledouble+0x20>
jmp    1337 <hpsort_doubledouble+0x237>
mov    -0x8(%rbp),%rax
shr    $0x1,%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
cmpq   $0x1,-0x38(%rbp)
jbe    117c <hpsort_doubledouble+0x7c>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
add    $0xffffffffffffffff,%rcx
mov    %rcx,-0x38(%rbp)
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    1214 <hpsort_doubledouble+0x114>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
sub    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
sub    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x28(%rbp)
cmp    $0x1,%rax
jne    120f <hpsort_doubledouble+0x10f>
movsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1337 <hpsort_doubledouble+0x237>
jmp    1214 <hpsort_doubledouble+0x114>
mov    -0x38(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x38(%rbp),%rax
add    -0x38(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     1306 <hpsort_doubledouble+0x206>
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    1278 <hpsort_doubledouble+0x178>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1278 <hpsort_doubledouble+0x178>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    12f5 <hpsort_doubledouble+0x1f5>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x30(%rbp),%rax
shl    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1301 <hpsort_doubledouble+0x201>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1228 <hpsort_doubledouble+0x128>
movsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1138 <hpsort_doubledouble+0x38>
pop    %rbp
ret

x86 assembly (clang, -O2):
<hpsort_doubledouble>:
test   %rdi,%rdi
je     11f6 <hpsort_doubledouble+0xf6>
mov    %rdi,%r8
shr    %r8
add    $0x1,%r8
jmp    112e <hpsort_doubledouble+0x2e>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm0,-0x8(%rsi,%r10,8)
movsd  %xmm1,-0x8(%rdx,%r10,8)
cmp    $0x2,%r8
jb     1150 <hpsort_doubledouble+0x50>
movsd  -0x10(%rsi,%r8,8),%xmm0
movsd  -0x10(%rdx,%r8,8),%xmm1
add    $0xffffffffffffffff,%r8
jmp    1183 <hpsort_doubledouble+0x83>
nopl   0x0(%rax,%rax,1)
lea    -0x1(%rdi),%rax
movsd  -0x8(%rsi,%rdi,8),%xmm0
movsd  -0x8(%rdx,%rdi,8),%xmm1
movsd  (%rsi),%xmm2
movsd  %xmm2,-0x8(%rsi,%rdi,8)
movsd  (%rdx),%xmm2
movsd  %xmm2,-0x8(%rdx,%rdi,8)
mov    $0x1,%r8d
cmp    $0x1,%rax
je     11ee <hpsort_doubledouble+0xee>
mov    %rax,%rdi
lea    (%r8,%r8,1),%rcx
mov    %r8,%r10
cmp    %rdi,%rcx
ja     1120 <hpsort_doubledouble+0x20>
lea    0x1(%rdi),%r9
mov    %r8,%r10
jmp    11c7 <hpsort_doubledouble+0xc7>
nopl   0x0(%rax,%rax,1)
movsd  %xmm2,-0x8(%rsi,%r10,8)
movsd  -0x8(%rdx,%rcx,8),%xmm2
movsd  %xmm2,-0x8(%rdx,%r10,8)
lea    (%rcx,%rcx,1),%rax
mov    %rcx,%r10
mov    %rax,%rcx
cmp    %rdi,%rax
ja     1120 <hpsort_doubledouble+0x20>
cmp    %rdi,%rcx
jae    11dd <hpsort_doubledouble+0xdd>
movsd  (%rsi,%rcx,8),%xmm2
ucomisd -0x8(%rsi,%rcx,8),%xmm2
jbe    11dd <hpsort_doubledouble+0xdd>
add    $0x1,%rcx
movsd  -0x8(%rsi,%rcx,8),%xmm2
ucomisd %xmm0,%xmm2
ja     11a0 <hpsort_doubledouble+0xa0>
mov    %r9,%rax
jmp    11bb <hpsort_doubledouble+0xbb>
movsd  %xmm0,(%rsi)
movsd  %xmm1,(%rdx)
ret

LLVM IR (-O0):
define void @hpsort_doubledouble(i64 %0, double* %1, double* %2) {
%4 = alloca i64
%5 = alloca double*
%6 = alloca double*
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
%11 = alloca double
%12 = alloca double
store i64 %0, i64* %4
store double* %1, double** %5
store double* %2, double** %6
%13 = load i64, i64* %4
%14 = icmp ult i64 %13, 1
br i1 %14, label %15, label %16
15: ; preds = %3
br label %144
16: ; preds = %3
%17 = load i64, i64* %4
%18 = lshr i64 %17, 1
%19 = add i64 %18, 1
store i64 %19, i64* %10
%20 = load i64, i64* %4
store i64 %20, i64* %8
br label %21
21: ; preds = %133, %16
%22 = load i64, i64* %10
%23 = icmp ugt i64 %22, 1
br i1 %23, label %24, label %36
24: ; preds = %21
%25 = load double*, double** %5
%26 = load i64, i64* %10
%27 = add i64 %26, -1
store i64 %27, i64* %10
%28 = sub i64 %27, 1
%29 = getelementptr inbounds double, double* %25, i64 %28
%30 = load double, double* %29
store double %30, double* %11
%31 = load double*, double** %6
%32 = load i64, i64* %10
%33 = sub i64 %32, 1
%34 = getelementptr inbounds double, double* %31, i64 %33
%35 = load double, double* %34
store double %35, double* %12
br label %72
36: ; preds = %21
%37 = load double*, double** %5
%38 = load i64, i64* %8
%39 = sub i64 %38, 1
%40 = getelementptr inbounds double, double* %37, i64 %39
%41 = load double, double* %40
store double %41, double* %11
%42 = load double*, double** %6
%43 = load i64, i64* %8
%44 = sub i64 %43, 1
%45 = getelementptr inbounds double, double* %42, i64 %44
%46 = load double, double* %45
store double %46, double* %12
%47 = load double*, double** %5
%48 = getelementptr inbounds double, double* %47, i64 0
%49 = load double, double* %48
%50 = load double*, double** %5
%51 = load i64, i64* %8
%52 = sub i64 %51, 1
%53 = getelementptr inbounds double, double* %50, i64 %52
store double %49, double* %53
%54 = load double*, double** %6
%55 = getelementptr inbounds double, double* %54, i64 0
%56 = load double, double* %55
%57 = load double*, double** %6
%58 = load i64, i64* %8
%59 = sub i64 %58, 1
%60 = getelementptr inbounds double, double* %57, i64 %59
store double %56, double* %60
%61 = load i64, i64* %8
%62 = add i64 %61, -1
store i64 %62, i64* %8
%63 = icmp eq i64 %62, 1
br i1 %63, label %64, label %71
64: ; preds = %36
%65 = load double, double* %11
%66 = load double*, double** %5
%67 = getelementptr inbounds double, double* %66, i64 0
store double %65, double* %67
%68 = load double, double* %12
%69 = load double*, double** %6
%70 = getelementptr inbounds double, double* %69, i64 0
store double %68, double* %70
br label %144
71: ; preds = %36
br label %72
72: ; preds = %71, %24
%73 = load i64, i64* %10
store i64 %73, i64* %7
%74 = load i64, i64* %10
%75 = load i64, i64* %10
%76 = add i64 %74, %75
store i64 %76, i64* %9
br label %77
77: ; preds = %132, %72
%78 = load i64, i64* %9
%79 = load i64, i64* %8
%80 = icmp ule i64 %78, %79
br i1 %80, label %81, label %133
81: ; preds = %77
%82 = load i64, i64* %9
%83 = load i64, i64* %8
%84 = icmp ult i64 %82, %83
br i1 %84, label %85, label %99
85: ; preds = %81
%86 = load double*, double** %5
%87 = load i64, i64* %9
%88 = sub i64 %87, 1
%89 = getelementptr inbounds double, double* %86, i64 %88
%90 = load double, double* %89
%91 = load double*, double** %5
%92 = load i64, i64* %9
%93 = getelementptr inbounds double, double* %91, i64 %92
%94 = load double, double* %93
%95 = fcmp olt double %90, %94
br i1 %95, label %96, label %99
96: ; preds = %85
%97 = load i64, i64* %9
%98 = add i64 %97, 1
store i64 %98, i64* %9
br label %99
99: ; preds = %96, %85, %81
%100 = load double, double* %11
%101 = load double*, double** %5
%102 = load i64, i64* %9
%103 = sub i64 %102, 1
%104 = getelementptr inbounds double, double* %101, i64 %103
%105 = load double, double* %104
%106 = fcmp olt double %100, %105
br i1 %106, label %107, label %129
107: ; preds = %99
%108 = load double*, double** %5
%109 = load i64, i64* %9
%110 = sub i64 %109, 1
%111 = getelementptr inbounds double, double* %108, i64 %110
%112 = load double, double* %111
%113 = load double*, double** %5
%114 = load i64, i64* %7
%115 = sub i64 %114, 1
%116 = getelementptr inbounds double, double* %113, i64 %115
store double %112, double* %116
%117 = load double*, double** %6
%118 = load i64, i64* %9
%119 = sub i64 %118, 1
%120 = getelementptr inbounds double, double* %117, i64 %119
%121 = load double, double* %120
%122 = load double*, double** %6
%123 = load i64, i64* %7
%124 = sub i64 %123, 1
%125 = getelementptr inbounds double, double* %122, i64 %124
store double %121, double* %125
%126 = load i64, i64* %9
store i64 %126, i64* %7
%127 = load i64, i64* %9
%128 = shl i64 %127, 1
store i64 %128, i64* %9
br label %132
129: ; preds = %99
%130 = load i64, i64* %8
%131 = add i64 %130, 1
store i64 %131, i64* %9
br label %132
132: ; preds = %129, %107
br label %77
133: ; preds = %77
%134 = load double, double* %11
%135 = load double*, double** %5
%136 = load i64, i64* %7
%137 = sub i64 %136, 1
%138 = getelementptr inbounds double, double* %135, i64 %137
store double %134, double* %138
%139 = load double, double* %12
%140 = load double*, double** %6
%141 = load i64, i64* %7
%142 = sub i64 %141, 1
%143 = getelementptr inbounds double, double* %140, i64 %142
store double %139, double* %143
br label %21
144: ; preds = %15, %64
ret void
}

LLVM IR (-O2):
define void @hpsort_doubledouble(i64 %0, double* nocapture %1, double* nocapture %2) {
%4 = icmp eq i64 %0, 0
br i1 %4, label %73, label %5
5: ; preds = %3
%6 = lshr i64 %0, 1
%7 = add nuw i64 %6, 1
br label %8
8: ; preds = %68, %5
%9 = phi i64 [ %0, %5 ], [ %30, %68 ]
%10 = phi i64 [ %7, %5 ], [ %31, %68 ]
%11 = icmp ugt i64 %10, 1
br i1 %11, label %12, label %19
12: ; preds = %8
%13 = add i64 %10, -1
%14 = add i64 %10, -2
%15 = getelementptr inbounds double, double* %1, i64 %14
%16 = load double, double* %15
%17 = getelementptr inbounds double, double* %2, i64 %14
%18 = load double, double* %17
br label %29
19: ; preds = %8
%20 = add i64 %9, -1
%21 = getelementptr inbounds double, double* %1, i64 %20
%22 = load double, double* %21
%23 = getelementptr inbounds double, double* %2, i64 %20
%24 = load double, double* %23
%25 = load double, double* %1
store double %25, double* %21
%26 = load double, double* %2
store double %26, double* %23
%27 = icmp eq i64 %20, 1
br i1 %27, label %28, label %29
28: ; preds = %19
store double %22, double* %1
store double %24, double* %2
br label %73
29: ; preds = %19, %12
%30 = phi i64 [ %9, %12 ], [ %20, %19 ]
%31 = phi i64 [ %13, %12 ], [ 1, %19 ]
%32 = phi double [ %16, %12 ], [ %22, %19 ]
%33 = phi double [ %18, %12 ], [ %24, %19 ]
%34 = shl i64 %31, 1
%35 = icmp ugt i64 %34, %30
br i1 %35, label %68, label %36
36: ; preds = %29
%37 = add i64 %30, 1
br label %38
38: ; preds = %36, %64
%39 = phi i64 [ %34, %36 ], [ %66, %64 ]
%40 = phi i64 [ %31, %36 ], [ %65, %64 ]
%41 = icmp ult i64 %39, %30
br i1 %41, label %42, label %51
42: ; preds = %38
%43 = add i64 %39, -1
%44 = getelementptr inbounds double, double* %1, i64 %43
%45 = load double, double* %44
%46 = getelementptr inbounds double, double* %1, i64 %39
%47 = load double, double* %46
%48 = fcmp olt double %45, %47
br i1 %48, label %49, label %51
49: ; preds = %42
%50 = add nuw i64 %39, 1
br label %51
51: ; preds = %49, %42, %38
%52 = phi i64 [ %50, %49 ], [ %39, %42 ], [ %39, %38 ]
%53 = add i64 %52, -1
%54 = getelementptr inbounds double, double* %1, i64 %53
%55 = load double, double* %54
%56 = fcmp olt double %32, %55
br i1 %56, label %57, label %64
57: ; preds = %51
%58 = add i64 %40, -1
%59 = getelementptr inbounds double, double* %1, i64 %58
store double %55, double* %59
%60 = getelementptr inbounds double, double* %2, i64 %53
%61 = load double, double* %60
%62 = getelementptr inbounds double, double* %2, i64 %58
store double %61, double* %62
%63 = shl i64 %52, 1
br label %64
64: ; preds = %51, %57
%65 = phi i64 [ %52, %57 ], [ %40, %51 ]
%66 = phi i64 [ %63, %57 ], [ %37, %51 ]
%67 = icmp ugt i64 %66, %30
br i1 %67, label %68, label %38
68: ; preds = %64, %29
%69 = phi i64 [ %31, %29 ], [ %65, %64 ]
%70 = add i64 %69, -1
%71 = getelementptr inbounds double, double* %1, i64 %70
store double %32, double* %71
%72 = getelementptr inbounds double, double* %2, i64 %70
store double %33, double* %72
br label %8
73: ; preds = %3, %28
ret void
}