user: 这段汇编代码用C语言表示是什么样的:

```asm
<hqrx>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xf8,%rsp
  mov    %edi,0x24(%rsp)
  mov    %esi,0x88(%rsp)
  mov    %edx,%r14d
  mov    %rcx,%r15
  mov    %r8,0xe0(%rsp)
  mov    %r9,0xe8(%rsp)
  movl   $0x0,(%r9)
  test   %edi,%edi
  jle    L689
  mov    %edi,%ebx
  mov    %esi,%r12d
  mov    %r8,%rbp
  lea    0x1(%rdi),%ecx
  mov    %ecx,%r11d
  movslq %edi,%rsi
  shl    $0x3,%rsi
  mov    $0x1,%edi
  mov    $0x0,%r8d
  mov    $0x0,%r9d
  pxor   %xmm5,%xmm5
  movsd  %xmm5,0x68(%rsp)
  mov    $0x1,%r10d
  movq   D149(%rip),%xmm1
  jmp    L452
L593:
  movslq %r9d,%rax
  movsd  (%r15,%rax,8),%xmm0
  movslq %r8d,%rax
  movsd  %xmm0,0x0(%rbp,%rax,8)
  lea    0x1(%r8),%eax
  cltq
  movq   $0x0,0x0(%rbp,%rax,8)
L445:
  add    %ecx,%r9d
  add    $0x2,%r8d
  add    $0x1,%rdi
  cmp    %r11,%rdi
  je     L552
L452:
  mov    %r10d,%eax
  mov    %edi,%r10d
  cmp    %eax,%ebx
  jl     L489
  lea    -0x1(%rax),%edx
  imul   %ebx,%edx
  movslq %edx,%rdx
  add    %rdi,%rdx
  lea    (%r15,%rdx,8),%rdx
L383:
  movsd  -0x8(%rdx),%xmm0
  andpd  %xmm1,%xmm0
  addsd  0x68(%rsp),%xmm0
  movsd  %xmm0,0x68(%rsp)
  add    $0x1,%eax
  add    %rsi,%rdx
  cmp    %ecx,%eax
  jne    L383
L489:
  cmp    %edi,%r12d
  jg     L593
  cmp    %edi,%r14d
  jl     L593
  jmp    L445
L689:
  pxor   %xmm4,%xmm4
  movsd  %xmm4,0x68(%rsp)
L552:
  mov    0x88(%rsp),%ecx
  cmp    %ecx,%r14d
  jl     L24
  mov    0x24(%rsp),%ebx
  mov    %ebx,%eax
  not    %eax
  cltq
  shl    $0x3,%rax
  mov    %rax,0x58(%rsp)
  movslq %ebx,%rdx
  lea    0x8(,%rdx,8),%rax
  mov    %rax,0x50(%rsp)
  lea    -0x1(%rcx),%eax
  imul   %ebx,%eax
  cltq
  movslq %ecx,%rcx
  add    %rcx,%rax
  lea    (%r15,%rax,8),%rax
  mov    %rax,0xd8(%rsp)
  lea    0x0(,%rdx,8),%rax
  mov    %rax,0x38(%rsp)
  pxor   %xmm7,%xmm7
  movsd  %xmm7,0xc8(%rsp)
  movsd  %xmm7,0x28(%rsp)
  movsd  %xmm7,0x8(%rsp)
  movsd  %xmm7,(%rsp)
  movl   $0x0,0x48(%rsp)
  movl   $0x0,0x64(%rsp)
  mov    %r14d,%r13d
  jmp    L121
L299:
  mov    %edx,0x64(%rsp)
  jmp    L381
L148:
  mov    %edx,0x64(%rsp)
L381:
  mov    0x80(%rsp),%rax
  movsd  (%rax),%xmm7
  movsd  %xmm7,0x10(%rsp)
  cmp    %ebx,0x64(%rsp)
  je     L341
  mov    0x90(%rsp),%rax
  movsd  (%rax),%xmm10
  mov    0x78(%rsp),%rax
  movsd  (%rax),%xmm0
  mov    0x98(%rsp),%rax
  mulsd  (%rax),%xmm0
  movsd  %xmm0,0x18(%rsp)
  cmp    %r14d,0x64(%rsp)
  je     L565
  cmpl   $0x1e,0x74(%rsp)
  je     L224
  mov    0x74(%rsp),%eax
  cmp    $0xa,%eax
  je     L385
  cmp    $0x14,%eax
  je     L385
L390:
  addl   $0x1,0x74(%rsp)
  mov    0x64(%rsp),%eax
  mov    0x8c(%rsp),%edx
  cmp    %edx,%eax
  jg     L618
  mov    %eax,%edi
  mov    0xc0(%rsp),%rsi
  mov    0xb8(%rsp),%rax
  mov    0xb0(%rsp),%rcx
  movq   D149(%rip),%xmm3
  movsd  D38(%rip),%xmm4
  mov    0x58(%rsp),%r8
  jmp    L65
L565:
  subsd  0x10(%rsp),%xmm10
  movapd %xmm10,%xmm0
  mulsd  D100(%rip),%xmm0
  movsd  %xmm0,(%rsp)
  mulsd  %xmm0,%xmm0
  addsd  0x18(%rsp),%xmm0
  movsd  %xmm0,0x8(%rsp)
  andpd  D149(%rip),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L103
  sqrtsd %xmm0,%xmm0
L376:
  movsd  0x10(%rsp),%xmm5
  addsd  0xc8(%rsp),%xmm5
  movsd  %xmm5,0x10(%rsp)
  pxor   %xmm1,%xmm1
  comisd 0x8(%rsp),%xmm1
  jbe    L513
  mov    0x8c(%rsp),%r13d
  lea    0x0(%r13,%r13,1),%eax
  movslq %eax,%rdx
  movsd  0x10(%rsp),%xmm1
  addsd  (%rsp),%xmm1
  mov    0xe0(%rsp),%rbx
  movsd  %xmm1,(%rbx,%rdx,8)
  movsd  %xmm1,0x10(%rbx,%rdx,8)
  lea    0x1(%rax),%edx
  movslq %edx,%rdx
  movsd  %xmm0,(%rbx,%rdx,8)
  add    $0x3,%eax
  cltq
  xorpd  D34(%rip),%xmm0
  movsd  %xmm0,(%rbx,%rax,8)
  mov    %r14d,0x64(%rsp)
  jmp    L509
L103:
  call   <sqrt@plt>
  jmp    L376
L224:
  mov    0xe8(%rsp),%rax
  mov    %ebx,(%rax)
L24:
  add    $0xf8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L385:
  movsd  0xc8(%rsp),%xmm5
  addsd  0x10(%rsp),%xmm5
  movsd  %xmm5,0xc8(%rsp)
  mov    0x88(%rsp),%edx
  cmp    %ebx,%edx
  jg     L541
  mov    0x70(%rsp),%eax
  lea    0x1(%rax),%ecx
  mov    0xd8(%rsp),%rax
L60:
  movsd  -0x8(%rax),%xmm0
  subsd  0x10(%rsp),%xmm0
  movsd  %xmm0,-0x8(%rax)
  add    $0x1,%edx
  add    0x50(%rsp),%rax
  cmp    %ecx,%edx
  jne    L60
L541:
  mov    0x78(%rsp),%rax
  movsd  (%rax),%xmm0
  movq   D149(%rip),%xmm2
  andpd  %xmm2,%xmm0
  mov    0xd0(%rsp),%rax
  movsd  (%rax),%xmm1
  andpd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm10
  mulsd  D120(%rip),%xmm10
  movsd  %xmm10,0x10(%rsp)
  movapd %xmm0,%xmm1
  mulsd  D165(%rip),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x18(%rsp)
  jmp    L390
L112:
  mov    %r10d,%edx
L65:
  movsd  -0x8(%rcx),%xmm5
  movsd  0x10(%rsp),%xmm8
  subsd  %xmm5,%xmm8
  movapd %xmm10,%xmm1
  subsd  %xmm5,%xmm1
  movapd %xmm8,%xmm0
  mulsd  %xmm1,%xmm0
  subsd  0x18(%rsp),%xmm0
  divsd  (%rcx),%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  (%rax),%xmm7
  movapd %xmm7,%xmm6
  subsd  %xmm5,%xmm6
  subsd  %xmm8,%xmm6
  subsd  %xmm1,%xmm6
  movsd  0x8(%rax),%xmm8
  movapd %xmm0,%xmm1
  andpd  %xmm3,%xmm1
  movapd %xmm6,%xmm9
  andpd  %xmm3,%xmm9
  addsd  %xmm9,%xmm1
  movapd %xmm8,%xmm9
  andpd  %xmm3,%xmm9
  addsd  %xmm9,%xmm1
  movapd %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  movsd  %xmm2,(%rsp)
  divsd  %xmm1,%xmm6
  movsd  %xmm6,0x8(%rsp)
  divsd  %xmm1,%xmm8
  cmp    %edi,%edx
  je     L668
  andpd  D149(%rip),%xmm6
  movapd %xmm6,%xmm1
  movapd %xmm8,%xmm0
  andpd  D149(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  0x8(%rsi),%xmm0
  andpd  D149(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  (%rsi),%xmm0
  andpd  D149(%rip),%xmm0
  andpd  D149(%rip),%xmm5
  addsd  %xmm5,%xmm0
  andpd  D149(%rip),%xmm7
  addsd  %xmm7,%xmm0
  andpd  D149(%rip),%xmm2
  mulsd  %xmm4,%xmm2
  mulsd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jae    L350
  lea    -0x1(%rdx),%r10d
  add    %r8,%rcx
  add    %r8,%rax
  add    %r8,%rsi
  cmp    %edi,%edx
  jne    L112
  mov    %edx,0x48(%rsp)
  movsd  %xmm8,0x28(%rsp)
  jmp    L618
L668:
  mov    %edx,0x48(%rsp)
  movsd  %xmm8,0x28(%rsp)
L618:
  mov    0x48(%rsp),%eax
  lea    0x2(%rax),%edi
  cmp    %ebx,%edi
  jg     L669
  mov    0x48(%rsp),%r10d
  lea    -0x1(%r10),%edx
  mov    0x24(%rsp),%esi
  mov    %esi,%eax
  imul   %edx,%eax
  lea    0x1(%r10,%rax,1),%ecx
  movslq %ecx,%rcx
  movq   $0x0,(%r15,%rcx,8)
  lea    0x1(%rsi),%r9d
  lea    0x2(%r10,%rax,1),%ecx
  add    %esi,%eax
  cltq
  movslq %edi,%rdx
  add    %rdx,%rax
  lea    (%r15,%rax,8),%rdx
  mov    0x70(%rsp),%eax
  lea    0x1(%rax),%r8d
  mov    %edi,%eax
  jmp    L661
L350:
  mov    %edx,0x48(%rsp)
  movsd  %xmm8,0x28(%rsp)
  jmp    L618
L432:
  movslq %ecx,%rsi
  movq   $0x0,(%r15,%rsi,8)
L682:
  add    %r9d,%ecx
  add    0x50(%rsp),%rdx
L661:
  add    $0x1,%eax
  cmp    %r8d,%eax
  je     L669
  movq   $0x0,(%rdx)
  cmp    %eax,%edi
  jne    L432
  jmp    L682
L669:
  mov    0x48(%rsp),%ebp
  cmp    %r14d,%ebp
  jg     L243
  mov    0x24(%rsp),%edi
  lea    0x1(%rdi),%ecx
  mov    %ecx,0x60(%rsp)
  lea    -0x1(%rbp),%r13d
  lea    -0x2(%rbp),%r12d
  imul   %edi,%r12d
  add    %r13d,%r12d
  imul   %edi,%r13d
  movslq %ebp,%rcx
  mov    %rcx,0x40(%rsp)
  lea    0x1(%rbp),%eax
  imul   %edi,%eax
  mov    %eax,0x4c(%rsp)
  jmp    L200
L278:
  movsd  (%rsp),%xmm0
  andpd  D149(%rip),%xmm0
  movsd  0x8(%rsp),%xmm1
  andpd  D149(%rip),%xmm1
  addsd  %xmm1,%xmm0
  movsd  0x28(%rsp),%xmm1
  andpd  D149(%rip),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%rsp)
  ucomisd D22(%rip),%xmm0
  jp     L631
  jne    L631
L688:
  add    $0x1,%ebp
  add    0x60(%rsp),%r12d
  mov    0x24(%rsp),%edi
  add    %edi,%r13d
  addq   $0x1,0x40(%rsp)
  add    %edi,0x4c(%rsp)
  cmp    %ebx,%ebp
  je     L243
L200:
  cmp    %ebp,0x48(%rsp)
  je     L639
  movslq %r12d,%rax
  movsd  (%r15,%rax,8),%xmm4
  movsd  %xmm4,(%rsp)
  lea    0x1(%r12),%eax
  cltq
  movsd  (%r15,%rax,8),%xmm7
  movsd  %xmm7,0x8(%rsp)
  pxor   %xmm2,%xmm2
  movsd  %xmm2,0x28(%rsp)
  cmp    %ebp,%r14d
  je     L278
  lea    0x2(%r12),%eax
  cltq
  movsd  (%r15,%rax,8),%xmm3
  movsd  %xmm3,0x28(%rsp)
  jmp    L278
L631:
  movsd  (%rsp),%xmm3
  movsd  0x10(%rsp),%xmm2
  divsd  %xmm2,%xmm3
  movsd  %xmm3,(%rsp)
  movsd  0x8(%rsp),%xmm5
  divsd  %xmm2,%xmm5
  movsd  %xmm5,0x8(%rsp)
  movsd  0x28(%rsp),%xmm4
  divsd  %xmm2,%xmm4
  movsd  %xmm4,0x28(%rsp)
L639:
  movsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm1
  mulsd  %xmm5,%xmm1
  movsd  (%rsp),%xmm4
  movapd %xmm4,%xmm0
  mulsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm0
  movsd  0x28(%rsp),%xmm7
  movapd %xmm7,%xmm1
  mulsd  %xmm7,%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm7,%xmm7
  ucomisd %xmm0,%xmm7
  ja     L319
  sqrtsd %xmm0,%xmm0
L358:
  movsd  (%rsp),%xmm1
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  cmp    %ebp,0x48(%rsp)
  je     L129
  movslq %r12d,%rax
  movapd %xmm0,%xmm1
  xorpd  D34(%rip),%xmm1
  mulsd  0x10(%rsp),%xmm1
  movsd  %xmm1,(%r15,%rax,8)
L174:
  movsd  (%rsp),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,(%rsp)
  movapd %xmm2,%xmm4
  divsd  %xmm0,%xmm4
  movsd  %xmm4,0x10(%rsp)
  movsd  0x8(%rsp),%xmm3
  movapd %xmm3,%xmm4
  divsd  %xmm0,%xmm4
  movsd  %xmm4,0x18(%rsp)
  movsd  0x28(%rsp),%xmm4
  movapd %xmm4,%xmm7
  divsd  %xmm0,%xmm7
  movsd  %xmm7,0x30(%rsp)
  divsd  %xmm2,%xmm3
  movsd  %xmm3,0x8(%rsp)
  divsd  %xmm2,%xmm4
  movapd %xmm4,%xmm2
  movsd  %xmm4,0x28(%rsp)
  cmp    %ebx,%ebp
  jg     L536
  movslq %r13d,%rax
  add    0x40(%rsp),%rax
  lea    (%r15,%rax,8),%rax
  mov    0x24(%rsp),%edi
  lea    0x2(%r12,%rdi,1),%esi
  mov    0x70(%rsp),%edi
  add    $0x1,%edi
  mov    %ebp,%ecx
  jmp    L63
L319:
  call   <sqrt@plt>
  jmp    L358
L129:
  mov    0x48(%rsp),%esi
  cmp    %esi,0x64(%rsp)
  je     L174
  movslq %r12d,%rax
  lea    (%r15,%rax,8),%rax
  movsd  (%rax),%xmm1
  xorpd  D34(%rip),%xmm1
  movsd  %xmm1,(%rax)
  jmp    L174
L340:
  movsd  (%rdx),%xmm0
  movsd  (%rsp),%xmm4
  movapd %xmm4,%xmm1
  mulsd  0x18(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx)
  movsd  -0x8(%rdx),%xmm0
  mulsd  0x10(%rsp),%xmm4
  subsd  %xmm4,%xmm0
  movsd  %xmm0,-0x8(%rdx)
  add    $0x1,%ecx
  add    0x38(%rsp),%rax
  add    0x24(%rsp),%esi
  cmp    %edi,%ecx
  je     L536
L63:
  mov    %rax,%rdx
  movsd  0x8(%rsp),%xmm0
  mulsd  (%rax),%xmm0
  movapd %xmm0,%xmm7
  addsd  -0x8(%rax),%xmm7
  movsd  %xmm7,(%rsp)
  cmp    %ebp,%r14d
  je     L340
  movslq %esi,%r8
  lea    (%r15,%r8,8),%r8
  movsd  (%r8),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm7,%xmm1
  movapd %xmm1,%xmm6
  movsd  %xmm1,(%rsp)
  movsd  0x30(%rsp),%xmm1
  mulsd  %xmm6,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%r8)
  jmp    L340
L536:
  lea    0x3(%rbp),%esi
  mov    %ebx,%edi
  mov    $0x0,%eax
  call   <imin@plt>
  mov    0x64(%rsp),%ecx
  cmp    %eax,%ecx
  jg     L688
  mov    0x24(%rsp),%esi
  lea    (%rsi,%r13,1),%r8d
  movslq %r8d,%r8
  movslq %ecx,%rdx
  lea    -0x1(%r8,%rdx,1),%rdx
  lea    (%r15,%rdx,8),%rdx
  mov    0x4c(%rsp),%esi
  lea    -0x1(%rcx,%rsi,1),%ecx
  lea    (%rax,%rsi,1),%edi
  movslq %r13d,%rsi
  sub    %r8,%rsi
  movsd  0x28(%rsp),%xmm2
  jmp    L44
L291:
  movsd  (%rax),%xmm0
  movsd  (%rsp),%xmm3
  movapd %xmm3,%xmm1
  mulsd  0x8(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  (%rax,%rsi,8),%xmm0
  subsd  %xmm3,%xmm0
  movsd  %xmm0,(%rax,%rsi,8)
  add    $0x8,%rdx
  add    $0x1,%ecx
  cmp    %edi,%ecx
  je     L688
L44:
  mov    %rdx,%rax
  movsd  0x10(%rsp),%xmm1
  mulsd  (%rdx,%rsi,8),%xmm1
  movsd  0x18(%rsp),%xmm0
  mulsd  (%rdx),%xmm0
  movapd %xmm1,%xmm5
  addsd  %xmm0,%xmm5
  movsd  %xmm5,(%rsp)
  cmp    %ebp,%r14d
  je     L291
  movslq %ecx,%r8
  lea    (%r15,%r8,8),%r8
  movsd  (%r8),%xmm0
  movsd  0x30(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm5,%xmm1
  movsd  %xmm1,(%rsp)
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%r8)
  jmp    L291
L513:
  movsd  (%rsp),%xmm1
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  addsd  (%rsp),%xmm1
  mov    0x8c(%rsp),%eax
  add    %eax,%eax
  movslq %eax,%rdx
  movsd  0x10(%rsp),%xmm0
  addsd  %xmm1,%xmm0
  mov    0xe0(%rsp),%rbx
  movsd  %xmm0,(%rbx,%rdx,8)
  lea    0x10(%rbx,%rdx,8),%rdx
  ucomisd D22(%rip),%xmm1
  jp     L201
  jne    L201
  movsd  %xmm0,(%rdx)
L21:
  lea    0x1(%rax),%edx
  movslq %edx,%rdx
  mov    0xe0(%rsp),%rbx
  movq   $0x0,(%rbx,%rdx,8)
  add    $0x3,%eax
  cltq
  movq   $0x0,(%rbx,%rax,8)
  mov    0x8c(%rsp),%r13d
  mov    %r14d,0x64(%rsp)
  jmp    L509
L201:
  movsd  0x18(%rsp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  0x10(%rsp),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rdx)
  jmp    L21
L596:
  mov    0x80(%rsp),%rax
  movsd  (%rax),%xmm4
  movsd  %xmm4,0x10(%rsp)
  mov    %ebx,0x64(%rsp)
L341:
  lea    (%r14,%r14,1),%eax
  movslq %eax,%rdx
  movsd  0x10(%rsp),%xmm0
  addsd  0xc8(%rsp),%xmm0
  mov    0xe0(%rsp),%rbx
  movsd  %xmm0,(%rbx,%rdx,8)
  add    $0x1,%eax
  cltq
  movq   $0x0,(%rbx,%rax,8)
  mov    %r14d,%r13d
L509:
  cmp    %r13d,0x88(%rsp)
  jg     L24
L121:
  lea    -0x1(%r13),%r14d
  lea    -0x2(%r13),%ebx
  mov    %ebx,0x8c(%rsp)
  mov    0x24(%rsp),%edi
  mov    %edi,%ecx
  imul   %r14d,%ecx
  lea    (%rcx,%r14,1),%eax
  cltq
  lea    (%r15,%rax,8),%rax
  mov    %rax,0x80(%rsp)
  mov    %ebx,%eax
  imul   %edi,%eax
  lea    (%rbx,%rax,1),%edx
  movslq %edx,%rdx
  lea    (%r15,%rdx,8),%rsi
  mov    %rsi,0x90(%rsp)
  lea    (%r14,%rax,1),%edx
  movslq %edx,%rdx
  lea    (%r15,%rdx,8),%rsi
  mov    %rsi,0x78(%rsp)
  lea    (%rbx,%rcx,1),%edx
  movslq %edx,%rdx
  lea    (%r15,%rdx,8),%rsi
  mov    %rsi,0x98(%rsp)
  mov    %r13d,0x70(%rsp)
  lea    -0x3(%r13),%edx
  imul   %edi,%edx
  lea    (%rbx,%rdx,1),%esi
  movslq %esi,%rsi
  lea    (%r15,%rsi,8),%rsi
  mov    %rsi,0xd0(%rsp)
  cltq
  movslq %r13d,%rsi
  add    %rax,%rsi
  lea    (%r15,%rsi,8),%rsi
  mov    %rsi,0xa0(%rsp)
  movslq %ecx,%rcx
  sub    %rax,%rcx
  mov    %rcx,0xa8(%rsp)
  movslq %ebx,%rcx
  movslq %edx,%rdx
  add    %rcx,%rdx
  lea    (%r15,%rdx,8),%rbx
  mov    %rbx,0xb0(%rsp)
  add    %rcx,%rax
  lea    (%r15,%rax,8),%rax
  mov    %rax,0xb8(%rsp)
  lea    -0x4(%r13),%eax
  imul   %edi,%eax
  cltq
  lea    -0x2(%rcx,%rax,1),%rax
  lea    (%r15,%rax,8),%rax
  mov    %rax,0xc0(%rsp)
  movl   $0x0,0x74(%rsp)
  mov    %r13d,%ebx
L243:
  mov    0x88(%rsp),%r9d
  cmp    %ebx,%r9d
  jg     L381
  je     L596
  mov    0xa0(%rsp),%rax
  mov    %ebx,%edx
  pxor   %xmm3,%xmm3
  movapd %xmm3,%xmm4
  movsd  D38(%rip),%xmm2
  mov    0x70(%rsp),%edi
  lea    0x1(%rdi,%r9,1),%esi
  lea    -0x1(%rbx),%ecx
  sub    %edi,%ecx
  mov    0x68(%rsp),%r10
  mov    0xa8(%rsp),%r8
L317:
  mov    %rax,%rdi
  movsd  -0x10(%rax),%xmm0
  andpd  D149(%rip),%xmm0
  movsd  -0x8(%rax,%r8,8),%xmm1
  andpd  D149(%rip),%xmm1
  addsd  %xmm1,%xmm0
  ucomisd %xmm3,%xmm0
  jp     L177
  ucomisd %xmm4,%xmm0
  movq   %xmm0,%r11
  cmove  %r10,%r11
  movq   %r11,%xmm0
L177:
  movsd  -0x8(%rdi),%xmm1
  andpd  D149(%rip),%xmm1
  mulsd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jae    L148
  mov    %esi,%edi
  sub    %edx,%edi
  cmp    %edi,%ebx
  jl     L299
  lea    (%rcx,%rdx,1),%edi
  mov    %edi,%edx
  add    0x58(%rsp),%rax
  cmp    %edi,%r9d
  jne    L317
  mov    %edi,0x64(%rsp)
  jmp    L381
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D149
type: qword[2]
name: extract_symbolic_data
arguments: data_label: D38
type: byte[8]
name: extract_symbolic_data
arguments: data_label: D100
type: f32[2]
name: extract_symbolic_data
arguments: data_label: D34
type: f32[4]
name: extract_symbolic_data
arguments: data_label: D120
type: f64
name: extract_symbolic_data
arguments: data_label: D165
type: f32[2]
name: extract_symbolic_data
arguments: data_label: D22
type: f32[2]

tool: {"D149": [0x7fffffffffffffff,0x0000000000000000]}

tool: {"D38": [0xbb,0xbd,0xd7,0xd9,0xdf,0x7c,0xdb,0x3d]}

tool: {"D100": [0.0,1.75]}

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

tool: D120: 0.75

tool: {"D165": [0.0,-1.71875]}

tool: {"D22": [0.0,0.0]}

assistant: ```c
void hqrx(n, low, igh, h, ev, ierr) int n, low, igh, *ierr;
double *h, *ev;
{
  int i, j, k, l = 0, m = 0, en, ll, mm, na, its, mp2, enm2;
  double p = 0.0, q = 0.0, r = 0.0, s, t, w, x, y, zz, norm, machep = 1.e-10;
  int notlas;
  *ierr = 0;
  norm = 0.0;
  k = 1;
  for (i = 1; i <= n; i++) {
    for (j = k; j <= n; j++)
      norm = norm + fabs(h[i - 1 + (j - 1) * n]);
    k = i;
    if ((i >= low) && (i <= igh))
      continue;
    ev[(i - 1) * 2] = h[i - 1 + (i - 1) * n];
    ev[1 + (i - 1) * 2] = 0.0;
  }
  en = igh;
  t = 0.0;
l60:
  if (en < low)
    return;
  its = 0;
  na = en - 1;
  enm2 = na - 1;
l70:
  for (ll = low; ll <= en; ll++) {
    l = en + low - ll;
    if (l == low)
      break;
    s = fabs(h[l - 2 + (l - 2) * n]) + fabs(h[l - 1 + (l - 1) * n]);
    if (s == 0.0)
      s = norm;
    if (fabs(h[l - 1 + (l - 2) * n]) <= machep * s)
      break;
  }
  x = h[en - 1 + (en - 1) * n];
  if (l == en)
    goto l270;
  y = h[na - 1 + (na - 1) * n];
  w = h[en - 1 + (na - 1) * n] * h[na - 1 + (en - 1) * n];
  if (l == na)
    goto l280;
  if (its == 30)
    goto l1000;
  if ((its != 10) && (its != 20))
    goto l130;
  t = t + x;
  for (i = low; i <= en; i++)
    h[i - 1 + (i - 1) * n] = h[i - 1 + (i - 1) * n] - x;
  s = fabs(h[en - 1 + (na - 1) * n]) + fabs(h[na - 1 + (enm2 - 1) * n]);
  x = 0.75 * s;
  y = x;
  w = -0.4375 * s * s;
l130:
  its++;
  for (mm = l; mm <= enm2; mm++) {
    m = enm2 + l - mm;
    zz = h[m - 1 + (m - 1) * n];
    r = x - zz;
    s = y - zz;
    p = (r * s - w) / h[m + (m - 1) * n] + h[m - 1 + m * n];
    q = h[m + m * n] - zz - r - s;
    r = h[m + 1 + m * n];
    s = fabs(p) + fabs(q) + fabs(r);
    p = p / s;
    q = q / s;
    r = r / s;
    if (m == l)
      break;
    if ((fabs(h[m - 1 + (m - 2) * n]) * (fabs(q) + fabs(r))) <= (machep * fabs(p) * (fabs(h[m - 2 + (m - 2) * n]) + fabs(zz) + fabs(h[m + m * n]))))
      break;
  }
  mp2 = m + 2;
  for (i = mp2; i <= en; i++) {
    h[i - 1 + (i - 3) * n] = 0.0;
    if (i == mp2)
      continue;
    h[i - 1 + (i - 4) * n] = 0.0;
  }
  for (k = m; k <= na; k++) {
    notlas = 0;
    if (k != na)
      notlas = 1;
    if (k == m)
      goto l170;
    p = h[k - 1 + (k - 2) * n];
    q = h[k + (k - 2) * n];
    r = 0.0;
    if (notlas)
      r = h[k + 1 + (k - 2) * n];
    x = fabs(p) + fabs(q) + fabs(r);
    if (x == 0.0)
      continue;
    p = p / x;
    q = q / x;
    r = r / x;
  l170:
    s = sign(sqrt(p * p + q * q + r * r), p);
    if (k != m)
      h[k - 1 + (k - 2) * n] = -s * x;
    else if (l != m)
      h[k - 1 + (k - 2) * n] = -h[k - 1 + (k - 2) * n];
    p = p + s;
    x = p / s;
    y = q / s;
    zz = r / s;
    q = q / p;
    r = r / p;
    for (j = k; j <= en; j++) {
      p = h[k - 1 + (j - 1) * n] + q * h[k + (j - 1) * n];
      if (notlas) {
        p = p + r * h[k + 1 + (j - 1) * n];
        h[k + 1 + (j - 1) * n] = h[k + 1 + (j - 1) * n] - p * zz;
      }
      h[k + (j - 1) * n] = h[k + (j - 1) * n] - p * y;
      h[k - 1 + (j - 1) * n] = h[k - 1 + (j - 1) * n] - p * x;
    }
    j = imin(en, k + 3);
    for (i = l; i <= j; i++) {
      p = x * h[i - 1 + (k - 1) * n] + y * h[i - 1 + k * n];
      if (notlas) {
        p = p + zz * h[i - 1 + (k + 1) * n];
        h[i - 1 + (k + 1) * n] = h[i - 1 + (k + 1) * n] - p * r;
      }
      h[i - 1 + k * n] = h[i - 1 + k * n] - p * q;
      h[i - 1 + (k - 1) * n] = h[i - 1 + (k - 1) * n] - p;
    }
  }
  goto l70;
l270:
  ev[(en - 1) * 2] = x + t;
  ev[1 + (en - 1) * 2] = 0.0;
  en = na;
  goto l60;
l280:
  p = (y - x) / 2.0;
  q = p * p + w;
  zz = sqrt(fabs(q));
  x = x + t;
  if (q < 0.0)
    goto l320;
  zz = p + sign(zz, p);
  ev[(na - 1) * 2] = x + zz;
  ev[(en - 1) * 2] = ev[(na - 1) * 2];
  if (zz != 0.0)
    ev[(en - 1) * 2] = x - w / zz;
  ev[1 + (na - 1) * 2] = 0.0;
  ev[1 + (en - 1) * 2] = 0.0;
  goto l330;
l320:
  ev[(na - 1) * 2] = x + p;
  ev[(en - 1) * 2] = x + p;
  ev[1 + (na - 1) * 2] = zz;
  ev[1 + (en - 1) * 2] = -zz;
l330:
  en = enm2;
  goto l60;
l1000:
  *ierr = en;
}
```