user: Please show me the C implementation that matches the assembly code listed below:

```asm
<packforblock13>:
  endbr64
  push   %r12
  mov    %rsi,%r9
  mov    %rdx,%r10
  mov    %rdi,%rcx
  push   %rbp
  push   %rbx
  mov    (%rsi),%rax
  mov    (%rdx),%r8
  mov    0x20(%rax),%esi
  mov    0x28(%rax),%edx
  mov    0x30(%rax),%ebx
  mov    0x58(%rax),%r11d
  sub    %edi,%esi
  sub    %edi,%edx
  mov    (%rax),%ebp
  mov    0xc0(%rax),%r12d
  sar    $0xc,%esi
  sub    %edi,%ebx
  add    %edx,%edx
  sub    %edi,%r11d
  shl    $0xe,%ebx
  or     %esi,%edx
  mov    0x48(%rax),%esi
  shl    $0xf,%r11d
  or     %edx,%ebx
  mov    0x50(%rax),%edx
  sub    %edi,%ebp
  sub    %edi,%esi
  sub    %edi,%edx
  sar    $0xb,%esi
  shl    $0x2,%edx
  or     %esi,%edx
  mov    0x70(%rax),%esi
  or     %edx,%r11d
  mov    0x78(%rax),%edx
  sub    %edi,%esi
  sub    %edi,%edx
  mov    0x80(%rax),%edi
  sar    $0xa,%esi
  shl    $0x3,%edx
  sub    %ecx,%edi
  or     %esi,%edx
  mov    0x98(%rax),%esi
  shl    $0x10,%edi
  or     %edx,%edi
  mov    0xa0(%rax),%edx
  sub    %ecx,%esi
  sar    $0x9,%esi
  sub    %ecx,%edx
  shl    $0x4,%edx
  or     %esi,%edx
  mov    0xa8(%rax),%esi
  sub    %ecx,%esi
  shl    $0x11,%esi
  or     %edx,%esi
  mov    0xc8(%rax),%edx
  sub    %ecx,%r12d
  sar    $0x8,%r12d
  sub    %ecx,%edx
  shl    $0x5,%edx
  or     %r12d,%edx
  mov    0xd0(%rax),%r12d
  sub    %ecx,%r12d
  shl    $0x12,%r12d
  or     %edx,%r12d
  mov    0xb0(%rax),%edx
  sub    %ecx,%edx
  shl    $0x1e,%edx
  or     %esi,%edx
  mov    0xd8(%rax),%esi
  movd   %edx,%xmm1
  mov    0xf0(%rax),%edx
  sub    %ecx,%esi
  shl    $0x1f,%esi
  sub    %ecx,%edx
  or     %r12d,%esi
  shl    $0x6,%edx
  mov    0x10(%rax),%r12d
  movd   %esi,%xmm3
  mov    0xe8(%rax),%esi
  sub    %ecx,%r12d
  punpckldq %xmm3,%xmm1
  sub    %ecx,%esi
  shl    $0x1a,%r12d
  sar    $0x7,%esi
  or     %esi,%edx
  mov    0xf8(%rax),%esi
  sub    %ecx,%esi
  shl    $0x13,%esi
  or     %edx,%esi
  mov    0x8(%rax),%edx
  sub    %ecx,%edx
  shl    $0xd,%edx
  or     %r12d,%edx
  or     %ebp,%edx
  movd   %edx,%xmm0
  mov    0x38(%rax),%edx
  sub    %ecx,%edx
  shl    $0x1b,%edx
  or     %edx,%ebx
  mov    0x60(%rax),%edx
  movd   %ebx,%xmm5
  sub    %ecx,%edx
  punpckldq %xmm5,%xmm0
  add    $0x100,%rax
  add    $0xd0,%r8
  shl    $0x1c,%edx
  or     %r11d,%edx
  mov    -0x78(%rax),%r11d
  movd   %edx,%xmm2
  sub    %ecx,%r11d
  mov    %r11d,%ecx
  shl    $0x1d,%ecx
  or     %edi,%ecx
  movd   %ecx,%xmm4
  punpckldq %xmm4,%xmm2
  punpcklqdq %xmm2,%xmm0
  movups %xmm0,-0xd0(%r8)
  movq   %xmm1,-0xc0(%r8)
  mov    %esi,-0xb8(%r8)
  pop    %rbx
  mov    %rax,(%r9)
  pop    %rbp
  mov    %r8,(%r10)
  pop    %r12
  ret
```

assistant: ```c
void packforblock13(const uint32_t base, const uint32_t **pin, uint8_t **pw) {
  uint64_t *pw64 = *(uint64_t **)pw;
  const uint32_t *in = *pin;
  uint64_t w0;
  uint64_t w1;
  uint64_t w2;
  uint64_t w3;
  uint64_t w4;
  uint64_t w5;
  uint64_t w6;
  w0 = (uint64_t)(in[0] - base);
  w0 |= (uint64_t)(in[1] - base) << 13;
  w0 |= (uint64_t)(in[2] - base) << 26;
  w0 |= (uint64_t)(in[3] - base) << 39;
  w0 |= (uint64_t)(in[4] - base) << 52;
  w1 = (uint64_t)(in[4] - base) >> 12;
  w1 |= (uint64_t)(in[5] - base) << 1;
  w1 |= (uint64_t)(in[6] - base) << 14;
  w1 |= (uint64_t)(in[7] - base) << 27;
  w1 |= (uint64_t)(in[8] - base) << 40;
  w1 |= (uint64_t)(in[9] - base) << 53;
  w2 = (uint64_t)(in[9] - base) >> 11;
  w2 |= (uint64_t)(in[10] - base) << 2;
  w2 |= (uint64_t)(in[11] - base) << 15;
  w2 |= (uint64_t)(in[12] - base) << 28;
  w2 |= (uint64_t)(in[13] - base) << 41;
  w2 |= (uint64_t)(in[14] - base) << 54;
  w3 = (uint64_t)(in[14] - base) >> 10;
  w3 |= (uint64_t)(in[15] - base) << 3;
  w3 |= (uint64_t)(in[16] - base) << 16;
  w3 |= (uint64_t)(in[17] - base) << 29;
  w3 |= (uint64_t)(in[18] - base) << 42;
  w3 |= (uint64_t)(in[19] - base) << 55;
  w4 = (uint64_t)(in[19] - base) >> 9;
  w4 |= (uint64_t)(in[20] - base) << 4;
  w4 |= (uint64_t)(in[21] - base) << 17;
  w4 |= (uint64_t)(in[22] - base) << 30;
  w4 |= (uint64_t)(in[23] - base) << 43;
  w4 |= (uint64_t)(in[24] - base) << 56;
  w5 = (uint64_t)(in[24] - base) >> 8;
  w5 |= (uint64_t)(in[25] - base) << 5;
  w5 |= (uint64_t)(in[26] - base) << 18;
  w5 |= (uint64_t)(in[27] - base) << 31;
  w5 |= (uint64_t)(in[28] - base) << 44;
  w5 |= (uint64_t)(in[29] - base) << 57;
  w6 = (uint64_t)(in[29] - base) >> 7;
  w6 |= (uint64_t)(in[30] - base) << 6;
  w6 |= (uint64_t)(in[31] - base) << 19;
  pw64[0] = w0;
  pw64[1] = w1;
  pw64[2] = w2;
  pw64[3] = w3;
  pw64[4] = w4;
  pw64[5] = w5;
  pw64[6] = w6;
  *pin += 32;
  *pw += 52;
}
```