void eigsort(double d[], double **v, int n) {
  int k, j, i;
  double p;
  for (i = 0; i < n - 1; i++) {
    p = d[k = i];
    for (j = i + 1; j < n; j++)
      if (d[j] >= p)
        p = d[k = j];
    if (k != i) {
      d[k] = d[i];
      d[i] = p;
      for (j = 0; j < n; j++) {
        p = v[j][i];
        v[j][i] = v[j][k];
        v[j][k] = p;
      }
    }
  }
}


x86 assembly (clang, -O0):
<eigsort>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1264 <eigsort+0x164>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    %ecx,-0x18(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1195 <eigsort+0x95>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jb     1182 <eigsort+0x82>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    %ecx,-0x18(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    1187 <eigsort+0x87>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1147 <eigsort+0x47>
mov    -0x18(%rbp),%eax
cmp    -0x20(%rbp),%eax
je     1251 <eigsort+0x151>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    124c <eigsort+0x14c>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11d4 <eigsort+0xd4>
jmp    1251 <eigsort+0x151>
jmp    1256 <eigsort+0x156>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1116 <eigsort+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<eigsort>:
push   %r15
push   %r14
push   %rbx
cmp    $0x2,%edx
jl     1227 <eigsort+0x127>
lea    -0x1(%rdx),%r8d
mov    %edx,%r14d
mov    %r14d,%r11d
and    $0xfffffffe,%r11d
mov    $0x1,%r9d
xor    %r10d,%r10d
jmp    113d <eigsort+0x3d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
cmp    %r8,%r10
je     1227 <eigsort+0x127>
mov    %r10,%r15
movsd  (%rdi,%r10,8),%xmm0
add    $0x1,%r10
cmp    %r14,%r10
jae    1190 <eigsort+0x90>
mov    %r9,%rcx
movapd %xmm0,%xmm2
mov    %r15d,%ebx
movapd %xmm0,%xmm1
nopl   (%rax)
movsd  (%rdi,%rcx,8),%xmm3
ucomisd %xmm2,%xmm3
cmovae %ecx,%ebx
cmpnlesd %xmm3,%xmm1
movapd %xmm1,%xmm4
andnpd %xmm3,%xmm4
andpd  %xmm2,%xmm1
orpd   %xmm4,%xmm1
add    $0x1,%rcx
movapd %xmm1,%xmm2
cmp    %rcx,%r14
jne    1160 <eigsort+0x60>
jmp    1197 <eigsort+0x97>
mov    %r15d,%ebx
movapd %xmm0,%xmm1
mov    %ebx,%ecx
cmp    %rcx,%r15
je     1130 <eigsort+0x30>
movslq %ebx,%rbx
movsd  %xmm0,(%rdi,%rbx,8)
movsd  %xmm1,(%rdi,%r15,8)
xor    %ecx,%ecx
cmp    $0x1,%edx
je     11fe <eigsort+0xfe>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    (%rsi,%rcx,8),%rax
movsd  (%rax,%r15,8),%xmm0
movsd  (%rax,%rbx,8),%xmm1
movsd  %xmm1,(%rax,%r15,8)
movsd  %xmm0,(%rax,%rbx,8)
mov    0x8(%rsi,%rcx,8),%rax
movsd  (%rax,%r15,8),%xmm0
movsd  (%rax,%rbx,8),%xmm1
movsd  %xmm1,(%rax,%r15,8)
movsd  %xmm0,(%rax,%rbx,8)
add    $0x2,%rcx
cmp    %rcx,%r11
jne    11c0 <eigsort+0xc0>
test   $0x1,%r14b
je     1130 <eigsort+0x30>
mov    (%rsi,%rcx,8),%rax
movsd  (%rax,%r15,8),%xmm0
movsd  (%rax,%rbx,8),%xmm1
movsd  %xmm1,(%rax,%r15,8)
movsd  %xmm0,(%rax,%rbx,8)
jmp    1130 <eigsort+0x30>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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