user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<transpose_64x64>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x68(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movl   $0x55555555,-0x40(%rbp)
  movl   $0xaaaaaaaa,-0x3c(%rbp)
  movl   $0x33333333,-0x38(%rbp)
  movl   $0xcccccccc,-0x34(%rbp)
  movl   $0xf0f0f0f,-0x30(%rbp)
  movl   $0xf0f0f0f0,-0x2c(%rbp)
  movl   $0xff00ff,-0x28(%rbp)
  movl   $0xff00ff00,-0x24(%rbp)
  movl   $0xffff,-0x20(%rbp)
  movl   $0xffff0000,-0x1c(%rbp)
  movl   $0xffffffff,-0x18(%rbp)
  movl   $0x0,-0x14(%rbp)
  movl   $0x0,-0x58(%rbp)
  jmp    L55
L21:
  mov    -0x58(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    -0x58(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,4),%rcx
  mov    -0x68(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  addl   $0x1,-0x58(%rbp)
L55:
  cmpl   $0x3f,-0x58(%rbp)
  jle    L21
  movl   $0x5,-0x50(%rbp)
  jmp    L56
L76:
  mov    -0x50(%rbp),%eax
  mov    $0x1,%edx
  mov    %eax,%ecx
  shl    %cl,%edx
  mov    %edx,%eax
  mov    %eax,-0x4c(%rbp)
  movl   $0x0,-0x58(%rbp)
  jmp    L69
L63:
  mov    -0x58(%rbp),%eax
  mov    %eax,-0x54(%rbp)
  jmp    L16
L35:
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x50(%rbp),%eax
  cltq
  mov    -0x40(%rbp,%rax,8),%eax
  mov    %edx,%esi
  and    %eax,%esi
  mov    -0x54(%rbp),%edx
  mov    -0x4c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x50(%rbp),%eax
  cltq
  mov    -0x40(%rbp,%rax,8),%eax
  and    %eax,%edx
  mov    -0x4c(%rbp),%eax
  mov    %eax,%ecx
  shl    %cl,%edx
  mov    %edx,%eax
  or     %esi,%eax
  mov    %eax,-0x48(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x50(%rbp),%eax
  cltq
  mov    -0x3c(%rbp,%rax,8),%eax
  and    %eax,%edx
  mov    -0x4c(%rbp),%eax
  mov    %eax,%ecx
  sar    %cl,%edx
  mov    %edx,%ecx
  mov    -0x54(%rbp),%edx
  mov    -0x4c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x50(%rbp),%eax
  cltq
  mov    -0x3c(%rbp,%rax,8),%eax
  and    %edx,%eax
  or     %ecx,%eax
  mov    %eax,-0x44(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rax,%rdx
  mov    -0x48(%rbp),%eax
  mov    %eax,(%rdx)
  mov    -0x54(%rbp),%edx
  mov    -0x4c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rax,%rdx
  mov    -0x44(%rbp),%eax
  mov    %eax,(%rdx)
  addl   $0x1,-0x54(%rbp)
L16:
  mov    -0x58(%rbp),%edx
  mov    -0x4c(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x54(%rbp)
  jl     L35
  mov    -0x4c(%rbp),%eax
  add    %eax,%eax
  add    %eax,-0x58(%rbp)
L69:
  cmpl   $0x3f,-0x58(%rbp)
  jle    L63
  subl   $0x1,-0x50(%rbp)
L56:
  cmpl   $0x0,-0x50(%rbp)
  jns    L76
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L44
  call   <__stack_chk_fail@plt>
L44:
  leave
  ret
```

assistant: ```c
void transpose_64x64(uint64_t *out, uint64_t *in) {
  int i, j, s, d;
  uint64_t x, y;
  uint64_t masks[6][2] = {{0x5555555555555555, 0xAAAAAAAAAAAAAAAA}, {0x3333333333333333, 0xCCCCCCCCCCCCCCCC}, {0x0F0F0F0F0F0F0F0F, 0xF0F0F0F0F0F0F0F0}, {0x00FF00FF00FF00FF, 0xFF00FF00FF00FF00}, {0x0000FFFF0000FFFF, 0xFFFF0000FFFF0000}, {0x00000000FFFFFFFF, 0xFFFFFFFF00000000}};
  for (i = 0; i < 64; i++)
    out[i] = in[i];
  for (d = 5; d >= 0; d--) {
    s = 1 << d;
    for (i = 0; i < 64; i += s * 2)
      for (j = i; j < i + s; j++) {
        x = (out[j] & masks[d][0]) | ((out[j + s] & masks[d][0]) << s);
        y = ((out[j] & masks[d][1]) >> s) | (out[j + s] & masks[d][1]);
        out[j + 0] = x;
        out[j + s] = y;
      }
  }
}
```