float medianf(float *a, int n) {
  int i, j, l, m;
  float x;
  int k = n & 1 ? n / 2 : n / 2 - 1;
  l = 0;
  m = n - 1;
  while (l < m) {
    x = a[k];
    i = l;
    j = m;
    do {
      while (a[i] < x)
        i++;
      while (x < a[j])
        j--;
      if (i <= j) {
        float t = a[i];
        a[i] = a[j];
        a[j] = t;
        i++;
        j--;
      }
    } while (i <= j);
    if (j < k)
      l = i;
    if (k < i)
      m = j;
  }
  return a[k];
}


x86 assembly (clang, -O0):
<medianf>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     112d <medianf+0x2d>
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x2c(%rbp)
jmp    113e <medianf+0x3e>
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1270 <medianf+0x170>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x10(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    1183 <medianf+0x83>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x20(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    11ac <medianf+0xac>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1183 <medianf+0x83>
jmp    11b1 <medianf+0xb1>
movss  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss %xmm1,%xmm0
jbe    11da <medianf+0xda>
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    11b1 <medianf+0xb1>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1236 <medianf+0x136>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    123b <medianf+0x13b>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jle    117e <medianf+0x7e>
mov    -0x14(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1259 <medianf+0x159>
mov    -0x10(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    126b <medianf+0x16b>
mov    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    1154 <medianf+0x54>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<medianf>:
push   %rbp
push   %r14
push   %rbx
mov    %esi,%eax
and    $0x1,%eax
mov    %esi,%ecx
shr    $0x1f,%ecx
add    %esi,%ecx
sar    %ecx
lea    (%rax,%rcx,1),%r9d
add    $0xffffffff,%r9d
cmp    $0x1,%esi
jle    11ca <medianf+0xca>
add    $0xffffffff,%esi
movslq %r9d,%r8
lea    0x4(%rdi),%r11
xor    %r10d,%r10d
jmp    1153 <medianf+0x53>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cmp    %r9d,%r14d
cmovl  %ebx,%r10d
cmp    %ebx,%r9d
cmovl  %r14d,%esi
cmp    %esi,%r10d
jge    11cd <medianf+0xcd>
movss  (%rdi,%r8,4),%xmm0
mov    %esi,%r14d
mov    %r10d,%ebx
jmp    1175 <medianf+0x75>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
cmp    %r14d,%ebx
jg     1140 <medianf+0x40>
movslq %ebx,%rax
lea    (%rdi,%rax,4),%rdx
mov    %ebx,%ecx
xchg   %ax,%ax
movss  (%rdx),%xmm1
add    $0x4,%rdx
add    $0x1,%ecx
ucomiss %xmm1,%xmm0
ja     1180 <medianf+0x80>
lea    -0x1(%rcx),%ebx
movslq %r14d,%rax
lea    (%r11,%rax,4),%rax
mov    %r14d,%ebp
nopl   (%rax)
movss  -0x4(%rax),%xmm2
add    $0xfffffffffffffffc,%rax
add    $0xffffffff,%ebp
ucomiss %xmm0,%xmm2
ja     11a0 <medianf+0xa0>
lea    0x1(%rbp),%r14d
cmp    %r14d,%ebx
jg     1170 <medianf+0x70>
movss  %xmm2,-0x4(%rdx)
movss  %xmm1,(%rax)
mov    %ebp,%r14d
mov    %ecx,%ebx
jmp    1170 <medianf+0x70>
movslq %r9d,%r8
movss  (%rdi,%r8,4),%xmm0
pop    %rbx
pop    %r14
pop    %rbp
ret

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

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