user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<radb3>:
  endbr64
  test   %esi,%esi
  jle    L44
  push   %r15
  mov    %rdx,%rax
  mov    %rcx,%rdx
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  lea    (%rsi,%rsi,1),%ebx
  cmp    $0x1,%edi
  jne    L17
  movslq %esi,%r11
  lea    -0x1(%rsi),%ebp
  lea    0x8(%rax),%rdi
  lea    (%r11,%r11,2),%r9
  lea    0x0(,%r11,8),%r10
  mov    %r11,%rbx
  mov    %ebp,-0x10(%rsp)
  lea    -0x8(%rax,%r9,8),%r13
  shl    $0x3,%r9
  shl    $0x4,%rbx
  lea    (%rcx,%r10,1),%rcx
  lea    (%rax,%r9,1),%r12
  add    %rdx,%r9
  lea    (%rdx,%rbx,1),%r14
  cmp    $0x8,%r10
  setne  %r15b
  cmp    $0x1,%ebp
  lea    (%rcx,%r10,1),%r8
  seta   %bpl
  and    %r15d,%ebp
  lea    0x8(%r10),%r15
  cmp    %rbx,%r15
  setne  %bl
  and    %ebp,%ebx
  cmp    %rcx,%rax
  setae  %bpl
  cmp    %r12,%rdx
  setae  %r15b
  or     %r15d,%ebp
  and    %ebx,%ebp
  cmp    %r14,%rax
  setae  %bl
  cmp    %rcx,%r12
  setbe  %r15b
  or     %r15d,%ebx
  and    %ebx,%ebp
  cmp    %r9,%rax
  setae  %bl
  cmp    %r12,%r8
  setae  %r12b
  or     %r12d,%ebx
  and    %ebp,%ebx
  cmp    %r14,%rdi
  setae  %bpl
  cmp    %r13,%rcx
  setae  %r12b
  or     %r12d,%ebp
  and    %ebp,%ebx
  cmp    %rdi,%r9
  setbe  %r9b
  cmp    %r13,%r8
  setae  %bpl
  or     %ebp,%r9d
  test   %r9b,%bl
  je     L89
  cmp    %rdi,%rcx
  setbe  %r9b
  cmp    %r13,%rdx
  setae  %bl
  or     %bl,%r9b
  je     L89
  mov    %esi,%r11d
  movapd D31(%rip),%xmm5
  mov    %rax,%r9
  mov    %rdx,%r10
  shr    %r11d
  movapd D21(%rip),%xmm4
  sub    $0x1,%r11d
  shl    $0x4,%r11
  lea    0x10(%rdx,%r11,1),%r11
  cs nopw 0x0(%rax,%rax,1)
L60:
  add    $0x10,%r10
  add    $0x30,%rdi
  add    $0x30,%r9
  add    $0x10,%rcx
  movupd -0x20(%rdi),%xmm1
  movupd -0x20(%r9),%xmm0
  add    $0x10,%r8
  movupd -0x20(%r9),%xmm3
  movlpd -0x30(%rdi),%xmm1
  movhpd -0x30(%r9),%xmm0
  movhps -0x8(%r9),%xmm0
  addpd  %xmm1,%xmm1
  movlpd -0x30(%r9),%xmm3
  mulpd  %xmm4,%xmm0
  movapd %xmm1,%xmm2
  addpd  %xmm3,%xmm1
  mulpd  %xmm5,%xmm2
  movups %xmm1,-0x10(%r10)
  addpd  %xmm3,%xmm2
  movapd %xmm2,%xmm1
  addpd  %xmm0,%xmm2
  subpd  %xmm0,%xmm1
  movups %xmm1,-0x10(%rcx)
  movups %xmm2,-0x10(%r8)
  cmp    %r10,%r11
  jne    L60
  mov    %esi,%ecx
  and    $0xfffffffe,%ecx
  test   $0x1,%sil
  je     L62
  lea    (%rcx,%rcx,2),%edi
  movsd  D31(%rip),%xmm1
  lea    0x1(%rdi),%r8d
  movslq %r8d,%r8
  movsd  (%rax,%r8,8),%xmm0
  movslq %edi,%r8
  add    $0x2,%edi
  movsd  (%rax,%r8,8),%xmm2
  movslq %edi,%rdi
  mov    %ecx,%r8d
  add    %esi,%ecx
  addsd  %xmm0,%xmm0
  mulsd  %xmm0,%xmm1
  addsd  %xmm2,%xmm0
  movsd  %xmm0,(%rdx,%r8,8)
  movsd  D21(%rip),%xmm0
  mulsd  (%rax,%rdi,8),%xmm0
  movslq %ecx,%rax
  add    %esi,%ecx
  addsd  %xmm2,%xmm1
  movslq %ecx,%rcx
  movapd %xmm1,%xmm2
  addsd  %xmm0,%xmm1
  subsd  %xmm0,%xmm2
  movsd  %xmm2,(%rdx,%rax,8)
  movsd  %xmm1,(%rdx,%rcx,8)
L62:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L17:
  mov    %esi,%r13d
  mov    %ebx,%r12d
  mov    %rcx,%r10
  movslq %edi,%r14
  imul   %edi,%r13d
  lea    (%rdi,%rdi,1),%ecx
  mov    %r14,-0x10(%rsp)
  xor    %r11d,%r11d
  imul   %edi,%r12d
  lea    (%rcx,%rdi,1),%ebp
  movsd  D31(%rip),%xmm4
  movslq %ecx,%rcx
  lea    0x0(,%r14,8),%r15
  movslq %ebp,%rbp
  lea    (%rax,%rcx,8),%rcx
  neg    %r14
  movsd  D21(%rip),%xmm3
  shl    $0x3,%rbp
  shl    $0x4,%r14
  movslq %r13d,%r13
  movslq %r12d,%r12
  nopw   0x0(%rax,%rax,1)
L50:
  movsd  -0x8(%rcx),%xmm1
  movsd  (%rcx,%r14,1),%xmm2
  add    $0x1,%r11d
  addsd  %xmm1,%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm2,%xmm1
  mulsd  %xmm4,%xmm0
  movsd  %xmm1,(%r10)
  movsd  (%rcx),%xmm1
  add    %rbp,%rcx
  mulsd  %xmm3,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm2
  subsd  %xmm1,%xmm2
  addsd  %xmm1,%xmm0
  movsd  %xmm2,(%r10,%r13,8)
  movsd  %xmm0,(%r10,%r12,8)
  add    %r15,%r10
  cmp    %r11d,%esi
  jne    L50
  cmp    $0x2,%edi
  jle    L62
  mov    -0x10(%rsp),%rcx
  imul   %edi,%ebx
  lea    (%rdi,%rdi,2),%r12d
  xor    %r15d,%r15d
  movslq %r12d,%r12
  movsd  D11(%rip),%xmm5
  mov    %rcx,%rbp
  shl    $0x3,%rcx
  shl    $0x4,%rbp
  mov    %rcx,%r13
  mov    %esi,%ecx
  shl    $0x3,%r12
  imul   %edi,%ecx
  add    %rax,%rbp
  movslq %ecx,%rcx
  lea    (%rdx,%rcx,8),%r14
  movslq %ebx,%rcx
  lea    (%rdx,%rcx,8),%rbx
  xchg   %ax,%ax
L29:
  lea    -0x18(%rbp),%r11
  mov    %rbp,%r10
  mov    $0x2,%ecx
  nopl   0x0(%rax)
L88:
  movsd  0x8(%r10),%xmm1
  addsd  (%r11),%xmm1
  add    $0x10,%r10
  sub    $0x10,%r11
  movsd  -0x8(%rax,%rcx,8),%xmm0
  movapd %xmm1,%xmm3
  mulsd  %xmm4,%xmm3
  addsd  %xmm0,%xmm3
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rdx,%rcx,8)
  movsd  (%r10),%xmm2
  subsd  0x18(%r11),%xmm2
  movapd %xmm3,%xmm6
  movsd  (%rax,%rcx,8),%xmm1
  movapd %xmm2,%xmm0
  mulsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rdx,%rcx,8)
  movsd  (%r10),%xmm1
  addsd  0x18(%r11),%xmm1
  movsd  -0x8(%r10),%xmm2
  subsd  0x10(%r11),%xmm2
  movsd  -0x8(%r8,%rcx,8),%xmm7
  mulsd  %xmm5,%xmm1
  mulsd  %xmm5,%xmm2
  subsd  %xmm1,%xmm6
  addsd  %xmm1,%xmm3
  movapd %xmm0,%xmm1
  addsd  %xmm2,%xmm1
  subsd  %xmm2,%xmm0
  movsd  -0x10(%r8,%rcx,8),%xmm2
  mulsd  %xmm6,%xmm2
  mulsd  %xmm1,%xmm7
  subsd  %xmm7,%xmm2
  movsd  %xmm2,-0x8(%r14,%rcx,8)
  mulsd  -0x10(%r8,%rcx,8),%xmm1
  mulsd  -0x8(%r8,%rcx,8),%xmm6
  addsd  %xmm6,%xmm1
  movsd  %xmm1,(%r14,%rcx,8)
  movsd  -0x10(%r9,%rcx,8),%xmm1
  movsd  -0x8(%r9,%rcx,8),%xmm2
  mulsd  %xmm3,%xmm1
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  %xmm1,-0x8(%rbx,%rcx,8)
  mulsd  -0x10(%r9,%rcx,8),%xmm0
  mulsd  -0x8(%r9,%rcx,8),%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%rbx,%rcx,8)
  add    $0x2,%rcx
  cmp    %ecx,%edi
  jg     L88
  add    $0x1,%r15d
  add    %r12,%rbp
  add    %r12,%rax
  add    %r13,%rdx
  add    %r13,%r14
  add    %r13,%rbx
  cmp    %r15d,%esi
  jne    L29
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L89:
  mov    -0x10(%rsp),%esi
  movsd  D31(%rip),%xmm4
  neg    %r10
  movsd  D21(%rip),%xmm3
  add    %r11,%rsi
  lea    0x8(%rdx,%rsi,8),%rdx
  nopl   (%rax)
L1:
  movsd  0x8(%rax),%xmm1
  movsd  (%rax),%xmm2
  add    $0x18,%rax
  addsd  %xmm1,%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm2,%xmm1
  mulsd  %xmm4,%xmm0
  movsd  %xmm1,(%rcx,%r10,1)
  movsd  -0x8(%rax),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm2
  subsd  %xmm1,%xmm2
  addsd  %xmm1,%xmm0
  movsd  %xmm2,(%rcx)
  movsd  %xmm0,(%rcx,%r11,8)
  add    $0x8,%rcx
  cmp    %rdx,%rcx
  jne    L1
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L44:
  ret
```

tool_calls: name: disassemble_data
arguments: data_tag: D31
type_guessed: byte[16]
name: disassemble_data
arguments: data_tag: D21
type_guessed: byte[16]
name: disassemble_data
arguments: data_tag: D11
type_guessed: f64

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

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

tool: D11: 0.866025403784439

assistant: ```c
void radb3(int ido, int l1, const double cc[], double ch[], const double wa1[], const double wa2[]) {
  static const double taur = -0.5;
  static const double taui = 0.866025403784439;
  int 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 + (3 * k + 1) * ido];
    cr2 = cc[3 * k * ido] + taur * tr2;
    ch[k * ido] = cc[3 * k * ido] + tr2;
    ci3 = 2 * taui * cc[(3 * k + 2) * ido];
    ch[(k + l1) * ido] = cr2 - ci3;
    ch[(k + 2 * l1) * ido] = 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 + (3 * k + 2) * ido] + cc[ic - 1 + (3 * k + 1) * ido];
      cr2 = cc[i - 1 + 3 * k * ido] + taur * tr2;
      ch[i - 1 + k * ido] = cc[i - 1 + 3 * k * ido] + tr2;
      ti2 = cc[i + (3 * k + 2) * ido] - cc[ic + (3 * k + 1) * ido];
      ci2 = cc[i + 3 * k * ido] + taur * ti2;
      ch[i + k * ido] = cc[i + 3 * k * ido] + ti2;
      cr3 = taui * (cc[i - 1 + (3 * k + 2) * ido] - cc[ic - 1 + (3 * k + 1) * ido]);
      ci3 = taui * (cc[i + (3 * k + 2) * ido] + cc[ic + (3 * k + 1) * ido]);
      dr2 = cr2 - ci3;
      dr3 = cr2 + ci3;
      di2 = ci2 + cr3;
      di3 = ci2 - cr3;
      ch[i - 1 + (k + l1) * ido] = wa1[i - 2] * dr2 - wa1[i - 1] * di2;
      ch[i + (k + l1) * ido] = wa1[i - 2] * di2 + wa1[i - 1] * dr2;
      ch[i - 1 + (k + 2 * l1) * ido] = wa2[i - 2] * dr3 - wa2[i - 1] * di3;
      ch[i + (k + 2 * l1) * ido] = wa2[i - 2] * di3 + wa2[i - 1] * dr3;
    }
  }
}
```