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


x86 assembly (clang, -O0):
<eigsrt_d>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1261 <eigsrt_d+0x161>
mov    -0x20(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%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
jg     1192 <eigsrt_d+0x92>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jb     117f <eigsrt_d+0x7f>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    1184 <eigsrt_d+0x84>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1143 <eigsrt_d+0x43>
mov    -0x18(%rbp),%eax
cmp    -0x20(%rbp),%eax
je     124e <eigsrt_d+0x14e>
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   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1249 <eigsrt_d+0x149>
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    11d1 <eigsrt_d+0xd1>
jmp    124e <eigsrt_d+0x14e>
jmp    1253 <eigsrt_d+0x153>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1116 <eigsrt_d+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<eigsrt_d>:
push   %r14
push   %rbx
cmp    $0x2,%edx
jl     122a <eigsrt_d+0x12a>
mov    %edx,%r8d
lea    0x1(%r8),%r14
mov    %r8d,%r11d
and    $0xfffffffe,%r11d
neg    %r11
mov    $0x2,%r10d
mov    $0x1,%r9d
jmp    113d <eigsrt_d+0x3d>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r10
cmp    %r8,%r9
je     122a <eigsrt_d+0x12a>
mov    %r9,%rax
movsd  (%rdi,%r9,8),%xmm0
add    $0x1,%r9
mov    %r10,%rbx
movapd %xmm0,%xmm1
mov    %eax,%edx
movapd %xmm0,%xmm2
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rbx,8),%xmm3
ucomisd %xmm1,%xmm3
cmovae %ebx,%edx
cmpnlesd %xmm3,%xmm2
movapd %xmm2,%xmm4
andnpd %xmm3,%xmm4
andpd  %xmm1,%xmm2
orpd   %xmm4,%xmm2
movapd %xmm2,%xmm1
add    $0x1,%rbx
cmp    %rbx,%r14
jne    1160 <eigsrt_d+0x60>
mov    %edx,%ebx
cmp    %rbx,%rax
je     1130 <eigsrt_d+0x30>
movslq %edx,%rdx
movsd  %xmm0,(%rdi,%rdx,8)
movsd  %xmm1,(%rdi,%rax,8)
mov    $0x1,%ebx
cmp    $0x2,%r14d
je     1203 <eigsrt_d+0x103>
mov    $0x1,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%rsi,%rbx,8),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  (%rcx,%rdx,8),%xmm1
movsd  %xmm1,(%rcx,%rax,8)
movsd  %xmm0,(%rcx,%rdx,8)
mov    0x8(%rsi,%rbx,8),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  (%rcx,%rdx,8),%xmm1
movsd  %xmm1,(%rcx,%rax,8)
movsd  %xmm0,(%rcx,%rdx,8)
lea    (%r11,%rbx,1),%rcx
add    $0x2,%rcx
add    $0x2,%rbx
cmp    $0x1,%rcx
jne    11c0 <eigsrt_d+0xc0>
test   $0x1,%r8b
je     1130 <eigsrt_d+0x30>
mov    (%rsi,%rbx,8),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  (%rcx,%rdx,8),%xmm1
movsd  %xmm1,(%rcx,%rax,8)
movsd  %xmm0,(%rcx,%rdx,8)
jmp    1130 <eigsrt_d+0x30>
pop    %rbx
pop    %r14
ret

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

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