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

```asm
<scbitrevR2>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%r11
  mov    %rdi,-0x40(%rsp)
  mov    %esi,%r8d
  mov    %rdx,%r13
  lea    0x3(%rsi),%eax
  mov    %eax,%ecx
  shr    $0x1f,%ecx
  add    %eax,%ecx
  sar    %ecx
  mov    $0x1,%edx
  mov    %rdx,%r10
  shl    %cl,%r10
  mov    %rdx,%rdi
  mov    %esi,%ecx
  shl    %cl,%rdi
  mov    %edi,%esi
  sub    %r10d,%esi
  mov    %r8d,%eax
  shr    $0x1f,%eax
  add    %r8d,%eax
  sar    %eax
  lea    -0x1(%rax),%ecx
  shl    %cl,%rdx
  lea    -0x1(%rdx),%r15d
  mov    %r15d,-0x34(%rsp)
  add    $0x1,%r8d
  mov    %r8d,%ecx
  shr    $0x1f,%ecx
  add    %r8d,%ecx
  sar    %ecx
  add    $0x1,%ecx
  mov    %ecx,-0x18(%rsp)
  mov    %r10d,%r10d
  shl    $0x2,%r10
  mov    %r10,-0x60(%rsp)
  mov    %r11,%rbx
  add    %r10,%rbx
  mov    %rbx,-0x8(%rsp)
  cmp    %rbx,%r11
  jae    L58
  lea    0x1(%rdi),%r9d
  lea    0x2(%rdi),%r8d
  lea    0x3(%rdi),%ecx
  movslq %r15d,%rbx
  mov    %rbx,-0x20(%rsp)
  lea    0x0(%r13,%rbx,2),%rbx
  mov    %rbx,-0x30(%rsp)
  mov    %edi,%edi
  lea    0x0(,%rdi,4),%r15
  mov    %r9d,%r9d
  lea    0x0(,%r9,4),%r14
  mov    %r8d,%r8d
  lea    0x0(,%r8,4),%rbp
  mov    %ecx,%ecx
  lea    0x0(,%rcx,4),%r12
  add    $0x1,%eax
  mov    $0x4,%edi
  mov    %eax,%ecx
  shl    %cl,%rdi
  mov    %rdi,-0x10(%rsp)
  mov    %esi,%eax
  mov    %rax,-0x28(%rsp)
  lea    -0x2(%rdx),%eax
  mov    %eax,-0x14(%rsp)
  mov    %rbp,%rbx
  mov    %r15,%rbp
  mov    %r13,%r15
  mov    %r12,%r13
L48:
  cmpl   $0x0,-0x34(%rsp)
  js     L56
  mov    -0x20(%rsp),%r12
L49:
  mov    %r12d,-0x38(%rsp)
  mov    %r12d,-0x54(%rsp)
  movswl (%r15,%r12,2),%eax
  shl    $0x2,%eax
  cltq
  add    -0x28(%rsp),%rax
  mov    -0x40(%rsp),%rsi
  lea    (%rsi,%rax,4),%r10
  mov    %r12d,%eax
  movzbl -0x18(%rsp),%ecx
  shl    %cl,%eax
  cltq
  lea    (%rsi,%rax,4),%r11
  mov    -0x30(%rsp),%rax
  movswl (%rax),%eax
  shl    $0x2,%eax
  cltq
  lea    (%r11,%rax,4),%rdx
  movss  (%r10),%xmm2
  movss  0x4(%r10),%xmm1
  movss  (%r10,%rbp,1),%xmm11
  movss  (%r10,%r14,1),%xmm10
  cmp    %r12d,-0x34(%rsp)
  jle    L11
  mov    -0x20(%rsp),%rcx
  mov    %r10,%rax
  mov    %r10,-0x50(%rsp)
  mov    %r12,-0x48(%rsp)
  mov    %r12d,%r10d
L44:
  movss  0x8(%rax),%xmm8
  movss  0xc(%rax),%xmm7
  mov    (%rax,%rbx,1),%r12d
  movss  (%rax,%r13,1),%xmm5
  movss  %xmm5,-0x58(%rsp)
  movss  (%rdx),%xmm6
  movss  0x4(%rdx),%xmm5
  lea    (%rdx,%rbp,1),%r9
  movss  (%r9),%xmm15
  lea    (%rdx,%r14,1),%r8
  movss  (%r8),%xmm14
  movss  0x8(%rdx),%xmm4
  movss  0xc(%rdx),%xmm3
  lea    (%rdx,%rbx,1),%rdi
  movss  (%rdi),%xmm13
  lea    (%rdx,%r13,1),%rsi
  movss  (%rsi),%xmm12
  movaps %xmm2,%xmm9
  addss  %xmm11,%xmm9
  mulss  %xmm0,%xmm9
  movss  %xmm9,(%rdx)
  movaps %xmm1,%xmm9
  addss  %xmm10,%xmm9
  mulss  %xmm0,%xmm9
  movss  %xmm9,0x4(%rdx)
  subss  %xmm11,%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,0x8(%rdx)
  subss  %xmm10,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,0xc(%rdx)
  movd   %r12d,%xmm1
  addss  %xmm8,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%r9)
  movss  -0x58(%rsp),%xmm2
  movaps %xmm2,%xmm1
  addss  %xmm7,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%r8)
  movd   %r12d,%xmm1
  subss  %xmm1,%xmm8
  mulss  %xmm0,%xmm8
  movss  %xmm8,(%rdi)
  subss  %xmm2,%xmm7
  mulss  %xmm0,%xmm7
  movss  %xmm7,(%rsi)
  movaps %xmm6,%xmm1
  addss  %xmm15,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%rax)
  movaps %xmm5,%xmm1
  addss  %xmm14,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,0x4(%rax)
  subss  %xmm15,%xmm6
  mulss  %xmm0,%xmm6
  movss  %xmm6,0x8(%rax)
  subss  %xmm14,%xmm5
  mulss  %xmm0,%xmm5
  movss  %xmm5,0xc(%rax)
  movaps %xmm4,%xmm1
  addss  %xmm13,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%rax,%rbp,1)
  movaps %xmm3,%xmm1
  addss  %xmm12,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%rax,%r14,1)
  subss  %xmm13,%xmm4
  mulss  %xmm0,%xmm4
  movss  %xmm4,(%rax,%rbx,1)
  subss  %xmm12,%xmm3
  mulss  %xmm0,%xmm3
  movss  %xmm3,(%rax,%r13,1)
  sub    -0x60(%rsp),%rax
  movss  (%rax),%xmm2
  movss  0x4(%rax),%xmm1
  movss  (%rax,%rbp,1),%xmm11
  movss  (%rax,%r14,1),%xmm10
  movswl -0x2(%r15,%rcx,2),%edx
  shl    $0x2,%edx
  movslq %edx,%rdx
  lea    (%r11,%rdx,4),%rdx
  sub    $0x1,%rcx
  cmp    %ecx,%r10d
  jl     L44
  mov    -0x50(%rsp),%r10
  mov    -0x48(%rsp),%r12
  mov    -0x14(%rsp),%edi
  mov    %edi,%eax
  sub    -0x38(%rsp),%eax
  mov    -0x54(%rsp),%ecx
  cmp    %ecx,%edi
  mov    $0x0,%ecx
  cmovl  %rcx,%rax
  mov    -0x60(%rsp),%rdi
  imul   %rdi,%rax
  sub    %rdi,%r10
  sub    %rax,%r10
L11:
  movss  0x8(%r10),%xmm4
  movss  0xc(%r10),%xmm3
  lea    (%r10,%rbx,1),%rdx
  movss  (%rdx),%xmm6
  lea    (%r10,%r13,1),%rax
  movss  (%rax),%xmm5
  movaps %xmm2,%xmm7
  addss  %xmm11,%xmm7
  mulss  %xmm0,%xmm7
  movss  %xmm7,(%r10)
  movaps %xmm1,%xmm7
  addss  %xmm10,%xmm7
  mulss  %xmm0,%xmm7
  movss  %xmm7,0x4(%r10)
  subss  %xmm11,%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,0x8(%r10)
  subss  %xmm10,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,0xc(%r10)
  movaps %xmm4,%xmm1
  addss  %xmm6,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%r10,%rbp,1)
  movaps %xmm3,%xmm1
  addss  %xmm5,%xmm1
  mulss  %xmm0,%xmm1
  movss  %xmm1,(%r10,%r14,1)
  subss  %xmm6,%xmm4
  mulss  %xmm0,%xmm4
  movss  %xmm4,(%rdx)
  subss  %xmm5,%xmm3
  mulss  %xmm0,%xmm3
  movss  %xmm3,(%rax)
  sub    $0x1,%r12
  test   %r12d,%r12d
  jns    L49
L56:
  mov    -0x10(%rsp),%rsi
  add    %rsi,-0x40(%rsp)
  mov    -0x40(%rsp),%rax
  cmp    %rax,-0x8(%rsp)
  ja     L48
L58:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void scbitrevR2(float *ioptr, int M, short *BRLow, float scale) {
  float f0r;
  float f0i;
  float f1r;
  float f1i;
  float f2r;
  float f2i;
  float f3r;
  float f3i;
  float f4r;
  float f4i;
  float f5r;
  float f5i;
  float f6r;
  float f6i;
  float f7r;
  float f7i;
  float t0r;
  float t0i;
  float t1r;
  float t1i;
  float *p0r;
  float *p1r;
  float *IOP;
  float *iolimit;
  int Colstart;
  int iCol;
  unsigned int posA;
  unsigned int posAi;
  unsigned int posB;
  unsigned int posBi;
  const unsigned int Nrems2 = ((unsigned long)1 << ((M + 3) / 2));
  const unsigned int Nroot_1_ColInc = ((unsigned long)1 << (M)) - Nrems2;
  const unsigned int Nroot_1 = ((unsigned long)1 << (M / 2 - 1)) - 1;
  const unsigned int ColstartShift = (M + 1) / 2 + 1;
  posA = ((unsigned long)1 << (M));
  posAi = posA + 1;
  posB = posA + 2;
  posBi = posB + 1;
  iolimit = ioptr + Nrems2;
  for (; ioptr < iolimit; ioptr += ((unsigned long)1 << (M / 2 + 1))) {
    for (Colstart = Nroot_1; Colstart >= 0; Colstart--) {
      iCol = Nroot_1;
      p0r = ioptr + Nroot_1_ColInc + BRLow[Colstart] * 4;
      IOP = ioptr + (Colstart << ColstartShift);
      p1r = IOP + BRLow[iCol] * 4;
      f0r = *(p0r);
      f0i = *(p0r + 1);
      f1r = *(p0r + posA);
      f1i = *(p0r + posAi);
      for (; iCol > Colstart;) {
        f2r = *(p0r + 2);
        f2i = *(p0r + (2 + 1));
        f3r = *(p0r + posB);
        f3i = *(p0r + posBi);
        f4r = *(p1r);
        f4i = *(p1r + 1);
        f5r = *(p1r + posA);
        f5i = *(p1r + posAi);
        f6r = *(p1r + 2);
        f6i = *(p1r + (2 + 1));
        f7r = *(p1r + posB);
        f7i = *(p1r + posBi);
        t0r = f0r + f1r;
        t0i = f0i + f1i;
        f1r = f0r - f1r;
        f1i = f0i - f1i;
        t1r = f2r + f3r;
        t1i = f2i + f3i;
        f3r = f2r - f3r;
        f3i = f2i - f3i;
        f0r = f4r + f5r;
        f0i = f4i + f5i;
        f5r = f4r - f5r;
        f5i = f4i - f5i;
        f2r = f6r + f7r;
        f2i = f6i + f7i;
        f7r = f6r - f7r;
        f7i = f6i - f7i;
        *(p1r) = scale * t0r;
        *(p1r + 1) = scale * t0i;
        *(p1r + 2) = scale * f1r;
        *(p1r + (2 + 1)) = scale * f1i;
        *(p1r + posA) = scale * t1r;
        *(p1r + posAi) = scale * t1i;
        *(p1r + posB) = scale * f3r;
        *(p1r + posBi) = scale * f3i;
        *(p0r) = scale * f0r;
        *(p0r + 1) = scale * f0i;
        *(p0r + 2) = scale * f5r;
        *(p0r + (2 + 1)) = scale * f5i;
        *(p0r + posA) = scale * f2r;
        *(p0r + posAi) = scale * f2i;
        *(p0r + posB) = scale * f7r;
        *(p0r + posBi) = scale * f7i;
        p0r -= Nrems2;
        f0r = *(p0r);
        f0i = *(p0r + 1);
        f1r = *(p0r + posA);
        f1i = *(p0r + posAi);
        iCol -= 1;
        p1r = IOP + BRLow[iCol] * 4;
      };
      f2r = *(p0r + 2);
      f2i = *(p0r + (2 + 1));
      f3r = *(p0r + posB);
      f3i = *(p0r + posBi);
      t0r = f0r + f1r;
      t0i = f0i + f1i;
      f1r = f0r - f1r;
      f1i = f0i - f1i;
      t1r = f2r + f3r;
      t1i = f2i + f3i;
      f3r = f2r - f3r;
      f3i = f2i - f3i;
      *(p0r) = scale * t0r;
      *(p0r + 1) = scale * t0i;
      *(p0r + 2) = scale * f1r;
      *(p0r + (2 + 1)) = scale * f1i;
      *(p0r + posA) = scale * t1r;
      *(p0r + posAi) = scale * t1i;
      *(p0r + posB) = scale * f3r;
      *(p0r + posBi) = scale * f3i;
    };
  };
}
```