user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<radb2>:
  endbr64
  push   %r13
  mov    %rsi,%r11
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  and    $0x1,%ebx
  test   %rsi,%rsi
  je     L12
  mov    %rcx,%rsi
  mov    %r11,%r12
  mov    %rdi,%rcx
  mov    %rdi,%r9
  imul   %rdi,%r12
  mov    %rdx,%r10
  mov    %rdx,%rax
  shl    $0x4,%rcx
  lea    0x0(,%rdi,8),%rbp
  mov    %rsi,%rdx
  xor    %edi,%edi
  nopl   0x0(%rax,%rax,1)
L5:
  movsd  (%rax),%xmm0
  addsd  -0x8(%rax,%rcx,1),%xmm0
  add    $0x1,%rdi
  movsd  %xmm0,(%rdx)
  movsd  (%rax),%xmm0
  subsd  -0x8(%rax,%rcx,1),%xmm0
  add    %rcx,%rax
  movsd  %xmm0,(%rdx,%r12,8)
  add    %rbp,%rdx
  cmp    %r11,%rdi
  jne    L5
  test   %rbx,%rbx
  jne    L41
  movsd  D5(%rip),%xmm1
  lea    (%r10,%rbp,1),%rdx
  lea    -0x8(%rsi,%r9,8),%rax
  nopl   0x0(%rax,%rax,1)
L53:
  movsd  -0x8(%rdx),%xmm0
  add    $0x1,%rbx
  addsd  %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  movsd  (%rdx),%xmm0
  add    %rcx,%rdx
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%r12,8)
  add    %rbp,%rax
  cmp    %r11,%rbx
  jne    L53
L41:
  cmp    $0x2,%r9
  jbe    L12
  lea    0x0(,%r9,8),%rbp
  mov    %r11,%rdi
  mov    %r9,%r12
  xor    %ebx,%ebx
  imul   %rbp,%rdi
  shl    $0x4,%r12
  add    $0x8,%r10
  lea    -0x20(%r12),%r13
  add    %rsi,%rdi
  nopl   0x0(%rax,%rax,1)
L61:
  add    $0x1,%rbx
  lea    0x0(%r13,%r10,1),%rcx
  mov    %r10,%rdx
  mov    $0x2,%eax
  nopl   0x0(%rax)
L2:
  movsd  (%rdx),%xmm0
  addsd  (%rcx),%xmm0
  add    $0x10,%rdx
  sub    $0x10,%rcx
  movsd  %xmm0,-0x8(%rsi,%rax,8)
  movsd  0x18(%rcx),%xmm3
  movsd  -0x8(%rdx),%xmm2
  movsd  -0x10(%rdx),%xmm0
  subsd  0x10(%rcx),%xmm0
  movapd %xmm2,%xmm1
  subsd  %xmm3,%xmm2
  addsd  %xmm3,%xmm1
  movsd  %xmm2,(%rsi,%rax,8)
  movsd  -0x10(%r8,%rax,8),%xmm2
  movsd  -0x8(%r8,%rax,8),%xmm3
  mulsd  %xmm1,%xmm2
  mulsd  %xmm0,%xmm3
  addsd  %xmm3,%xmm2
  movsd  %xmm2,(%rdi,%rax,8)
  mulsd  -0x10(%r8,%rax,8),%xmm0
  mulsd  -0x8(%r8,%rax,8),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rdi,%rax,8)
  add    $0x2,%rax
  cmp    %rax,%r9
  ja     L2
  add    %r12,%r10
  add    %rbp,%rsi
  add    %rbp,%rdi
  cmp    %r11,%rbx
  jb     L61
L12:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D5
data_type: byte[8]

tool: {"D5": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xc0]}

assistant: ```c
void radb2(size_t ido, size_t l1, const double *cc, double *ch, const double *wa) {
  const size_t cdim = 2;
  size_t i, k, ic;
  double ti2, tr2;
  for (k = 0; k < l1; k++) {
    ch[(0) + ido * ((k) + l1 * (0))] = cc[(0) + ido * ((0) + cdim * (k))] + cc[(ido - 1) + ido * ((1) + cdim * (k))];
    ch[(0) + ido * ((k) + l1 * (1))] = cc[(0) + ido * ((0) + cdim * (k))] - cc[(ido - 1) + ido * ((1) + cdim * (k))];
  }
  if ((ido & 1) == 0)
    for (k = 0; k < l1; k++) {
      ch[(ido - 1) + ido * ((k) + l1 * (0))] = 2 * cc[(ido - 1) + ido * ((0) + cdim * (k))];
      ch[(ido - 1) + ido * ((k) + l1 * (1))] = -2 * cc[(0) + ido * ((1) + cdim * (k))];
    }
  if (ido <= 2)
    return;
  for (k = 0; k < l1; ++k)
    for (i = 2; i < ido; i += 2) {
      ic = ido - i;
      {
        ch[(i - 1) + ido * ((k) + l1 * (0))] = cc[(i - 1) + ido * ((0) + cdim * (k))] + cc[(ic - 1) + ido * ((1) + cdim * (k))];
        tr2 = cc[(i - 1) + ido * ((0) + cdim * (k))] - cc[(ic - 1) + ido * ((1) + cdim * (k))];
      }
      {
        ti2 = cc[(i) + ido * ((0) + cdim * (k))] + cc[(ic) + ido * ((1) + cdim * (k))];
        ch[(i) + ido * ((k) + l1 * (0))] = cc[(i) + ido * ((0) + cdim * (k))] - cc[(ic) + ido * ((1) + cdim * (k))];
      }
      {
        ch[(i) + ido * ((k) + l1 * (1))] = wa[(i - 2) + (0) * ido] * ti2 + wa[(i - 1) + (0) * ido] * tr2;
        ch[(i - 1) + ido * ((k) + l1 * (1))] = wa[(i - 2) + (0) * ido] * tr2 - wa[(i - 1) + (0) * ido] * ti2;
      }
    }
}
```