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

```asm
<dradbg>:
  endbr64
  push   %r15
  mov    %rsi,%rax
  pxor   %xmm1,%xmm1
  mov    %rsi,%r15
  push   %r14
  cvtsi2sd %rax,%xmm1
  push   %r13
  imul   %rdi,%r15
  mov    %r9,%r13
  push   %r12
  mov    %rdi,%r12
  mov    %rdx,%rdi
  push   %rbp
  imul   %r12,%rdi
  push   %rbx
  mov    %r15,%rbp
  mov    %r8,%rbx
  mov    %rax,%r15
  sub    $0x168,%rsp
  movsd  D2(%rip),%xmm0
  mov    %rsi,0xc8(%rsp)
  lea    0x150(%rsp),%rsi
  mov    0x1a8(%rsp),%r14
  divsd  %xmm1,%xmm0
  mov    %rdi,0x78(%rsp)
  lea    0x158(%rsp),%rdi
  mov    %rdx,0x70(%rsp)
  mov    %rcx,0x130(%rsp)
  call   <sincos@plt>
  lea    -0x1(%r12),%rax
  mov    0x70(%rsp),%rdx
  movsd  0x150(%rsp),%xmm1
  mov    %rax,0x10(%rsp)
  sar    %rax
  movsd  0x158(%rsp),%xmm2
  mov    %rax,0x138(%rsp)
  lea    0x1(%r15),%rax
  mov    0x1b0(%rsp),%r11
  sar    %rax
  cmp    %rdx,%r12
  mov    %rax,0x38(%rsp)
  jl     L202
  test   %rdx,%rdx
  jle    L517
  test   %r12,%r12
  jle    L955
  mov    %r12,%rax
  lea    0x8(%rbx),%rcx
  mov    %r12,%rsi
  mov    %r13,0x18(%rsp)
  sub    %rbp,%rax
  mov    %rcx,0x8(%rsp)
  shr    %rsi
  mov    %r12,%r9
  shl    $0x3,%rax
  shl    $0x4,%rsi
  and    $0xfffffffffffffffe,%r9
  mov    %r14,%r8
  xor    %edi,%edi
  xor    %r10d,%r10d
  xor    %r15d,%r15d
  mov    %rax,%r11
L1155:
  mov    0x8(%rsp),%rax
  lea    0x0(,%rdi,8),%rdx
  lea    (%r14,%r10,8),%rcx
  add    %rdx,%rax
  cmp    %rax,%rcx
  je     L841
  cmpq   $0x2,0x10(%rsp)
  seta   %r13b
  add    %rbx,%rdx
  xor    %eax,%eax
  test   %r13b,%r13b
  je     L841
  nopl   0x0(%rax,%rax,1)
L626:
  movupd (%rdx,%rax,1),%xmm6
  movups %xmm6,(%rcx,%rax,1)
  add    $0x10,%rax
  cmp    %rsi,%rax
  jne    L626
  lea    (%r10,%r9,1),%rax
  lea    (%rdi,%r9,1),%rdx
  cmp    %r12,%r9
  je     L24
  movsd  (%rbx,%rdx,8),%xmm0
  movsd  %xmm0,(%r14,%rax,8)
L24:
  add    %r12,%r10
  add    %rbp,%rdi
  add    $0x1,%r15
  add    %r11,%r8
  cmp    %r15,0x70(%rsp)
  jne    L1155
  mov    0x18(%rsp),%r13
  mov    0x1b0(%rsp),%r11
L955:
  mov    0xc8(%rsp),%rax
  imul   0x78(%rsp),%rax
  cmpq   $0x1,0x38(%rsp)
  mov    %rax,0x98(%rsp)
  lea    (%r12,%r12,1),%rax
  mov    %rax,0x18(%rsp)
  jle    L11
L1068:
  cmp    $0x1,%r12
  jne    L1196
  mov    0x78(%rsp),%r8
  mov    %rbp,%rsi
  mov    0x98(%rsp),%rdx
  movq   $0x1,0x8(%rsp)
  movq   $0x0,0x30(%rsp)
  mov    %r8,%rax
  sub    %r8,%rdx
  lea    (%r14,%r8,8),%r15
  mov    %r12,0x90(%rsp)
  neg    %rax
  lea    0x0(,%rdx,8),%r10
  mov    %r15,0x88(%rsp)
  mov    %rbp,%r15
  lea    0x0(,%rax,8),%rcx
  sub    %r8,%rax
  mov    %r11,0x1b0(%rsp)
  mov    %rcx,0x50(%rsp)
  neg    %rcx
  mov    %rcx,%r9
  lea    0x8(%rbx),%rcx
  mov    %rax,0x48(%rsp)
  mov    %rcx,0x10(%rsp)
  mov    %rdx,%rcx
  sub    %r8,%rcx
  mov    %r13,0xa0(%rsp)
  mov    %rcx,0x18(%rsp)
  lea    -0x1(%r8),%rcx
  mov    0x18(%rsp),%r13
  mov    %rcx,0x40(%rsp)
  imul   %rbp,%rcx
  shl    $0x3,%rsi
  mov    %rsi,%rdx
  mov    $0x8,%esi
  lea    0x8(,%rcx,8),%rax
  cmovns %rax,%rsi
  sub    $0x8,%rax
  xor    %ecx,%ecx
  test   %rdx,%rdx
  cmovs  %rax,%rcx
  mov    %r8,%rax
  and    $0xfffffffffffffffe,%r8
  shr    %rax
  imul   %r8,%rbp
  sub    $0x8,%rsi
  mov    %r8,0x68(%rsp)
  mov    %rax,%rdi
  mov    0x98(%rsp),%rax
  mov    %rsi,0x58(%rsp)
  lea    -0x8(%rcx),%rsi
  mov    %rsi,0x60(%rsp)
  shl    $0x4,%r15
  shl    $0x4,%rdi
  mov    %rax,0x28(%rsp)
  mov    0x88(%rsp),%rax
  mov    %rbp,0x80(%rsp)
  mov    $0x2,%ebp
  sub    $0x8,%rax
  mov    %rax,0x20(%rsp)
L120:
  mov    0x58(%rsp),%rax
  mov    0x78(%rsp),%rcx
  lea    0x8(%r9),%r11
  lea    (%r14,%r10,1),%rsi
  add    %rcx,0x30(%rsp)
  sub    %rcx,0x28(%rsp)
  lea    0x0(,%rbp,8),%rcx
  lea    (%rax,%rcx,1),%r8
  add    0x60(%rsp),%rcx
  lea    (%r14,%r9,1),%rax
  add    %rbx,%r8
  add    %rbx,%rcx
  cmp    %r10,%r11
  setne  %r12b
  cmpq   $0x5,0x40(%rsp)
  seta   %r11b
  and    %r11d,%r12d
  mov    0x20(%rsp),%r11
  add    %r9,%r11
  cmp    %r11,%rcx
  seta   %r11b
  cmp    %rax,%r8
  setb   0x18(%rsp)
  or     0x18(%rsp),%r11b
  test   %r11b,%r12b
  je     L674
  mov    0x20(%rsp),%r11
  add    %r10,%r11
  cmp    %r11,%rcx
  seta   %cl
  cmp    %rsi,%r8
  setb   %r8b
  or     %cl,%r8b
  je     L674
  mov    0x10(%rsp),%rcx
  xor    %r8d,%r8d
  nopl   0x0(%rax)
L1161:
  movsd  (%rcx),%xmm0
  movhpd (%rcx,%rdx,1),%xmm0
  addpd  %xmm0,%xmm0
  movups %xmm0,(%rax,%r8,1)
  movsd  0x8(%rcx),%xmm0
  movhpd 0x8(%rdx,%rcx,1),%xmm0
  add    %r15,%rcx
  addpd  %xmm0,%xmm0
  movups %xmm0,(%rsi,%r8,1)
  add    $0x10,%r8
  cmp    %rdi,%r8
  jne    L1161
  mov    0x68(%rsp),%r11
  mov    0x30(%rsp),%rsi
  mov    0x28(%rsp),%rcx
  mov    0x80(%rsp),%rax
  mov    0x78(%rsp),%r8
  add    %r11,%rsi
  add    %r11,%rcx
  add    %rbp,%rax
  cmp    %r8,%r11
  je     L817
  movsd  -0x8(%rbx,%rax,8),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,(%r14,%rsi,8)
  movsd  (%rbx,%rax,8),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,(%r14,%rcx,8)
L817:
  mov    0x50(%rsp),%rcx
  addq   $0x1,0x8(%rsp)
  add    $0x2,%rbp
  addq   $0x10,0x10(%rsp)
  mov    0x8(%rsp),%rax
  add    %rcx,%r10
  sub    %rcx,%r9
  add    0x48(%rsp),%r13
  cmp    %rax,0x38(%rsp)
  jne    L120
  mov    0x90(%rsp),%r12
  mov    0xa0(%rsp),%r13
  mov    0x1b0(%rsp),%r11
L1087:
  mov    0x130(%rsp),%rax
  pxor   %xmm5,%xmm5
  mov    0xc8(%rsp),%r15
  movq   $0x0,0xa8(%rsp)
  mov    0x1a0(%rsp),%rdx
  movsd  D0(%rip),%xmm6
  movq   $0x1,0xa0(%rsp)
  imul   %rax,%r15
  mov    %rax,%rdi
  neg    %rdi
  lea    0x0(,%rdi,8),%rbp
  mov    %rdi,%rsi
  mov    %rax,%rdi
  mov    %r15,%r8
  mov    %r15,0xb0(%rsp)
  mov    %rdi,%rbx
  mov    %rbp,%rcx
  sub    %rax,%r8
  lea    (%rax,%rax,1),%rax
  neg    %rcx
  lea    0x10(,%r8,8),%r15
  sub    %r8,%rbx
  mov    %rcx,0x48(%rsp)
  lea    0x0(,%r8,8),%r10
  mov    %r15,0x58(%rsp)
  add    %r10,%rdx
  mov    0x1a0(%rsp),%rcx
  mov    %rax,0xc0(%rsp)
  mov    $0x10,%eax
  sub    %rbp,%rax
  mov    %r10,0x40(%rsp)
  sub    %rbp,%rcx
  mov    %rax,0x60(%rsp)
  mov    %rax,%r9
  mov    %r11,%rax
  add    0x1a0(%rsp),%r9
  sub    %rbp,%rax
  mov    %r9,0x50(%rsp)
  lea    0x8(%rax),%r15
  mov    %rax,0xe8(%rsp)
  mov    %rax,%r9
  mov    %rdi,%rax
  mov    %r15,0xd0(%rsp)
  lea    (%r11,%r10,1),%r15
  shr    %rax
  sub    %rbp,%r9
  mov    %r15,0xd8(%rsp)
  lea    0x10(%r11,%r10,1),%r15
  mov    %r15,0x128(%rsp)
  lea    0x8(%r11,%r10,1),%r15
  mov    %r15,0xe0(%rsp)
  mov    %rax,%r15
  mov    %rdi,%rax
  and    $0xfffffffffffffffe,%rdi
  lea    (%r11,%rdi,8),%r10
  mov    %rdi,0x8(%rsp)
  shl    $0x4,%r15
  mov    %r10,0x80(%rsp)
  mov    %rax,%r10
  add    %rdi,%rax
  lea    (%r11,%rax,8),%rax
  mov    %r8,0x120(%rsp)
  sub    %r10,%rsi
  mov    %rax,0x88(%rsp)
  mov    %r8,%rax
  sub    %r10,%r8
  add    %rdi,%rax
  mov    %r9,0x108(%rsp)
  lea    (%r11,%rax,8),%rdi
  mov    %r10,%rax
  shl    $0x4,%rax
  mov    %rdi,0x90(%rsp)
  lea    0x8(%r11,%rax,1),%rdi
  mov    %r8,%rax
  sub    0xc0(%rsp),%rax
  mov    %rax,0x110(%rsp)
  lea    0x8(%r11),%rax
  mov    %rax,0xf0(%rsp)
  lea    -0x1(%r10),%rax
  mov    %rax,0xf8(%rsp)
  mov    $0xfffffffffffffff8,%rax
  mov    %rdi,0x100(%rsp)
  sub    %rbp,%rax
  lea    (%r11,%r8,8),%rdi
  mov    %rdi,0x118(%rsp)
  mov    %rsi,0x68(%rsp)
  mov    %rax,0xb8(%rsp)
  mov    %r12,0x140(%rsp)
  mov    %r10,%r12
  mov    %r13,0x148(%rsp)
  mov    %rbx,%r13
  mov    %r14,0x1a8(%rsp)
  mov    %rbp,%r14
  nopw   0x0(%rax,%rax,1)
L1185:
  movapd %xmm6,%xmm0
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm6
  add    %r12,0xa8(%rsp)
  mulsd  %xmm0,%xmm6
  sub    %r12,0xb0(%rsp)
  mulsd  %xmm5,%xmm3
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm5
  subsd  %xmm3,%xmm6
  addsd  %xmm0,%xmm5
  test   %r12,%r12
  jle    L688
  mov    0x48(%rsp),%rax
  mov    0xd0(%rsp),%rdi
  mov    0xf0(%rsp),%rbx
  add    $0x8,%rax
  cmp    0x40(%rsp),%rax
  setne  %sil
  cmp    %rdi,%rcx
  setne  %al
  and    %esi,%eax
  cmpq   $0x1,0xf8(%rsp)
  seta   %sil
  and    %eax,%esi
  cmp    %rbx,%rcx
  setne  %al
  and    %esi,%eax
  cmp    %rdx,%rdi
  setne  %sil
  and    %eax,%esi
  cmp    0xe0(%rsp),%rdx
  setne  %al
  and    %esi,%eax
  cmp    %rdx,%rbx
  setne  %sil
  test   %sil,%al
  je     L295
  mov    0x50(%rsp),%rax
  mov    0xd8(%rsp),%rdi
  cmp    %rax,%rdi
  setae  %sil
  cmp    0x128(%rsp),%rcx
  setae  %al
  or     %sil,%al
  je     L295
  movapd %xmm6,%xmm4
  movapd %xmm5,%xmm3
  mov    %rdi,%rsi
  xor    %eax,%eax
  mov    0xe8(%rsp),%rdi
  unpcklpd %xmm4,%xmm4
  unpcklpd %xmm3,%xmm3
  nopl   0x0(%rax)
L189:
  movupd (%rdi,%rax,1),%xmm0
  movupd (%r11,%rax,1),%xmm7
  mulpd  %xmm4,%xmm0
  addpd  %xmm7,%xmm0
  movups %xmm0,(%rcx,%rax,1)
  movupd (%rsi,%rax,1),%xmm0
  mulpd  %xmm3,%xmm0
  movups %xmm0,(%rdx,%rax,1)
  add    $0x10,%rax
  cmp    %r15,%rax
  jne    L189
  mov    0x8(%rsp),%rdi
  mov    0xa8(%rsp),%rsi
  mov    0xb0(%rsp),%rax
  add    %rdi,%rsi
  add    %rdi,%rax
  cmp    %rdi,%r12
  je     L1388
  mov    0x88(%rsp),%rdi
  movsd  (%rdi),%xmm0
  mov    0x80(%rsp),%rdi
  mulsd  %xmm6,%xmm0
  addsd  (%rdi),%xmm0
  mov    0x1a0(%rsp),%rdi
  movsd  %xmm0,(%rdi,%rsi,8)
  mov    0x90(%rsp),%rdi
  movsd  (%rdi),%xmm0
  mov    0x1a0(%rsp),%rdi
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,(%rdi,%rax,8)
L1388:
  cmpq   $0x2,0x38(%rsp)
  jle    L339
  mov    0x8(%rsp),%rbx
  mov    %r12,0x20(%rsp)
  movapd %xmm6,%xmm4
  movapd %xmm5,%xmm3
  mov    0xa8(%rsp),%rax
  mov    0x1a0(%rsp),%rdi
  movq   $0x2,0x10(%rsp)
  mov    %r13,0x18(%rsp)
  mov    0x110(%rsp),%rbp
  add    %rbx,%rax
  mov    0x118(%rsp),%r8
  mov    0x108(%rsp),%r10
  lea    (%rdi,%rax,8),%rdi
  mov    0xb0(%rsp),%rax
  add    %rbx,%rax
  mov    0x1a0(%rsp),%rbx
  lea    (%rbx,%rax,8),%rsi
  mov    0x120(%rsp),%rax
  mov    0x100(%rsp),%rbx
  mov    %rax,0x28(%rsp)
  nopw   0x0(%rax,%rax,1)
L1270:
  mov    0x40(%rsp),%r9
  add    %r12,0x20(%rsp)
  movapd %xmm4,%xmm0
  movapd %xmm5,%xmm7
  sub    %r12,0x28(%rsp)
  cmp    %r9,0x60(%rsp)
  movapd %xmm6,%xmm4
  mov    0x48(%rsp),%r13
  mulsd  %xmm0,%xmm4
  setle  %r9b
  cmp    %r13,0x58(%rsp)
  mulsd  %xmm3,%xmm7
  setle  %al
  mulsd  %xmm6,%xmm3
  or     %r9d,%eax
  cmp    $0x1,%r12
  mulsd  %xmm5,%xmm0
  setne  %r9b
  and    %eax,%r9d
  cmp    %rcx,%rbx
  setne  %al
  and    %r9d,%eax
  cmp    %rdx,%rbx
  subsd  %xmm7,%xmm4
  setne  %r9b
  and    %r9d,%eax
  lea    0x8(%r8),%r9
  addsd  %xmm0,%xmm3
  cmp    %r9,%rdx
  setne  %r9b
  test   %r9b,%al
  je     L773
  lea    0x10(%r8),%rax
  movapd %xmm4,%xmm8
  movapd %xmm3,%xmm7
  cmp    %rax,%rcx
  unpcklpd %xmm8,%xmm8
  unpcklpd %xmm7,%xmm7
  setae  0x30(%rsp)
  cmp    0x50(%rsp),%r8
  movzbl 0x30(%rsp),%r13d
  setae  %r9b
  xor    %eax,%eax
  or     %r13b,%r9b
  mov    0x18(%rsp),%r13
  je     L773
  nopl   0x0(%rax,%rax,1)
L1258:
  movupd (%r10,%rax,1),%xmm0
  movupd (%rcx,%rax,1),%xmm9
  mulpd  %xmm8,%xmm0
  addpd  %xmm9,%xmm0
  movups %xmm0,(%rcx,%rax,1)
  movupd (%r8,%rax,1),%xmm0
  movupd (%rdx,%rax,1),%xmm10
  mulpd  %xmm7,%xmm0
  addpd  %xmm10,%xmm0
  movups %xmm0,(%rdx,%rax,1)
  add    $0x10,%rax
  cmp    %r15,%rax
  jne    L1258
  mov    0x8(%rsp),%rax
  mov    0x20(%rsp),%r9
  mov    %r13,0x18(%rsp)
  mov    0x28(%rsp),%r13
  add    %rax,%r9
  add    %rax,%r13
  cmp    0x8(%rsp),%r12
  je     L564
  movsd  (%r11,%r9,8),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  (%rdi),%xmm0
  movsd  %xmm0,(%rdi)
  movsd  (%r11,%r13,8),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  (%rsi),%xmm0
  movsd  %xmm0,(%rsi)
L564:
  addq   $0x1,0x10(%rsp)
  sub    %r14,%rbx
  sub    %r14,%r10
  add    %r14,%r8
  mov    0x10(%rsp),%rax
  add    0x68(%rsp),%rbp
  cmp    %rax,0x38(%rsp)
  jne    L1270
  mov    0x18(%rsp),%r13
L584:
  addq   $0x1,0xa0(%rsp)
  add    %r14,%rdx
  sub    %r14,%rcx
  mov    0xa0(%rsp),%rax
  sub    %r14,0x60(%rsp)
  add    0xc0(%rsp),%r13
  sub    %r14,0x48(%rsp)
  add    %r14,0x58(%rsp)
  sub    %r14,0x50(%rsp)
  add    %r14,0x40(%rsp)
  cmp    %rax,0x38(%rsp)
  jg     L1185
L339:
  mov    0x130(%rsp),%rbx
  mov    %r14,%rbp
  mov    0x140(%rsp),%r12
  mov    0x148(%rsp),%r13
  mov    0x1a8(%rsp),%r14
  test   %rbx,%rbx
  jle    L838
  mov    $0x8,%eax
  mov    %r12,0x10(%rsp)
  mov    %r11,%rdi
  xor    %ecx,%ecx
  mov    0x8(%rsp),%r12
  sub    %r11,%rax
  mov    %r13,0x8(%rsp)
  sub    %rbp,%rdi
  mov    0xe8(%rsp),%rdx
  mov    0x80(%rsp),%r10
  mov    $0x1,%esi
  mov    %rax,%r8
  mov    0x38(%rsp),%r13
  lea    -0x1(%rbx),%r9
  nopl   0x0(%rax)
L1042:
  mov    %r8,%rax
  add    %rbx,%rcx
  add    %rdx,%rax
  je     L1273
  xor    %eax,%eax
  cmp    $0x1,%r9
  jbe    L1273
  nopl   0x0(%rax,%rax,1)
L423:
  movupd (%rdx,%rax,1),%xmm0
  movupd (%r11,%rax,1),%xmm6
  addpd  %xmm6,%xmm0
  movups %xmm0,(%r11,%rax,1)
  add    $0x10,%rax
  cmp    %r15,%rax
  jne    L423
  lea    (%rcx,%r12,1),%rax
  cmp    %r12,%rbx
  je     L609
  movsd  (%r10),%xmm0
  addsd  (%r11,%rax,8),%xmm0
  movsd  %xmm0,(%r10)
L609:
  add    $0x1,%rsi
  sub    %rbp,%rdx
  cmp    %rsi,%r13
  jg     L1042
  mov    0x10(%rsp),%r12
  mov    0x8(%rsp),%r13
L838:
  cmpq   $0x0,0x70(%rsp)
  jle    L689
  cmp    $0x1,%r12
  jne    L759
  mov    0x78(%rsp),%rdi
  mov    0x98(%rsp),%r10
  xor    %ebx,%ebx
  mov    %rdi,%rax
  sub    %rdi,%r10
  mov    %rdi,%r11
  mov    %rdi,%r15
  shl    $0x3,%r10
  neg    %rax
  shr    %r11
  and    $0xfffffffffffffffe,%r15
  lea    0x0(,%rdi,8),%rbp
  shl    $0x3,%rax
  shl    $0x4,%r11
  mov    %rax,0x8(%rsp)
  lea    (%r14,%rbp,1),%r9
  lea    0x0(%r13,%rbp,1),%rdx
  mov    %rbp,0x10(%rsp)
  lea    0x0(%r13,%r10,1),%rax
  mov    %rdi,%rbp
  add    %r14,%r10
L621:
  lea    0x8(%rdx),%rcx
  lea    0x8(%rax),%rsi
  sub    %rbp,0x98(%rsp)
  add    %rbp,%rbx
  cmp    %rcx,%r10
  lea    0x10(%r9),%rdi
  setne  %cl
  cmp    %rsi,%r10
  setne  %sil
  and    %ecx,%esi
  cmp    $0x1,%rbp
  setne  %cl
  and    %ecx,%esi
  lea    0x10(%rdx),%rcx
  cmp    %rcx,%r9
  setae  %cl
  cmp    %rdi,%rdx
  setae  %r8b
  or     %r8d,%ecx
  and    %esi,%ecx
  lea    0x10(%rax),%rsi
  cmp    %rsi,%r9
  setae  %sil
  cmp    %rdi,%rax
  setae  %dil
  or     %edi,%esi
  test   %sil,%cl
  je     L1002
  mov    %rax,%rsi
  lea    0x8(,%rbx,8),%rcx
  sub    %r13,%rsi
  cmp    %rsi,%rcx
  je     L1002
  xor    %ecx,%ecx
  nopl   0x0(%rax)
L1209:
  movupd (%rax,%rcx,1),%xmm6
  movupd (%rdx,%rcx,1),%xmm0
  subpd  %xmm6,%xmm0
  movups %xmm0,(%r9,%rcx,1)
  movupd (%rdx,%rcx,1),%xmm0
  movupd (%rax,%rcx,1),%xmm6
  addpd  %xmm6,%xmm0
  movups %xmm0,(%r10,%rcx,1)
  add    $0x10,%rcx
  cmp    %r11,%rcx
  jne    L1209
  mov    0x98(%rsp),%rdi
  lea    (%rbx,%r15,1),%r8
  add    %r15,%rdi
  cmp    %rbp,%r15
  je     L861
  lea    0x0(%r13,%r8,8),%rsi
  lea    0x0(%r13,%rdi,8),%rcx
  movsd  (%rsi),%xmm0
  subsd  (%rcx),%xmm0
  movsd  %xmm0,(%r14,%r8,8)
  movsd  (%rsi),%xmm0
  addsd  (%rcx),%xmm0
  movsd  %xmm0,(%r14,%rdi,8)
L861:
  mov    0x10(%rsp),%rdi
  add    $0x1,%r12
  add    %rdi,%r9
  add    %rdi,%rdx
  mov    0x8(%rsp),%rdi
  add    %rdi,%rax
  add    %rdi,%r10
  cmp    %r12,0x38(%rsp)
  jg     L621
L1327:
  add    $0x168,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L202:
  test   %r12,%r12
  jle    L986
  cmpq   $0x0,0x70(%rsp)
  jle    L517
  cmp    $0x1,%r12
  jne    L402
  mov    0x78(%rsp),%rdi
  lea    0x0(,%rbp,8),%rsi
  mov    %rbx,%rdx
  mov    %r14,%rax
  lea    (%r14,%rdi,8),%rcx
  xchg   %ax,%ax
L140:
  movsd  (%rdx),%xmm0
  add    $0x8,%rax
  add    %rsi,%rdx
  movsd  %xmm0,-0x8(%rax)
  cmp    %rcx,%rax
  jne    L140
  jmp    L955
  nopw   0x0(%rax,%rax,1)
L773:
  mov    0xb8(%rsp),%r9
  mov    0x40(%rsp),%rax
  mov    %rbx,0x30(%rsp)
  add    0x1a0(%rsp),%rax
  lea    (%r9,%rbx,1),%r13
  mov    %r10,%r9
  mov    %r13,%rbx
  mov    0x18(%rsp),%r13
  nopl   0x0(%rax)
L1116:
  movsd  (%r9),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  (%rax,%r13,8),%xmm0
  movsd  %xmm0,(%rax,%r13,8)
  movsd  (%r9,%rbp,8),%xmm0
  add    $0x8,%r9
  add    $0x8,%rax
  mulsd  %xmm3,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %r9,%rbx
  jne    L1116
  mov    %r13,0x18(%rsp)
  mov    0x30(%rsp),%rbx
  jmp    L564
L295:
  mov    0x1a0(%rsp),%rax
  mov    0xa8(%rsp),%rdi
  mov    0x1a0(%rsp),%rbx
  mov    0x40(%rsp),%rsi
  mov    0xd8(%rsp),%r8
  mov    0xe8(%rsp),%r9
  lea    (%rax,%rdi,8),%rdi
  xor    %eax,%eax
  add    %rbx,%rsi
  nopl   0x0(%rax,%rax,1)
L197:
  add    $0x1,%rax
  movsd  -0x8(%r9,%rax,8),%xmm0
  mulsd  %xmm6,%xmm0
  addsd  -0x8(%r11,%rax,8),%xmm0
  movsd  %xmm0,-0x8(%rdi,%rax,8)
  movsd  -0x8(%r8,%rax,8),%xmm0
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,-0x8(%rsi,%rax,8)
  cmp    %rax,%r12
  jne    L197
  jmp    L1388
L674:
  mov    0x88(%rsp),%rcx
  lea    (%r9,%rcx,1),%rsi
  mov    0x10(%rsp),%rcx
  nopl   0x0(%rax,%rax,1)
L1139:
  movsd  (%rcx),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x8(%rcx),%xmm0
  add    %rdx,%rcx
  addsd  %xmm0,%xmm0
  movsd  %xmm0,(%rax,%r13,8)
  add    $0x8,%rax
  cmp    %rsi,%rax
  jne    L1139
  jmp    L817
L1273:
  mov    %r11,%rax
  nop
L409:
  movsd  (%rax),%xmm0
  addsd  (%rax,%rcx,8),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rdi
  jne    L409
  jmp    L609
L1002:
  mov    0x98(%rsp),%rsi
  mov    %rbx,%rcx
  lea    (%rbx,%rbp,1),%rdi
  nopl   0x0(%rax,%rax,1)
L982:
  movsd  0x0(%r13,%rcx,8),%xmm0
  subsd  0x0(%r13,%rsi,8),%xmm0
  movsd  %xmm0,(%r14,%rcx,8)
  movsd  0x0(%r13,%rcx,8),%xmm0
  addsd  0x0(%r13,%rsi,8),%xmm0
  add    $0x1,%rcx
  movsd  %xmm0,(%r14,%rsi,8)
  add    $0x1,%rsi
  cmp    %rdi,%rcx
  jne    L982
  jmp    L861
L688:
  cmpq   $0x2,0x38(%rsp)
  jg     L584
  jmp    L339
L841:
  mov    %rdi,%rax
  lea    (%r12,%rdi,1),%rdx
  nopw   0x0(%rax,%rax,1)
L13:
  movsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,(%r8,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L13
  jmp    L24
L689:
  cmp    $0x1,%r12
  je     L1327
L1291:
  mov    0x138(%rsp),%rdi
  cmp    %rdi,0x70(%rsp)
  jg     L984
  cmpq   $0x1,0x38(%rsp)
  jle    L175
  cmpq   $0x0,0x70(%rsp)
  jle    L175
  cmp    $0x2,%r12
  jle    L175
  mov    0x78(%rsp),%rdi
  mov    0x98(%rsp),%rbx
  movq   $0x0,0x88(%rsp)
  lea    0x0(,%r12,8),%r15
  movq   $0x1,0x90(%rsp)
  mov    %rdi,%rax
  mov    %r11,0x1b0(%rsp)
  neg    %rax
  mov    %r12,0x40(%rsp)
  sub    %rdi,%rax
  mov    %rax,0xb0(%rsp)
  mov    %rbx,%rax
  sub    %rdi,%rax
  sub    %rdi,%rax
  mov    %rax,%rcx
  mov    %rdi,%rax
  lea    0x0(,%rdi,8),%rdi
  lea    0x8(,%rax,8),%rax
  mov    %rdi,0xb8(%rsp)
  lea    -0x3(%r12),%rdi
  mov    %rcx,%r12
  mov    %rax,0xa0(%rsp)
  mov    %rdi,%rax
  shr    $0x2,%rax
  mov    %rdi,0x58(%rsp)
  and    $0xfffffffffffffffe,%rdi
  mov    %rax,%r10
  lea    0x0(,%rax,4),%rsi
  lea    0x4(,%rax,4),%rax
  mov    %rdi,0xa8(%rsp)
  mov    %rax,0x68(%rsp)
  lea    0x20(%r14),%rax
  shl    $0x5,%r10
  mov    %rax,0x48(%rsp)
  lea    0x8(%r13),%rax
  mov    %r10,%rbp
  mov    %rax,0x50(%rsp)
  lea    0x20(%r13),%rax
  mov    %rax,0x60(%rsp)
  lea    0x8(,%rbx,8),%rax
  mov    %rsi,0x80(%rsp)
  mov    %rax,0xc0(%rsp)
L1122:
  mov    0x78(%rsp),%rax
  add    %rax,0x88(%rsp)
  sub    %rax,0x98(%rsp)
  mov    0x88(%rsp),%rbx
  mov    0x98(%rsp),%rcx
  mov    0xa0(%rsp),%r8
  movq   $0x0,0x8(%rsp)
  mov    0xc0(%rsp),%rdi
  mov    %rbx,%rdx
  mov    %rbx,0x18(%rsp)
  mov    %rcx,0x30(%rsp)
  neg    %rdx
  lea    (%r14,%r8,1),%rax
  lea    (%rdi,%rdx,8),%rdi
  nopw   0x0(%rax,%rax,1)
L1129:
  mov    0x48(%rsp),%r10
  mov    0x60(%rsp),%rsi
  lea    (%r8,%r13,1),%rdx
  lea    (%rdi,%r13,1),%rcx
  lea    (%r10,%rdi,1),%r9
  lea    (%r10,%r8,1),%rbx
  mov    %r9,0x20(%rsp)
  lea    0x20(%rdi),%r9
  lea    (%rsi,%r8,1),%r11
  add    %rdi,%rsi
  lea    0x20(%r8),%r10
  cmp    %r9,%r8
  mov    %rsi,0x10(%rsp)
  lea    (%r14,%rdi,1),%rsi
  setge  %r9b
  cmp    %rdi,%r10
  setle  %r10b
  or     %r9d,%r10d
  cmpq   $0x3,0x58(%rsp)
  seta   %r9b
  and    %r9d,%r10d
  cmp    %r11,%rax
  setae  %r9b
  cmp    %rdx,%rbx
  setbe  0x28(%rsp)
  or     0x28(%rsp),%r9b
  and    %r10d,%r9d
  cmp    0x10(%rsp),%rax
  setae  %r10b
  cmp    %rcx,%rbx
  setbe  %bl
  or     %ebx,%r10d
  mov    0x20(%rsp),%rbx
  and    %r9d,%r10d
  cmp    %rdx,%rbx
  setbe  %r9b
  cmp    %rsi,%r11
  setbe  %r11b
  or     %r11d,%r9d
  test   %r9b,%r10b
  je     L1251
  cmp    %rsi,0x10(%rsp)
  mov    0x50(%rsp),%r10
  setbe  %r11b
  cmp    %rcx,%rbx
  setbe  %bl
  add    %r8,%r10
  xor    %r9d,%r9d
  or     %r11b,%bl
  je     L1251
  nopl   0x0(%rax,%rax,1)
L1188:
  movupd (%rdx,%r9,1),%xmm0
  movupd 0x10(%rcx,%r9,1),%xmm1
  movlpd 0x8(%rcx,%r9,1),%xmm1
  movupd (%rcx,%r9,1),%xmm3
  movhpd 0x10(%rdx,%r9,1),%xmm0
  movapd %xmm0,%xmm2
  addpd  %xmm1,%xmm0
  movhpd 0x10(%rcx,%r9,1),%xmm3
  subpd  %xmm1,%xmm2
  movupd 0x10(%rdx,%r9,1),%xmm1
  movlpd 0x8(%rdx,%r9,1),%xmm1
  addpd  %xmm3,%xmm1
  movapd %xmm2,%xmm4
  unpcklpd %xmm1,%xmm4
  unpckhpd %xmm1,%xmm2
  movups %xmm2,0x10(%rax,%r9,1)
  movapd %xmm0,%xmm2
  movups %xmm4,(%rax,%r9,1)
  movupd (%r10,%r9,1),%xmm1
  movhpd 0x10(%r10,%r9,1),%xmm1
  subpd  %xmm3,%xmm1
  unpcklpd %xmm1,%xmm2
  unpckhpd %xmm1,%xmm0
  movups %xmm2,(%rsi,%r9,1)
  movups %xmm0,0x10(%rsi,%r9,1)
  add    $0x20,%r9
  cmp    %r9,%rbp
  jne    L1188
  mov    0x18(%rsp),%rcx
  mov    0x80(%rsp),%rbx
  mov    0x30(%rsp),%rsi
  lea    0x2(%rcx,%rbx,1),%r11
  lea    0x2(%rsi,%rbx,1),%rsi
  lea    0x0(,%r11,8),%rcx
  lea    0x0(,%rsi,8),%rdx
  lea    -0x8(%r13,%rcx,1),%r10
  lea    0x0(%r13,%rdx,1),%r9
  movsd  (%r10),%xmm0
  subsd  (%r9),%xmm0
  movsd  %xmm0,-0x8(%r14,%rcx,1)
  movsd  (%r10),%xmm0
  addsd  (%r9),%xmm0
  lea    0x0(%r13,%rcx,1),%r10
  lea    -0x8(%r13,%rdx,1),%r9
  movsd  %xmm0,-0x8(%r14,%rdx,1)
  movsd  (%r10),%xmm0
  addsd  (%r9),%xmm0
  movsd  %xmm0,(%r14,%r11,8)
  movsd  (%r10),%xmm0
  subsd  (%r9),%xmm0
  movsd  %xmm0,(%r14,%rsi,8)
  mov    0x40(%rsp),%rsi
  cmp    %rsi,0x68(%rsp)
  jge    L1356
  lea    0x8(%r13,%rcx,1),%r9
  lea    0x10(%r13,%rdx,1),%rsi
  movsd  (%r9),%xmm0
  subsd  (%rsi),%xmm0
  movsd  %xmm0,0x8(%r14,%rcx,1)
  movsd  (%r9),%xmm0
  addsd  (%rsi),%xmm0
  lea    0x10(%r13,%rcx,1),%r9
  lea    0x8(%r13,%rdx,1),%rsi
  movsd  %xmm0,0x8(%r14,%rdx,1)
  movsd  (%r9),%xmm0
  addsd  (%rsi),%xmm0
  movsd  %xmm0,0x10(%r14,%rcx,1)
  movsd  (%r9),%xmm0
  subsd  (%rsi),%xmm0
  movsd  %xmm0,0x10(%r14,%rdx,1)
L1356:
  mov    0x40(%rsp),%rsi
  add    %r15,%rax
  add    %rsi,0x18(%rsp)
  add    %r15,%rdi
  addq   $0x1,0x8(%rsp)
  add    %r15,%r8
  mov    0x8(%rsp),%rbx
  add    %rsi,0x30(%rsp)
  cmp    %rbx,0x70(%rsp)
  jne    L1129
  mov    0xb8(%rsp),%rbx
  add    0xb0(%rsp),%r12
  addq   $0x1,0x90(%rsp)
  mov    0x90(%rsp),%rax
  add    %rbx,0xa0(%rsp)
  cmp    %rax,0x38(%rsp)
  jne    L1122
  mov    0x1b0(%rsp),%r11
  mov    %rsi,%r12
L175:
  cmpq   $0x0,0x130(%rsp)
  jle    L0
  lea    0x8(%r11),%rax
  cmp    %rax,0x1a0(%rsp)
  je     L498
  mov    0x130(%rsp),%rdi
  lea    -0x1(%rdi),%rax
  cmp    $0x2,%rax
  jbe    L498
  shr    %rdi
  mov    0x1a0(%rsp),%rcx
  xor    %eax,%eax
  mov    %rdi,%rdx
  shl    $0x4,%rdx
L363:
  movupd (%r11,%rax,1),%xmm6
  movups %xmm6,(%rcx,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%rdx
  jne    L363
  mov    0x130(%rsp),%rdi
  mov    %rdi,%rax
  and    $0xfffffffffffffffe,%rax
  and    $0x1,%edi
  je     L0
  movsd  (%r11,%rax,8),%xmm0
  mov    0x1a0(%rsp),%rdi
  movsd  %xmm0,(%rdi,%rax,8)
L0:
  cmpq   $0x1,0xc8(%rsp)
  jle    L1327
  mov    0x70(%rsp),%rdi
  test   %rdi,%rdi
  jle    L1327
  mov    0x78(%rsp),%r8
  mov    0xc8(%rsp),%r9
  xor    %esi,%esi
  mov    $0x1,%ecx
L1381:
  add    %r8,%rsi
  xor    %edx,%edx
  mov    %rsi,%rax
  nopl   0x0(%rax,%rax,1)
L101:
  movsd  (%r14,%rax,8),%xmm0
  add    $0x1,%rdx
  movsd  %xmm0,0x0(%r13,%rax,8)
  add    %r12,%rax
  cmp    %rdx,%rdi
  jne    L101
  add    $0x1,%rcx
  cmp    %rcx,%r9
  jne    L1381
  mov    0x138(%rsp),%rdi
  cmp    %rdi,0x70(%rsp)
  jl     L159
  mov    %r12,%rax
  not    %rax
  mov    %rax,0x50(%rsp)
  cmp    $0x2,%r12
  jle    L1327
  mov    0x78(%rsp),%rdi
  lea    0x0(,%r12,8),%rdx
  mov    0x70(%rsp),%rbx
  movq   $0x0,0x58(%rsp)
  lea    0xf(%rdx),%rcx
  mov    0x1b8(%rsp),%r15
  movq   $0x1,0x60(%rsp)
  mov    %rdi,%rax
  mov    %r12,0x98(%rsp)
  sub    %r12,%rax
  shl    $0x3,%rax
  mov    %rax,0x88(%rsp)
  lea    0x0(,%rdi,8),%rax
  mov    %rax,0x80(%rsp)
  sub    0x1b8(%rsp),%rax
  add    $0x8,%rax
  cmp    $0x1e,%rcx
  mov    %rax,0x8(%rsp)
  lea    0x8(,%rdi,8),%rax
  seta   %sil
  mov    %rax,0x68(%rsp)
  lea    -0x1(%rbx),%rax
  cmp    $0xa,%rax
  seta   %cl
  imul   %r12,%rax
  and    %ecx,%esi
  xor    %ecx,%ecx
  mov    %sil,0x30(%rsp)
  shl    $0x3,%rax
  test   %rdx,%rdx
  cmovs  %rax,%rcx
  add    $0x8,%rax
  test   %rdx,%rdx
  mov    %rcx,0x28(%rsp)
  mov    $0x8,%ecx
  cmovns %rax,%rcx
  mov    $0x0,%eax
  cmovns %rdx,%rax
  mov    %rcx,0x20(%rsp)
  mov    %r12,%rcx
  mov    %rax,%rsi
  mov    $0x0,%eax
  cmovle %rdx,%rax
  add    $0x8,%r15
  add    $0x8,%rsi
  shl    $0x4,%rcx
  mov    %r15,0x90(%rsp)
  mov    %rax,0x18(%rsp)
  mov    %rbx,%rax
  shr    %rbx
  mov    %rbx,%rdi
  mov    %rax,%rbx
  and    $0xfffffffffffffffe,%rax
  mov    %rsi,0x38(%rsp)
  mov    %rax,0x48(%rsp)
  imul   %r12,%rax
  mov    %rax,0x40(%rsp)
  lea    -0x3(%r12),%rax
  mov    %rbx,%r12
  and    $0xfffffffffffffffe,%rax
  add    $0x2,%rax
  mov    %rax,0x70(%rsp)
L347:
  mov    0x98(%rsp),%rbx
  mov    0x70(%rsp),%r15
  add    %rbx,0x50(%rsp)
  mov    0x78(%rsp),%rbx
  add    %rbx,0x58(%rsp)
  mov    0x58(%rsp),%rsi
  mov    0x90(%rsp),%rbx
  mov    0x50(%rsp),%rax
  add    %rsi,%r15
  mov    %rsi,%r11
  lea    (%rbx,%rax,8),%rax
  mov    %r15,0x10(%rsp)
  mov    0x68(%rsp),%rbx
  lea    (%r14,%rbx,1),%rbp
  add    %r13,%rbx
  nopl   0x0(%rax)
L407:
  mov    0x8(%rsp),%rsi
  mov    0x38(%rsp),%r15
  add    $0x2,%r11
  add    %rax,%rsi
  lea    (%r15,%rsi,1),%r8
  mov    0x18(%rsp),%r15
  lea    (%r15,%rsi,1),%r9
  mov    0x28(%rsp),%r15
  lea    (%r15,%rsi,1),%r10
  lea    0x8(%rax),%r15
  add    %r13,%r10
  cmp    %r15,%r10
  seta   %r10b
  add    0x20(%rsp),%rsi
  add    %r13,%rsi
  cmp    %rsi,%rax
  seta   %sil
  or     %esi,%r10d
  test   %r10b,0x30(%rsp)
  je     L942
  lea    0x0(%r13,%r8,1),%rsi
  lea    (%r14,%r9,1),%r10
  cmp    %rsi,%r10
  seta   %sil
  add    %r14,%r8
  add    %r13,%r9
  cmp    %r8,%r9
  seta   %r8b
  or     %sil,%r8b
  je     L942
  movsd  (%rax),%xmm4
  movsd  0x8(%rax),%xmm3
  mov    %rbp,%r8
  mov    %rbx,%rsi
  xor    %r9d,%r9d
  unpcklpd %xmm4,%xmm4
  unpcklpd %xmm3,%xmm3
  nopw   0x0(%rax,%rax,1)
L20:
  movsd  (%r8),%xmm0
  movsd  0x8(%r8),%xmm1
  add    $0x1,%r9
  movhpd (%r8,%rdx,1),%xmm0
  movhpd 0x8(%rdx,%r8,1),%xmm1
  add    %rcx,%r8
  movapd %xmm0,%xmm2
  movapd %xmm1,%xmm5
  mulpd  %xmm4,%xmm2
  mulpd  %xmm3,%xmm5
  mulpd  %xmm3,%xmm0
  mulpd  %xmm4,%xmm1
  subpd  %xmm5,%xmm2
  addpd  %xmm1,%xmm0
  movlpd %xmm2,(%rsi)
  movhpd %xmm2,(%rsi,%rdx,1)
  movlpd %xmm0,0x8(%rsi)
  movhpd %xmm0,0x8(%rdx,%rsi,1)
  add    %rcx,%rsi
  cmp    %r9,%rdi
  jne    L20
  mov    0x40(%rsp),%rsi
  lea    (%r11,%rsi,1),%r10
  cmp    %r12,0x48(%rsp)
  je     L178
  lea    0x0(,%r10,8),%rsi
  lea    -0x8(%r14,%rsi,1),%r9
  lea    (%r14,%rsi,1),%r8
  movsd  (%r9),%xmm0
  movsd  (%r8),%xmm1
  mulsd  (%rax),%xmm0
  mulsd  0x8(%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r13,%rsi,1)
  movsd  (%r9),%xmm0
  movsd  (%r8),%xmm1
  mulsd  0x8(%rax),%xmm0
  mulsd  (%rax),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%r10,8)
L178:
  add    $0x10,%rax
  add    $0x10,%rbp
  add    $0x10,%rbx
  cmp    0x10(%rsp),%r11
  jne    L407
  mov    0x88(%rsp),%rsi
  addq   $0x1,0x60(%rsp)
  add    %rsi,0x8(%rsp)
  mov    0x60(%rsp),%rax
  mov    0x80(%rsp),%rsi
  add    %rsi,0x68(%rsp)
  cmp    %rax,0xc8(%rsp)
  jne    L347
  jmp    L1327
  nop
L942:
  lea    0x8(%rbp),%rsi
  lea    0x8(%rbx),%r8
  xor    %r9d,%r9d
  nopl   0x0(%rax,%rax,1)
L1067:
  movsd  -0x8(%rsi),%xmm0
  movsd  (%rsi),%xmm1
  add    $0x1,%r9
  mulsd  0x8(%rax),%xmm1
  mulsd  (%rax),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r8)
  movsd  -0x8(%rsi),%xmm0
  movsd  (%rsi),%xmm1
  add    %rdx,%rsi
  mulsd  0x8(%rax),%xmm0
  mulsd  (%rax),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%r8)
  add    %rdx,%r8
  cmp    %r9,%r12
  jne    L1067
  jmp    L178
L1251:
  mov    0x50(%rsp),%rbx
  mov    0xa8(%rsp),%r9
  lea    0x8(%rax),%rdx
  add    0x18(%rsp),%r9
  lea    (%r8,%rbx,1),%rsi
  lea    (%rdi,%rbx,1),%rcx
  mov    0x48(%rsp),%rbx
  lea    (%rbx,%r9,8),%r9
  nopw   0x0(%rax,%rax,1)
L15:
  movsd  -0x8(%rsi),%xmm0
  subsd  (%rcx),%xmm0
  add    $0x10,%rsi
  add    $0x10,%rcx
  movsd  %xmm0,-0x8(%rdx)
  movsd  -0x18(%rsi),%xmm0
  addsd  -0x10(%rcx),%xmm0
  movsd  %xmm0,-0x8(%rdx,%r12,8)
  movsd  -0x10(%rsi),%xmm0
  addsd  -0x18(%rcx),%xmm0
  movsd  %xmm0,(%rdx)
  movsd  -0x10(%rsi),%xmm0
  subsd  -0x18(%rcx),%xmm0
  movsd  %xmm0,(%rdx,%r12,8)
  add    $0x10,%rdx
  cmp    %rdx,%r9
  jne    L15
  jmp    L1356
L984:
  cmpq   $0x1,0x38(%rsp)
  jle    L175
  cmp    $0x2,%r12
  jle    L175
  mov    0x70(%rsp),%rsi
  test   %rsi,%rsi
  jle    L175
  mov    0x78(%rsp),%rbx
  xor    %r9d,%r9d
  mov    0x98(%rsp),%rdi
  mov    $0x1,%r10d
  mov    %rbx,%rax
  sub    %rbx,%rdi
  neg    %rax
  sub    %rbx,%rdi
  sub    %rbx,%rax
  mov    %rax,%rbp
  lea    -0x3(%r12),%rax
  and    $0xfffffffffffffffe,%rax
  lea    0x2(%rax),%rbx
L985:
  add    0x78(%rsp),%r9
  lea    (%rbx,%r9,1),%r15
  mov    %r9,%r8
L1400:
  add    $0x2,%r8
  xor    %ecx,%ecx
  lea    (%r8,%rdi,1),%rdx
  mov    %r8,%rax
  nopl   (%rax)
L741:
  movsd  -0x8(%r13,%rax,8),%xmm0
  subsd  0x0(%r13,%rdx,8),%xmm0
  add    $0x1,%rcx
  movsd  %xmm0,-0x8(%r14,%rax,8)
  movsd  -0x8(%r13,%rax,8),%xmm0
  addsd  0x0(%r13,%rdx,8),%xmm0
  movsd  %xmm0,-0x8(%r14,%rdx,8)
  movsd  0x0(%r13,%rax,8),%xmm0
  addsd  -0x8(%r13,%rdx,8),%xmm0
  movsd  %xmm0,(%r14,%rax,8)
  movsd  0x0(%r13,%rax,8),%xmm0
  add    %r12,%rax
  subsd  -0x8(%r13,%rdx,8),%xmm0
  movsd  %xmm0,(%r14,%rdx,8)
  add    %r12,%rdx
  cmp    %rcx,%rsi
  jne    L741
  cmp    %r15,%r8
  jne    L1400
  add    $0x1,%r10
  add    %rbp,%rdi
  cmp    %r10,0x38(%rsp)
  jne    L985
  jmp    L175
L159:
  mov    %r12,%rax
  not    %rax
  mov    %rax,0x40(%rsp)
  cmp    $0x2,%r12
  jle    L1327
  mov    0x78(%rsp),%rax
  lea    -0x3(%r12),%rbx
  movq   $0x0,0x48(%rsp)
  lea    0x0(,%r12,8),%rbp
  mov    %rbx,0x20(%rsp)
  mov    0x1b8(%rsp),%r9
  lea    0x0(,%rax,8),%rdi
  mov    %r12,0x18(%rsp)
  mov    %rdi,0x80(%rsp)
  lea    0x28(%r13,%rax,8),%rdi
  lea    0x8(%r14,%rax,8),%rax
  mov    %rax,0x60(%rsp)
  mov    0x1b8(%rsp),%rax
  mov    %rdi,0x68(%rsp)
  movq   $0x1,0x50(%rsp)
  add    $0x20,%rax
  mov    %rax,0x10(%rsp)
  mov    %rbx,%rax
  shr    %rax
  lea    0x1(%rax),%rdi
  mov    %rdi,%rax
  mov    %rdi,0x38(%rsp)
  and    $0xfffffffffffffffe,%rdi
  shr    %rax
  mov    %rdi,0x30(%rsp)
  shl    $0x5,%rax
  mov    %rax,%r15
  lea    (%rdi,%rdi,1),%rax
  mov    %rax,0x28(%rsp)
  mov    %rbx,%rax
  and    $0xfffffffffffffffe,%rax
  add    $0x2,%rax
  mov    %rax,0x58(%rsp)
L102:
  mov    0x18(%rsp),%rdi
  mov    0x28(%rsp),%rcx
  movq   $0x0,0x8(%rsp)
  add    %rdi,0x40(%rsp)
  mov    0x40(%rsp),%rax
  mov    0x68(%rsp),%rbx
  mov    0x78(%rsp),%rdi
  lea    0x2(%rax,%rcx,1),%r10
  add    %rdi,0x48(%rsp)
  mov    0x60(%rsp),%r8
  mov    0x1b8(%rsp),%rax
  shl    $0x3,%r10
  mov    0x48(%rsp),%r12
  lea    -0x20(%rbx),%rdi
  lea    -0x8(%rax,%r10,1),%r11
  add    %rax,%r10
  xchg   %ax,%ax
L495:
  cmp    %rbx,%r9
  setae  %dl
  cmp    0x10(%rsp),%rdi
  setae  %al
  or     %edx,%eax
  cmpq   $0x1,0x20(%rsp)
  seta   %dl
  test   %dl,%al
  je     L739
  lea    0x20(%r8),%rax
  cmp    %rax,%rdi
  setae  %cl
  cmp    %rbx,%r8
  setae  %dl
  xor    %eax,%eax
  or     %cl,%dl
  je     L739
  nopl   0x0(%rax,%rax,1)
L97:
  movupd (%r8,%rax,1),%xmm0
  movupd (%r9,%rax,1),%xmm4
  movupd 0x10(%r8,%rax,1),%xmm1
  movlpd 0x8(%r8,%rax,1),%xmm1
  movhpd 0x10(%r8,%rax,1),%xmm0
  movupd 0x10(%r9,%rax,1),%xmm3
  movlpd 0x8(%r9,%rax,1),%xmm3
  movhpd 0x10(%r9,%rax,1),%xmm4
  movapd %xmm0,%xmm2
  movapd %xmm1,%xmm5
  mulpd  %xmm4,%xmm2
  mulpd  %xmm3,%xmm5
  mulpd  %xmm4,%xmm1
  mulpd  %xmm3,%xmm0
  subpd  %xmm5,%xmm2
  addpd  %xmm1,%xmm0
  movapd %xmm2,%xmm1
  unpcklpd %xmm0,%xmm1
  unpckhpd %xmm0,%xmm2
  movups %xmm1,(%rdi,%rax,1)
  movups %xmm2,0x10(%rdi,%rax,1)
  add    $0x20,%rax
  cmp    %r15,%rax
  jne    L97
  mov    0x28(%rsp),%rax
  mov    0x38(%rsp),%rsi
  add    %r12,%rax
  cmp    %rsi,0x30(%rsp)
  je     L459
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  lea    -0x8(%r14,%rdx,1),%rsi
  lea    (%r14,%rdx,1),%rcx
  movsd  (%rsi),%xmm0
  movsd  (%rcx),%xmm1
  mulsd  (%r11),%xmm0
  mulsd  (%r10),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r13,%rdx,1)
  movsd  (%rsi),%xmm0
  movsd  (%rcx),%xmm1
  mulsd  (%r10),%xmm0
  mulsd  (%r11),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rax,8)
L459:
  addq   $0x1,0x8(%rsp)
  add    %rbp,%rdi
  add    %rbp,%rbx
  add    %rbp,%r8
  mov    0x8(%rsp),%rax
  add    0x18(%rsp),%r12
  cmp    %rax,0x70(%rsp)
  jne    L495
  mov    0x80(%rsp),%rbx
  addq   $0x1,0x50(%rsp)
  add    %rbp,%r9
  add    %rbx,0x68(%rsp)
  mov    0x50(%rsp),%rax
  add    %rbx,0x60(%rsp)
  add    %rbp,0x10(%rsp)
  cmp    %rax,0xc8(%rsp)
  jne    L102
  jmp    L1327
L517:
  mov    0xc8(%rsp),%rax
  imul   0x78(%rsp),%rax
  cmpq   $0x1,0x38(%rsp)
  mov    %rax,0x98(%rsp)
  lea    (%r12,%r12,1),%rax
  mov    %rax,0x18(%rsp)
  jle    L11
L1290:
  cmp    $0x1,%r12
  je     L1087
L1:
  mov    0x138(%rsp),%rdi
  cmp    %rdi,0x70(%rsp)
  jle    L1166
  cmp    $0x2,%r12
  jle    L1087
  mov    0x70(%rsp),%r9
  test   %r9,%r9
  jle    L1087
  mov    0x78(%rsp),%rdi
  mov    0x18(%rsp),%rsi
  movq   $0x1,0x8(%rsp)
  lea    0x10(%r14),%r15
  mov    %r12,0x50(%rsp)
  lea    0x0(,%r12,8),%r8
  mov    %rdi,%rax
  mov    %r13,0x58(%rsp)
  neg    %rax
  mov    %r14,0x1a8(%rsp)
  mov    %rax,%rcx
  mov    %rdi,%rax
  mov    %r15,0x48(%rsp)
  sub    %rdi,%rcx
  mov    0x98(%rsp),%rdi
  mov    %r11,0x1b0(%rsp)
  mov    %rcx,0x20(%rsp)
  mov    %rsi,%rcx
  add    %rax,%rsi
  sub    %rax,%rdi
  sub    %rax,%rcx
  mov    %rsi,0x30(%rsp)
  sub    %rax,%rdi
  lea    -0x3(%r12),%rax
  mov    %rsi,0x18(%rsp)
  and    $0xfffffffffffffffe,%rax
  mov    %rcx,0x28(%rsp)
  add    $0x2,%rax
  mov    %rcx,0x10(%rsp)
  xor    %ecx,%ecx
  mov    %rax,0x40(%rsp)
  mov    %rcx,%r15
L396:
  mov    0x48(%rsp),%rax
  add    0x78(%rsp),%r15
  mov    %r15,%r10
  lea    (%rax,%r15,8),%r11
  mov    0x40(%rsp),%rax
  lea    (%rax,%r15,1),%r14
  mov    0x10(%rsp),%rax
  lea    0x2(%rax),%r13
  mov    0x18(%rsp),%rax
  lea    -0x2(%rax),%r12
  nopl   0x0(%rax,%rax,1)
L754:
  mov    %r10,%rax
  mov    %r12,%rdx
  add    $0x2,%r10
  xor    %esi,%esi
  lea    0x0(%r13,%rax,1),%rcx
  sub    %rax,%rdx
  mov    %r11,%rax
  nop
L66:
  movsd  -0x8(%rbx,%rcx,8),%xmm0
  addsd  -0x8(%rbx,%rdx,8),%xmm0
  add    $0x1,%rsi
  movsd  %xmm0,-0x8(%rax)
  movsd  -0x8(%rbx,%rcx,8),%xmm0
  subsd  -0x8(%rbx,%rdx,8),%xmm0
  movsd  %xmm0,-0x8(%rax,%rdi,8)
  movsd  (%rbx,%rcx,8),%xmm0
  subsd  (%rbx,%rdx,8),%xmm0
  movsd  %xmm0,(%rax)
  movsd  (%rbx,%rcx,8),%xmm0
  addsd  (%rbx,%rdx,8),%xmm0
  add    %rbp,%rcx
  add    %rbp,%rdx
  movsd  %xmm0,(%rax,%rdi,8)
  add    %r8,%rax
  cmp    %rsi,%r9
  jne    L66
  add    $0x10,%r11
  cmp    %r10,%r14
  jne    L754
  mov    0x28(%rsp),%rsi
  addq   $0x1,0x8(%rsp)
  add    %rsi,0x10(%rsp)
  mov    0x8(%rsp),%rax
  mov    0x30(%rsp),%rsi
  add    0x20(%rsp),%rdi
  add    %rsi,0x18(%rsp)
  cmp    %rax,0x38(%rsp)
  jne    L396
  mov    0x50(%rsp),%r12
  mov    0x58(%rsp),%r13
  mov    0x1a8(%rsp),%r14
  mov    0x1b0(%rsp),%r11
  jmp    L1087
L986:
  mov    0xc8(%rsp),%rax
  imul   0x78(%rsp),%rax
  cmpq   $0x1,0x38(%rsp)
  mov    %rax,0x98(%rsp)
  lea    (%r12,%r12,1),%rax
  mov    %rax,0x18(%rsp)
  jle    L11
  cmpq   $0x0,0x70(%rsp)
  jg     L1068
  jmp    L1290
  nopl   (%rax)
L739:
  mov    0x58(%rsp),%rax
  lea    0x8(%r9),%rdx
  lea    (%rax,%r12,1),%rcx
  mov    %r12,%rax
  nopl   0x0(%rax,%rax,1)
L497:
  add    $0x2,%rax
  add    $0x10,%rdx
  movsd  -0x8(%r14,%rax,8),%xmm0
  movsd  (%r14,%rax,8),%xmm1
  mulsd  -0x18(%rdx),%xmm0
  mulsd  -0x10(%rdx),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r13,%rax,8)
  movsd  -0x8(%r14,%rax,8),%xmm0
  movsd  (%r14,%rax,8),%xmm1
  mulsd  -0x10(%rdx),%xmm0
  mulsd  -0x18(%rdx),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rax,8)
  cmp    %rax,%rcx
  jne    L497
  jmp    L459
L1166:
  cmpq   $0x1,0x38(%rsp)
  jle    L175
  mov    0x70(%rsp),%rcx
  test   %rcx,%rcx
  jle    L1087
  cmp    $0x2,%r12
  jle    L1087
  mov    0x78(%rsp),%rdi
  xor    %edx,%edx
  mov    %r13,0x30(%rsp)
  lea    0x10(%r14),%r15
  mov    0x98(%rsp),%rsi
  mov    %r14,0x1a8(%rsp)
  shl    $0x3,%rbp
  mov    %rdi,%rax
  mov    %r11,0x1b0(%rsp)
  mov    %rcx,%r11
  movq   $0x0,0x8(%rsp)
  neg    %rax
  sub    %rdi,%rsi
  movq   $0x1,0x10(%rsp)
  sub    %rdi,%rax
  sub    %rdi,%rsi
  lea    0x20(%r14),%rdi
  mov    %rax,0x20(%rsp)
  lea    0x10(%rbx),%rax
  mov    %rdx,%r14
  mov    %rdi,%r13
  mov    %rax,0x28(%rsp)
  lea    -0x3(%r12),%rax
  and    $0xfffffffffffffffe,%rax
  mov    %rax,%rbx
L304:
  mov    0x78(%rsp),%rax
  add    %rax,0x8(%rsp)
  xor    %r10d,%r10d
  mov    0x28(%rsp),%rax
  add    0x18(%rsp),%r14
  mov    0x8(%rsp),%r8
  lea    (%rax,%r14,8),%r9
  xchg   %ax,%ax
L443:
  lea    (%rbx,%r8,1),%rdx
  lea    -0x20(%r9),%rcx
  lea    0x0(%r13,%rdx,8),%rdi
  lea    (%r15,%r8,8),%rax
  mov    %r9,%rdx
  nopl   0x0(%rax)
L795:
  movsd  -0x8(%rdx),%xmm0
  addsd  -0x8(%rcx),%xmm0
  add    $0x10,%rdx
  sub    $0x10,%rcx
  movsd  %xmm0,-0x8(%rax)
  movsd  -0x18(%rdx),%xmm0
  subsd  0x8(%rcx),%xmm0
  movsd  %xmm0,-0x8(%rax,%rsi,8)
  movsd  -0x10(%rdx),%xmm0
  subsd  0x10(%rcx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x10(%rdx),%xmm0
  addsd  0x10(%rcx),%xmm0
  movsd  %xmm0,(%rax,%rsi,8)
  add    $0x10,%rax
  cmp    %rdi,%rax
  jne    L795
  add    $0x1,%r10
  add    %r12,%r8
  add    %rbp,%r9
  cmp    %r10,%r11
  jne    L443
  addq   $0x1,0x10(%rsp)
  add    0x20(%rsp),%rsi
  mov    0x10(%rsp),%rax
  cmp    %rax,0x38(%rsp)
  jne    L304
  mov    0x30(%rsp),%r13
  mov    0x1a8(%rsp),%r14
  mov    0x1b0(%rsp),%r11
  jmp    L1087
L1196:
  mov    0x78(%rsp),%rdi
  mov    0x98(%rsp),%rsi
  xor    %r10d,%r10d
  mov    %rbp,0x20(%rsp)
  mov    0x18(%rsp),%r9
  lea    0x0(,%r12,8),%r8
  mov    $0x1,%r15d
  mov    %rdi,%rax
  sub    %rdi,%rsi
  neg    %rax
  sub    %rdi,%rsi
  sub    %rdi,%rax
  lea    0x0(,%rbp,8),%rdi
  mov    0x70(%rsp),%rbp
  mov    %rax,0x8(%rsp)
  lea    -0x8(%rbx),%rax
  mov    %rax,0x10(%rsp)
L342:
  mov    0x10(%rsp),%rax
  add    0x78(%rsp),%r10
  xor    %ecx,%ecx
  lea    (%r14,%r10,8),%rdx
  lea    (%rax,%r9,8),%rax
  nop
L394:
  movsd  (%rax),%xmm0
  add    $0x1,%rcx
  addsd  %xmm0,%xmm0
  movsd  %xmm0,(%rdx)
  movsd  0x8(%rax),%xmm0
  add    %rdi,%rax
  addsd  %xmm0,%xmm0
  movsd  %xmm0,(%rdx,%rsi,8)
  add    %r8,%rdx
  cmp    %rcx,%rbp
  jne    L394
  add    0x18(%rsp),%r9
  add    $0x1,%r15
  add    0x8(%rsp),%rsi
  cmp    %r15,0x38(%rsp)
  jne    L342
  mov    0x20(%rsp),%rbp
  jmp    L1
L759:
  mov    0x98(%rsp),%rsi
  mov    0x78(%rsp),%r10
  xor    %edi,%edi
  mov    $0x1,%r8d
  mov    0x70(%rsp),%r9
L1058:
  add    %r10,%rdi
  sub    %r10,%rsi
  xor    %ecx,%ecx
  mov    %rsi,%rdx
  mov    %rdi,%rax
L297:
  movsd  0x0(%r13,%rax,8),%xmm0
  subsd  0x0(%r13,%rdx,8),%xmm0
  add    $0x1,%rcx
  movsd  %xmm0,(%r14,%rax,8)
  movsd  0x0(%r13,%rax,8),%xmm0
  add    %r12,%rax
  addsd  0x0(%r13,%rdx,8),%xmm0
  movsd  %xmm0,(%r14,%rdx,8)
  add    %r12,%rdx
  cmp    %rcx,%r9
  jne    L297
  add    $0x1,%r8
  cmp    %r8,0x38(%rsp)
  jg     L1058
  jmp    L1291
L498:
  mov    0x1a0(%rsp),%rdx
  xor    %eax,%eax
L1236:
  movsd  (%r11,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    0x130(%rsp),%rax
  jne    L1236
  jmp    L0
L11:
  cmp    $0x1,%r12
  jne    L175
  jmp    L1327
L402:
  mov    %rdx,%r15
  mov    %rbx,%r8
  mov    %r14,%r9
  xor    %r10d,%r10d
  lea    0x0(,%rbp,8),%rdi
  lea    0x0(,%r12,8),%rsi
L1260:
  mov    %r9,%rcx
  mov    %r8,%rdx
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L515:
  movsd  (%rdx),%xmm0
  add    $0x1,%rax
  add    %rdi,%rdx
  movsd  %xmm0,(%rcx)
  add    %rsi,%rcx
  cmp    %rax,%r15
  jne    L515
  add    $0x1,%r10
  add    $0x8,%r8
  add    $0x8,%r9
  cmp    %r10,%r12
  jne    L1260
  jmp    L955
```

tool_calls: name: unpack_data
arguments: data_label: D2
datatype: f64
name: unpack_data
arguments: data_label: D0
datatype: f64

tool: D2: 6.283185307179586

tool: D0: 1

assistant: ```c
void dradbg(long ido, long ip, long l1, long idl1, const double *cc, double *c1, double *c2, double *ch, double *ch2, const double *wa) {
  static double tpi = 6.28318530717958647692528676655900577;
  long idij, ipph, i, j, k, l, ik, is, t0, t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12;
  double dc2, ai1, ai2, ar1, ar2, ds2;
  long nbd;
  double dcp, arg, dsp, ar1h, ar2h;
  long ipp2;
  t10 = ip * ido;
  t0 = l1 * ido;
  arg = tpi / (double)ip;
  dcp = cos(arg);
  dsp = sin(arg);
  nbd = (ido - 1) >> 1;
  ipp2 = ip;
  ipph = (ip + 1) >> 1;
  if (ido < l1)
    goto L103;
  t1 = 0;
  t2 = 0;
  for (k = 0; k < l1; k++) {
    t3 = t1;
    t4 = t2;
    for (i = 0; i < ido; i++) {
      ch[t3] = cc[t4];
      t3++;
      t4++;
    }
    t1 += ido;
    t2 += t10;
  }
  goto L106;
L103:
  t1 = 0;
  for (i = 0; i < ido; i++) {
    t2 = t1;
    t3 = t1;
    for (k = 0; k < l1; k++) {
      ch[t2] = cc[t3];
      t2 += ido;
      t3 += t10;
    }
    t1++;
  }
L106:
  t1 = 0;
  t2 = ipp2 * t0;
  t7 = (t5 = ido << 1);
  for (j = 1; j < ipph; j++) {
    t1 += t0;
    t2 -= t0;
    t3 = t1;
    t4 = t2;
    t6 = t5;
    for (k = 0; k < l1; k++) {
      ch[t3] = cc[t6 - 1] + cc[t6 - 1];
      ch[t4] = cc[t6] + cc[t6];
      t3 += ido;
      t4 += ido;
      t6 += t10;
    }
    t5 += t7;
  }
  if (ido == 1)
    goto L116;
  if (nbd < l1)
    goto L112;
  t1 = 0;
  t2 = ipp2 * t0;
  t7 = 0;
  for (j = 1; j < ipph; j++) {
    t1 += t0;
    t2 -= t0;
    t3 = t1;
    t4 = t2;
    t7 += (ido << 1);
    t8 = t7;
    for (k = 0; k < l1; k++) {
      t5 = t3;
      t6 = t4;
      t9 = t8;
      t11 = t8;
      for (i = 2; i < ido; i += 2) {
        t5 += 2;
        t6 += 2;
        t9 += 2;
        t11 -= 2;
        ch[t5 - 1] = cc[t9 - 1] + cc[t11 - 1];
        ch[t6 - 1] = cc[t9 - 1] - cc[t11 - 1];
        ch[t5] = cc[t9] - cc[t11];
        ch[t6] = cc[t9] + cc[t11];
      }
      t3 += ido;
      t4 += ido;
      t8 += t10;
    }
  }
  goto L116;
L112:
  t1 = 0;
  t2 = ipp2 * t0;
  t7 = 0;
  for (j = 1; j < ipph; j++) {
    t1 += t0;
    t2 -= t0;
    t3 = t1;
    t4 = t2;
    t7 += (ido << 1);
    t8 = t7;
    t9 = t7;
    for (i = 2; i < ido; i += 2) {
      t3 += 2;
      t4 += 2;
      t8 += 2;
      t9 -= 2;
      t5 = t3;
      t6 = t4;
      t11 = t8;
      t12 = t9;
      for (k = 0; k < l1; k++) {
        ch[t5 - 1] = cc[t11 - 1] + cc[t12 - 1];
        ch[t6 - 1] = cc[t11 - 1] - cc[t12 - 1];
        ch[t5] = cc[t11] - cc[t12];
        ch[t6] = cc[t11] + cc[t12];
        t5 += ido;
        t6 += ido;
        t11 += t10;
        t12 += t10;
      }
    }
  }
L116:
  ar1 = 1.;
  ai1 = 0.;
  t1 = 0;
  t9 = (t2 = ipp2 * idl1);
  t3 = (ip - 1) * idl1;
  for (l = 1; l < ipph; l++) {
    t1 += idl1;
    t2 -= idl1;
    ar1h = dcp * ar1 - dsp * ai1;
    ai1 = dcp * ai1 + dsp * ar1;
    ar1 = ar1h;
    t4 = t1;
    t5 = t2;
    t6 = 0;
    t7 = idl1;
    t8 = t3;
    for (ik = 0; ik < idl1; ik++) {
      c2[t4++] = ch2[t6++] + ar1 * ch2[t7++];
      c2[t5++] = ai1 * ch2[t8++];
    }
    dc2 = ar1;
    ds2 = ai1;
    ar2 = ar1;
    ai2 = ai1;
    t6 = idl1;
    t7 = t9 - idl1;
    for (j = 2; j < ipph; j++) {
      t6 += idl1;
      t7 -= idl1;
      ar2h = dc2 * ar2 - ds2 * ai2;
      ai2 = dc2 * ai2 + ds2 * ar2;
      ar2 = ar2h;
      t4 = t1;
      t5 = t2;
      t11 = t6;
      t12 = t7;
      for (ik = 0; ik < idl1; ik++) {
        c2[t4++] += ar2 * ch2[t11++];
        c2[t5++] += ai2 * ch2[t12++];
      }
    }
  }
  t1 = 0;
  for (j = 1; j < ipph; j++) {
    t1 += idl1;
    t2 = t1;
    for (ik = 0; ik < idl1; ik++)
      ch2[ik] += ch2[t2++];
  }
  t1 = 0;
  t2 = ipp2 * t0;
  for (j = 1; j < ipph; j++) {
    t1 += t0;
    t2 -= t0;
    t3 = t1;
    t4 = t2;
    for (k = 0; k < l1; k++) {
      ch[t3] = c1[t3] - c1[t4];
      ch[t4] = c1[t3] + c1[t4];
      t3 += ido;
      t4 += ido;
    }
  }
  if (ido == 1)
    goto L132;
  if (nbd < l1)
    goto L128;
  t1 = 0;
  t2 = ipp2 * t0;
  for (j = 1; j < ipph; j++) {
    t1 += t0;
    t2 -= t0;
    t3 = t1;
    t4 = t2;
    for (k = 0; k < l1; k++) {
      t5 = t3;
      t6 = t4;
      for (i = 2; i < ido; i += 2) {
        t5 += 2;
        t6 += 2;
        ch[t5 - 1] = c1[t5 - 1] - c1[t6];
        ch[t6 - 1] = c1[t5 - 1] + c1[t6];
        ch[t5] = c1[t5] + c1[t6 - 1];
        ch[t6] = c1[t5] - c1[t6 - 1];
      }
      t3 += ido;
      t4 += ido;
    }
  }
  goto L132;
L128:
  t1 = 0;
  t2 = ipp2 * t0;
  for (j = 1; j < ipph; j++) {
    t1 += t0;
    t2 -= t0;
    t3 = t1;
    t4 = t2;
    for (i = 2; i < ido; i += 2) {
      t3 += 2;
      t4 += 2;
      t5 = t3;
      t6 = t4;
      for (k = 0; k < l1; k++) {
        ch[t5 - 1] = c1[t5 - 1] - c1[t6];
        ch[t6 - 1] = c1[t5 - 1] + c1[t6];
        ch[t5] = c1[t5] + c1[t6 - 1];
        ch[t6] = c1[t5] - c1[t6 - 1];
        t5 += ido;
        t6 += ido;
      }
    }
  }
L132:
  if (ido == 1)
    return;
  for (ik = 0; ik < idl1; ik++)
    c2[ik] = ch2[ik];
  t1 = 0;
  for (j = 1; j < ip; j++) {
    t2 = (t1 += t0);
    for (k = 0; k < l1; k++) {
      c1[t2] = ch[t2];
      t2 += ido;
    }
  }
  if (nbd > l1)
    goto L139;
  is = -ido - 1;
  t1 = 0;
  for (j = 1; j < ip; j++) {
    is += ido;
    t1 += t0;
    idij = is;
    t2 = t1;
    for (i = 2; i < ido; i += 2) {
      t2 += 2;
      idij += 2;
      t3 = t2;
      for (k = 0; k < l1; k++) {
        c1[t3 - 1] = wa[idij - 1] * ch[t3 - 1] - wa[idij] * ch[t3];
        c1[t3] = wa[idij - 1] * ch[t3] + wa[idij] * ch[t3 - 1];
        t3 += ido;
      }
    }
  }
  return;
L139:
  is = -ido - 1;
  t1 = 0;
  for (j = 1; j < ip; j++) {
    is += ido;
    t1 += t0;
    t2 = t1;
    for (k = 0; k < l1; k++) {
      idij = is;
      t3 = t2;
      for (i = 2; i < ido; i += 2) {
        idij += 2;
        t3 += 2;
        c1[t3 - 1] = wa[idij - 1] * ch[t3 - 1] - wa[idij] * ch[t3];
        c1[t3] = wa[idij - 1] * ch[t3] + wa[idij] * ch[t3 - 1];
      }
      t2 += ido;
    }
  }
}
```