user: What is the c source code of the assembly code below?

```asm
<radb3>:
  endbr64
  test   %rsi,%rsi
  je     L28
  push   %r15
  mov    %rcx,%r9
  lea    (%rsi,%rsi,1),%r10
  push   %r14
  push   %r13
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  cmp    $0x1,%rdi
  jne    L96
  lea    0x0(,%rsi,8),%r11
  lea    (%r10,%rsi,1),%r8
  lea    (%rcx,%r11,1),%rax
  mov    %rsi,%rcx
  lea    (%rdx,%r8,8),%rbx
  shl    $0x4,%rcx
  lea    (%rax,%r11,1),%rsi
  cmp    $0x8,%rcx
  setne  %dil
  cmp    $0x8,%r11
  setne  %bpl
  and    %ebp,%edi
  cmp    $0x1,%r12
  setne  %bpl
  and    %ebp,%edi
  mov    %r12,%rbp
  not    %rbp
  shl    $0x3,%rbp
  setne  %bpl
  and    %ebp,%edi
  cmp    %rsi,%rdx
  setae  %bpl
  cmp    %rbx,%rax
  setae  %r13b
  add    %rsi,%r11
  or     %r13d,%ebp
  and    %ebp,%edi
  cmp    %r11,%rdx
  setae  %r11b
  cmp    %rbx,%rsi
  setae  %bpl
  or     %ebp,%r11d
  test   %r11b,%dil
  je     L74
  cmp    %rax,%rdx
  setae  %dil
  cmp    %rbx,%r9
  setae  %r11b
  or     %r11b,%dil
  je     L74
  mov    %r12,%r8
  movapd D22(%rip),%xmm7
  mov    %r9,%rdi
  mov    %rdx,%rcx
  shr    %r8
  movapd D59(%rip),%xmm6
  shl    $0x4,%r8
  add    %rax,%r8
  nopl   0x0(%rax,%rax,1)
L14:
  movupd (%rcx),%xmm4
  add    $0x10,%rax
  add    $0x10,%rsi
  add    $0x10,%rdi
  movupd 0x20(%rcx),%xmm1
  movupd 0x10(%rcx),%xmm3
  add    $0x30,%rcx
  movapd %xmm4,%xmm0
  shufpd $0x1,%xmm4,%xmm0
  unpcklpd %xmm1,%xmm0
  movapd %xmm3,%xmm5
  unpcklpd %xmm4,%xmm3
  addpd  %xmm0,%xmm0
  movsd  %xmm4,%xmm5
  movsd  %xmm3,%xmm1
  movapd %xmm0,%xmm2
  addpd  %xmm5,%xmm0
  mulpd  %xmm7,%xmm2
  movups %xmm0,-0x10(%rdi)
  movapd %xmm1,%xmm0
  mulpd  %xmm6,%xmm0
  addpd  %xmm5,%xmm2
  movapd %xmm0,%xmm1
  addpd  %xmm2,%xmm1
  subpd  %xmm0,%xmm2
  movups %xmm1,-0x10(%rsi)
  movups %xmm2,-0x10(%rax)
  cmp    %r8,%rax
  jne    L14
  mov    %r12,%rax
  and    $0xfffffffffffffffe,%rax
  test   $0x1,%r12b
  je     L55
  lea    (%rax,%rax,2),%rcx
  movsd  D22(%rip),%xmm0
  shl    $0x3,%rcx
  movsd  0x8(%rdx,%rcx,1),%xmm1
  movsd  (%rdx,%rcx,1),%xmm2
  addsd  %xmm1,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm2
  movsd  D59(%rip),%xmm1
  movsd  %xmm2,(%r9,%rax,8)
  mulsd  0x10(%rdx,%rcx,1),%xmm1
  add    %r10,%rax
  movapd %xmm1,%xmm2
  addsd  %xmm0,%xmm2
  subsd  %xmm1,%xmm0
  movsd  %xmm2,(%r9,%rax,8)
  sub    %r12,%rax
  movsd  %xmm0,(%r9,%rax,8)
L55:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L96:
  mov    %rsi,%rbp
  mov    %r10,%r13
  lea    (%rdi,%rdi,2),%r11
  shl    $0x4,%rdi
  imul   %rbx,%r13
  lea    (%rdx,%rdi,1),%rcx
  shl    $0x3,%r11
  mov    %r9,%rax
  imul   %rbx,%rbp
  neg    %rdi
  xor    %esi,%esi
  movsd  D22(%rip),%xmm5
  movsd  D59(%rip),%xmm3
  lea    0x0(,%rbx,8),%r14
  nop
L10:
  movsd  -0x8(%rcx),%xmm1
  movsd  (%rcx,%rdi,1),%xmm2
  add    $0x1,%rsi
  addsd  %xmm1,%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm2,%xmm1
  mulsd  %xmm5,%xmm0
  movsd  %xmm1,(%rax)
  movsd  (%rcx),%xmm1
  add    %r11,%rcx
  mulsd  %xmm3,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm1,%xmm2
  subsd  %xmm1,%xmm0
  movsd  %xmm2,(%rax,%r13,8)
  movsd  %xmm0,(%rax,%rbp,8)
  add    %r14,%rax
  cmp    %rsi,%r12
  jne    L10
  cmp    $0x2,%rbx
  jbe    L55
  mov    %rbx,%rax
  mov    %r12,%r11
  lea    (%rbx,%rbx,2),%r15
  xor    %r14d,%r14d
  shl    $0x4,%rax
  lea    0x0(,%rbx,8),%r13
  movsd  D9(%rip),%xmm6
  shl    $0x3,%r15
  lea    -0x18(%rdx,%rax,1),%rbp
  mov    $0x28,%edx
  lea    (%r8,%r13,1),%rdi
  imul   %r13,%r11
  sub    %rax,%rdx
  mov    %rdi,-0x18(%rsp)
  imul   %r13,%r10
  mov    %rdx,-0x10(%rsp)
  add    %r9,%r11
  add    %r9,%r10
  nopl   0x0(%rax,%rax,1)
L57:
  mov    -0x10(%rsp),%rax
  mov    -0x18(%rsp),%rcx
  lea    0x20(%rbp),%rsi
  mov    %rbp,%rdx
  lea    (%rax,%rbp,1),%rdi
  mov    $0x2,%eax
  nopw   0x0(%rax,%rax,1)
L43:
  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(%r9,%rax,8)
  addsd  -0x10(%rdi),%xmm2
  movsd  %xmm2,(%r9,%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,(%r11,%rax,8)
  mulsd  -0x10(%r8,%rax,8),%xmm0
  mulsd  -0x8(%r8,%rax,8),%xmm4
  subsd  %xmm4,%xmm0
  movsd  %xmm0,-0x8(%r11,%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     L43
  add    $0x1,%r14
  add    %r15,%rbp
  add    %r13,%r9
  add    %r13,%r11
  add    %r13,%r10
  cmp    %r14,%r12
  jne    L57
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L74:
  lea    (%r9,%rcx,1),%rax
  mov    %r12,%rcx
  lea    (%r9,%r8,8),%rsi
  sub    %r10,%r12
  neg    %rcx
  movsd  D22(%rip),%xmm5
  movsd  D59(%rip),%xmm3
  shl    $0x4,%rcx
  nopw   0x0(%rax,%rax,1)
L94:
  movsd  0x8(%rdx),%xmm1
  movsd  (%rdx),%xmm2
  add    $0x18,%rdx
  addsd  %xmm1,%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm2,%xmm1
  mulsd  %xmm5,%xmm0
  movsd  %xmm1,(%rax,%rcx,1)
  movsd  -0x8(%rdx),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm1,%xmm2
  addsd  %xmm0,%xmm2
  subsd  %xmm1,%xmm0
  movsd  %xmm2,(%rax)
  movsd  %xmm0,(%rax,%r12,8)
  add    $0x8,%rax
  cmp    %rsi,%rax
  jne    L94
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L28:
  ret
```

tool_calls: name: disassemble_data
arguments: label: D22
data_type: byte[16]
name: disassemble_data
arguments: label: D59
data_type: byte[16]
name: disassemble_data
arguments: label: D9
data_type: f64

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

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

tool: D9: 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;
      }
    }
}
```