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

```asm
<BGBBTJ_JPG_TransIDCT_Vert2>:
  endbr64
  mov    %rdi,%rax
  push   %r15
  mov    %rsi,%rdi
  mov    %edx,%r9d
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x40(%rax),%esi
  imul   $0x5b,(%rax),%r10d
  mov    0xc0(%rax),%r11d
  imul   $0x76,%esi,%ebp
  mov    0x20(%rax),%r13d
  imul   $0x31,%esi,%esi
  imul   $0x5b,0x80(%rax),%ebx
  mov    %r10d,%ecx
  imul   $0x76,%r11d,%r12d
  lea    (%r10,%rbp,1),%r8d
  sub    %esi,%ecx
  imul   $0x31,%r11d,%r11d
  lea    (%r10,%rsi,1),%edx
  sub    %ebp,%r10d
  sub    %ebx,%edx
  sub    %ebx,%ecx
  lea    (%r10,%rbx,1),%esi
  add    %ebx,%r8d
  sub    %r12d,%edx
  imul   $0x7e,%r13d,%r10d
  add    %r12d,%ecx
  mov    0x60(%rax),%r12d
  add    %r11d,%r8d
  sub    %r11d,%esi
  mov    0xa0(%rax),%ebx
  mov    0xe0(%rax),%eax
  imul   $0x6a,%r12d,%r11d
  mov    %r12d,%r15d
  imul   $0xffffffe7,%r12d,%ebp
  mov    %ebx,%r14d
  shl    $0x6,%r15d
  shl    $0x6,%r14d
  add    %r11d,%r10d
  imul   $0x47,%ebx,%r11d
  add    %r11d,%r10d
  lea    (%rax,%rax,4),%r11d
  lea    (%r11,%r11,4),%r11d
  add    %r11d,%r10d
  imul   $0x6a,%r13d,%r11d
  add    %ebp,%r11d
  mov    %ebx,%ebp
  sub    %r14d,%ebp
  mov    %r12d,%r14d
  lea    (%r11,%rbp,2),%r11d
  sub    %r15d,%r14d
  imul   $0xffffffb9,%eax,%ebp
  add    %ebp,%r11d
  imul   $0x47,%r13d,%ebp
  lea    0x0(%r13,%r13,4),%r13d
  lea    0x0(%r13,%r13,4),%r13d
  lea    0x0(%rbp,%r14,2),%ebp
  lea    (%rbx,%rbx,4),%r14d
  lea    (%r14,%r14,4),%r14d
  add    %r14d,%ebp
  imul   $0x6a,%eax,%r14d
  imul   $0xffffffb9,%r12d,%r12d
  imul   $0x6a,%ebx,%ebx
  add    %r14d,%ebp
  add    %r13d,%r12d
  add    %r12d,%ebx
  mov    %eax,%r12d
  shl    $0x6,%r12d
  sub    %r12d,%eax
  lea    (%rbx,%rax,2),%ebx
  lea    (%r8,%r10,1),%eax
  sar    $0x10,%eax
  add    %r9d,%eax
  mov    %eax,(%rdi)
  lea    (%rdx,%r11,1),%eax
  sub    %r11d,%edx
  sar    $0x10,%eax
  sar    $0x10,%edx
  add    %r9d,%eax
  add    %r9d,%edx
  mov    %eax,0x20(%rdi)
  lea    (%rcx,%rbp,1),%eax
  sub    %ebp,%ecx
  sar    $0x10,%eax
  sar    $0x10,%ecx
  mov    %edx,0xc0(%rdi)
  add    %r9d,%eax
  add    %r9d,%ecx
  mov    %eax,0x40(%rdi)
  lea    (%rsi,%rbx,1),%eax
  sub    %ebx,%esi
  pop    %rbx
  sar    $0x10,%eax
  sar    $0x10,%esi
  pop    %rbp
  mov    %ecx,0xa0(%rdi)
  add    %r9d,%eax
  add    %r9d,%esi
  pop    %r12
  pop    %r13
  mov    %eax,0x60(%rdi)
  mov    %r8d,%eax
  pop    %r14
  sub    %r10d,%eax
  mov    %esi,0x80(%rdi)
  pop    %r15
  sar    $0x10,%eax
  add    %r9d,%eax
  mov    %eax,0xe0(%rdi)
  ret
```

assistant: ```c
void BGBBTJ_JPG_TransIDCT_Vert2(int *iblk, int *oblk, int dcbias) {
  int a, b, c, d, e, f;
  int i, j, k, l;
  int m, n, o, p;
  a = iblk[0] * 91;
  b = iblk[32] * 91;
  c = iblk[16] * 118;
  d = iblk[16] * 49;
  e = iblk[48] * 118;
  f = iblk[48] * 49;
  i = a + c + b + f;
  j = a + d - b - e;
  k = a - d - b + e;
  l = a - c + b - f;
  m = iblk[8] * 126 + iblk[24] * 106 + iblk[40] * 71 + iblk[56] * 25;
  n = iblk[8] * 106 - iblk[24] * 25 - iblk[40] * 126 - iblk[56] * 71;
  o = iblk[8] * 71 - iblk[24] * 126 + iblk[40] * 25 + iblk[56] * 106;
  p = iblk[8] * 25 - iblk[24] * 71 + iblk[40] * 106 - iblk[56] * 126;
  oblk[0] = ((i + m) >> 16) + dcbias;
  oblk[8] = ((j + n) >> 16) + dcbias;
  oblk[16] = ((k + o) >> 16) + dcbias;
  oblk[24] = ((l + p) >> 16) + dcbias;
  oblk[32] = ((l - p) >> 16) + dcbias;
  oblk[40] = ((k - o) >> 16) + dcbias;
  oblk[48] = ((j - n) >> 16) + dcbias;
  oblk[56] = ((i - m) >> 16) + dcbias;
}
```