void eigsrt(d, v, n) float d[], **v;
int n;
{
  int k, j, i;
  float 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>:
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    1260 <eigsrt+0x160>
mov    -0x20(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1191 <eigsrt+0x91>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x24(%rbp),%xmm0
jb     117e <eigsrt+0x7e>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
jmp    1183 <eigsrt+0x83>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1143 <eigsrt+0x43>
mov    -0x18(%rbp),%eax
cmp    -0x20(%rbp),%eax
je     124d <eigsrt+0x14d>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1248 <eigsrt+0x148>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x24(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11d0 <eigsrt+0xd0>
jmp    124d <eigsrt+0x14d>
jmp    1252 <eigsrt+0x152>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1116 <eigsrt+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<eigsrt>:
push   %r14
push   %rbx
cmp    $0x2,%edx
jl     121a <eigsrt+0x11a>
mov    %edx,%r8d
lea    0x1(%r8),%r14
mov    %r8d,%r11d
and    $0xfffffffe,%r11d
neg    %r11
mov    $0x2,%r10d
mov    $0x1,%r9d
jmp    113d <eigsrt+0x3d>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r10
cmp    %r8,%r9
je     121a <eigsrt+0x11a>
mov    %r9,%rax
movss  (%rdi,%r9,4),%xmm0
add    $0x1,%r9
mov    %r10,%rbx
movaps %xmm0,%xmm1
mov    %eax,%edx
movaps %xmm0,%xmm2
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rdi,%rbx,4),%xmm3
ucomiss %xmm1,%xmm3
cmovae %ebx,%edx
cmpnless %xmm3,%xmm2
movaps %xmm2,%xmm4
andnps %xmm3,%xmm4
andps  %xmm1,%xmm2
orps   %xmm4,%xmm2
movaps %xmm2,%xmm1
add    $0x1,%rbx
cmp    %rbx,%r14
jne    1160 <eigsrt+0x60>
mov    %edx,%ebx
cmp    %rbx,%rax
je     1130 <eigsrt+0x30>
movslq %edx,%rdx
movss  %xmm0,(%rdi,%rdx,4)
movss  %xmm1,(%rdi,%rax,4)
mov    $0x1,%ebx
cmp    $0x2,%r14d
je     11f3 <eigsrt+0xf3>
mov    $0x1,%ebx
nopl   0x0(%rax)
mov    (%rsi,%rbx,8),%rcx
movss  (%rcx,%rax,4),%xmm0
movss  (%rcx,%rdx,4),%xmm1
movss  %xmm1,(%rcx,%rax,4)
movss  %xmm0,(%rcx,%rdx,4)
mov    0x8(%rsi,%rbx,8),%rcx
movss  (%rcx,%rax,4),%xmm0
movss  (%rcx,%rdx,4),%xmm1
movss  %xmm1,(%rcx,%rax,4)
movss  %xmm0,(%rcx,%rdx,4)
lea    (%r11,%rbx,1),%rcx
add    $0x2,%rcx
add    $0x2,%rbx
cmp    $0x1,%rcx
jne    11b0 <eigsrt+0xb0>
test   $0x1,%r8b
je     1130 <eigsrt+0x30>
mov    (%rsi,%rbx,8),%rcx
movss  (%rcx,%rax,4),%xmm0
movss  (%rcx,%rdx,4),%xmm1
movss  %xmm1,(%rcx,%rax,4)
movss  %xmm0,(%rcx,%rdx,4)
jmp    1130 <eigsrt+0x30>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @eigsrt(float* %0, float** %1, i32 %2) {
%4 = alloca float*
%5 = alloca float**
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca float
store float* %0, float** %4
store float** %1, float*** %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 float*, float** %4
%18 = load i32, i32* %9
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds float, float* %17, i64 %19
%21 = load float, float* %20
store float %21, float* %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 float*, float** %4
%30 = load i32, i32* %8
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds float, float* %29, i64 %31
%33 = load float, float* %32
%34 = load float, float* %10
%35 = fcmp oge float %33, %34
br i1 %35, label %36, label %43
36: ; preds = %28
%37 = load i32, i32* %8
store i32 %37, i32* %7
%38 = load float*, float** %4
%39 = load i32, i32* %8
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %38, i64 %40
%42 = load float, float* %41
store float %42, float* %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 float*, float** %4
%53 = load i32, i32* %9
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %52, i64 %54
%56 = load float, float* %55
%57 = load float*, float** %4
%58 = load i32, i32* %7
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %57, i64 %59
store float %56, float* %60
%61 = load float, float* %10
%62 = load float*, float** %4
%63 = load i32, i32* %9
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds float, float* %62, i64 %64
store float %61, float* %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 float**, float*** %5
%72 = load i32, i32* %8
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float*, float** %71, i64 %73
%75 = load float*, float** %74
%76 = load i32, i32* %9
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float, float* %75, i64 %77
%79 = load float, float* %78
store float %79, float* %10
%80 = load float**, float*** %5
%81 = load i32, i32* %8
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds float*, float** %80, i64 %82
%84 = load float*, float** %83
%85 = load i32, i32* %7
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds float, float* %84, i64 %86
%88 = load float, float* %87
%89 = load float**, float*** %5
%90 = load i32, i32* %8
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds float*, float** %89, i64 %91
%93 = load float*, float** %92
%94 = load i32, i32* %9
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds float, float* %93, i64 %95
store float %88, float* %96
%97 = load float, float* %10
%98 = load float**, float*** %5
%99 = load i32, i32* %8
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float*, float** %98, i64 %100
%102 = load float*, float** %101
%103 = load i32, i32* %7
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds float, float* %102, i64 %104
store float %97, float* %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(float* nocapture %0, float** 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 float, float* %0, i64 %15
%18 = load float, float* %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 float [ %18, %14 ], [ %30, %21 ]
%24 = phi i32 [ %20, %14 ], [ %29, %21 ]
%25 = getelementptr inbounds float, float* %0, i64 %22
%26 = load float, float* %25
%27 = fcmp ult float %26, %23
%28 = trunc i64 %22 to i32
%29 = select i1 %27, i32 %24, i32 %28
%30 = select i1 %27, float %23, float %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 float, float* %0, i64 %37
store float %18, float* %38
store float %30, float* %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 float*, float** %1, i64 %40
%43 = load float*, float** %42
%44 = getelementptr inbounds float, float* %43, i64 %15
%45 = load float, float* %44
%46 = getelementptr inbounds float, float* %43, i64 %37
%47 = load float, float* %46
store float %47, float* %44
store float %45, float* %46
%48 = add nuw nsw i64 %40, 1
%49 = getelementptr inbounds float*, float** %1, i64 %48
%50 = load float*, float** %49
%51 = getelementptr inbounds float, float* %50, i64 %15
%52 = load float, float* %51
%53 = getelementptr inbounds float, float* %50, i64 %37
%54 = load float, float* %53
store float %54, float* %51
store float %52, float* %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 float*, float** %1, i64 %59
%62 = load float*, float** %61
%63 = getelementptr inbounds float, float* %62, i64 %15
%64 = load float, float* %63
%65 = getelementptr inbounds float, float* %62, i64 %37
%66 = load float, float* %65
store float %66, float* %63
store float %64, float* %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
}