double median(int a[], int a_size, int b[], int b_size) {
  int a_idx = 0;
  int b_idx = 0;
  int current = 0;
  int previous = 0;
  int index = 0;
  while (1) {
    previous = current;
    if (a_idx < a_size && b_idx < b_size) {
      current = a[a_idx] < b[b_idx] ? a[a_idx++] : b[b_idx++];
      index++;
    } else if (a_idx < a_size) {
      current = a[a_idx++];
      index++;
    } else if (b_idx < b_size) {
      current = b[b_idx++];
      index++;
    } else {
      break;
    }
    if (index == (a_size + b_size) / 2 + 1)
      break;
  }
  return (a_size + b_size) % 2 == 0 ? (current + previous) / 2.0 : current;
}


x86 assembly (clang, -O0):
<median>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11b8 <median+0xb8>
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11b8 <median+0xb8>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x18(%rbp),%rcx
movslq -0x24(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jge    118c <median+0x8c>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x20(%rbp)
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x34(%rbp)
jmp    11a4 <median+0xa4>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x24(%rbp)
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    122b <median+0x12b>
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11ea <median+0xea>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x20(%rbp)
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x28(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1226 <median+0x126>
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    121c <median+0x11c>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x24(%rbp)
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x28(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1221 <median+0x121>
jmp    1259 <median+0x159>
jmp    1226 <median+0x126>
jmp    122b <median+0x12b>
mov    -0x30(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
add    -0x1c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x38(%rbp),%eax
add    $0x1,%ecx
cmp    %ecx,%eax
jne    1254 <median+0x154>
jmp    1259 <median+0x159>
jmp    1135 <median+0x35>
mov    -0xc(%rbp),%eax
add    -0x1c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1290 <median+0x190>
mov    -0x28(%rbp),%eax
add    -0x2c(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xd7e(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    129a <median+0x19a>
cvtsi2sdl -0x28(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<median>:
push   %rbp
push   %r15
push   %r14
push   %rbx
lea    (%rcx,%rsi,1),%r8d
mov    %r8d,%eax
shr    $0x1f,%eax
add    %r8d,%eax
sar    %eax
add    $0x1,%eax
xor    %r10d,%r10d
xor    %r11d,%r11d
xor    %r15d,%r15d
jmp    1140 <median+0x40>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movslq %r10d,%rbx
add    $0x1,%r10d
mov    (%rdi,%rbx,4),%r15d
add    $0xffffffff,%eax
je     119a <median+0x9a>
mov    %r15d,%r9d
cmp    %esi,%r10d
jge    1180 <median+0x80>
cmp    %ecx,%r11d
jge    1180 <median+0x80>
movslq %r10d,%rbx
mov    (%rdi,%rbx,4),%r14d
movslq %r11d,%rbx
mov    (%rdx,%rbx,4),%r15d
xor    %ebp,%ebp
xor    %ebx,%ebx
cmp    %r15d,%r14d
setge  %bpl
setl   %bl
cmovl  %r14d,%r15d
add    %ebx,%r10d
add    %ebp,%r11d
jmp    113b <median+0x3b>
cs nopw 0x0(%rax,%rax,1)
nop
cmp    %esi,%r10d
jl     1130 <median+0x30>
cmp    %ecx,%r11d
jge    1197 <median+0x97>
movslq %r11d,%rbx
add    $0x1,%r11d
mov    (%rdx,%rbx,4),%r15d
jmp    113b <median+0x3b>
mov    %r9d,%r15d
test   $0x1,%r8b
je     11a7 <median+0xa7>
cvtsi2sd %r15d,%xmm0
jmp    11b7 <median+0xb7>
add    %r9d,%r15d
cvtsi2sd %r15d,%xmm0
mulsd  0xe49(%rip),%xmm0
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

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