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

```asm
<MM_matrix_transpose_r2c>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %edx,-0x54(%rbp)
  mov    %ecx,-0x58(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  mov    -0x58(%rbp),%eax
  lea    0x3(%rax),%edx
  test   %eax,%eax
  cmovs  %edx,%eax
  sar    $0x2,%eax
  mov    %eax,-0x34(%rbp)
  movl   $0x0,-0x38(%rbp)
  jmp    L1
L27:
  mov    -0x34(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpacklo_epi32@plt>
  mov    %eax,-0x30(%rbp)
  mov    -0x34(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%ecx
  mov    -0x34(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpacklo_epi32@plt>
  mov    %eax,-0x2c(%rbp)
  mov    -0x34(%rbp),%eax
  lea    (%rax,%rax,1),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpackhi_epi32@plt>
  mov    %eax,-0x28(%rbp)
  mov    -0x34(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%ecx
  mov    -0x34(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpackhi_epi32@plt>
  mov    %eax,-0x24(%rbp)
  mov    -0x2c(%rbp),%edx
  mov    -0x30(%rbp),%eax
  mov    -0x38(%rbp),%ecx
  imul   -0x54(%rbp),%ecx
  movslq %ecx,%rcx
  lea    0x0(,%rcx,4),%rsi
  mov    -0x50(%rbp),%rcx
  lea    (%rsi,%rcx,1),%rbx
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpacklo_epi32@plt>
  mov    %eax,(%rbx)
  mov    -0x2c(%rbp),%edx
  mov    -0x30(%rbp),%eax
  mov    -0x38(%rbp),%ecx
  imul   -0x54(%rbp),%ecx
  movslq %ecx,%rcx
  add    $0x1,%rcx
  lea    0x0(,%rcx,4),%rsi
  mov    -0x50(%rbp),%rcx
  lea    (%rsi,%rcx,1),%rbx
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpackhi_epi32@plt>
  mov    %eax,(%rbx)
  mov    -0x24(%rbp),%edx
  mov    -0x28(%rbp),%eax
  mov    -0x38(%rbp),%ecx
  imul   -0x54(%rbp),%ecx
  movslq %ecx,%rcx
  add    $0x2,%rcx
  lea    0x0(,%rcx,4),%rsi
  mov    -0x50(%rbp),%rcx
  lea    (%rsi,%rcx,1),%rbx
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpacklo_epi32@plt>
  mov    %eax,(%rbx)
  mov    -0x24(%rbp),%edx
  mov    -0x28(%rbp),%eax
  mov    -0x38(%rbp),%ecx
  imul   -0x54(%rbp),%ecx
  movslq %ecx,%rcx
  add    $0x3,%rcx
  lea    0x0(,%rcx,4),%rsi
  mov    -0x50(%rbp),%rcx
  lea    (%rsi,%rcx,1),%rbx
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_unpackhi_epi32@plt>
  mov    %eax,(%rbx)
  addl   $0x1,-0x38(%rbp)
L1:
  mov    -0x38(%rbp),%eax
  cmp    -0x34(%rbp),%eax
  jl     L27
  nop
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L24
  call   <__stack_chk_fail@plt>
L24:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
void MM_matrix_transpose_r2c(__m128i *inbuf, __m128i *outbuf, uint32_t rows, uint32_t colums) {
  __m128i block[4];
  uint32_t col = colums / 4;
  uint32_t i;
  for (i = 0; i < col; i++) {
    block[0] = _mm_unpacklo_epi32(inbuf[i], inbuf[i + 2 * col]);
    block[1] = _mm_unpacklo_epi32(inbuf[i + col], inbuf[i + 3 * col]);
    block[2] = _mm_unpackhi_epi32(inbuf[i], inbuf[i + 2 * col]);
    block[3] = _mm_unpackhi_epi32(inbuf[i + col], inbuf[i + 3 * col]);
    outbuf[i * rows + 0] = _mm_unpacklo_epi32(block[0], block[1]);
    outbuf[i * rows + 1] = _mm_unpackhi_epi32(block[0], block[1]);
    outbuf[i * rows + 2] = _mm_unpacklo_epi32(block[2], block[3]);
    outbuf[i * rows + 3] = _mm_unpackhi_epi32(block[2], block[3]);
  }
}
```