user: 请从给定的汇编指令推导出C源代码：

```asm
<col221111_16>:
  endbr64
  push   %r15
  movslq %edx,%rax
  add    %edx,%edx
  push   %r14
  add    %rsi,%rax
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  lea    0x100(%rdi),%rbx
  mov    %rax,-0x20(%rsp)
  lea    0x400(%rdi),%rax
  mov    %rax,-0x28(%rsp)
  movslq %edx,%rax
  mov    %rax,-0x30(%rsp)
  shl    $0x2,%rax
  mov    %rbx,-0x38(%rsp)
  lea    0x500(%rdi),%rbx
  mov    %rsi,-0x18(%rsp)
  mov    %rbx,-0x8(%rsp)
  mov    %rax,-0x10(%rsp)
L30:
  mov    -0x38(%rsp),%rax
  mov    -0x28(%rsp),%r15
  mov    -0x20(%rsp),%rbp
  mov    -0x18(%rsp),%rbx
  lea    -0x100(%rax),%r11
L231:
  xor    %eax,%eax
  jmp    L46
  nopw   0x0(%rax,%rax,1)
L102:
  shl    $0x8,%edx
  and    $0xf800,%edx
L130:
  mov    %r9d,%r14d
  sub    %r8d,%r14d
  lea    0x3(%r14),%r13d
  cmp    $0xff,%r13d
  ja     L77
  shl    $0x3,%r13d
  and    $0x7e0,%r13d
  or     %r13d,%edx
L201:
  add    %esi,%r9d
  lea    0x7(%r9),%r13d
  cmp    $0xff,%r13d
  ja     L272
  sar    $0x3,%r13d
  or     %r13d,%edx
L195:
  mov    %dh,(%rbx,%rax,4)
  mov    %dl,0x1(%rbx,%rax,4)
  mov    -0x40(%rsp),%rdx
  mov    0x4(%r11,%rdx,4),%edx
  lea    (%rdi,%rdx,1),%ecx
  lea    0x1(%rcx),%r9d
  cmp    $0xff,%r9d
  ja     L122
  shl    $0x8,%r9d
  and    $0xf800,%r9d
L76:
  mov    %edx,%ecx
  sub    %r8d,%ecx
  cmp    $0xff,%ecx
  ja     L170
  shl    $0x3,%ecx
  and    $0x7e0,%ecx
  or     %ecx,%r9d
L27:
  add    %esi,%edx
  lea    0x1(%rdx),%ecx
  cmp    $0xff,%ecx
  ja     L17
  sar    $0x3,%ecx
  or     %ecx,%r9d
L144:
  mov    %r9d,%edx
  mov    %r9b,0x3(%rbx,%rax,4)
  mov    %dh,0x2(%rbx,%rax,4)
  lea    0x1(,%r12,8),%edx
  lea    (%r10,%rdx,4),%edx
  add    %edx,%edx
  movslq %edx,%r10
  mov    (%r11,%r10,4),%edx
  lea    (%rdi,%rdx,1),%ecx
  lea    0x3(%rcx),%r9d
  cmp    $0xff,%r9d
  ja     L11
  shl    $0x8,%r9d
  and    $0xf800,%r9d
L166:
  mov    %edx,%r12d
  sub    %r8d,%r12d
  lea    0x1(%r12),%ecx
  cmp    $0xff,%ecx
  ja     L62
  shl    $0x3,%ecx
  and    $0x7e0,%ecx
  or     %ecx,%r9d
L120:
  add    %esi,%edx
  lea    0x3(%rdx),%ecx
  cmp    $0xff,%ecx
  ja     L235
  sar    $0x3,%ecx
  or     %ecx,%r9d
L174:
  mov    %r9d,%ecx
  mov    %r9b,0x1(%rbp,%rax,4)
  mov    %ch,0x0(%rbp,%rax,4)
  mov    0x4(%r11,%r10,4),%r9d
  add    %r9d,%edi
  lea    0x5(%rdi),%edx
  cmp    $0xff,%edx
  ja     L225
  shl    $0x8,%edx
  and    $0xf800,%edx
L121:
  mov    %r9d,%edi
  sub    %r8d,%edi
  lea    0x2(%rdi),%ecx
  cmp    $0xff,%ecx
  ja     L135
  shl    $0x3,%ecx
  and    $0x7e0,%ecx
  or     %ecx,%edx
L82:
  add    %r9d,%esi
  lea    0x5(%rsi),%ecx
  cmp    $0xff,%ecx
  ja     L265
  sar    $0x3,%ecx
  or     %ecx,%edx
L0:
  mov    %dh,0x2(%rbp,%rax,4)
  mov    %dl,0x3(%rbp,%rax,4)
  add    $0x1,%rax
  cmp    $0x8,%rax
  je     L100
L46:
  mov    0x100(%r15,%rax,4),%edi
  mov    (%r15,%rax,4),%esi
  mov    %eax,%r12d
  mov    %eax,%r10d
  sar    $0x2,%r12d
  and    $0x3,%r10d
  imul   $0x32,%esi,%ecx
  imul   $0x82,%edi,%edx
  lea    0x80(%rcx,%rdx,1),%r8d
  mov    %r12d,%ecx
  shl    $0x5,%ecx
  sar    $0x8,%r8d
  add    %r10d,%ecx
  add    %ecx,%ecx
  movslq %ecx,%rdx
  mov    (%r11,%rdx,4),%r9d
  mov    %rdx,-0x40(%rsp)
  lea    (%rdi,%r9,1),%r13d
  lea    0x7(%r13),%edx
  cmp    $0xff,%edx
  jbe    L102
  xor    %edx,%edx
  mov    $0xf800,%ecx
  cmp    $0xfffffff8,%r13d
  cmovg  %ecx,%edx
  jmp    L130
  nopl   (%rax)
L265:
  mov    %edx,%ecx
  or     $0x1f,%ecx
  cmp    $0xfffffffa,%esi
  cmovg  %ecx,%edx
  jmp    L0
L135:
  mov    %edx,%ecx
  or     $0x7e0,%ecx
  cmp    $0xfffffffd,%edi
  cmovg  %ecx,%edx
  jmp    L82
  nopl   0x0(%rax,%rax,1)
L225:
  xor    %edx,%edx
  cmp    $0xfffffffa,%edi
  mov    $0xf800,%edi
  cmovg  %edi,%edx
  jmp    L121
  nopw   0x0(%rax,%rax,1)
L235:
  mov    %r9d,%ecx
  or     $0x1f,%ecx
  cmp    $0xfffffffc,%edx
  cmovg  %ecx,%r9d
  jmp    L174
  nopw   0x0(%rax,%rax,1)
L62:
  mov    %r9d,%ecx
  or     $0x7e0,%ecx
  cmp    $0xfffffffe,%r12d
  cmovg  %ecx,%r9d
  jmp    L120
  xchg   %ax,%ax
L11:
  xor    %r9d,%r9d
  cmp    $0xfffffffc,%ecx
  mov    $0xf800,%ecx
  cmovg  %ecx,%r9d
  jmp    L166
  nopl   0x0(%rax)
L17:
  mov    %r9d,%ecx
  or     $0x1f,%ecx
  cmp    $0xfffffffe,%edx
  cmovg  %ecx,%r9d
  jmp    L144
  nopw   0x0(%rax,%rax,1)
L170:
  mov    %r9d,%r13d
  or     $0x7e0,%r13d
  test   %ecx,%ecx
  cmovns %r13d,%r9d
  jmp    L27
  nopl   (%rax)
L122:
  xor    %r9d,%r9d
  cmp    $0xfffffffe,%ecx
  mov    $0xf800,%ecx
  cmovg  %ecx,%r9d
  jmp    L76
  nopl   0x0(%rax)
L272:
  mov    %edx,%r13d
  or     $0x1f,%r13d
  cmp    $0xfffffff8,%r9d
  cmovg  %r13d,%edx
  jmp    L195
  nopl   0x0(%rax)
L77:
  mov    %edx,%r13d
  or     $0x7e0,%r13d
  cmp    $0xfffffffc,%r14d
  cmovg  %r13d,%edx
  jmp    L201
  nop
L100:
  mov    -0x30(%rsp),%rax
  add    $0x20,%r15
  add    $0x40,%r11
  add    %rax,%rbx
  add    %rax,%rbp
  cmp    -0x38(%rsp),%r11
  jne    L231
  lea    0x200(%r11),%rax
  mov    -0x10(%rsp),%rbx
  subq   $0xffffffffffffff80,-0x28(%rsp)
  add    %rbx,-0x18(%rsp)
  add    %rbx,-0x20(%rsp)
  mov    %rax,-0x38(%rsp)
  cmp    -0x8(%rsp),%rax
  jne    L30
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void col221111_16(out, pic, width) int *out;
unsigned char *pic;
int width;
{
  int i, j, k;
  unsigned char *pic0, *pic1;
  int *outy, *outc;
  int cr, cg, cb, y;
  pic0 = pic;
  pic1 = pic + width;
  outy = out;
  outc = out + 64 * 4;
  for (i = 2; i > 0; i--) {
    for (j = 4; j > 0; j--) {
      for (k = 0; k < 8; k++) {
        ((cb = outc[0 + 0 * 8 + k], cr = outc[64 + 0 * 8 + k], cg = (50 * cb + 130 * cr + 128) >> 8), (y = outy[(k / 4 * 8 + 0) * 8 + (k & 3) * 2 + 0], y = ((((unsigned int)(y + cr + 3 * 2 + 1) >= 256 ? ((y + cr + 3 * 2 + 1) < 0 ? 0 : 255) : (y + cr + 3 * 2 + 1)) & 0xf8) << 8) | ((((unsigned int)(y - cg + 3) >= 256 ? ((y - cg + 3) < 0 ? 0 : 255) : (y - cg + 3)) & 0xfc) << 3) | ((((unsigned int)(y + cb + 3 * 2 + 1) >= 256 ? ((y + cb + 3 * 2 + 1) < 0 ? 0 : 255) : (y + cb + 3 * 2 + 1))) >> 3), pic0[(k * 2 + 0) * 2 + 0] = y >> 8, pic0[(k * 2 + 0) * 2 + 1] = y & 0xff), (y = outy[(k / 4 * 8 + 0) * 8 + (k & 3) * 2 + 1], y = ((((unsigned int)(y + cr + 0 * 2 + 1) >= 256 ? ((y + cr + 0 * 2 + 1) < 0 ? 0 : 255) : (y + cr + 0 * 2 + 1)) & 0xf8) << 8) | ((((unsigned int)(y - cg + 0) >= 256 ? ((y - cg + 0) < 0 ? 0 : 255) : (y - cg + 0)) & 0xfc) << 3) | ((((unsigned int)(y + cb + 0 * 2 + 1) >= 256 ? ((y + cb + 0 * 2 + 1) < 0 ? 0 : 255) : (y + cb + 0 * 2 + 1))) >> 3), pic0[(k * 2 + 1) * 2 + 0] = y >> 8, pic0[(k * 2 + 1) * 2 + 1] = y & 0xff), (y = outy[(k / 4 * 8 + 1) * 8 + (k & 3) * 2 + 0], y = ((((unsigned int)(y + cr + 1 * 2 + 1) >= 256 ? ((y + cr + 1 * 2 + 1) < 0 ? 0 : 255) : (y + cr + 1 * 2 + 1)) & 0xf8) << 8) | ((((unsigned int)(y - cg + 1) >= 256 ? ((y - cg + 1) < 0 ? 0 : 255) : (y - cg + 1)) & 0xfc) << 3) | ((((unsigned int)(y + cb + 1 * 2 + 1) >= 256 ? ((y + cb + 1 * 2 + 1) < 0 ? 0 : 255) : (y + cb + 1 * 2 + 1))) >> 3), pic1[(k * 2 + 0) * 2 + 0] = y >> 8, pic1[(k * 2 + 0) * 2 + 1] = y & 0xff), (y = outy[(k / 4 * 8 + 1) * 8 + (k & 3) * 2 + 1], y = ((((unsigned int)(y + cr + 2 * 2 + 1) >= 256 ? ((y + cr + 2 * 2 + 1) < 0 ? 0 : 255) : (y + cr + 2 * 2 + 1)) & 0xf8) << 8) | ((((unsigned int)(y - cg + 2) >= 256 ? ((y - cg + 2) < 0 ? 0 : 255) : (y - cg + 2)) & 0xfc) << 3) | ((((unsigned int)(y + cb + 2 * 2 + 1) >= 256 ? ((y + cb + 2 * 2 + 1) < 0 ? 0 : 255) : (y + cb + 2 * 2 + 1))) >> 3), pic1[(k * 2 + 1) * 2 + 0] = y >> 8, pic1[(k * 2 + 1) * 2 + 1] = y & 0xff));
      }
      outc += 8;
      outy += 16;
      pic0 += 2 * width;
      pic1 += 2 * width;
    }
    outy += 64 * 2 - 16 * 4;
  }
}
```