user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<median_select>:
  endbr64
  push   %rbx
  mov    %edi,%r10d
  mov    %rdx,%rcx
  lea    -0x1(%rsi),%r11d
  cmp    $0x1,%r11d
  jbe    L15
  mov    $0x0,%ebx
  mov    $0x1,%edi
  jmp    L19
L108:
  mov    %edi,%ebx
L58:
  lea    0x1(%rbx),%edi
  cmp    %r11d,%edi
  jae    L3
L19:
  lea    (%r11,%rbx,1),%eax
  shr    %eax
  mov    %eax,%eax
  lea    (%rcx,%rax,4),%rdx
  movss  (%rdx),%xmm0
  mov    %edi,%eax
  lea    (%rcx,%rax,4),%rax
  movss  (%rax),%xmm1
  movss  %xmm1,(%rdx)
  movss  %xmm0,(%rax)
  mov    %r11d,%edx
  lea    (%rcx,%rdx,4),%rdx
  movss  (%rdx),%xmm1
  comiss %xmm1,%xmm0
  jbe    L91
  movss  %xmm1,(%rax)
  movss  %xmm0,(%rdx)
L91:
  mov    %ebx,%esi
  lea    (%rcx,%rsi,4),%r9
  movss  (%r9),%xmm0
  movss  (%rdx),%xmm1
  comiss %xmm1,%xmm0
  jbe    L126
  movss  %xmm1,(%r9)
  movss  %xmm0,(%rdx)
L126:
  movss  (%rax),%xmm0
  movss  (%r9),%xmm1
  comiss %xmm1,%xmm0
  jbe    L69
  movss  %xmm1,(%rax)
  movss  %xmm0,(%r9)
L69:
  movss  (%r9),%xmm1
  mov    %r11d,%eax
L57:
  add    $0x1,%edi
  mov    %edi,%edx
  lea    (%rcx,%rdx,4),%r8
  movss  (%r8),%xmm2
  comiss %xmm2,%xmm1
  ja     L57
L104:
  mov    %eax,%esi
  sub    $0x1,%eax
  mov    %eax,%edx
  lea    (%rcx,%rdx,4),%rdx
  movss  (%rdx),%xmm0
  comiss %xmm1,%xmm0
  ja     L104
  cmp    %eax,%edi
  ja     L93
  movss  %xmm0,(%r8)
  movss  %xmm2,(%rdx)
  jmp    L57
L15:
  mov    $0x0,%ebx
  mov    $0x1,%edi
L3:
  cmp    %r11d,%edi
  je     L130
L120:
  mov    %r10d,%r10d
  movss  (%rcx,%r10,4),%xmm0
  pop    %rbx
  ret
L130:
  mov    %r11d,%r11d
  lea    (%rcx,%r11,4),%rax
  movss  (%rax),%xmm1
  mov    %ebx,%ebx
  lea    (%rcx,%rbx,4),%rdx
  movss  (%rdx),%xmm0
  comiss %xmm1,%xmm0
  jbe    L120
  movss  %xmm1,(%rdx)
  movss  %xmm0,(%rax)
  jmp    L120
L93:
  movss  %xmm0,(%r9)
  movss  %xmm1,(%rdx)
  cmp    %r10d,%eax
  jb     L108
  lea    -0x2(%rsi),%r11d
  cmove  %edi,%ebx
  jmp    L58
```

assistant: ```c
float median_select(k, n, arr)
float arr[];
unsigned int k, n;
{
  unsigned int i, ir, j, l, mid;
  float a, temp;
  l = 0;
  ir = n - 1;
  for (;;) {
    if (ir <= l + 1) {
      if (ir == l + 1 && arr[ir] < arr[l]) {
        temp = (arr[l]);
        (arr[l]) = (arr[ir]);
        (arr[ir]) = temp;
      }
      return arr[k];
    } else {
      mid = (l + ir) >> 1;
      temp = (arr[mid]);
      (arr[mid]) = (arr[l + 1]);
      (arr[l + 1]) = temp;
      if (arr[l + 1] > arr[ir]) {
        temp = (arr[l + 1]);
        (arr[l + 1]) = (arr[ir]);
        (arr[ir]) = temp;
      }
      if (arr[l] > arr[ir]) {
        temp = (arr[l]);
        (arr[l]) = (arr[ir]);
        (arr[ir]) = temp;
      }
      if (arr[l + 1] > arr[l]) {
        temp = (arr[l + 1]);
        (arr[l + 1]) = (arr[l]);
        (arr[l]) = temp;
      }
      i = l + 1;
      j = ir;
      a = arr[l];
      for (;;) {
        do
          i++;
        while (arr[i] < a);
        do
          j--;
        while (arr[j] > a);
        if (j < i)
          break;
        temp = (arr[i]);
        (arr[i]) = (arr[j]);
        (arr[j]) = temp;
      }
      arr[l] = arr[j];
      arr[j] = a;
      if (j >= k)
        ir = j - 1;
      if (j <= k)
        l = i;
    }
  }
}
```