user: What is the c source code of the assembly code below?

```asm
<col221111>:
  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,-0x10(%rsp)
  mov    %rax,-0x8(%rsp)
L226:
  mov    -0x38(%rsp),%rax
  mov    -0x28(%rsp),%r13
  mov    -0x20(%rsp),%r15
  mov    -0x18(%rsp),%r14
  lea    -0x100(%rax),%rbx
L68:
  mov    %r15,%rdi
  mov    %r14,%rsi
  xor    %r9d,%r9d
  jmp    L196
  xchg   %ax,%ax
L152:
  mov    %r10d,%r12d
  sub    %eax,%r12d
  mov    %r12b,0x1(%rsi)
  lea    (%r8,%r10,1),%r12d
  cmp    $0xff,%r12d
  ja     L170
L121:
  add    %r8d,%r10d
  mov    %r10b,0x2(%rsi)
L145:
  mov    0x4(%rbx,%rcx,4),%ecx
  lea    (%rdx,%rcx,1),%r10d
  cmp    $0xff,%r10d
  ja     L141
  lea    (%rcx,%rdx,1),%r10d
  mov    %r10b,0x3(%rsi)
L274:
  mov    %ecx,%r10d
  sub    %eax,%r10d
  cmp    $0xff,%r10d
  ja     L48
  mov    %ecx,%r10d
  sub    %eax,%r10d
  mov    %r10b,0x4(%rsi)
L157:
  lea    (%r8,%rcx,1),%r10d
  cmp    $0xff,%r10d
  ja     L168
  add    %r8d,%ecx
  mov    %cl,0x5(%rsi)
L268:
  lea    0x1(,%rbp,8),%ecx
  lea    (%r11,%rcx,4),%r10d
  add    %r10d,%r10d
  movslq %r10d,%r10
  mov    (%rbx,%r10,4),%ecx
  lea    (%rdx,%rcx,1),%r11d
  cmp    $0xff,%r11d
  ja     L267
  lea    (%rcx,%rdx,1),%r11d
  mov    %r11b,(%rdi)
L80:
  mov    %ecx,%r11d
  sub    %eax,%r11d
  cmp    $0xff,%r11d
  ja     L217
  mov    %ecx,%r11d
  sub    %eax,%r11d
  mov    %r11b,0x1(%rdi)
L16:
  lea    (%r8,%rcx,1),%r11d
  cmp    $0xff,%r11d
  ja     L235
  add    %r8d,%ecx
  mov    %cl,0x2(%rdi)
L208:
  mov    0x4(%rbx,%r10,4),%ecx
  lea    (%rdx,%rcx,1),%r10d
  cmp    $0xff,%r10d
  ja     L203
  add    %ecx,%edx
  mov    %dl,0x3(%rdi)
L8:
  mov    %ecx,%edx
  sub    %eax,%edx
  cmp    $0xff,%edx
  ja     L78
  mov    %ecx,%edx
  sub    %eax,%edx
  mov    %dl,0x4(%rdi)
L187:
  lea    (%r8,%rcx,1),%eax
  cmp    $0xff,%eax
  ja     L151
  lea    (%rcx,%r8,1),%eax
  mov    %al,0x5(%rdi)
L98:
  add    $0x1,%r9
  add    $0x6,%rdi
  add    $0x6,%rsi
  cmp    $0x8,%r9
  je     L17
L196:
  mov    0x100(%r13,%r9,4),%edx
  mov    0x0(%r13,%r9,4),%r8d
  mov    %r9d,%ebp
  mov    %r9d,%r11d
  sar    $0x2,%ebp
  and    $0x3,%r11d
  imul   $0x32,%r8d,%ecx
  imul   $0x82,%edx,%eax
  lea    0x80(%rcx,%rax,1),%eax
  mov    %ebp,%ecx
  shl    $0x5,%ecx
  sar    $0x8,%eax
  add    %r11d,%ecx
  add    %ecx,%ecx
  movslq %ecx,%rcx
  mov    (%rbx,%rcx,4),%r10d
  lea    (%rdx,%r10,1),%r12d
  cmp    $0xff,%r12d
  ja     L234
  lea    (%rdx,%r10,1),%r12d
  mov    %r12b,(%rsi)
L185:
  mov    %r10d,%r12d
  sub    %eax,%r12d
  cmp    $0xff,%r12d
  jbe    L152
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,0x1(%rsi)
  lea    (%r8,%r10,1),%r12d
  cmp    $0xff,%r12d
  jbe    L121
L170:
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,0x2(%rsi)
  jmp    L145
  nopl   0x0(%rax)
L234:
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,(%rsi)
  jmp    L185
  nopl   0x0(%rax)
L151:
  not    %eax
  sar    $0x1f,%eax
  mov    %al,0x5(%rdi)
  jmp    L98
  nopl   (%rax)
L78:
  not    %edx
  sar    $0x1f,%edx
  mov    %dl,0x4(%rdi)
  jmp    L187
  nopl   (%rax)
L203:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x3(%rdi)
  jmp    L8
  nopl   0x0(%rax,%rax,1)
L235:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,0x2(%rdi)
  jmp    L208
L217:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,0x1(%rdi)
  jmp    L16
L267:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,(%rdi)
  jmp    L80
  nop
L168:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x5(%rsi)
  jmp    L268
L48:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x4(%rsi)
  jmp    L157
L141:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x3(%rsi)
  jmp    L274
L17:
  mov    -0x30(%rsp),%rax
  add    $0x20,%r13
  add    $0x40,%rbx
  add    %rax,%r14
  add    %rax,%r15
  cmp    -0x38(%rsp),%rbx
  jne    L68
  lea    0x200(%rbx),%rax
  mov    -0x8(%rsp),%rsi
  subq   $0xffffffffffffff80,-0x28(%rsp)
  add    %rsi,-0x18(%rsp)
  add    %rsi,-0x20(%rsp)
  mov    %rax,-0x38(%rsp)
  cmp    %rax,-0x10(%rsp)
  jne    L226
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void col221111(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], ((pic0[(k * 2 + 0) * 3 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic0[(k * 2 + 0) * 3 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 0) * 3 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic0[(k * 2 + 0) * 3 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 0) * 3 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic0[(k * 2 + 0) * 3 + 2]) = (y + cb) < 0 ? 0 : 255) : 0)), (y = outy[(k / 4 * 8 + 0) * 8 + (k & 3) * 2 + 1], ((pic0[(k * 2 + 1) * 3 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic0[(k * 2 + 1) * 3 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 1) * 3 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic0[(k * 2 + 1) * 3 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 1) * 3 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic0[(k * 2 + 1) * 3 + 2]) = (y + cb) < 0 ? 0 : 255) : 0)), (y = outy[(k / 4 * 8 + 1) * 8 + (k & 3) * 2 + 0], ((pic1[(k * 2 + 0) * 3 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic1[(k * 2 + 0) * 3 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 0) * 3 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic1[(k * 2 + 0) * 3 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 0) * 3 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic1[(k * 2 + 0) * 3 + 2]) = (y + cb) < 0 ? 0 : 255) : 0)), (y = outy[(k / 4 * 8 + 1) * 8 + (k & 3) * 2 + 1], ((pic1[(k * 2 + 1) * 3 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic1[(k * 2 + 1) * 3 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 1) * 3 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic1[(k * 2 + 1) * 3 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 1) * 3 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic1[(k * 2 + 1) * 3 + 2]) = (y + cb) < 0 ? 0 : 255) : 0)));
      outc += 8;
      outy += 16;
      pic0 += 2 * width;
      pic1 += 2 * width;
    }
    outy += 64 * 2 - 16 * 4;
  }
}
```