void praxis_SortDV(double *d, double **v, unsigned n) {
  unsigned i, j, k;
  double s;
  for (i = 0; i < n - 1; i++) {
    k = i;
    s = d[i];
    for (j = i + 1; j < n; j++) {
      if (d[j] > s) {
        k = j;
        s = d[j];
      }
    }
    if (k > i) {
      d[k] = d[i];
      d[i] = s;
      for (j = 0; j < n; j++) {
        s = v[j][i];
        v[j][i] = v[j][k];
        v[j][k] = s;
      }
    }
  }
}


x86 assembly (clang, -O0):
<praxis_SortDV>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jae    1258 <praxis_SortDV+0x158>
mov    -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    1194 <praxis_SortDV+0x94>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    1181 <praxis_SortDV+0x81>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    1186 <praxis_SortDV+0x86>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1147 <praxis_SortDV+0x47>
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jbe    1245 <praxis_SortDV+0x145>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    1240 <praxis_SortDV+0x140>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    (%rax,%rcx,8),%rax
mov    -0x18(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    (%rax,%rcx,8),%rax
mov    -0x18(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11d0 <praxis_SortDV+0xd0>
jmp    1245 <praxis_SortDV+0x145>
jmp    124a <praxis_SortDV+0x14a>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1116 <praxis_SortDV+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<praxis_SortDV>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edx,%eax
add    $0xffffffff,%eax
mov    %eax,-0x8(%rsp)
je     12c7 <praxis_SortDV+0x1c7>
mov    %edx,%r11d
mov    -0x8(%rsp),%r10d
lea    -0x2(%rdx),%eax
mov    %eax,-0x4(%rsp)
mov    %r11d,%r15d
and    $0xfffffffe,%r15d
lea    0x3(%rdx),%r9d
mov    $0x1,%r12d
xor    %r8d,%r8d
jmp    1151 <praxis_SortDV+0x51>
xchg   %ax,%ax
add    $0x1,%r12
add    $0x3,%r9b
cmp    %r10,%r8
je     12c7 <praxis_SortDV+0x1c7>
mov    %r8,%r14
movzbl %r9b,%r9d
and    $0x3,%r9d
movsd  (%rdi,%r8,8),%xmm0
add    $0x1,%r8
cmp    %r11,%r8
jae    11c0 <praxis_SortDV+0xc0>
mov    -0x4(%rsp),%r13d
sub    %r14d,%r13d
mov    -0x8(%rsp),%eax
sub    %r14d,%eax
test   $0x3,%al
je     11c9 <praxis_SortDV+0xc9>
mov    %r9d,%eax
mov    %r12,%rbx
movapd %xmm0,%xmm2
mov    %r14d,%ebp
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rbx,8),%xmm1
ucomisd %xmm2,%xmm1
cmova  %ebx,%ebp
maxsd  %xmm2,%xmm1
add    $0x1,%rbx
movapd %xmm1,%xmm2
add    $0xffffffff,%eax
jne    1190 <praxis_SortDV+0x90>
cmp    $0x3,%r13d
jae    11e0 <praxis_SortDV+0xe0>
jmp    1234 <praxis_SortDV+0x134>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r14d,%ebp
movapd %xmm0,%xmm1
jmp    1234 <praxis_SortDV+0x134>
mov    %r12,%rbx
movapd %xmm0,%xmm1
mov    %r14d,%ebp
cmp    $0x3,%r13d
jb     1234 <praxis_SortDV+0x134>
nopl   0x0(%rax)
movsd  (%rdi,%rbx,8),%xmm2
ucomisd %xmm1,%xmm2
cmova  %ebx,%ebp
movsd  0x8(%rdi,%rbx,8),%xmm3
maxsd  %xmm1,%xmm2
ucomisd %xmm2,%xmm3
lea    0x1(%rbx),%eax
cmovbe %ebp,%eax
maxsd  %xmm2,%xmm3
movsd  0x10(%rdi,%rbx,8),%xmm2
ucomisd %xmm3,%xmm2
lea    0x2(%rbx),%ecx
cmovbe %eax,%ecx
maxsd  %xmm3,%xmm2
movsd  0x18(%rdi,%rbx,8),%xmm1
ucomisd %xmm2,%xmm1
lea    0x3(%rbx),%ebp
cmovbe %ecx,%ebp
maxsd  %xmm2,%xmm1
add    $0x4,%rbx
cmp    %edx,%ebx
jne    11e0 <praxis_SortDV+0xe0>
mov    %ebp,%ebx
cmp    %rbx,%r14
jae    1140 <praxis_SortDV+0x40>
movsd  %xmm0,(%rdi,%rbx,8)
movsd  %xmm1,(%rdi,%r14,8)
test   %edx,%edx
je     1140 <praxis_SortDV+0x40>
xor    %eax,%eax
cmp    $0x1,%edx
je     129e <praxis_SortDV+0x19e>
nopl   0x0(%rax)
mov    (%rsi,%rax,8),%rcx
movsd  (%rcx,%r14,8),%xmm0
movsd  (%rcx,%rbx,8),%xmm1
movsd  %xmm1,(%rcx,%r14,8)
movsd  %xmm0,(%rcx,%rbx,8)
mov    0x8(%rsi,%rax,8),%rcx
movsd  (%rcx,%r14,8),%xmm0
movsd  (%rcx,%rbx,8),%xmm1
movsd  %xmm1,(%rcx,%r14,8)
movsd  %xmm0,(%rcx,%rbx,8)
add    $0x2,%rax
cmp    %rax,%r15
jne    1260 <praxis_SortDV+0x160>
test   $0x1,%r11b
je     1140 <praxis_SortDV+0x40>
mov    (%rsi,%rax,8),%rax
movsd  (%rax,%r14,8),%xmm0
movsd  (%rax,%rbx,8),%xmm1
movsd  %xmm1,(%rax,%r14,8)
movsd  %xmm0,(%rax,%rbx,8)
jmp    1140 <praxis_SortDV+0x40>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @praxis_SortDV(double* nocapture %0, double** nocapture %1, i32 %2) {
%4 = add i32 %2, -1
%5 = icmp eq i32 %4, 0
br i1 %5, label %122, label %6
6: ; preds = %3
%7 = zext i32 %2 to i64
%8 = zext i32 %4 to i64
%9 = add i32 %2, -2
%10 = and i64 %7, 1
%11 = icmp eq i32 %2, 1
%12 = and i64 %7, 4294967294
%13 = icmp eq i64 %10, 0
br label %14
14: ; preds = %6, %119
%15 = phi i64 [ 0, %6 ], [ %21, %119 ]
%16 = phi i64 [ 1, %6 ], [ %120, %119 ]
%17 = trunc i64 %15 to i32
%18 = sub i32 %9, %17
%19 = getelementptr inbounds double, double* %0, i64 %15
%20 = load double, double* %19
%21 = add nuw nsw i64 %15, 1
%22 = icmp ult i64 %21, %7
%23 = trunc i64 %15 to i32
br i1 %22, label %24, label %84
24: ; preds = %14
%25 = trunc i64 %15 to i32
%26 = sub i32 %4, %25
%27 = and i32 %26, 3
%28 = icmp eq i32 %27, 0
br i1 %28, label %43, label %29
29: ; preds = %24, %29
%30 = phi i64 [ %40, %29 ], [ %16, %24 ]
%31 = phi double [ %39, %29 ], [ %20, %24 ]
%32 = phi i32 [ %38, %29 ], [ %23, %24 ]
%33 = phi i32 [ %41, %29 ], [ 0, %24 ]
%34 = getelementptr inbounds double, double* %0, i64 %30
%35 = load double, double* %34
%36 = fcmp ogt double %35, %31
%37 = trunc i64 %30 to i32
%38 = select i1 %36, i32 %37, i32 %32
%39 = select i1 %36, double %35, double %31
%40 = add nuw nsw i64 %30, 1
%41 = add i32 %33, 1
%42 = icmp eq i32 %41, %27
br i1 %42, label %43, label %29
43: ; preds = %29, %24
%44 = phi i32 [ undef, %24 ], [ %38, %29 ]
%45 = phi double [ undef, %24 ], [ %39, %29 ]
%46 = phi i64 [ %16, %24 ], [ %40, %29 ]
%47 = phi double [ %20, %24 ], [ %39, %29 ]
%48 = phi i32 [ %23, %24 ], [ %38, %29 ]
%49 = icmp ult i32 %18, 3
br i1 %49, label %84, label %50
50: ; preds = %43, %50
%51 = phi i64 [ %81, %50 ], [ %46, %43 ]
%52 = phi double [ %80, %50 ], [ %47, %43 ]
%53 = phi i32 [ %79, %50 ], [ %48, %43 ]
%54 = getelementptr inbounds double, double* %0, i64 %51
%55 = load double, double* %54
%56 = fcmp ogt double %55, %52
%57 = trunc i64 %51 to i32
%58 = select i1 %56, i32 %57, i32 %53
%59 = select i1 %56, double %55, double %52
%60 = add nuw nsw i64 %51, 1
%61 = getelementptr inbounds double, double* %0, i64 %60
%62 = load double, double* %61
%63 = fcmp ogt double %62, %59
%64 = trunc i64 %60 to i32
%65 = select i1 %63, i32 %64, i32 %58
%66 = select i1 %63, double %62, double %59
%67 = add nuw nsw i64 %51, 2
%68 = getelementptr inbounds double, double* %0, i64 %67
%69 = load double, double* %68
%70 = fcmp ogt double %69, %66
%71 = trunc i64 %67 to i32
%72 = select i1 %70, i32 %71, i32 %65
%73 = select i1 %70, double %69, double %66
%74 = add nuw nsw i64 %51, 3
%75 = getelementptr inbounds double, double* %0, i64 %74
%76 = load double, double* %75
%77 = fcmp ogt double %76, %73
%78 = trunc i64 %74 to i32
%79 = select i1 %77, i32 %78, i32 %72
%80 = select i1 %77, double %76, double %73
%81 = add nuw nsw i64 %51, 4
%82 = trunc i64 %81 to i32
%83 = icmp eq i32 %82, %2
br i1 %83, label %84, label %50
84: ; preds = %43, %50, %14
%85 = phi i32 [ %23, %14 ], [ %44, %43 ], [ %79, %50 ]
%86 = phi double [ %20, %14 ], [ %45, %43 ], [ %80, %50 ]
%87 = zext i32 %85 to i64
%88 = icmp ult i64 %15, %87
br i1 %88, label %89, label %119
89: ; preds = %84
%90 = getelementptr inbounds double, double* %0, i64 %87
store double %20, double* %90
store double %86, double* %19
switch i32 %2, label %91 [
i32 0, label %119
i32 1, label %110
]
91: ; preds = %89, %91
%92 = phi i64 [ %107, %91 ], [ 0, %89 ]
%93 = phi i64 [ %108, %91 ], [ 0, %89 ]
%94 = getelementptr inbounds double*, double** %1, i64 %92
%95 = load double*, double** %94
%96 = getelementptr inbounds double, double* %95, i64 %15
%97 = load double, double* %96
%98 = getelementptr inbounds double, double* %95, i64 %87
%99 = load double, double* %98
store double %99, double* %96
store double %97, double* %98
%100 = or i64 %92, 1
%101 = getelementptr inbounds double*, double** %1, i64 %100
%102 = load double*, double** %101
%103 = getelementptr inbounds double, double* %102, i64 %15
%104 = load double, double* %103
%105 = getelementptr inbounds double, double* %102, i64 %87
%106 = load double, double* %105
store double %106, double* %103
store double %104, double* %105
%107 = add nuw nsw i64 %92, 2
%108 = add i64 %93, 2
%109 = icmp eq i64 %108, %12
br i1 %109, label %110, label %91
110: ; preds = %91, %89
%111 = phi i64 [ 0, %89 ], [ %107, %91 ]
br i1 %13, label %119, label %112
112: ; preds = %110
%113 = getelementptr inbounds double*, double** %1, i64 %111
%114 = load double*, double** %113
%115 = getelementptr inbounds double, double* %114, i64 %15
%116 = load double, double* %115
%117 = getelementptr inbounds double, double* %114, i64 %87
%118 = load double, double* %117
store double %118, double* %115
store double %116, double* %117
br label %119
119: ; preds = %112, %110, %89, %84
%120 = add nuw nsw i64 %16, 1
%121 = icmp eq i64 %21, %8
br i1 %121, label %122, label %14
122: ; preds = %119, %3
ret void
}