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

```asm
<mmx_blur16>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x18,%rsp
  mov    %rdi,(%rsp)
  mov    %esi,%r12d
  mov    %edx,%r13d
  mov    %ecx,0xc(%rsp)
  mov    %r8d,%r14d
  test   %edx,%edx
  jle    L98
  movslq %ecx,%rax
  lea    (%rax,%rax,1),%r15
  mov    %rdi,%rbp
  mov    $0x0,%ebx
L50:
  mov    %r14d,%edx
  mov    %r12d,%esi
  mov    %rbp,%rdi
  call   <mmx16_blurrow@plt>
  add    $0x1,%ebx
  add    %r15,%rbp
  cmp    %ebx,%r13d
  jne    L50
L98:
  test   %r12d,%r12d
  je     L86
  testb  $0x2,(%rsp)
  jne    L125
L82:
  cmp    $0x1,%r12d
  jle    L99
  testb  $0x3,(%rsp)
  jne    L152
  cmp    $0x3,%r12d
  jg     L42
L154:
  mov    0xc(%rsp),%eax
  lea    (%rax,%rax,1),%edx
  mov    %r14d,%ecx
  mov    %r13d,%esi
  mov    (%rsp),%rbx
  mov    %rbx,%rdi
  call   <mmx16_blurcol_2@plt>
  mov    %rbx,%rax
  add    $0x4,%rax
  mov    %rax,(%rsp)
  sub    $0x2,%r12d
  jmp    L86
L125:
  mov    (%rsp),%rdi
  lea    0x2(%rdi),%rbx
  mov    0xc(%rsp),%eax
  lea    (%rax,%rax,1),%edx
  mov    %r14d,%ecx
  mov    %r13d,%esi
  call   <mmx16_blurcol@plt>
  sub    $0x1,%r12d
  mov    %rbx,(%rsp)
  jmp    L82
L152:
  mov    0xc(%rsp),%eax
  lea    (%rax,%rax,1),%edx
  mov    %r14d,%ecx
  mov    %r13d,%esi
  mov    (%rsp),%rbx
  mov    %rbx,%rdi
  call   <mmx16_blurcol_2@plt>
  mov    %rbx,%rax
  add    $0x4,%rax
  mov    %rax,(%rsp)
  sub    $0x2,%r12d
  cmp    $0x3,%r12d
  jg     L42
  jmp    L119
L39:
  mov    0xc(%rsp),%edx
  add    %edx,%edx
  mov    %r14d,%ecx
  mov    %r13d,%esi
  mov    (%rsp),%rdi
  call   <mmx16_blurcol@plt>
  jmp    L19
L99:
  cmp    $0x3,%r12d
  jle    L86
L42:
  mov    0xc(%rsp),%eax
  lea    (%rax,%rax,1),%r15d
  mov    %r12d,%ebp
  mov    (%rsp),%rbx
L113:
  mov    %r14d,%ecx
  mov    %r15d,%edx
  mov    %r13d,%esi
  mov    %rbx,%rdi
  call   <mmx16_blurcol_4@plt>
  add    $0x8,%rbx
  sub    $0x4,%ebp
  cmp    $0x3,%ebp
  jg     L113
  lea    -0x4(%r12),%ecx
  mov    %ecx,%eax
  shr    $0x2,%eax
  lea    0x1(%rax),%edx
  mov    %edx,%edx
  shl    $0x3,%rdx
  cmp    $0x3,%r12d
  mov    $0x8,%esi
  cmovle %rsi,%rdx
  add    %rdx,(%rsp)
  neg    %eax
  shl    $0x2,%eax
  cmp    $0x3,%r12d
  mov    $0x0,%edx
  cmovle %edx,%eax
  lea    (%rcx,%rax,1),%r12d
L119:
  cmp    $0x1,%r12d
  jg     L154
L86:
  test   %r12d,%r12d
  jne    L39
L19:
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void mmx_blur16(unsigned short *img, int width, int height, int step_width, int alpha32) {
  int row, col;
  int ii;
  for (row = 0; row < height; ++row)
    mmx16_blurrow(img + row * step_width, width, alpha32);
  if (width && (unsigned int)img & 0x2) {
    mmx16_blurcol(img++, height, step_width * 2, alpha32);
    width--;
  }
  if (width >= 2 && (unsigned int)img & 0x3) {
    mmx16_blurcol_2(img, height, step_width * 2, alpha32);
    img += 2;
    width -= 2;
  }
  while (width >= 4) {
    mmx16_blurcol_4(img, height, step_width * 2, alpha32);
    img += 4;
    width -= 4;
  }
  if (width >= 2) {
    mmx16_blurcol_2(img, height, step_width * 2, alpha32);
    img += 2;
    width -= 2;
  }
  if (width) {
    mmx16_blurcol(img, height, step_width * 2, alpha32);
    img += 1;
    width -= 1;
  }
}
```