user: Please show me the C implementation that matches the assembly code listed below:

```asm
<cxpose>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  mov    %rsi,%rax
  mov    %rsi,-0x10(%rsp)
  mov    %rdx,%r14
  mov    %rdx,-0x20(%rsp)
  mov    %rcx,-0x8(%rsp)
  mov    %r8,%rbp
  mov    %r8,-0x28(%rsp)
  mov    %r9,%r12
  lea    (%rsi,%rsi,1),%rdi
  lea    0x0(,%rsi,4),%rsi
  lea    (%rdi,%rsi,1),%r8
  shl    $0x3,%rax
  lea    0x3(%rbp),%rdx
  test   %rbp,%rbp
  cmovns %rbp,%rdx
  mov    %rdx,%r15
  sar    $0x2,%r15
  mov    %r15,-0x18(%rsp)
  cmp    $0x3,%rbp
  jle    L47
  lea    0x1(%rdi),%r9
  lea    0x1(%rsi),%r11
  lea    0x1(%r8),%r10
  shl    $0x2,%rdi
  shl    $0x2,%r9
  shl    $0x2,%rsi
  shl    $0x2,%r11
  shl    $0x2,%r8
  shl    $0x2,%r10
  lea    0x0(,%rcx,8),%rbp
  shl    $0x2,%rax
  mov    %rax,-0x30(%rsp)
  mov    %rbx,%r13
  jmp    L67
L73:
  movss  (%rax),%xmm7
  movss  0x4(%rax),%xmm6
  movss  (%rax,%rdi,1),%xmm5
  movss  (%rax,%r9,1),%xmm4
  movss  (%rax,%rsi,1),%xmm3
  movss  (%rax,%r11,1),%xmm2
  movss  (%rax,%r8,1),%xmm1
  movss  (%rax,%r10,1),%xmm0
  movss  %xmm7,(%rdx)
  movss  %xmm6,0x4(%rdx)
  movss  %xmm5,0x8(%rdx)
  movss  %xmm4,0xc(%rdx)
  movss  %xmm3,0x10(%rdx)
  movss  %xmm2,0x14(%rdx)
  movss  %xmm1,0x18(%rdx)
  movss  %xmm0,0x1c(%rdx)
  add    $0x8,%rax
  add    %rbp,%rdx
  sub    $0x1,%rcx
  jne    L73
L20:
  add    -0x30(%rsp),%r13
  add    $0x20,%r14
  sub    $0x1,%r15
  test   %r15,%r15
  jle    L104
L67:
  mov    %r12,%rcx
  mov    %r14,%rdx
  mov    %r13,%rax
  test   %r12,%r12
  jg     L73
  jmp    L20
L104:
  mov    -0x28(%rsp),%rcx
  cmp    $0x3,%rcx
  mov    $0x1,%eax
  mov    -0x18(%rsp),%rdi
  cmovg  %rdi,%rax
  imul   -0x30(%rsp),%rax
  add    %rax,%rbx
  mov    %rdi,%rax
  shl    $0x5,%rax
  cmp    $0x3,%rcx
  mov    $0x20,%edx
  cmovle %rdx,%rax
  add    %rax,-0x20(%rsp)
L47:
  mov    -0x28(%rsp),%rdi
  test   $0x3,%dil
  je     L75
  test   %r12,%r12
  jle    L75
  mov    %rdi,%rax
  sar    $0x3f,%rax
  shr    $0x3e,%rax
  add    %rax,%rdi
  and    $0x3,%edi
  sub    %rax,%rdi
  mov    -0x10(%rsp),%rsi
  shl    $0x3,%rsi
  mov    -0x8(%rsp),%r8
  shl    $0x3,%r8
  mov    -0x20(%rsp),%rax
  lea    (%rax,%rdi,8),%rcx
  mov    %rdi,%r9
  neg    %r9
  shl    $0x3,%r9
L86:
  lea    (%rcx,%r9,1),%rax
  test   %rdi,%rdi
  jle    L28
  mov    %rbx,%rdx
L105:
  movss  (%rdx),%xmm1
  movss  0x4(%rdx),%xmm0
  movss  %xmm1,(%rax)
  movss  %xmm0,0x4(%rax)
  add    $0x8,%rax
  add    %rsi,%rdx
  cmp    %rcx,%rax
  jne    L105
L28:
  add    $0x8,%rbx
  add    %r8,%rcx
  sub    $0x1,%r12
  jne    L86
L75:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void cxpose(float *indata, long iRsiz, float *outdata, long oRsiz, long Nrows, long Ncols) {
  float *irow;
  float *ocol;
  float *idata;
  float *odata;
  long RowCnt;
  long ColCnt;
  float T0r;
  float T0i;
  float T1r;
  float T1i;
  float T2r;
  float T2i;
  float T3r;
  float T3i;
  const long inRsizd1 = 2 * iRsiz;
  const long inRsizd1i = 2 * iRsiz + 1;
  const long inRsizd2 = 4 * iRsiz;
  const long inRsizd2i = 4 * iRsiz + 1;
  const long inRsizd3 = inRsizd2 + inRsizd1;
  const long inRsizd3i = inRsizd2 + inRsizd1 + 1;
  const long inRsizd4 = 8 * iRsiz;
  ocol = outdata;
  irow = indata;
  for (RowCnt = Nrows / 4; RowCnt > 0; RowCnt--) {
    idata = irow;
    odata = ocol;
    for (ColCnt = Ncols; ColCnt > 0; ColCnt--) {
      T0r = *idata;
      T0i = *(idata + 1);
      T1r = *(idata + inRsizd1);
      T1i = *(idata + inRsizd1i);
      T2r = *(idata + inRsizd2);
      T2i = *(idata + inRsizd2i);
      T3r = *(idata + inRsizd3);
      T3i = *(idata + inRsizd3i);
      *odata = T0r;
      *(odata + 1) = T0i;
      *(odata + 2) = T1r;
      *(odata + 3) = T1i;
      *(odata + 4) = T2r;
      *(odata + 5) = T2i;
      *(odata + 6) = T3r;
      *(odata + 7) = T3i;
      idata += 2;
      odata += 2 * oRsiz;
    }
    irow += inRsizd4;
    ocol += 8;
  }
  if (Nrows % 4 != 0) {
    for (ColCnt = Ncols; ColCnt > 0; ColCnt--) {
      idata = irow;
      odata = ocol;
      for (RowCnt = Nrows % 4; RowCnt > 0; RowCnt--) {
        T0r = *idata;
        T0i = *(idata + 1);
        *odata = T0r;
        *(odata + 1) = T0i;
        odata += 2;
        idata += 2 * iRsiz;
      }
      irow += 2;
      ocol += 2 * oRsiz;
    }
  }
}
```