void qsd(double *a, int n) {
  int i, j;
  double temp, p;
  p = *(a + (n >> 1));
  i = 0;
  j = n;
  do {
    while (*(a + i) < p)
      i++;
    while (*(a + j) > p)
      j--;
    if (i <= j) {
      temp = *(a + i);
      *(a + i) = *(a + j);
      *(a + j) = temp;
      i++;
      j--;
    }
  } while (i <= j);
  if (j > 0)
    qsd(a, j);
  if (n > i)
    qsd(a + i, n - i);
}


x86 assembly (clang, -O0):
<qsd>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sar    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    1148 <qsd+0x38>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1172 <qsd+0x62>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1148 <qsd+0x38>
jmp    1177 <qsd+0x67>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    119d <qsd+0x8d>
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    1177 <qsd+0x67>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11f9 <qsd+0xe9>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    11fe <qsd+0xee>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jle    1143 <qsd+0x33>
cmpl   $0x0,-0x14(%rbp)
jle    1220 <qsd+0x110>
mov    -0x8(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1030 <qsd@plt>
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jle    1246 <qsd+0x136>
mov    -0x8(%rbp),%rdi
movslq -0x10(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0xc(%rbp),%esi
sub    -0x10(%rbp),%esi
call   1030 <qsd@plt>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qsd>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %esi,%r14d
mov    %rdi,%r15
nopl   (%rax)
mov    %r14d,%eax
sar    %eax
cltq
movsd  (%r15,%rax,8),%xmm0
lea    0x8(%r15),%rax
xor    %ebp,%ebp
mov    %r14d,%esi
jmp    1144 <qsd+0x34>
nopl   0x0(%rax,%rax,1)
cmp    %esi,%ebp
jg     11a0 <qsd+0x90>
movslq %ebp,%rcx
lea    (%r15,%rcx,8),%rdx
mov    %ebp,%ecx
nopl   (%rax)
movsd  (%rdx),%xmm1
add    $0x8,%rdx
add    $0x1,%ecx
ucomisd %xmm1,%xmm0
ja     1150 <qsd+0x40>
lea    -0x1(%rcx),%ebp
movslq %esi,%rdi
lea    (%rax,%rdi,8),%rdi
mov    %esi,%ebx
nopl   (%rax)
movsd  -0x8(%rdi),%xmm2
add    $0xfffffffffffffff8,%rdi
add    $0xffffffff,%ebx
ucomisd %xmm0,%xmm2
ja     1170 <qsd+0x60>
lea    0x1(%rbx),%esi
cmp    %esi,%ebp
jg     1140 <qsd+0x30>
movsd  %xmm2,-0x8(%rdx)
movsd  %xmm1,(%rdi)
mov    %ebx,%esi
mov    %ecx,%ebp
jmp    1140 <qsd+0x30>
nopl   0x0(%rax,%rax,1)
test   %esi,%esi
jle    11ac <qsd+0x9c>
mov    %r15,%rdi
call   1030 <qsd@plt>
sub    %ebp,%r14d
jle    11bd <qsd+0xad>
movslq %ebp,%rax
lea    (%r15,%rax,8),%r15
jmp    1120 <qsd+0x10>
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

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