user: 请从给定的汇编指令推导出C源代码：

```asm
<radb3>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdx,-0x10(%rsp)
  test   %rsi,%rsi
  je     L40
  lea    (%rsi,%rsi,1),%rax
  mov    %rsi,%r9
  mov    %rdi,%rbx
  mov    %rdi,%rsi
  lea    (%rdi,%rdi,2),%r15
  mov    %rax,-0x8(%rsp)
  shl    $0x4,%rsi
  mov    %rcx,%r10
  imul   %rdi,%rax
  mov    %rsi,%r14
  mov    %rcx,%r11
  shl    $0x3,%r15
  lea    0x0(,%rdi,8),%r12
  mov    %r9,%rdi
  add    %rsi,%rdx
  neg    %r14
  imul   %rbx,%rdi
  movsd  D5(%rip),%xmm5
  movsd  D27(%rip),%xmm3
  mov    %rax,%rbp
  mov    %rcx,%rax
  xor    %ecx,%ecx
  nopl   0x0(%rax)
L63:
  movsd  -0x8(%rdx),%xmm1
  movsd  (%rdx,%r14,1),%xmm2
  mov    %rcx,%r13
  lea    0x1(%rcx),%rcx
  addsd  %xmm1,%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm2,%xmm1
  mulsd  %xmm5,%xmm0
  movsd  %xmm1,(%rax)
  movsd  (%rdx),%xmm1
  add    %r15,%rdx
  mulsd  %xmm3,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm1,%xmm2
  subsd  %xmm1,%xmm0
  movsd  %xmm2,(%rax,%rbp,8)
  movsd  %xmm0,(%rax,%rdi,8)
  add    %r12,%rax
  cmp    %rcx,%r9
  jne    L63
  cmp    $0x1,%rbx
  je     L40
  mov    -0x10(%rsp),%rax
  imul   %r12,%r9
  movsd  D19(%rip),%xmm6
  xor    %r14d,%r14d
  lea    -0x18(%rax,%rsi,1),%rbp
  mov    -0x8(%rsp),%rax
  add    %r10,%r9
  imul   %r12,%rax
  add    %rax,%r10
  lea    (%r8,%r12,1),%rax
  mov    %rax,-0x10(%rsp)
  mov    $0x28,%eax
  sub    %rsi,%rax
  mov    %rax,-0x8(%rsp)
  nopl   0x0(%rax,%rax,1)
L32:
  cmp    $0x2,%rbx
  jbe    L28
  mov    -0x8(%rsp),%rax
  mov    -0x10(%rsp),%rcx
  lea    0x20(%rbp),%rsi
  mov    %rbp,%rdx
  lea    (%rax,%rbp,1),%rdi
  mov    $0x2,%eax
  nopl   0x0(%rax)
L38:
  movsd  (%rsi),%xmm4
  addsd  (%rdx),%xmm4
  add    $0x10,%rsi
  sub    $0x10,%rdx
  movsd  -0x8(%rdi),%xmm3
  movsd  -0x8(%rsi),%xmm2
  add    $0x10,%rcx
  add    $0x10,%rdi
  subsd  0x18(%rdx),%xmm2
  movapd %xmm4,%xmm0
  mulsd  %xmm5,%xmm0
  movapd %xmm2,%xmm1
  mulsd  %xmm5,%xmm1
  addsd  -0x10(%rdi),%xmm1
  addsd  %xmm3,%xmm0
  addsd  %xmm4,%xmm3
  movsd  %xmm3,-0x8(%r11,%rax,8)
  addsd  -0x10(%rdi),%xmm2
  movsd  %xmm2,(%r11,%rax,8)
  movsd  -0x8(%rsi),%xmm4
  addsd  0x18(%rdx),%xmm4
  movapd %xmm0,%xmm2
  movsd  -0x10(%rsi),%xmm3
  subsd  0x10(%rdx),%xmm3
  movsd  -0x8(%r8,%rax,8),%xmm7
  mulsd  %xmm6,%xmm4
  mulsd  %xmm6,%xmm3
  addsd  %xmm4,%xmm2
  subsd  %xmm4,%xmm0
  movapd %xmm1,%xmm4
  addsd  %xmm3,%xmm4
  subsd  %xmm3,%xmm1
  movsd  -0x10(%r8,%rax,8),%xmm3
  mulsd  %xmm0,%xmm7
  mulsd  %xmm4,%xmm3
  addsd  %xmm7,%xmm3
  movsd  %xmm3,(%r9,%rax,8)
  mulsd  -0x10(%r8,%rax,8),%xmm0
  mulsd  -0x8(%r8,%rax,8),%xmm4
  subsd  %xmm4,%xmm0
  movsd  %xmm0,-0x8(%r9,%rax,8)
  movsd  -0x10(%rcx),%xmm0
  movsd  -0x8(%rcx),%xmm3
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%r10,%rax,8)
  mulsd  -0x10(%rcx),%xmm2
  mulsd  -0x8(%rcx),%xmm1
  subsd  %xmm1,%xmm2
  movsd  %xmm2,-0x8(%r10,%rax,8)
  add    $0x2,%rax
  cmp    %rax,%rbx
  ja     L38
L28:
  lea    0x1(%r14),%rax
  add    %r15,%rbp
  add    %r12,%r11
  add    %r12,%r9
  add    %r12,%r10
  cmp    %r13,%r14
  je     L40
  mov    %rax,%r14
  jmp    L32
  xchg   %ax,%ax
L40:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: label: D5
data_type: byte[8]
name: read_data
arguments: label: D27
data_type: byte[8]
name: read_data
arguments: label: D19
data_type: f64

tool: {"D5": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0xbf]}

tool: {"D27": [0xaa,0x4c,0x58,0xe8,0x7a,0xb6,0xfb,0x3f]}

tool: D19: 0.8660254037844386

assistant: ```c
void radb3(size_t ido, size_t l1, const double *cc, double *ch, const double *wa) {
  const size_t cdim = 3;
  static const double taur = -0.5, taui = 0.86602540378443864676;
  size_t i, k, ic;
  double ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
  for (k = 0; k < l1; k++) {
    tr2 = 2 * cc[(ido - 1) + ido * ((1) + cdim * (k))];
    cr2 = cc[(0) + ido * ((0) + cdim * (k))] + taur * tr2;
    ch[(0) + ido * ((k) + l1 * (0))] = cc[(0) + ido * ((0) + cdim * (k))] + tr2;
    ci3 = 2 * taui * cc[(0) + ido * ((2) + cdim * (k))];
    {
      ch[(0) + ido * ((k) + l1 * (2))] = cr2 + ci3;
      ch[(0) + ido * ((k) + l1 * (1))] = cr2 - ci3;
    };
  }
  if (ido == 1)
    return;
  for (k = 0; k < l1; k++)
    for (i = 2; i < ido; i += 2) {
      ic = ido - i;
      tr2 = cc[(i - 1) + ido * ((2) + cdim * (k))] + cc[(ic - 1) + ido * ((1) + cdim * (k))];
      ti2 = cc[(i) + ido * ((2) + cdim * (k))] - cc[(ic) + ido * ((1) + cdim * (k))];
      cr2 = cc[(i - 1) + ido * ((0) + cdim * (k))] + taur * tr2;
      ci2 = cc[(i) + ido * ((0) + cdim * (k))] + taur * ti2;
      ch[(i - 1) + ido * ((k) + l1 * (0))] = cc[(i - 1) + ido * ((0) + cdim * (k))] + tr2;
      ch[(i) + ido * ((k) + l1 * (0))] = cc[(i) + ido * ((0) + cdim * (k))] + ti2;
      cr3 = taui * (cc[(i - 1) + ido * ((2) + cdim * (k))] - cc[(ic - 1) + ido * ((1) + cdim * (k))]);
      ci3 = taui * (cc[(i) + ido * ((2) + cdim * (k))] + cc[(ic) + ido * ((1) + cdim * (k))]);
      {
        dr3 = cr2 + ci3;
        dr2 = cr2 - ci3;
      }
      {
        di2 = ci2 + cr3;
        di3 = ci2 - cr3;
      }
      {
        ch[(i) + ido * ((k) + l1 * (1))] = wa[(i - 2) + (0) * ido] * di2 + wa[(i - 1) + (0) * ido] * dr2;
        ch[(i - 1) + ido * ((k) + l1 * (1))] = wa[(i - 2) + (0) * ido] * dr2 - wa[(i - 1) + (0) * ido] * di2;
      }
      {
        ch[(i) + ido * ((k) + l1 * (2))] = wa[(i - 2) + (1) * ido] * di3 + wa[(i - 1) + (1) * ido] * dr3;
        ch[(i - 1) + ido * ((k) + l1 * (2))] = wa[(i - 2) + (1) * ido] * dr3 - wa[(i - 1) + (1) * ido] * di3;
      }
    }
}
```