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

```asm
<radfg>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0xd8,%rsp
  mov    %rdi,-0xb8(%rbp)
  mov    %rsi,-0xc0(%rbp)
  mov    %rdx,-0xc8(%rbp)
  mov    %rcx,-0xd0(%rbp)
  mov    %r8,-0xd8(%rbp)
  mov    %r9,-0xe0(%rbp)
  mov    -0xc0(%rbp),%rax
  mov    %rax,-0x60(%rbp)
  mov    -0xc0(%rbp),%rax
  add    $0x1,%rax
  shr    %rax
  mov    %rax,-0x58(%rbp)
  cmpq   $0x1,-0xb8(%rbp)
  je     L124
  mov    -0xd0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rcx
  mov    -0xe0(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
  movq   $0x1,-0x98(%rbp)
  jmp    L368
L521:
  movq   $0x0,-0x90(%rbp)
  jmp    L114
L190:
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rax,%rdx
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rcx
  mov    -0x90(%rbp),%rax
  add    %rcx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0xe0(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x98(%rbp),%rax
  sub    $0x1,%rax
  imul   -0xb8(%rbp),%rax
  add    $0x1,%rax
  mov    %rax,-0xa8(%rbp)
  movq   $0x2,-0xa0(%rbp)
  jmp    L146
L169:
  mov    -0xa8(%rbp),%rax
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0xa8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0xa8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x2,-0xa0(%rbp)
  addq   $0x2,-0xa8(%rbp)
L146:
  mov    -0xa0(%rbp),%rax
  cmp    -0xb8(%rbp),%rax
  jb     L169
  addq   $0x1,-0x90(%rbp)
L114:
  mov    -0x90(%rbp),%rax
  cmp    -0xc8(%rbp),%rax
  jb     L190
  addq   $0x1,-0x98(%rbp)
L368:
  mov    -0x98(%rbp),%rax
  cmp    -0xc0(%rbp),%rax
  jb     L521
  movq   $0x1,-0x98(%rbp)
  mov    -0xc0(%rbp),%rax
  sub    $0x1,%rax
  mov    %rax,-0x80(%rbp)
  jmp    L89
L451:
  movq   $0x0,-0x90(%rbp)
  jmp    L36
L254:
  movq   $0x2,-0xa0(%rbp)
  jmp    L35
L324:
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x2,-0xa0(%rbp)
L35:
  mov    -0xa0(%rbp),%rax
  cmp    -0xb8(%rbp),%rax
  jb     L324
  addq   $0x1,-0x90(%rbp)
L36:
  mov    -0x90(%rbp),%rax
  cmp    -0xc8(%rbp),%rax
  jb     L254
  addq   $0x1,-0x98(%rbp)
  subq   $0x1,-0x80(%rbp)
L89:
  mov    -0x98(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L451
  jmp    L76
L124:
  mov    -0xd0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rcx
  mov    -0xd8(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
L76:
  movq   $0x1,-0x98(%rbp)
  mov    -0xc0(%rbp),%rax
  sub    $0x1,%rax
  mov    %rax,-0x80(%rbp)
  jmp    L495
L421:
  movq   $0x0,-0x90(%rbp)
  jmp    L50
L270:
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x90(%rbp)
L50:
  mov    -0x90(%rbp),%rax
  cmp    -0xc8(%rbp),%rax
  jb     L270
  addq   $0x1,-0x98(%rbp)
  subq   $0x1,-0x80(%rbp)
L495:
  mov    -0x98(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L421
  mov    -0xc0(%rbp),%rax
  shl    $0x4,%eax
  mov    %eax,%edi
  call   <util_malloc_@plt>
  mov    %rax,-0x50(%rbp)
  movsd  D29(%rip),%xmm1
  mov    -0xc0(%rbp),%rax
  test   %rax,%rax
  js     L63
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L448
L63:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L448:
  divsd  %xmm0,%xmm1
  movsd  %xmm1,-0x48(%rbp)
  mov    -0x50(%rbp),%rax
  movsd  D18(%rip),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  add    $0x8,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x48(%rbp),%xmm0
  comisd D62(%rip),%xmm0
  jae    L524
  movsd  -0x48(%rbp),%xmm0
  cvttsd2si %xmm0,%rax
  jmp    L256
L524:
  movsd  -0x48(%rbp),%xmm0
  movsd  D62(%rip),%xmm1
  subsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%rax
  movabs $0x8000000000000000,%rdx
  xor    %rdx,%rax
L256:
  mov    -0xc0(%rbp),%rdx
  shl    $0x4,%rdx
  lea    -0x10(%rdx),%rcx
  mov    -0x50(%rbp),%rdx
  lea    (%rcx,%rdx,1),%rbx
  mov    %rax,%rdi
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0x50(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rbx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x48(%rbp),%xmm0
  comisd D62(%rip),%xmm0
  jae    L562
  movsd  -0x48(%rbp),%xmm0
  cvttsd2si %xmm0,%rax
  jmp    L358
L562:
  movsd  -0x48(%rbp),%xmm0
  movsd  D62(%rip),%xmm1
  subsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%rax
  movabs $0x8000000000000000,%rdx
  xor    %rdx,%rax
L358:
  mov    -0x50(%rbp),%rdx
  lea    0x18(%rdx),%rbx
  mov    %rax,%rdi
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0x50(%rbp),%rax
  add    $0x18,%rax
  movsd  (%rax),%xmm0
  mov    -0xc0(%rbp),%rax
  shl    $0x4,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movq   D13(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movq   $0x2,-0xa0(%rbp)
  jmp    L108
L231:
  mov    -0xa0(%rbp),%rax
  test   %rax,%rax
  js     L144
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L59
L144:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L59:
  mulsd  -0x48(%rbp),%xmm0
  comisd D62(%rip),%xmm0
  jae    L582
  cvttsd2si %xmm0,%rax
  jmp    L563
L582:
  movsd  D62(%rip),%xmm1
  subsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%rax
  movabs $0x8000000000000000,%rdx
  xor    %rdx,%rax
L563:
  mov    -0xc0(%rbp),%rdx
  sub    -0xa0(%rbp),%rdx
  mov    %rdx,%rcx
  shl    $0x4,%rcx
  mov    -0x50(%rbp),%rdx
  lea    (%rcx,%rdx,1),%rbx
  mov    %rax,%rdi
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0xa0(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rbx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa0(%rbp),%rax
  test   %rax,%rax
  js     L133
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L172
L133:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L172:
  mulsd  -0x48(%rbp),%xmm0
  comisd D62(%rip),%xmm0
  jae    L49
  cvttsd2si %xmm0,%rax
  jmp    L111
L49:
  movsd  D62(%rip),%xmm1
  subsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%rax
  movabs $0x8000000000000000,%rdx
  xor    %rdx,%rax
L111:
  mov    -0xa0(%rbp),%rdx
  shl    $0x4,%rdx
  lea    0x8(%rdx),%rcx
  mov    -0x50(%rbp),%rdx
  lea    (%rcx,%rdx,1),%rbx
  mov    %rax,%rdi
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0xa0(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc0(%rbp),%rax
  sub    -0xa0(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movq   D13(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0xa0(%rbp)
L108:
  mov    -0xc0(%rbp),%rax
  shr    %rax
  cmp    %rax,-0xa0(%rbp)
  jbe    L231
  movq   $0x1,-0x88(%rbp)
  mov    -0xc0(%rbp),%rax
  sub    $0x1,%rax
  mov    %rax,-0x78(%rbp)
  jmp    L593
L379:
  mov    -0x88(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x88(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movq   $0x0,-0x70(%rbp)
  jmp    L355
L198:
  mov    -0x70(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x70(%rbp),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  mov    -0xd0(%rbp),%rax
  imul   -0x88(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc0(%rbp),%rax
  sub    $0x1,%rax
  imul   -0xd0(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xd0(%rbp),%rax
  imul   -0x78(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  mulsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x70(%rbp)
L355:
  mov    -0x70(%rbp),%rax
  cmp    -0xd0(%rbp),%rax
  jb     L198
  mov    -0x88(%rbp),%rax
  add    %rax,%rax
  mov    %rax,-0x68(%rbp)
  movq   $0x2,-0x98(%rbp)
  mov    -0xc0(%rbp),%rax
  sub    $0x2,%rax
  mov    %rax,-0x80(%rbp)
  jmp    L54
L262:
  mov    -0x88(%rbp),%rax
  add    %rax,%rax
  add    %rax,-0x68(%rbp)
  mov    -0xc0(%rbp),%rax
  add    %rax,%rax
  cmp    %rax,-0x68(%rbp)
  jb     L314
  mov    -0xc0(%rbp),%rax
  add    %rax,%rax
  sub    %rax,-0x68(%rbp)
L314:
  mov    -0x68(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x68(%rbp),%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movq   $0x0,-0x70(%rbp)
  jmp    L498
L223:
  mov    -0xd0(%rbp),%rax
  imul   -0x88(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xd0(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  mov    -0xd0(%rbp),%rax
  imul   -0x88(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xd0(%rbp),%rax
  imul   -0x78(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xd0(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  mov    -0xd0(%rbp),%rax
  imul   -0x78(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x70(%rbp)
L498:
  mov    -0x70(%rbp),%rax
  cmp    -0xd0(%rbp),%rax
  jb     L223
  addq   $0x1,-0x98(%rbp)
  subq   $0x1,-0x80(%rbp)
L54:
  mov    -0x98(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L262
  addq   $0x1,-0x88(%rbp)
  subq   $0x1,-0x78(%rbp)
L593:
  mov    -0x88(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L379
  mov    -0x50(%rbp),%rax
  mov    %rax,%rdi
  call   <util_free_@plt>
  movq   $0x0,-0x50(%rbp)
  movq   $0x1,-0x98(%rbp)
  jmp    L564
L192:
  movq   $0x0,-0x70(%rbp)
  jmp    L435
L183:
  mov    -0x70(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xd0(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x70(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x70(%rbp)
L435:
  mov    -0x70(%rbp),%rax
  cmp    -0xd0(%rbp),%rax
  jb     L183
  addq   $0x1,-0x98(%rbp)
L564:
  mov    -0x98(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L192
  movq   $0x0,-0x90(%rbp)
  jmp    L356
L43:
  mov    -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  imul   -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0xe0(%rbp),%rax
  add    %rax,%rcx
  mov    -0x60(%rbp),%rax
  imul   -0x90(%rbp),%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rsi
  mov    -0xd8(%rbp),%rax
  add    %rsi,%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
  addq   $0x1,-0x90(%rbp)
L356:
  mov    -0x90(%rbp),%rax
  cmp    -0xc8(%rbp),%rax
  jb     L43
  movq   $0x1,-0x98(%rbp)
  jmp    L79
L472:
  mov    -0xc0(%rbp),%rax
  sub    -0x98(%rbp),%rax
  mov    %rax,-0x80(%rbp)
  mov    -0x98(%rbp),%rax
  add    %rax,%rax
  mov    %rax,-0x40(%rbp)
  movq   $0x0,-0x90(%rbp)
  jmp    L473
L442:
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rax,%rdx
  mov    -0x60(%rbp),%rax
  imul   -0x90(%rbp),%rax
  mov    %rax,%rcx
  mov    -0x40(%rbp),%rax
  add    %rcx,%rax
  imul   -0xb8(%rbp),%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0xd8(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rax,%rdx
  mov    -0x60(%rbp),%rax
  imul   -0x90(%rbp),%rax
  mov    %rax,%rcx
  mov    -0x40(%rbp),%rax
  add    %rcx,%rax
  imul   -0xb8(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0xd8(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x90(%rbp)
L473:
  mov    -0x90(%rbp),%rax
  cmp    -0xc8(%rbp),%rax
  jb     L442
  addq   $0x1,-0x98(%rbp)
L79:
  mov    -0x98(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L472
  cmpq   $0x1,-0xb8(%rbp)
  je     L318
  movq   $0x1,-0x98(%rbp)
  jmp    L496
L110:
  mov    -0xc0(%rbp),%rax
  sub    -0x98(%rbp),%rax
  mov    %rax,-0x80(%rbp)
  mov    -0x98(%rbp),%rax
  add    %rax,%rax
  mov    %rax,-0x40(%rbp)
  movq   $0x0,-0x90(%rbp)
  jmp    L484
L230:
  movq   $0x2,-0xa0(%rbp)
  jmp    L203
L307:
  mov    -0xb8(%rbp),%rax
  sub    -0xa0(%rbp),%rax
  mov    %rax,-0x38(%rbp)
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x60(%rbp),%rax
  imul   -0x90(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x60(%rbp),%rax
  imul   -0x90(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x60(%rbp),%rax
  imul   -0x90(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc8(%rbp),%rax
  imul   -0x80(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xc8(%rbp),%rax
  imul   -0x98(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x60(%rbp),%rax
  imul   -0x90(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  sub    $0x1,%rax
  imul   -0xb8(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x2,-0xa0(%rbp)
L203:
  mov    -0xa0(%rbp),%rax
  cmp    -0xb8(%rbp),%rax
  jb     L307
  addq   $0x1,-0x90(%rbp)
L484:
  mov    -0x90(%rbp),%rax
  cmp    -0xc8(%rbp),%rax
  jb     L230
  addq   $0x1,-0x98(%rbp)
L496:
  mov    -0x98(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L110
  jmp    L480
L318:
  nop
L480:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_label: D29
type_guessed: f64
name: parse_data
arguments: data_label: D18
type_guessed: f64
name: parse_data
arguments: data_label: D62
type_guessed: f32[4]
name: parse_data
arguments: data_label: D13
type_guessed: f32[2]

tool: D29: 6.283185307179586

tool: D18: 1

tool: {"D62": [0.0,448.0,0.0,0.0]}

tool: {"D13": [0.0,-0.0]}

assistant: ```c
void radfg(size_t ido, size_t ip, size_t l1, size_t idl1, double *cc, double *ch, const double *wa) {
  const size_t cdim = ip;
  static const double twopi = 6.28318530717958647692;
  size_t idij, ipph, i, j, k, l, j2, ic, jc, lc, ik;
  double ai1, ai2, ar1, ar2, arg;
  double *csarr;
  size_t aidx;
  ipph = (ip + 1) / 2;
  if (ido != 1) {
    memcpy(ch, cc, idl1 * sizeof(double));
    for (j = 1; j < ip; j++)
      for (k = 0; k < l1; k++) {
        ch[(0) + ido * ((k) + l1 * (j))] = cc[(0) + ido * ((k) + l1 * (j))];
        idij = (j - 1) * ido + 1;
        for (i = 2; i < ido; i += 2, idij += 2) {
          ch[(i - 1) + ido * ((k) + l1 * (j))] = wa[idij - 1] * cc[(i - 1) + ido * ((k) + l1 * (j))] + wa[idij] * cc[(i) + ido * ((k) + l1 * (j))];
          ch[(i) + ido * ((k) + l1 * (j))] = wa[idij - 1] * cc[(i) + ido * ((k) + l1 * (j))] - wa[idij] * cc[(i - 1) + ido * ((k) + l1 * (j))];
        }
      }
    for (j = 1, jc = ip - 1; j < ipph; j++, jc--)
      for (k = 0; k < l1; k++)
        for (i = 2; i < ido; i += 2) {
          {
            cc[(i - 1) + ido * ((k) + l1 * (j))] = ch[(i - 1) + ido * ((k) + l1 * (jc))] + ch[(i - 1) + ido * ((k) + l1 * (j))];
            cc[(i) + ido * ((k) + l1 * (jc))] = ch[(i - 1) + ido * ((k) + l1 * (jc))] - ch[(i - 1) + ido * ((k) + l1 * (j))];
          }
          {
            cc[(i) + ido * ((k) + l1 * (j))] = ch[(i) + ido * ((k) + l1 * (j))] + ch[(i) + ido * ((k) + l1 * (jc))];
            cc[(i - 1) + ido * ((k) + l1 * (jc))] = ch[(i) + ido * ((k) + l1 * (j))] - ch[(i) + ido * ((k) + l1 * (jc))];
          }
        }
  } else
    memcpy(cc, ch, idl1 * sizeof(double));
  for (j = 1, jc = ip - 1; j < ipph; j++, jc--)
    for (k = 0; k < l1; k++) {
      cc[(0) + ido * ((k) + l1 * (j))] = ch[(0) + ido * ((k) + l1 * (jc))] + ch[(0) + ido * ((k) + l1 * (j))];
      cc[(0) + ido * ((k) + l1 * (jc))] = ch[(0) + ido * ((k) + l1 * (jc))] - ch[(0) + ido * ((k) + l1 * (j))];
    }
  csarr = ((double *)util_malloc_((2 * ip) * sizeof(double)));
  arg = twopi / ip;
  csarr[0] = 1.;
  csarr[1] = 0.;
  csarr[2] = csarr[2 * ip - 2] = cos(arg);
  csarr[3] = sin(arg);
  csarr[2 * ip - 1] = -csarr[3];
  for (i = 2; i <= ip / 2; ++i) {
    csarr[2 * i] = csarr[2 * ip - 2 * i] = cos(i * arg);
    csarr[2 * i + 1] = sin(i * arg);
    csarr[2 * ip - 2 * i + 1] = -csarr[2 * i + 1];
  }
  for (l = 1, lc = ip - 1; l < ipph; l++, lc--) {
    ar1 = csarr[2 * l];
    ai1 = csarr[2 * l + 1];
    for (ik = 0; ik < idl1; ik++) {
      ch[(ik) + idl1 * (l)] = cc[(ik) + idl1 * (0)] + ar1 * cc[(ik) + idl1 * (1)];
      ch[(ik) + idl1 * (lc)] = ai1 * cc[(ik) + idl1 * (ip - 1)];
    }
    aidx = 2 * l;
    for (j = 2, jc = ip - 2; j < ipph; j++, jc--) {
      aidx += 2 * l;
      if (aidx >= 2 * ip)
        aidx -= 2 * ip;
      ar2 = csarr[aidx];
      ai2 = csarr[aidx + 1];
      for (ik = 0; ik < idl1; ik++) {
        ch[(ik) + idl1 * (l)] += ar2 * cc[(ik) + idl1 * (j)];
        ch[(ik) + idl1 * (lc)] += ai2 * cc[(ik) + idl1 * (jc)];
      }
    }
  }
  do {
    util_free_(csarr);
    (csarr) = NULL;
  } while (0);
  for (j = 1; j < ipph; j++)
    for (ik = 0; ik < idl1; ik++)
      ch[(ik) + idl1 * (0)] += cc[(ik) + idl1 * (j)];
  for (k = 0; k < l1; k++)
    memcpy(&cc[(0) + ido * ((0) + cdim * (k))], &ch[(0) + ido * ((k) + l1 * (0))], ido * sizeof(double));
  for (j = 1; j < ipph; j++) {
    jc = ip - j;
    j2 = 2 * j;
    for (k = 0; k < l1; k++) {
      cc[(ido - 1) + ido * ((j2 - 1) + cdim * (k))] = ch[(0) + ido * ((k) + l1 * (j))];
      cc[(0) + ido * ((j2) + cdim * (k))] = ch[(0) + ido * ((k) + l1 * (jc))];
    }
  }
  if (ido == 1)
    return;
  for (j = 1; j < ipph; j++) {
    jc = ip - j;
    j2 = 2 * j;
    for (k = 0; k < l1; k++)
      for (i = 2; i < ido; i += 2) {
        ic = ido - i;
        {
          cc[(i - 1) + ido * ((j2) + cdim * (k))] = ch[(i - 1) + ido * ((k) + l1 * (j))] + ch[(i - 1) + ido * ((k) + l1 * (jc))];
          cc[(ic - 1) + ido * ((j2 - 1) + cdim * (k))] = ch[(i - 1) + ido * ((k) + l1 * (j))] - ch[(i - 1) + ido * ((k) + l1 * (jc))];
        }
        {
          cc[(i) + ido * ((j2) + cdim * (k))] = ch[(i) + ido * ((k) + l1 * (jc))] + ch[(i) + ido * ((k) + l1 * (j))];
          cc[(ic) + ido * ((j2 - 1) + cdim * (k))] = ch[(i) + ido * ((k) + l1 * (jc))] - ch[(i) + ido * ((k) + l1 * (j))];
        }
      }
  }
}
```