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

```asm
<bitrv2conj>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rsi,%r9
  mov    %rdx,%rax
  movl   $0x0,(%rsi)
  cmp    $0x8,%edi
  jle    L32
  mov    $0x1,%r10d
  lea    0x4(%rsi),%r11
  jmp    L31
L80:
  mov    %ebx,%r10d
L31:
  sar    %edi
  test   %r10d,%r10d
  jle    L122
  mov    %r9,%rdx
  lea    -0x1(%r10),%ecx
  lea    (%r11,%rcx,4),%r8
  movslq %r10d,%rsi
L33:
  mov    %edi,%ecx
  add    (%rdx),%ecx
  mov    %ecx,(%rdx,%rsi,4)
  add    $0x4,%rdx
  cmp    %r8,%rdx
  jne    L33
L122:
  lea    (%r10,%r10,1),%ebx
  shl    $0x4,%r10d
  cmp    %edi,%r10d
  jl     L80
  lea    (%rbx,%rbx,1),%esi
  je     L174
L59:
  movsd  0x8(%rax),%xmm0
  movq   D0(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rax)
  movslq %esi,%rdx
  lea    0x8(%rax,%rdx,8),%rdx
  movsd  (%rdx),%xmm0
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx)
  cmp    $0x1,%ebx
  jle    L95
  lea    0x4(%r9),%r10
  add    %ebx,%ebx
  mov    $0x2,%ebp
  mov    $0x1,%r8d
  movapd %xmm1,%xmm0
  movapd %xmm1,%xmm5
  jmp    L170
L174:
  test   %ebx,%ebx
  jle    L95
L132:
  lea    0x0(,%rbx,4),%r10d
  mov    %r9,%r12
  mov    %ebx,%ebx
  mov    %rbx,-0x8(%rsp)
  movl   $0x0,-0xc(%rsp)
  mov    $0x1,%r11d
  movq   D0(%rip),%xmm3
  movapd %xmm3,%xmm0
  jmp    L39
L47:
  mov    $0x1,%ebx
  mov    $0x2,%esi
  jmp    L132
L198:
  lea    (%r11,%r11,1),%ebp
  mov    $0x0,%edi
L86:
  mov    0x4(%rbx),%edx
  lea    (%rdx,%rdi,2),%ecx
  mov    %ebp,%edx
  add    (%r9,%rdi,4),%edx
  movslq %ecx,%r8
  shl    $0x3,%r8
  lea    (%rax,%r8,1),%r14
  movsd  (%r14),%xmm4
  lea    0x8(%rax,%r8,1),%r15
  movsd  (%r15),%xmm1
  xorpd  %xmm0,%xmm1
  movslq %edx,%r8
  shl    $0x3,%r8
  lea    (%rax,%r8,1),%r13
  movsd  0x0(%r13),%xmm5
  lea    0x8(%rax,%r8,1),%r8
  movsd  (%r8),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  %xmm5,(%r14)
  movsd  %xmm2,(%r15)
  movsd  %xmm4,0x0(%r13)
  movsd  %xmm1,(%r8)
  add    %esi,%ecx
  add    %r10d,%edx
  movslq %ecx,%r8
  shl    $0x3,%r8
  lea    (%rax,%r8,1),%r14
  movsd  (%r14),%xmm4
  lea    0x8(%rax,%r8,1),%r15
  movsd  (%r15),%xmm1
  xorpd  %xmm0,%xmm1
  movslq %edx,%r8
  shl    $0x3,%r8
  lea    (%rax,%r8,1),%r13
  movsd  0x0(%r13),%xmm5
  lea    0x8(%rax,%r8,1),%r8
  movsd  (%r8),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  %xmm5,(%r14)
  movsd  %xmm2,(%r15)
  movsd  %xmm4,0x0(%r13)
  movsd  %xmm1,(%r8)
  add    %esi,%ecx
  sub    %esi,%edx
  movslq %ecx,%r8
  shl    $0x3,%r8
  lea    (%rax,%r8,1),%r14
  movsd  (%r14),%xmm4
  lea    0x8(%rax,%r8,1),%r15
  movsd  (%r15),%xmm1
  xorpd  %xmm0,%xmm1
  movslq %edx,%r8
  shl    $0x3,%r8
  lea    (%rax,%r8,1),%r13
  movsd  0x0(%r13),%xmm5
  lea    0x8(%rax,%r8,1),%r8
  movsd  (%r8),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  %xmm5,(%r14)
  movsd  %xmm2,(%r15)
  movsd  %xmm4,0x0(%r13)
  movsd  %xmm1,(%r8)
  add    %esi,%ecx
  movslq %ecx,%rcx
  shl    $0x3,%rcx
  lea    (%rax,%rcx,1),%r8
  movsd  (%r8),%xmm4
  lea    0x8(%rax,%rcx,1),%rcx
  movsd  (%rcx),%xmm1
  xorpd  %xmm0,%xmm1
  add    %r10d,%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  lea    (%rax,%rdx,1),%r13
  movsd  0x0(%r13),%xmm5
  lea    0x8(%rax,%rdx,1),%rdx
  movsd  (%rdx),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  %xmm5,(%r8)
  movsd  %xmm2,(%rcx)
  movsd  %xmm4,0x0(%r13)
  movsd  %xmm1,(%rdx)
  add    $0x1,%rdi
  cmp    %r11,%rdi
  jne    L86
L121:
  add    $0x1,%r11
  add    $0x4,%r12
  addl   $0x2,-0xc(%rsp)
L39:
  mov    %r12,%rbx
  mov    -0xc(%rsp),%edx
  add    (%r12),%edx
  movslq %edx,%rcx
  lea    0x8(%rax,%rcx,8),%rcx
  movsd  (%rcx),%xmm1
  xorpd  %xmm3,%xmm1
  movsd  %xmm1,(%rcx)
  add    %esi,%edx
  lea    (%rsi,%rdx,1),%ecx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  lea    (%rax,%rdx,1),%r8
  movsd  (%r8),%xmm4
  lea    0x8(%rax,%rdx,1),%rbp
  movsd  0x0(%rbp),%xmm1
  xorpd  %xmm3,%xmm1
  movslq %ecx,%rdx
  shl    $0x3,%rdx
  lea    (%rax,%rdx,1),%rdi
  movsd  (%rdi),%xmm5
  lea    0x8(%rax,%rdx,1),%rdx
  movsd  (%rdx),%xmm2
  xorpd  %xmm3,%xmm2
  movsd  %xmm5,(%r8)
  movsd  %xmm2,0x0(%rbp)
  movsd  %xmm4,(%rdi)
  movsd  %xmm1,(%rdx)
  add    %esi,%ecx
  movslq %ecx,%rcx
  lea    0x8(%rax,%rcx,8),%rdx
  movsd  (%rdx),%xmm1
  xorpd  %xmm3,%xmm1
  movsd  %xmm1,(%rdx)
  cmp    -0x8(%rsp),%r11
  je     L95
  test   %r11d,%r11d
  jg     L198
  jmp    L121
L46:
  mov    (%r10),%edx
  lea    (%rdx,%rdi,2),%ecx
  mov    %r11d,%edx
  add    (%r9,%rdi,4),%edx
  movslq %ecx,%r12
  shl    $0x3,%r12
  lea    (%rax,%r12,1),%r14
  movsd  (%r14),%xmm3
  lea    0x8(%rax,%r12,1),%r15
  movsd  (%r15),%xmm1
  xorpd  %xmm0,%xmm1
  movslq %edx,%r12
  shl    $0x3,%r12
  lea    (%rax,%r12,1),%r13
  movsd  0x0(%r13),%xmm4
  lea    0x8(%rax,%r12,1),%r12
  movsd  (%r12),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  %xmm4,(%r14)
  movsd  %xmm2,(%r15)
  movsd  %xmm3,0x0(%r13)
  movsd  %xmm1,(%r12)
  add    %esi,%ecx
  movslq %ecx,%rcx
  shl    $0x3,%rcx
  lea    (%rax,%rcx,1),%r12
  movsd  (%r12),%xmm3
  lea    0x8(%rax,%rcx,1),%rcx
  movsd  (%rcx),%xmm1
  xorpd  %xmm0,%xmm1
  add    %esi,%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  lea    (%rax,%rdx,1),%r13
  movsd  0x0(%r13),%xmm4
  lea    0x8(%rax,%rdx,1),%rdx
  movsd  (%rdx),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  %xmm4,(%r12)
  movsd  %xmm2,(%rcx)
  movsd  %xmm3,0x0(%r13)
  movsd  %xmm1,(%rdx)
  add    $0x1,%rdi
  cmp    %r8,%rdi
  jne    L46
L126:
  mov    %ebp,%edx
  add    (%r10),%edx
  movslq %edx,%rcx
  lea    0x8(%rax,%rcx,8),%rcx
  movsd  (%rcx),%xmm1
  xorpd  %xmm5,%xmm1
  movsd  %xmm1,(%rcx)
  add    %esi,%edx
  movslq %edx,%rdx
  lea    0x8(%rax,%rdx,8),%rdx
  movsd  (%rdx),%xmm1
  xorpd  %xmm5,%xmm1
  movsd  %xmm1,(%rdx)
  add    $0x1,%r8
  add    $0x4,%r10
  add    $0x2,%ebp
  cmp    %ebx,%ebp
  je     L95
L170:
  mov    %ebp,%r11d
  mov    $0x0,%edi
  test   %r8d,%r8d
  jg     L46
  jmp    L126
L32:
  je     L47
  mov    $0x1,%ebx
  mov    $0x2,%esi
  jmp    L59
L95:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D0
data_type: byte[8]

tool: {"D0": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

assistant: ```c
void bitrv2conj(int n, int *ip, double *a) {
  int j, j1, k, k1, l, m, m2;
  double xr, xi, yr, yi;
  ip[0] = 0;
  l = n;
  m = 1;
  while ((m << 3) < l) {
    l >>= 1;
    for (j = 0; j < m; j++) {
      ip[m + j] = ip[j] + l;
    }
    m <<= 1;
  }
  m2 = 2 * m;
  if ((m << 3) == l) {
    for (k = 0; k < m; k++) {
      for (j = 0; j < k; j++) {
        j1 = 2 * j + ip[k];
        k1 = 2 * k + ip[j];
        xr = a[j1];
        xi = -a[j1 + 1];
        yr = a[k1];
        yi = -a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += m2;
        k1 += 2 * m2;
        xr = a[j1];
        xi = -a[j1 + 1];
        yr = a[k1];
        yi = -a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += m2;
        k1 -= m2;
        xr = a[j1];
        xi = -a[j1 + 1];
        yr = a[k1];
        yi = -a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += m2;
        k1 += 2 * m2;
        xr = a[j1];
        xi = -a[j1 + 1];
        yr = a[k1];
        yi = -a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
      }
      k1 = 2 * k + ip[k];
      a[k1 + 1] = -a[k1 + 1];
      j1 = k1 + m2;
      k1 = j1 + m2;
      xr = a[j1];
      xi = -a[j1 + 1];
      yr = a[k1];
      yi = -a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
      k1 += m2;
      a[k1 + 1] = -a[k1 + 1];
    }
  } else {
    a[1] = -a[1];
    a[m2 + 1] = -a[m2 + 1];
    for (k = 1; k < m; k++) {
      for (j = 0; j < k; j++) {
        j1 = 2 * j + ip[k];
        k1 = 2 * k + ip[j];
        xr = a[j1];
        xi = -a[j1 + 1];
        yr = a[k1];
        yi = -a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += m2;
        k1 += m2;
        xr = a[j1];
        xi = -a[j1 + 1];
        yr = a[k1];
        yi = -a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
      }
      k1 = 2 * k + ip[k];
      a[k1 + 1] = -a[k1 + 1];
      a[k1 + m2 + 1] = -a[k1 + m2 + 1];
    }
  }
}
```