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

```asm
<jpegdct>:
  endbr64
  push   %r15
  push   %r14
  mov    %rsi,%r14
  push   %r13
  mov    $0xfffffffffffffffe,%r13
  push   %r12
  mov    $0x2,%r12d
  push   %rbp
  push   %rbx
  sub    $0xe8,%rsp
  mov    %rdi,0x58(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0xd8(%rsp)
  xor    %eax,%eax
  lea    0x70(%rsp),%rax
  movq   $0x6,0x10(%rsp)
  mov    %rax,0x50(%rsp)
  lea    0x90(%rsp),%rax
  movq   $0xa,0x18(%rsp)
  movq   $0xe,0x8(%rsp)
  movq   $0xc,0x20(%rsp)
  movq   $0x80,0x68(%rsp)
  movq   $0x0,0x60(%rsp)
  mov    %rax,0x48(%rsp)
L11:
  mov    0x60(%rsp),%rdi
  lea    0x0(,%rdi,8),%rax
  mov    %rax,%rdx
  sub    %rdi,%rdx
  lea    0x8(,%rdx,8),%rbx
  mov    %rbx,0x28(%rsp)
  lea    0x4(,%rdx,4),%rbx
  mov    %rdi,%rdx
  sub    %rax,%rdx
  mov    %rbx,0x30(%rsp)
  lea    0x10(%rdx,%rdx,1),%rax
  mov    %rax,0x38(%rsp)
  lea    (%rdi,%rdi,2),%eax
  mov    0x58(%rsp),%rdi
  mov    %eax,0x44(%rsp)
  lea    0xd(%rax),%ebp
  mov    $0x8,%eax
  mov    %ax,0x42(%rsp)
L3:
  mov    0x8(%rsp),%rax
  mov    0x50(%rsp),%rsi
  mov    %rdi,%r9
  mov    0x48(%rsp),%r10
  lea    (%rax,%rdi,1),%r8
  mov    %rsi,%rdx
L40:
  movswl (%r9),%eax
  movswl (%r8),%r11d
  add    $0x8,%rdx
  add    %r12,%r9
  add    %r13,%r8
  sub    $0x8,%rsi
  lea    (%rax,%r11,1),%ecx
  sub    %r11d,%eax
  movslq %ecx,%rcx
  cltq
  mov    %rcx,-0x8(%rdx)
  mov    %rax,0x40(%rsi)
  cmp    %rdx,%r10
  jne    L40
  mov    0x70(%rsp),%rsi
  mov    0x78(%rsp),%rdx
  mov    0x88(%rsp),%rax
  mov    0x80(%rsp),%r9
  mov    0x44(%rsp),%ebx
  add    %rdx,%r9
  add    %rsi,%rax
  sub    0x80(%rsp),%rdx
  sub    0x88(%rsp),%rsi
  mov    %ebx,%ecx
  lea    (%rax,%r9,1),%r8
  mov    0x28(%rsp),%rbx
  sub    %r9,%rax
  sar    %cl,%r8
  sar    %cl,%rax
  mov    %r8w,(%rdi)
  mov    %ax,(%rdi,%rbx,1)
  lea    (%rdx,%rsi,1),%eax
  movswl 0x14(%r14),%ecx
  cwtl
  mov    0x30(%rsp),%rbx
  imul   %ecx,%eax
  movswq 0x12(%r14),%rcx
  imul   %rcx,%rsi
  mov    %ebp,%ecx
  sar    %cl,%rsi
  add    %eax,%esi
  mov    %si,(%rdi,%rbx,1)
  movswq 0x16(%r14),%rcx
  mov    0x20(%rsp),%rbx
  mov    0x98(%rsp),%rsi
  imul   %rcx,%rdx
  mov    %ebp,%ecx
  sar    %cl,%rdx
  mov    0x90(%rsp),%rcx
  add    %edx,%eax
  mov    %ax,(%rdi,%rbx,1)
  mov    0xa8(%rsp),%rax
  movswl 0x4(%r14),%edx
  movswl (%r14),%r8d
  lea    (%rax,%rcx,1),%r11d
  lea    (%rax,%rsi,1),%r15
  movswl 0x10(%r14),%r9d
  movswl %r11w,%r11d
  imul   %edx,%r11d
  movslq %r11d,%rdx
  mov    %rdx,0x70(%rsp)
  mov    0xa0(%rsp),%rdx
  lea    (%rdx,%rsi,1),%ebx
  movswl %bx,%ebx
  imul   %r8d,%ebx
  movslq %ebx,%r8
  mov    %r8,0x78(%rsp)
  lea    (%rcx,%rdx,1),%r8
  imul   0x8(%r14),%cx
  lea    (%r8,%r15,1),%r10d
  movswl %r8w,%r8d
  movswl %r15w,%r15d
  movswl %r10w,%r10d
  imul   %r9d,%r10d
  movswl 0x2(%r14),%r9d
  imul   %r8d,%r9d
  movswl 0x6(%r14),%r8d
  movslq %r10d,%r10
  imul   %r15d,%r8d
  movslq %r9d,%r9
  add    %r10,%r9
  movslq %r8d,%r8
  mov    %r9,0x80(%rsp)
  add    %r10,%r8
  lea    (%rcx,%r11,1),%r10d
  mov    %ebp,%ecx
  add    %r9d,%r10d
  mov    %r8,0x88(%rsp)
  movswl %r10w,%r10d
  sar    %cl,%r10d
  mov    0x8(%rsp),%rcx
  mov    %r10w,(%rdi,%rcx,1)
  imul   0xc(%r14),%si
  mov    %ebp,%ecx
  add    %ebx,%esi
  add    %r8d,%esi
  movswl %si,%esi
  sar    %cl,%esi
  mov    0x18(%rsp),%rcx
  mov    %si,(%rdi,%rcx,1)
  imul   0xa(%r14),%dx
  mov    %ebp,%ecx
  add    %ebx,%edx
  mov    0x10(%rsp),%rbx
  add    %r9d,%edx
  movswl %dx,%edx
  sar    %cl,%edx
  mov    %dx,(%rdi,%rbx,1)
  imul   0xe(%r14),%ax
  add    %r11d,%eax
  add    %r8d,%eax
  cwtl
  sar    %cl,%eax
  mov    %ax,(%rdi,%r12,1)
  add    0x38(%rsp),%rdi
  subw   $0x1,0x42(%rsp)
  jne    L3
  mov    0x68(%rsp),%rax
  mov    0x58(%rsp),%rdi
  add    $0xe,%r12
  sub    $0xe,%r13
  addq   $0x54,0x20(%rsp)
  lea    -0x80(%rdi,%rax,1),%rdi
  addq   $0x62,0x8(%rsp)
  sub    $0x70,%rax
  addq   $0x46,0x18(%rsp)
  addq   $0x2a,0x10(%rsp)
  cmpq   $0x1,0x60(%rsp)
  mov    %rdi,0x58(%rsp)
  mov    %rax,0x68(%rsp)
  je     L12
  movq   $0x1,0x60(%rsp)
  jmp    L11
  nopl   (%rax)
L12:
  mov    0xd8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L49
  add    $0xe8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L49:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void jpegdct(short *d, short *r) {
  long int t[12];
  short i, j, k, m, n, p;
  for (k = 1, m = 0, n = 13, p = 8; k <= 8; k += 7, m += 3, n += 3, p -= 7, d -= 64) {
    for (i = 0; i < 8; i++, d += p) {
      for (j = 0; j < 4; j++) {
        t[j] = d[k * j] + d[k * (7 - j)];
        t[7 - j] = d[k * j] - d[k * (7 - j)];
      }
      t[8] = t[0] + t[3];
      t[9] = t[0] - t[3];
      t[10] = t[1] + t[2];
      t[11] = t[1] - t[2];
      d[0] = (t[8] + t[10]) >> m;
      d[4 * k] = (t[8] - t[10]) >> m;
      t[8] = (short)(t[11] + t[9]) * r[10];
      d[2 * k] = t[8] + (short)((t[9] * r[9]) >> n);
      d[6 * k] = t[8] + (short)((t[11] * r[11]) >> n);
      t[0] = (short)(t[4] + t[7]) * r[2];
      t[1] = (short)(t[5] + t[6]) * r[0];
      t[2] = t[4] + t[6];
      t[3] = t[5] + t[7];
      t[8] = (short)(t[2] + t[3]) * r[8];
      t[2] = (short)t[2] * r[1] + t[8];
      t[3] = (short)t[3] * r[3] + t[8];
      d[7 * k] = (short)(t[4] * r[4] + t[0] + t[2]) >> n;
      d[5 * k] = (short)(t[5] * r[6] + t[1] + t[3]) >> n;
      d[3 * k] = (short)(t[6] * r[5] + t[1] + t[2]) >> n;
      d[1 * k] = (short)(t[7] * r[7] + t[0] + t[3]) >> n;
    }
  }
}
```