user: 下面的汇编指令对应的C源代码是什么：

```asm
<filter_c_3d_median>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %esi,-0x5c(%rbp)
  mov    %edx,-0x60(%rbp)
  mov    %ecx,-0x64(%rbp)
  mov    %r8,-0x70(%rbp)
  mov    %r9d,-0x68(%rbp)
  mov    0x20(%rbp),%eax
  imul   %eax,%eax
  mov    0x20(%rbp),%edx
  imul   %edx,%eax
  mov    %eax,-0x30(%rbp)
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x2,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0x8(%rbp)
  mov    0x20(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x2c(%rbp)
  mov    0x20(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x28(%rbp)
  mov    0x20(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x24(%rbp)
  mov    -0x30(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x20(%rbp)
  mov    -0x60(%rbp),%eax
  imul   -0x64(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x24(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  jmp    L69
L35:
  mov    -0x44(%rbp),%eax
  imul   -0x1c(%rbp),%eax
  mov    %eax,-0x18(%rbp)
  mov    -0x28(%rbp),%eax
  mov    %eax,-0x48(%rbp)
  jmp    L42
L104:
  mov    -0x48(%rbp),%eax
  imul   -0x60(%rbp),%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x2c(%rbp),%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L62
L61:
  movl   $0x0,-0x34(%rbp)
  movl   $0x0,-0x38(%rbp)
  jmp    L26
L110:
  mov    -0x44(%rbp),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  sub    -0x24(%rbp),%eax
  mov    -0x1c(%rbp),%edx
  imul   %edx,%eax
  mov    %eax,-0x10(%rbp)
  movl   $0x0,-0x3c(%rbp)
  jmp    L66
L83:
  mov    -0x48(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  sub    -0x28(%rbp),%eax
  imul   -0x60(%rbp),%eax
  mov    %eax,%edx
  mov    -0x10(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0xc(%rbp)
  movl   $0x0,-0x40(%rbp)
  jmp    L97
L74:
  mov    -0xc(%rbp),%edx
  mov    -0x4c(%rbp),%eax
  add    %eax,%edx
  mov    -0x40(%rbp),%eax
  add    %edx,%eax
  sub    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rax,%rdx
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x8(%rbp),%rax
  add    %rcx,%rax
  movss  (%rdx),%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x34(%rbp)
  addl   $0x1,-0x40(%rbp)
L97:
  mov    -0x40(%rbp),%eax
  cmp    0x20(%rbp),%eax
  jl     L74
  addl   $0x1,-0x3c(%rbp)
L66:
  mov    -0x3c(%rbp),%eax
  cmp    0x20(%rbp),%eax
  jl     L83
  addl   $0x1,-0x38(%rbp)
L26:
  mov    -0x38(%rbp),%eax
  cmp    0x20(%rbp),%eax
  jl     L110
  mov    -0x30(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x8(%rbp),%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <inkernel_quicksort@plt>
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x8(%rbp),%rax
  add    %rax,%rdx
  mov    -0x14(%rbp),%ecx
  mov    -0x4c(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x70(%rbp),%rax
  add    %rcx,%rax
  movss  (%rdx),%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x4c(%rbp)
L62:
  mov    -0x64(%rbp),%eax
  sub    -0x2c(%rbp),%eax
  cmp    %eax,-0x4c(%rbp)
  jl     L61
  addl   $0x1,-0x48(%rbp)
L42:
  mov    -0x60(%rbp),%eax
  sub    -0x28(%rbp),%eax
  cmp    %eax,-0x48(%rbp)
  jl     L104
  addl   $0x1,-0x44(%rbp)
L69:
  mov    -0x5c(%rbp),%eax
  sub    -0x24(%rbp),%eax
  cmp    %eax,-0x44(%rbp)
  jl     L35
  nop
  nop
  leave
  ret
```

assistant: ```c
void filter_c_3d_median(float *im, int nz, int ny, int nx, float *res, int nzr, int nyr, int nxr, int w) {
  int nwin = w * w * w;
  float *win = (float *)malloc(nwin * sizeof(float));
  int edgex = w / 2;
  int edgey = w / 2;
  int edgez = w / 2;
  int mpos = nwin / 2;
  int step = ny * nx;
  int x, y, z, wx, wy, wz, ind, indy, indz, indw;
  int nwa;
  for (z = edgez; z < (nz - edgez); ++z) {
    indz = z * step;
    for (y = edgey; y < (ny - edgey); ++y) {
      ind = indz + y * ny;
      for (x = edgex; x < (nx - edgex); ++x) {
        nwa = 0;
        for (wz = 0; wz < w; ++wz) {
          indw = step * (z + wz - edgez);
          for (wy = 0; wy < w; ++wy) {
            indy = indw + ny * (y + wy - edgey);
            for (wx = 0; wx < w; ++wx) {
              win[nwa] = im[indy + x + wx - edgex];
              ++nwa;
            }
          }
        }
        inkernel_quicksort(win, 0, nwin - 1);
        res[ind + x] = win[mpos];
      }
    }
  }
}
```