user: What would the C source look like for the assembly code provided here:

```asm
<fdct_int32>:
  endbr64
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x118,%rsp
  mov    %rdi,%r9
  mov    %fs:0x28,%rax
  mov    %rax,0x108(%rsp)
  xor    %eax,%eax
  mov    %rdi,%r8
  mov    %rsp,%rax
  lea    0x100(%rsp),%r11
L17:
  movswl (%r8),%edx
  movswl 0xe(%r8),%ecx
  lea    (%rdx,%rcx,1),%r10d
  sub    %ecx,%edx
  movswl 0x2(%r8),%ecx
  movswl 0xc(%r8),%esi
  lea    (%rcx,%rsi,1),%ebp
  sub    %esi,%ecx
  movswl 0x4(%r8),%esi
  movswl 0xa(%r8),%edi
  lea    (%rsi,%rdi,1),%r12d
  sub    %edi,%esi
  movswl 0x6(%r8),%edi
  movswl 0x8(%r8),%ebx
  lea    (%rdi,%rbx,1),%r13d
  sub    %ebx,%edi
  lea    (%r10,%r13,1),%ebx
  sub    %r13d,%r10d
  lea    0x0(%rbp,%r12,1),%r13d
  sub    %r12d,%ebp
  lea    (%rbx,%r13,1),%r12d
  shl    $0x2,%r12d
  mov    %r12d,(%rax)
  sub    %r13d,%ebx
  shl    $0x2,%ebx
  mov    %ebx,0x10(%rax)
  lea    (%r10,%rbp,1),%ebx
  imul   $0x1151,%ebx,%ebx
  imul   $0x187e,%r10d,%r10d
  lea    0x400(%rbx,%r10,1),%r10d
  sar    $0xb,%r10d
  mov    %r10d,0x8(%rax)
  imul   $0xffffc4df,%ebp,%ebp
  lea    0x400(%rbx,%rbp,1),%r10d
  sar    $0xb,%r10d
  mov    %r10d,0x18(%rax)
  lea    (%rdx,%rdi,1),%ebx
  lea    (%rcx,%rsi,1),%r12d
  lea    (%rcx,%rdi,1),%ebp
  lea    (%rdx,%rsi,1),%r13d
  lea    0x0(%rbp,%r13,1),%r10d
  imul   $0x25a1,%r10d,%r10d
  imul   $0xffffe333,%ebx,%ebx
  imul   $0xffffadfd,%r12d,%r12d
  imul   $0xffffc13b,%ebp,%ebp
  imul   $0xfffff384,%r13d,%r13d
  add    %r10d,%ebp
  add    %r13d,%r10d
  imul   $0x98e,%edi,%edi
  add    %ebx,%edi
  lea    0x400(%rbp,%rdi,1),%edi
  sar    $0xb,%edi
  mov    %edi,0x1c(%rax)
  imul   $0x41b3,%esi,%esi
  add    %r12d,%esi
  lea    0x400(%r10,%rsi,1),%esi
  sar    $0xb,%esi
  mov    %esi,0x14(%rax)
  imul   $0x6254,%ecx,%ecx
  add    %r12d,%ecx
  lea    0x400(%rbp,%rcx,1),%ecx
  sar    $0xb,%ecx
  mov    %ecx,0xc(%rax)
  imul   $0x300b,%edx,%edx
  add    %ebx,%edx
  lea    0x400(%r10,%rdx,1),%edx
  sar    $0xb,%edx
  mov    %edx,0x4(%rax)
  add    $0x20,%rax
  add    $0x10,%r8
  cmp    %r11,%rax
  jne    L17
  mov    %rsp,%rax
  lea    0x20(%rsp),%r10
L21:
  mov    (%rax),%edx
  mov    %edx,%ebx
  add    0xe0(%rax),%ebx
  sub    0xe0(%rax),%edx
  mov    0x20(%rax),%ecx
  mov    %ecx,%r11d
  add    0xc0(%rax),%r11d
  sub    0xc0(%rax),%ecx
  mov    0x40(%rax),%esi
  mov    %esi,%ebp
  add    0xa0(%rax),%ebp
  sub    0xa0(%rax),%esi
  mov    0x60(%rax),%edi
  mov    %edi,%r12d
  add    0x80(%rax),%r12d
  sub    0x80(%rax),%edi
  lea    (%rbx,%r12,1),%r8d
  sub    %r12d,%ebx
  lea    (%r11,%rbp,1),%r12d
  sub    %ebp,%r11d
  lea    0x2(%r8,%r12,1),%ebp
  sar    $0x2,%ebp
  mov    %ebp,(%rax)
  sub    %r12d,%r8d
  add    $0x2,%r8d
  sar    $0x2,%r8d
  mov    %r8d,0x80(%rax)
  lea    (%rbx,%r11,1),%r8d
  imul   $0x1151,%r8d,%r8d
  imul   $0x187e,%ebx,%ebx
  lea    0x4000(%r8,%rbx,1),%ebx
  sar    $0xf,%ebx
  mov    %ebx,0x40(%rax)
  imul   $0xffffc4df,%r11d,%r11d
  lea    0x4000(%r8,%r11,1),%r8d
  sar    $0xf,%r8d
  mov    %r8d,0xc0(%rax)
  lea    (%rdx,%rdi,1),%r11d
  lea    (%rcx,%rsi,1),%ebp
  lea    (%rcx,%rdi,1),%ebx
  lea    (%rdx,%rsi,1),%r12d
  lea    (%rbx,%r12,1),%r8d
  imul   $0x25a1,%r8d,%r8d
  imul   $0xffffe333,%r11d,%r11d
  imul   $0xffffadfd,%ebp,%ebp
  imul   $0xffffc13b,%ebx,%ebx
  imul   $0xfffff384,%r12d,%r12d
  add    %r8d,%ebx
  add    %r12d,%r8d
  imul   $0x98e,%edi,%edi
  add    %r11d,%edi
  lea    0x4000(%rbx,%rdi,1),%edi
  sar    $0xf,%edi
  mov    %edi,0xe0(%rax)
  imul   $0x41b3,%esi,%esi
  add    %ebp,%esi
  lea    0x4000(%r8,%rsi,1),%esi
  sar    $0xf,%esi
  mov    %esi,0xa0(%rax)
  imul   $0x6254,%ecx,%ecx
  add    %ebp,%ecx
  lea    0x4000(%rbx,%rcx,1),%ecx
  sar    $0xf,%ecx
  mov    %ecx,0x60(%rax)
  imul   $0x300b,%edx,%edx
  add    %r11d,%edx
  lea    0x4000(%r8,%rdx,1),%edx
  sar    $0xf,%edx
  mov    %edx,0x20(%rax)
  add    $0x4,%rax
  cmp    %r10,%rax
  jne    L21
  mov    $0x0,%edx
L26:
  mov    (%rsp,%rdx,4),%eax
  add    $0x4,%eax
  sar    $0x3,%eax
  mov    %ax,(%r9,%rdx,2)
  add    $0x1,%rdx
  cmp    $0x40,%rdx
  jne    L26
  mov    0x108(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L19
  add    $0x118,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
L19:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void fdct_int32(short *const block) {
  int tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
  int tmp10, tmp11, tmp12, tmp13;
  int z1, z2, z3, z4, z5;
  short *blkptr;
  int *dataptr;
  int data[64];
  int i;
  dataptr = data;
  blkptr = block;
  for (i = 0; i < 8; i++) {
    tmp0 = blkptr[0] + blkptr[7];
    tmp7 = blkptr[0] - blkptr[7];
    tmp1 = blkptr[1] + blkptr[6];
    tmp6 = blkptr[1] - blkptr[6];
    tmp2 = blkptr[2] + blkptr[5];
    tmp5 = blkptr[2] - blkptr[5];
    tmp3 = blkptr[3] + blkptr[4];
    tmp4 = blkptr[3] - blkptr[4];
    tmp10 = tmp0 + tmp3;
    tmp13 = tmp0 - tmp3;
    tmp11 = tmp1 + tmp2;
    tmp12 = tmp1 - tmp2;
    dataptr[0] = (tmp10 + tmp11) << 2;
    dataptr[4] = (tmp10 - tmp11) << 2;
    z1 = (tmp12 + tmp13) * ((int)4433);
    dataptr[2] = (((z1 + tmp13 * ((int)6270)) + (((int)1) << ((13 - 2) - 1))) >> (13 - 2));
    dataptr[6] = (((z1 + tmp12 * (-((int)15137))) + (((int)1) << ((13 - 2) - 1))) >> (13 - 2));
    z1 = tmp4 + tmp7;
    z2 = tmp5 + tmp6;
    z3 = tmp4 + tmp6;
    z4 = tmp5 + tmp7;
    z5 = (z3 + z4) * ((int)9633);
    tmp4 *= ((int)2446);
    tmp5 *= ((int)16819);
    tmp6 *= ((int)25172);
    tmp7 *= ((int)12299);
    z1 *= -((int)7373);
    z2 *= -((int)20995);
    z3 *= -((int)16069);
    z4 *= -((int)3196);
    z3 += z5;
    z4 += z5;
    dataptr[7] = (((tmp4 + z1 + z3) + (((int)1) << ((13 - 2) - 1))) >> (13 - 2));
    dataptr[5] = (((tmp5 + z2 + z4) + (((int)1) << ((13 - 2) - 1))) >> (13 - 2));
    dataptr[3] = (((tmp6 + z2 + z3) + (((int)1) << ((13 - 2) - 1))) >> (13 - 2));
    dataptr[1] = (((tmp7 + z1 + z4) + (((int)1) << ((13 - 2) - 1))) >> (13 - 2));
    dataptr += 8;
    blkptr += 8;
  }
  dataptr = data;
  for (i = 0; i < 8; i++) {
    tmp0 = dataptr[0] + dataptr[56];
    tmp7 = dataptr[0] - dataptr[56];
    tmp1 = dataptr[8] + dataptr[48];
    tmp6 = dataptr[8] - dataptr[48];
    tmp2 = dataptr[16] + dataptr[40];
    tmp5 = dataptr[16] - dataptr[40];
    tmp3 = dataptr[24] + dataptr[32];
    tmp4 = dataptr[24] - dataptr[32];
    tmp10 = tmp0 + tmp3;
    tmp13 = tmp0 - tmp3;
    tmp11 = tmp1 + tmp2;
    tmp12 = tmp1 - tmp2;
    dataptr[0] = (((tmp10 + tmp11) + (((int)1) << ((2) - 1))) >> (2));
    dataptr[32] = (((tmp10 - tmp11) + (((int)1) << ((2) - 1))) >> (2));
    z1 = (tmp12 + tmp13) * ((int)4433);
    dataptr[16] = (((z1 + tmp13 * ((int)6270)) + (((int)1) << ((13 + 2) - 1))) >> (13 + 2));
    dataptr[48] = (((z1 + tmp12 * (-((int)15137))) + (((int)1) << ((13 + 2) - 1))) >> (13 + 2));
    z1 = tmp4 + tmp7;
    z2 = tmp5 + tmp6;
    z3 = tmp4 + tmp6;
    z4 = tmp5 + tmp7;
    z5 = (z3 + z4) * ((int)9633);
    tmp4 *= ((int)2446);
    tmp5 *= ((int)16819);
    tmp6 *= ((int)25172);
    tmp7 *= ((int)12299);
    z1 *= -((int)7373);
    z2 *= -((int)20995);
    z3 *= -((int)16069);
    z4 *= -((int)3196);
    z3 += z5;
    z4 += z5;
    dataptr[56] = (((tmp4 + z1 + z3) + (((int)1) << ((13 + 2) - 1))) >> (13 + 2));
    dataptr[40] = (((tmp5 + z2 + z4) + (((int)1) << ((13 + 2) - 1))) >> (13 + 2));
    dataptr[24] = (((tmp6 + z2 + z3) + (((int)1) << ((13 + 2) - 1))) >> (13 + 2));
    dataptr[8] = (((tmp7 + z1 + z4) + (((int)1) << ((13 + 2) - 1))) >> (13 + 2));
    dataptr++;
  }
  for (i = 0; i < 64; i++)
    block[i] = (short int)(((data[i]) + (((int)1) << ((3) - 1))) >> (3));
}
```