void sortd(double *x, int n) {
  int top, n1;
  double x0;
  if (n < 2)
    return;
  top = n / 2;
  n1 = n - 1;
  while (1) {
    if (top > 0) {
      top--;
      x0 = x[top];
    } else {
      x0 = x[n1];
      x[n1] = x[0];
      n1--;
      if (!n1) {
        x[0] = x0;
        return;
      }
    }
    {
      int parent, child;
      parent = top;
      child = top * 2 + 1;
      while (child <= n1) {
        if (child < n1 && x[child] < x[child + 1])
          child++;
        if (x0 < x[child]) {
          x[parent] = x[child];
          parent = child;
          child += (parent + 1);
        } else
          break;
      }
      x[parent] = x0;
    }
  }
}


x86 assembly (clang, -O0):
<sortd>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
cmpl   $0x2,-0xc(%rbp)
jge    111a <sortd+0x1a>
jmp    1276 <sortd+0x176>
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x10(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0x10(%rbp)
jle    115b <sortd+0x5b>
mov    -0x10(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x10(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11ac <sortd+0xac>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jne    11a7 <sortd+0xa7>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1276 <sortd+0x176>
jmp    11ac <sortd+0xac>
mov    -0x10(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     125f <sortd+0x15f>
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1208 <sortd+0x108>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1208 <sortd+0x108>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1255 <sortd+0x155>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    125a <sortd+0x15a>
jmp    125f <sortd+0x15f>
jmp    11be <sortd+0xbe>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1131 <sortd+0x31>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sortd>:
cmp    $0x2,%esi
jl     11d3 <sortd+0xd3>
mov    %esi,%r8d
shr    %r8d
add    $0xffffffff,%esi
jmp    112b <sortd+0x2b>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r9d,%ecx
movslq %ecx,%rax
movsd  %xmm0,(%rdi,%rax,8)
test   %r8d,%r8d
jle    1140 <sortd+0x40>
add    $0xffffffff,%r8d
movsd  (%rdi,%r8,8),%xmm0
jmp    1159 <sortd+0x59>
nopl   0x0(%rax)
movslq %esi,%rsi
movsd  (%rdi,%rsi,8),%xmm0
movsd  (%rdi),%xmm1
movsd  %xmm1,(%rdi,%rsi,8)
xor    %r8d,%r8d
add    $0xffffffff,%esi
je     11cf <sortd+0xcf>
lea    (%r8,%r8,1),%edx
add    $0x1,%edx
mov    %r8d,%ecx
cmp    %esi,%edx
jg     1123 <sortd+0x23>
lea    (%r8,%r8,1),%eax
mov    %r8d,%r9d
xchg   %ax,%ax
movslq %edx,%r10
cmp    %esi,%edx
jge    11a0 <sortd+0xa0>
add    $0x2,%eax
movslq %eax,%r11
movsd  (%rdi,%r11,8),%xmm1
ucomisd (%rdi,%r10,8),%xmm1
mov    %edx,%ecx
jbe    11a2 <sortd+0xa2>
mov    %r11,%r10
mov    %eax,%ecx
jmp    11a2 <sortd+0xa2>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %edx,%ecx
movsd  (%rdi,%r10,8),%xmm1
ucomisd %xmm0,%xmm1
jbe    1120 <sortd+0x20>
movslq %r9d,%rax
movsd  %xmm1,(%rdi,%rax,8)
lea    (%rcx,%rcx,1),%eax
lea    (%rcx,%rcx,1),%edx
add    $0x1,%edx
mov    %ecx,%r9d
cmp    %esi,%edx
jle    1170 <sortd+0x70>
jmp    1123 <sortd+0x23>
movsd  %xmm0,(%rdi)
ret

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

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