void mos_heapsort_d(int n, double a[]) {
  int i, j, k;
  double x;
  for (k = n / 2 - 1; k >= 0; k--) {
    i = k;
    x = a[i];
    j = 2 * i + 1;
    while (j <= n - 1) {
      if (j < n - 1 && a[j] < a[j + 1])
        j++;
      if (x >= a[j])
        break;
      a[i] = a[j];
      i = j;
      j = 2 * i + 1;
    }
    a[i] = x;
  }
  while (n > 0) {
    x = a[n - 1];
    a[n - 1] = a[0];
    n--;
    i = 0;
    j = 2 * i + 1;
    while (j <= n - 1) {
      if (j < n - 1 && a[j] < a[j + 1])
        j++;
      if (x >= a[j])
        break;
      a[i] = a[j];
      i = j;
      j = 2 * i + 1;
    }
    a[i] = x;
  }
}


x86 assembly (clang, -O0):
<mos_heapsort_d>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jl     120c <mos_heapsort_d+0x10c>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     11ec <mos_heapsort_d+0xec>
mov    -0x18(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    119e <mos_heapsort_d+0x9e>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    119e <mos_heapsort_d+0x9e>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     11bb <mos_heapsort_d+0xbb>
jmp    11ec <mos_heapsort_d+0xec>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    114a <mos_heapsort_d+0x4a>
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    111c <mos_heapsort_d+0x1c>
jmp    1211 <mos_heapsort_d+0x111>
cmpl   $0x0,-0x4(%rbp)
jle    1321 <mos_heapsort_d+0x221>
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     130a <mos_heapsort_d+0x20a>
mov    -0x18(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12bc <mos_heapsort_d+0x1bc>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    12bc <mos_heapsort_d+0x1bc>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     12d9 <mos_heapsort_d+0x1d9>
jmp    130a <mos_heapsort_d+0x20a>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1268 <mos_heapsort_d+0x168>
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1211 <mos_heapsort_d+0x111>
pop    %rbp
ret

x86 assembly (clang, -O2):
<mos_heapsort_d>:
push   %r14
push   %rbx
cmp    $0x2,%edi
jl     11c8 <mos_heapsort_d+0xc8>
mov    %edi,%r9d
shr    %r9d
lea    -0x1(%rdi),%r10d
jmp    1135 <mos_heapsort_d+0x35>
nopl   0x0(%rax,%rax,1)
mov    %r9d,%edx
movslq %edx,%rax
movsd  %xmm0,(%rsi,%rax,8)
cmp    $0x1,%r8
jle    11c8 <mos_heapsort_d+0xc8>
mov    %r9,%r8
add    $0xffffffffffffffff,%r9
mov    %r9d,%eax
movsd  (%rsi,%rax,8),%xmm0
lea    (%r9,%r9,1),%eax
add    $0x1,%eax
cmp    %edi,%eax
jge    1120 <mos_heapsort_d+0x20>
lea    (%r9,%r9,1),%ecx
mov    %r9d,%r11d
cs nopw 0x0(%rax,%rax,1)
movslq %eax,%rbx
cmp    %r10d,%eax
jge    1190 <mos_heapsort_d+0x90>
add    $0x2,%ecx
movslq %ecx,%r14
movsd  (%rsi,%r14,8),%xmm1
ucomisd (%rsi,%rbx,8),%xmm1
mov    %eax,%edx
jbe    1192 <mos_heapsort_d+0x92>
mov    %r14,%rbx
mov    %ecx,%edx
jmp    1192 <mos_heapsort_d+0x92>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %eax,%edx
movsd  (%rsi,%rbx,8),%xmm1
ucomisd %xmm1,%xmm0
jae    11c0 <mos_heapsort_d+0xc0>
movslq %r11d,%rax
movsd  %xmm1,(%rsi,%rax,8)
lea    (%rdx,%rdx,1),%ecx
lea    (%rdx,%rdx,1),%eax
add    $0x1,%eax
mov    %edx,%r11d
cmp    %edi,%eax
jl     1160 <mos_heapsort_d+0x60>
jmp    1123 <mos_heapsort_d+0x23>
nopw   0x0(%rax,%rax,1)
mov    %r11d,%edx
jmp    1123 <mos_heapsort_d+0x23>
test   %edi,%edi
jle    1267 <mos_heapsort_d+0x167>
mov    %edi,%r8d
jmp    11ef <mos_heapsort_d+0xef>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %edi,%eax
cltq
movsd  %xmm0,(%rsi,%rax,8)
cmp    $0x1,%r9
jle    1267 <mos_heapsort_d+0x167>
mov    %r8,%r9
add    $0xffffffffffffffff,%r8
mov    %r8d,%eax
movsd  (%rsi,%rax,8),%xmm0
movsd  (%rsi),%xmm1
movsd  %xmm1,(%rsi,%rax,8)
xor    %eax,%eax
cmp    $0x3,%r9
jb     11e2 <mos_heapsort_d+0xe2>
lea    -0x2(%r9),%r10d
xor    %edi,%edi
mov    $0x1,%edx
nopw   0x0(%rax,%rax,1)
movslq %edx,%rcx
cmp    %r10d,%edx
jge    1240 <mos_heapsort_d+0x140>
movsd  0x8(%rsi,%rcx,8),%xmm1
ucomisd (%rsi,%rcx,8),%xmm1
jbe    1240 <mos_heapsort_d+0x140>
lea    0x1(%rcx),%eax
add    $0x1,%rcx
jmp    1242 <mos_heapsort_d+0x142>
xchg   %ax,%ax
mov    %edx,%eax
movsd  (%rsi,%rcx,8),%xmm1
ucomisd %xmm1,%xmm0
jae    11e0 <mos_heapsort_d+0xe0>
movslq %edi,%rcx
movsd  %xmm1,(%rsi,%rcx,8)
lea    (%rax,%rax,1),%edx
add    $0x1,%edx
mov    %eax,%edi
cmp    %r10d,%edx
jle    1220 <mos_heapsort_d+0x120>
jmp    11e2 <mos_heapsort_d+0xe2>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @mos_heapsort_d(i32 %0, double* %1) {
%3 = alloca i32
%4 = alloca double*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca double
store i32 %0, i32* %3
store double* %1, double** %4
%9 = load i32, i32* %3
%10 = sdiv i32 %9, 2
%11 = sub nsw i32 %10, 1
store i32 %11, i32* %7
br label %12
12: ; preds = %80, %2
%13 = load i32, i32* %7
%14 = icmp sge i32 %13, 0
br i1 %14, label %15, label %83
15: ; preds = %12
%16 = load i32, i32* %7
store i32 %16, i32* %5
%17 = load double*, double** %4
%18 = load i32, i32* %5
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds double, double* %17, i64 %19
%21 = load double, double* %20
store double %21, double* %8
%22 = load i32, i32* %5
%23 = mul nsw i32 2, %22
%24 = add nsw i32 %23, 1
store i32 %24, i32* %6
br label %25
25: ; preds = %60, %15
%26 = load i32, i32* %6
%27 = load i32, i32* %3
%28 = sub nsw i32 %27, 1
%29 = icmp sle i32 %26, %28
br i1 %29, label %30, label %74
30: ; preds = %25
%31 = load i32, i32* %6
%32 = load i32, i32* %3
%33 = sub nsw i32 %32, 1
%34 = icmp slt i32 %31, %33
br i1 %34, label %35, label %51
35: ; preds = %30
%36 = load double*, double** %4
%37 = load i32, i32* %6
%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* %6
%43 = add nsw i32 %42, 1
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %41, i64 %44
%46 = load double, double* %45
%47 = fcmp olt double %40, %46
br i1 %47, label %48, label %51
48: ; preds = %35
%49 = load i32, i32* %6
%50 = add nsw i32 %49, 1
store i32 %50, i32* %6
br label %51
51: ; preds = %48, %35, %30
%52 = load double, double* %8
%53 = load double*, double** %4
%54 = load i32, i32* %6
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
%57 = load double, double* %56
%58 = fcmp oge double %52, %57
br i1 %58, label %59, label %60
59: ; preds = %51
br label %74
60: ; preds = %51
%61 = load double*, double** %4
%62 = load i32, i32* %6
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %61, i64 %63
%65 = load double, double* %64
%66 = load double*, double** %4
%67 = load i32, i32* %5
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %66, i64 %68
store double %65, double* %69
%70 = load i32, i32* %6
store i32 %70, i32* %5
%71 = load i32, i32* %5
%72 = mul nsw i32 2, %71
%73 = add nsw i32 %72, 1
store i32 %73, i32* %6
br label %25
74: ; preds = %59, %25
%75 = load double, double* %8
%76 = load double*, double** %4
%77 = load i32, i32* %5
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %76, i64 %78
store double %75, double* %79
br label %80
80: ; preds = %74
%81 = load i32, i32* %7
%82 = add nsw i32 %81, -1
store i32 %82, i32* %7
br label %12
83: ; preds = %12
br label %84
84: ; preds = %156, %83
%85 = load i32, i32* %3
%86 = icmp sgt i32 %85, 0
br i1 %86, label %87, label %162
87: ; preds = %84
%88 = load double*, double** %4
%89 = load i32, i32* %3
%90 = sub nsw i32 %89, 1
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %88, i64 %91
%93 = load double, double* %92
store double %93, double* %8
%94 = load double*, double** %4
%95 = getelementptr inbounds double, double* %94, i64 0
%96 = load double, double* %95
%97 = load double*, double** %4
%98 = load i32, i32* %3
%99 = sub nsw i32 %98, 1
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %97, i64 %100
store double %96, double* %101
%102 = load i32, i32* %3
%103 = add nsw i32 %102, -1
store i32 %103, i32* %3
store i32 0, i32* %5
%104 = load i32, i32* %5
%105 = mul nsw i32 2, %104
%106 = add nsw i32 %105, 1
store i32 %106, i32* %6
br label %107
107: ; preds = %142, %87
%108 = load i32, i32* %6
%109 = load i32, i32* %3
%110 = sub nsw i32 %109, 1
%111 = icmp sle i32 %108, %110
br i1 %111, label %112, label %156
112: ; preds = %107
%113 = load i32, i32* %6
%114 = load i32, i32* %3
%115 = sub nsw i32 %114, 1
%116 = icmp slt i32 %113, %115
br i1 %116, label %117, label %133
117: ; preds = %112
%118 = load double*, double** %4
%119 = load i32, i32* %6
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %118, i64 %120
%122 = load double, double* %121
%123 = load double*, double** %4
%124 = load i32, i32* %6
%125 = add nsw i32 %124, 1
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds double, double* %123, i64 %126
%128 = load double, double* %127
%129 = fcmp olt double %122, %128
br i1 %129, label %130, label %133
130: ; preds = %117
%131 = load i32, i32* %6
%132 = add nsw i32 %131, 1
store i32 %132, i32* %6
br label %133
133: ; preds = %130, %117, %112
%134 = load double, double* %8
%135 = load double*, double** %4
%136 = load i32, i32* %6
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %135, i64 %137
%139 = load double, double* %138
%140 = fcmp oge double %134, %139
br i1 %140, label %141, label %142
141: ; preds = %133
br label %156
142: ; preds = %133
%143 = load double*, double** %4
%144 = load i32, i32* %6
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds double, double* %143, i64 %145
%147 = load double, double* %146
%148 = load double*, double** %4
%149 = load i32, i32* %5
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %148, i64 %150
store double %147, double* %151
%152 = load i32, i32* %6
store i32 %152, i32* %5
%153 = load i32, i32* %5
%154 = mul nsw i32 2, %153
%155 = add nsw i32 %154, 1
store i32 %155, i32* %6
br label %107
156: ; preds = %141, %107
%157 = load double, double* %8
%158 = load double*, double** %4
%159 = load i32, i32* %5
%160 = sext i32 %159 to i64
%161 = getelementptr inbounds double, double* %158, i64 %160
store double %157, double* %161
br label %84
162: ; preds = %84
ret void
}

LLVM IR (-O2):
define void @mos_heapsort_d(i32 %0, double* nocapture %1) {
%3 = icmp sgt i32 %0, 1
br i1 %3, label %4, label %8
4: ; preds = %2
%5 = lshr i32 %0, 1
%6 = add nsw i32 %0, -1
%7 = zext i32 %5 to i64
br label %12
8: ; preds = %49, %2
%9 = icmp sgt i32 %0, 0
br i1 %9, label %10, label %95
10: ; preds = %8
%11 = zext i32 %0 to i64
br label %54
12: ; preds = %4, %49
%13 = phi i64 [ %7, %4 ], [ %14, %49 ]
%14 = add nsw i64 %13, -1
%15 = trunc i64 %14 to i32
%16 = and i64 %14, 4294967295
%17 = getelementptr inbounds double, double* %1, i64 %16
%18 = load double, double* %17
%19 = shl nsw i32 %15, 1
%20 = or i32 %19, 1
%21 = icmp slt i32 %20, %0
br i1 %21, label %22, label %49
22: ; preds = %12, %43
%23 = phi i32 [ %47, %43 ], [ %20, %12 ]
%24 = phi i32 [ %46, %43 ], [ %19, %12 ]
%25 = phi i32 [ %39, %43 ], [ %15, %12 ]
%26 = icmp slt i32 %23, %6
%27 = sext i32 %23 to i64
br i1 %26, label %28, label %37
28: ; preds = %22
%29 = getelementptr inbounds double, double* %1, i64 %27
%30 = load double, double* %29
%31 = add i32 %24, 2
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %1, i64 %32
%34 = load double, double* %33
%35 = fcmp olt double %30, %34
br i1 %35, label %36, label %37
36: ; preds = %28
br label %37
37: ; preds = %22, %36, %28
%38 = phi i64 [ %32, %36 ], [ %27, %28 ], [ %27, %22 ]
%39 = phi i32 [ %31, %36 ], [ %23, %28 ], [ %23, %22 ]
%40 = getelementptr inbounds double, double* %1, i64 %38
%41 = load double, double* %40
%42 = fcmp ult double %18, %41
br i1 %42, label %43, label %49
43: ; preds = %37
%44 = sext i32 %25 to i64
%45 = getelementptr inbounds double, double* %1, i64 %44
store double %41, double* %45
%46 = shl nsw i32 %39, 1
%47 = or i32 %46, 1
%48 = icmp slt i32 %47, %0
br i1 %48, label %22, label %49
49: ; preds = %43, %37, %12
%50 = phi i32 [ %15, %12 ], [ %25, %37 ], [ %39, %43 ]
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %1, i64 %51
store double %18, double* %52
%53 = icmp sgt i64 %13, 1
br i1 %53, label %12, label %8
54: ; preds = %10, %90
%55 = phi i64 [ %11, %10 ], [ %57, %90 ]
%56 = trunc i64 %55 to i32
%57 = add nsw i64 %55, -1
%58 = and i64 %57, 4294967295
%59 = getelementptr inbounds double, double* %1, i64 %58
%60 = load double, double* %59
%61 = load double, double* %1
store double %61, double* %59
%62 = add nsw i32 %56, -2
%63 = icmp ult i64 %55, 3
br i1 %63, label %90, label %64
64: ; preds = %54, %84
%65 = phi i32 [ %88, %84 ], [ 1, %54 ]
%66 = phi i32 [ %80, %84 ], [ 0, %54 ]
%67 = icmp slt i32 %65, %62
%68 = sext i32 %65 to i64
br i1 %67, label %69, label %78
69: ; preds = %64
%70 = getelementptr inbounds double, double* %1, i64 %68
%71 = load double, double* %70
%72 = add nsw i32 %65, 1
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %1, i64 %73
%75 = load double, double* %74
%76 = fcmp olt double %71, %75
br i1 %76, label %77, label %78
77: ; preds = %69
br label %78
78: ; preds = %64, %77, %69
%79 = phi i64 [ %73, %77 ], [ %68, %69 ], [ %68, %64 ]
%80 = phi i32 [ %72, %77 ], [ %65, %69 ], [ %65, %64 ]
%81 = getelementptr inbounds double, double* %1, i64 %79
%82 = load double, double* %81
%83 = fcmp ult double %60, %82
br i1 %83, label %84, label %90
84: ; preds = %78
%85 = sext i32 %66 to i64
%86 = getelementptr inbounds double, double* %1, i64 %85
store double %82, double* %86
%87 = shl nsw i32 %80, 1
%88 = or i32 %87, 1
%89 = icmp sgt i32 %88, %62
br i1 %89, label %90, label %64
90: ; preds = %84, %78, %54
%91 = phi i32 [ 0, %54 ], [ %66, %78 ], [ %80, %84 ]
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double, double* %1, i64 %92
store double %60, double* %93
%94 = icmp sgt i64 %55, 1
br i1 %94, label %54, label %95
95: ; preds = %90, %8
ret void
}