double median(int n, double x[3]) {
  double temp;
  int i, j;
  for (i = 0; i < n - 1; i++) {
    for (j = i + 1; j < n; j++) {
      if (x[j] < x[i]) {
        temp = x[i];
        x[i] = x[j];
        x[j] = temp;
      }
    }
  }
  if (n % 2 == 0) {
    return ((x[n / 2] + x[n / 2 - 1]) / 2.0);
  } else {
    return x[n / 2];
  }
}


x86 assembly (clang, -O0):
<median>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11c0 <median+0xc0>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11ad <median+0xad>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    119a <median+0x9a>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    119f <median+0x9f>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    112c <median+0x2c>
jmp    11b2 <median+0xb2>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1112 <median+0x12>
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    122f <median+0x12f>
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x30(%rbp),%rax
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xddf(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1255 <median+0x155>
mov    -0x18(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<median>:
cmp    $0x2,%edi
jl     1165 <median+0x65>
lea    -0x1(%rdi),%r8d
mov    %edi,%r9d
xor    %r10d,%r10d
mov    %r9,%r11
mov    %rsi,%rax
jmp    112d <median+0x2d>
nopw   0x0(%rax,%rax,1)
add    $0x8,%rax
add    $0xffffffffffffffff,%r11
cmp    %r8,%r10
je     1165 <median+0x65>
mov    %r10,%rcx
add    $0x1,%r10
cmp    %r9,%r10
jae    1120 <median+0x20>
mov    $0x1,%edx
jmp    1149 <median+0x49>
add    $0x1,%rdx
cmp    %rdx,%r11
je     1120 <median+0x20>
movsd  (%rax,%rdx,8),%xmm0
movsd  (%rsi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jbe    1140 <median+0x40>
movsd  %xmm0,(%rsi,%rcx,8)
movsd  %xmm1,(%rax,%rdx,8)
jmp    1140 <median+0x40>
mov    %edi,%eax
shr    $0x1f,%eax
add    %edi,%eax
sar    %eax
cltq
movsd  (%rsi,%rax,8),%xmm0
test   $0x1,%dil
jne    1189 <median+0x89>
addsd  -0x8(%rsi,%rax,8),%xmm0
mulsd  0xe77(%rip),%xmm0
ret

LLVM IR (-O0):
define double @median(i32 %0, double* %1) {
%3 = alloca double
%4 = alloca i32
%5 = alloca double*
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
store i32 %0, i32* %4
store double* %1, double** %5
store i32 0, i32* %7
br label %9
9: ; preds = %58, %2
%10 = load i32, i32* %7
%11 = load i32, i32* %4
%12 = sub nsw i32 %11, 1
%13 = icmp slt i32 %10, %12
br i1 %13, label %14, label %61
14: ; preds = %9
%15 = load i32, i32* %7
%16 = add nsw i32 %15, 1
store i32 %16, i32* %8
br label %17
17: ; preds = %54, %14
%18 = load i32, i32* %8
%19 = load i32, i32* %4
%20 = icmp slt i32 %18, %19
br i1 %20, label %21, label %57
21: ; preds = %17
%22 = load double*, double** %5
%23 = load i32, i32* %8
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %22, i64 %24
%26 = load double, double* %25
%27 = load double*, double** %5
%28 = load i32, i32* %7
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
%31 = load double, double* %30
%32 = fcmp olt double %26, %31
br i1 %32, label %33, label %53
33: ; preds = %21
%34 = load double*, double** %5
%35 = load i32, i32* %7
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
store double %38, double* %6
%39 = load double*, double** %5
%40 = load i32, i32* %8
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %39, i64 %41
%43 = load double, double* %42
%44 = load double*, double** %5
%45 = load i32, i32* %7
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %44, i64 %46
store double %43, double* %47
%48 = load double, double* %6
%49 = load double*, double** %5
%50 = load i32, i32* %8
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
store double %48, double* %52
br label %53
53: ; preds = %33, %21
br label %54
54: ; preds = %53
%55 = load i32, i32* %8
%56 = add nsw i32 %55, 1
store i32 %56, i32* %8
br label %17
57: ; preds = %17
br label %58
58: ; preds = %57
%59 = load i32, i32* %7
%60 = add nsw i32 %59, 1
store i32 %60, i32* %7
br label %9
61: ; preds = %9
%62 = load i32, i32* %4
%63 = srem i32 %62, 2
%64 = icmp eq i32 %63, 0
br i1 %64, label %65, label %81
65: ; preds = %61
%66 = load double*, double** %5
%67 = load i32, i32* %4
%68 = sdiv i32 %67, 2
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %66, i64 %69
%71 = load double, double* %70
%72 = load double*, double** %5
%73 = load i32, i32* %4
%74 = sdiv i32 %73, 2
%75 = sub nsw i32 %74, 1
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %72, i64 %76
%78 = load double, double* %77
%79 = fadd double %71, %78
%80 = fdiv double %79, 2.000000e+00
store double %80, double* %3
br label %88
81: ; preds = %61
%82 = load double*, double** %5
%83 = load i32, i32* %4
%84 = sdiv i32 %83, 2
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %82, i64 %85
%87 = load double, double* %86
store double %87, double* %3
br label %88
88: ; preds = %81, %65
%89 = load double, double* %3
ret double %89
}

LLVM IR (-O2):
define double @median(i32 %0, double* nocapture %1) {
%3 = icmp sgt i32 %0, 1
br i1 %3, label %4, label %29
4: ; preds = %2
%5 = add nsw i32 %0, -1
%6 = zext i32 %0 to i64
%7 = zext i32 %5 to i64
%8 = zext i32 %0 to i64
br label %12
9: ; preds = %26, %12
%10 = add nuw nsw i64 %14, 1
%11 = icmp eq i64 %15, %7
br i1 %11, label %29, label %12
12: ; preds = %4, %9
%13 = phi i64 [ 0, %4 ], [ %15, %9 ]
%14 = phi i64 [ 1, %4 ], [ %10, %9 ]
%15 = add nuw nsw i64 %13, 1
%16 = icmp ult i64 %15, %6
br i1 %16, label %17, label %9
17: ; preds = %12
%18 = getelementptr inbounds double, double* %1, i64 %13
br label %19
19: ; preds = %17, %26
%20 = phi i64 [ %14, %17 ], [ %27, %26 ]
%21 = getelementptr inbounds double, double* %1, i64 %20
%22 = load double, double* %21
%23 = load double, double* %18
%24 = fcmp olt double %22, %23
br i1 %24, label %25, label %26
25: ; preds = %19
store double %22, double* %18
store double %23, double* %21
br label %26
26: ; preds = %19, %25
%27 = add nuw nsw i64 %20, 1
%28 = icmp eq i64 %27, %8
br i1 %28, label %9, label %19
29: ; preds = %9, %2
%30 = and i32 %0, 1
%31 = icmp eq i32 %30, 0
%32 = sdiv i32 %0, 2
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %1, i64 %33
%35 = load double, double* %34
br i1 %31, label %36, label %43
36: ; preds = %29
%37 = add nsw i32 %32, -1
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %1, i64 %38
%40 = load double, double* %39
%41 = fadd double %35, %40
%42 = fmul double %41, 5.000000e-01
br label %43
43: ; preds = %29, %36
%44 = phi double [ %42, %36 ], [ %35, %29 ]
ret double %44
}