user: 请从给定的汇编指令推导出C源代码：

```asm
<CINTdmat_transpose>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %edx,-0x44(%rbp)
  mov    %ecx,-0x48(%rbp)
  movl   $0x0,-0x20(%rbp)
  jmp    L47
L142:
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  movl   $0x0,-0x24(%rbp)
  mov    -0x20(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  jmp    L88
L40:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x38(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x18(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x10(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x3,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x8(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x24(%rbp)
  mov    -0x48(%rbp),%eax
  add    %eax,-0x1c(%rbp)
L88:
  mov    -0x24(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L40
  mov    -0x44(%rbp),%eax
  shl    $0x2,%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x38(%rbp)
  addl   $0x4,-0x20(%rbp)
L47:
  mov    -0x48(%rbp),%eax
  sub    $0x3,%eax
  cmp    %eax,-0x20(%rbp)
  jl     L142
  mov    -0x48(%rbp),%eax
  sub    -0x20(%rbp),%eax
  cmp    $0x3,%eax
  je     L3
  cmp    $0x3,%eax
  jg     L20
  cmp    $0x1,%eax
  je     L76
  cmp    $0x2,%eax
  je     L121
  jmp    L20
L76:
  movl   $0x0,-0x24(%rbp)
  mov    -0x20(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  jmp    L10
L14:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x38(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x24(%rbp)
  mov    -0x48(%rbp),%eax
  add    %eax,-0x1c(%rbp)
L10:
  mov    -0x24(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L14
  jmp    L20
L121:
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x18(%rbp)
  movl   $0x0,-0x24(%rbp)
  mov    -0x20(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  jmp    L86
L166:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x38(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x18(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x24(%rbp)
  mov    -0x48(%rbp),%eax
  add    %eax,-0x1c(%rbp)
L86:
  mov    -0x24(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L166
  jmp    L20
L3:
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x10(%rbp)
  movl   $0x0,-0x24(%rbp)
  mov    -0x20(%rbp),%eax
  mov    %eax,-0x1c(%rbp)
  jmp    L134
L120:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x38(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x18(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rax,%rdx
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x10(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x24(%rbp)
  mov    -0x48(%rbp),%eax
  add    %eax,-0x1c(%rbp)
L134:
  mov    -0x24(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L120
  nop
L20:
  nop
  pop    %rbp
  ret
```

assistant: ```c
void CINTdmat_transpose(double *a_t, const double *a, const int m, const int n) {
  int i, j, k;
  double *pa1, *pa2, *pa3;
  for (j = 0; j < n - 3; j += 4) {
    pa1 = a_t + m;
    pa2 = pa1 + m;
    pa3 = pa2 + m;
    for (i = 0, k = j; i < m; i++, k += n) {
      a_t[i] = a[k + 0];
      pa1[i] = a[k + 1];
      pa2[i] = a[k + 2];
      pa3[i] = a[k + 3];
    }
    a_t += m * 4;
  }
  switch (n - j) {
  case 1:
    for (i = 0, k = j; i < m; i++, k += n) {
      a_t[i] = a[k];
    }
    break;
  case 2:
    pa1 = a_t + m;
    for (i = 0, k = j; i < m; i++, k += n) {
      a_t[i] = a[k + 0];
      pa1[i] = a[k + 1];
    }
    break;
  case 3:
    pa1 = a_t + m;
    pa2 = pa1 + m;
    for (i = 0, k = j; i < m; i++, k += n) {
      a_t[i] = a[k + 0];
      pa1[i] = a[k + 1];
      pa2[i] = a[k + 2];
    }
    break;
  }
}
```