void eigsrt(d, v, n) float d[], **v;
int n;
{
  int k, j, i;
  for (i = 0; i < n - 1; i++) {
    float 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):
<eigsrt>:
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    1263 <eigsrt+0x163>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    %ecx,-0x18(%rbp)
movslq %ecx,%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
jge    1194 <eigsrt+0x94>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x24(%rbp),%xmm0
jb     1181 <eigsrt+0x81>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    %ecx,-0x18(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
jmp    1186 <eigsrt+0x86>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1147 <eigsrt+0x47>
mov    -0x18(%rbp),%eax
cmp    -0x20(%rbp),%eax
je     1250 <eigsrt+0x150>
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   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    124b <eigsrt+0x14b>
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    11d3 <eigsrt+0xd3>
jmp    1250 <eigsrt+0x150>
jmp    1255 <eigsrt+0x155>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1116 <eigsrt+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<eigsrt>:
push   %r15
push   %r14
push   %rbx
cmp    $0x2,%edx
jl     1227 <eigsrt+0x127>
lea    -0x1(%rdx),%r8d
mov    %edx,%r14d
mov    %r14d,%r11d
and    $0xfffffffe,%r11d
mov    $0x1,%r9d
xor    %r10d,%r10d
jmp    113d <eigsrt+0x3d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
cmp    %r8,%r10
je     1227 <eigsrt+0x127>
mov    %r10,%r15
movss  (%rdi,%r10,4),%xmm0
add    $0x1,%r10
cmp    %r14,%r10
jae    1190 <eigsrt+0x90>
mov    %r9,%rcx
movaps %xmm0,%xmm2
mov    %r15d,%ebx
movaps %xmm0,%xmm1
nopl   0x0(%rax,%rax,1)
movss  (%rdi,%rcx,4),%xmm3
ucomiss %xmm2,%xmm3
cmovae %ecx,%ebx
cmpnless %xmm3,%xmm1
movaps %xmm1,%xmm4
andnps %xmm3,%xmm4
andps  %xmm2,%xmm1
orps   %xmm4,%xmm1
add    $0x1,%rcx
movaps %xmm1,%xmm2
cmp    %rcx,%r14
jne    1160 <eigsrt+0x60>
jmp    1196 <eigsrt+0x96>
nopw   0x0(%rax,%rax,1)
mov    %r15d,%ebx
movaps %xmm0,%xmm1
mov    %ebx,%ecx
cmp    %rcx,%r15
je     1130 <eigsrt+0x30>
movslq %ebx,%rbx
movss  %xmm0,(%rdi,%rbx,4)
movss  %xmm1,(%rdi,%r15,4)
xor    %ecx,%ecx
cmp    $0x1,%edx
je     11fe <eigsrt+0xfe>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%rsi,%rcx,8),%rax
movss  (%rax,%r15,4),%xmm0
movss  (%rax,%rbx,4),%xmm1
movss  %xmm1,(%rax,%r15,4)
movss  %xmm0,(%rax,%rbx,4)
mov    0x8(%rsi,%rcx,8),%rax
movss  (%rax,%r15,4),%xmm0
movss  (%rax,%rbx,4),%xmm1
movss  %xmm1,(%rax,%r15,4)
movss  %xmm0,(%rax,%rbx,4)
add    $0x2,%rcx
cmp    %rcx,%r11
jne    11c0 <eigsrt+0xc0>
test   $0x1,%r14b
je     1130 <eigsrt+0x30>
mov    (%rsi,%rcx,8),%rax
movss  (%rax,%r15,4),%xmm0
movss  (%rax,%rbx,4),%xmm1
movss  %xmm1,(%rax,%r15,4)
movss  %xmm0,(%rax,%rbx,4)
jmp    1130 <eigsrt+0x30>
pop    %rbx
pop    %r14
pop    %r15
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 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 float*, float** %4
%18 = load i32, i32* %9
store i32 %18, i32* %7
%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 = %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 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 %42
36: ; preds = %28
%37 = load float*, float** %4
%38 = load i32, i32* %8
store i32 %38, i32* %7
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds float, float* %37, i64 %39
%41 = load float, float* %40
store float %41, float* %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 float*, float** %4
%52 = load i32, i32* %9
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds float, float* %51, i64 %53
%55 = load float, float* %54
%56 = load float*, float** %4
%57 = load i32, i32* %7
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %56, i64 %58
store float %55, float* %59
%60 = load float, float* %10
%61 = load float*, float** %4
%62 = load i32, i32* %9
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds float, float* %61, i64 %63
store float %60, float* %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 float**, float*** %5
%71 = load i32, i32* %8
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds float*, float** %70, i64 %72
%74 = load float*, float** %73
%75 = load i32, i32* %9
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds float, float* %74, i64 %76
%78 = load float, float* %77
store float %78, float* %10
%79 = load float**, float*** %5
%80 = load i32, i32* %8
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float*, float** %79, i64 %81
%83 = load float*, float** %82
%84 = load i32, i32* %7
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds float, float* %83, i64 %85
%87 = load float, float* %86
%88 = load float**, float*** %5
%89 = load i32, i32* %8
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds float*, float** %88, i64 %90
%92 = load float*, float** %91
%93 = load i32, i32* %9
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds float, float* %92, i64 %94
store float %87, float* %95
%96 = load float, float* %10
%97 = load float**, float*** %5
%98 = load i32, i32* %8
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds float*, float** %97, i64 %99
%101 = load float*, float** %100
%102 = load i32, i32* %7
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds float, float* %101, i64 %103
store float %96, float* %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 @eigsrt(float* nocapture %0, float** 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 float, float* %0, i64 %15
%18 = load float, float* %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 float [ %31, %22 ], [ %18, %14 ]
%25 = phi i32 [ %30, %22 ], [ %21, %14 ]
%26 = getelementptr inbounds float, float* %0, i64 %23
%27 = load float, float* %26
%28 = fcmp ult float %27, %24
%29 = trunc i64 %23 to i32
%30 = select i1 %28, i32 %25, i32 %29
%31 = select i1 %28, float %24, float %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 float [ %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 float, float* %0, i64 %40
store float %18, float* %41
store float %36, float* %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 float*, float** %1, i64 %43
%46 = load float*, float** %45
%47 = getelementptr inbounds float, float* %46, i64 %15
%48 = load float, float* %47
%49 = getelementptr inbounds float, float* %46, i64 %40
%50 = load float, float* %49
store float %50, float* %47
store float %48, float* %49
%51 = or i64 %43, 1
%52 = getelementptr inbounds float*, float** %1, i64 %51
%53 = load float*, float** %52
%54 = getelementptr inbounds float, float* %53, i64 %15
%55 = load float, float* %54
%56 = getelementptr inbounds float, float* %53, i64 %40
%57 = load float, float* %56
store float %57, float* %54
store float %55, float* %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 float*, float** %1, i64 %62
%65 = load float*, float** %64
%66 = getelementptr inbounds float, float* %65, i64 %15
%67 = load float, float* %66
%68 = getelementptr inbounds float, float* %65, i64 %40
%69 = load float, float* %68
store float %69, float* %66
store float %67, float* %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
}