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

```asm
<hqrx>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xc0,%rsp
  mov    %edi,-0x94(%rbp)
  mov    %esi,-0x98(%rbp)
  mov    %edx,-0x9c(%rbp)
  mov    %rcx,-0xa8(%rbp)
  mov    %r8,-0xb0(%rbp)
  mov    %r9,-0xb8(%rbp)
  movl   $0x0,-0x80(%rbp)
  movl   $0x0,-0x7c(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  D96(%rip),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0xb8(%rbp),%rax
  movl   $0x0,(%rax)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x1,-0x84(%rbp)
  movl   $0x1,-0x8c(%rbp)
  jmp    L162
L348:
  mov    -0x84(%rbp),%eax
  mov    %eax,-0x88(%rbp)
  jmp    L190
L328:
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  addl   $0x1,-0x88(%rbp)
L190:
  mov    -0x88(%rbp),%eax
  cmp    -0x94(%rbp),%eax
  jle    L328
  mov    -0x8c(%rbp),%eax
  mov    %eax,-0x84(%rbp)
  mov    -0x8c(%rbp),%eax
  cmp    -0x98(%rbp),%eax
  jl     L275
  mov    -0x8c(%rbp),%eax
  cmp    -0x9c(%rbp),%eax
  jle    L141
L275:
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x8c(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rax,%rdx
  mov    -0x8c(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0xb0(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x8c(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  add    $0x1,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L345
L141:
  nop
L345:
  addl   $0x1,-0x8c(%rbp)
L162:
  mov    -0x8c(%rbp),%eax
  cmp    -0x94(%rbp),%eax
  jle    L348
  mov    -0x9c(%rbp),%eax
  mov    %eax,-0x78(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x38(%rbp)
L388:
  mov    -0x78(%rbp),%eax
  cmp    -0x98(%rbp),%eax
  jl     L120
  movl   $0x0,-0x6c(%rbp)
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x64(%rbp)
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x60(%rbp)
L319:
  mov    -0x98(%rbp),%eax
  mov    %eax,-0x74(%rbp)
  jmp    L248
L27:
  mov    -0x78(%rbp),%edx
  mov    -0x98(%rbp),%eax
  add    %edx,%eax
  sub    -0x74(%rbp),%eax
  mov    %eax,-0x80(%rbp)
  mov    -0x80(%rbp),%eax
  cmp    -0x98(%rbp),%eax
  je     L553
  mov    -0x80(%rbp),%eax
  lea    -0x2(%rax),%edx
  mov    -0x80(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm0,%xmm1
  mov    -0x80(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x80(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x40(%rbp),%xmm0
  jp     L262
  pxor   %xmm0,%xmm0
  ucomisd -0x40(%rbp),%xmm0
  jne    L262
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
L262:
  mov    -0x80(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x80(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jae    L584
  addl   $0x1,-0x74(%rbp)
L248:
  mov    -0x74(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jle    L27
  jmp    L66
L553:
  nop
  jmp    L66
L584:
  nop
L66:
  mov    -0x78(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0x80(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  je     L563
  mov    -0x64(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x78(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x64(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x80(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  je     L340
  cmpl   $0x1e,-0x6c(%rbp)
  je     L407
  cmpl   $0xa,-0x6c(%rbp)
  je     L249
  cmpl   $0x14,-0x6c(%rbp)
  jne    L456
L249:
  movsd  -0x38(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x98(%rbp),%eax
  mov    %eax,-0x8c(%rbp)
  jmp    L246
L130:
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x8c(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x8c(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  subsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x8c(%rbp)
L246:
  mov    -0x8c(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jle    L130
  mov    -0x78(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm0,%xmm1
  mov    -0x64(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x60(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D32(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D243(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x40(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  jmp    L217
L456:
  nop
L217:
  addl   $0x1,-0x6c(%rbp)
  mov    -0x80(%rbp),%eax
  mov    %eax,-0x70(%rbp)
  jmp    L513
L18:
  mov    -0x60(%rbp),%edx
  mov    -0x80(%rbp),%eax
  add    %edx,%eax
  sub    -0x70(%rbp),%eax
  mov    %eax,-0x7c(%rbp)
  mov    -0x7c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x7c(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x28(%rbp),%xmm0
  subsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x48(%rbp),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  subsd  -0x30(%rbp),%xmm0
  mov    -0x7c(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  mov    %eax,%edx
  mov    -0x7c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  mov    -0x7c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x7c(%rbp),%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  mov    -0x7c(%rbp),%eax
  imul   -0x94(%rbp),%eax
  mov    %eax,%edx
  mov    -0x7c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  subsd  -0x8(%rbp),%xmm0
  subsd  -0x48(%rbp),%xmm0
  subsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x7c(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x7c(%rbp),%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x58(%rbp),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x58(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x50(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x48(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x7c(%rbp),%eax
  cmp    -0x80(%rbp),%eax
  je     L236
  mov    -0x7c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x7c(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm0,%xmm2
  movsd  -0x48(%rbp),%xmm0
  movq   D136(%rip),%xmm3
  andpd  %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x58(%rbp),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm2,%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  mov    -0x7c(%rbp),%eax
  lea    -0x2(%rax),%edx
  mov    -0x7c(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm3
  andpd  %xmm0,%xmm3
  movsd  -0x8(%rbp),%xmm0
  movq   D136(%rip),%xmm4
  andpd  %xmm4,%xmm0
  addsd  %xmm0,%xmm3
  mov    -0x7c(%rbp),%eax
  imul   -0x94(%rbp),%eax
  mov    %eax,%edx
  mov    -0x7c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D136(%rip),%xmm4
  andpd  %xmm4,%xmm0
  addsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jae    L416
  addl   $0x1,-0x70(%rbp)
L513:
  mov    -0x70(%rbp),%eax
  cmp    -0x60(%rbp),%eax
  jle    L18
  jmp    L128
L236:
  nop
  jmp    L128
L416:
  nop
L128:
  mov    -0x7c(%rbp),%eax
  add    $0x2,%eax
  mov    %eax,-0x5c(%rbp)
  mov    -0x5c(%rbp),%eax
  mov    %eax,-0x8c(%rbp)
  jmp    L278
L472:
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x8c(%rbp),%eax
  sub    $0x3,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x8c(%rbp),%eax
  cmp    -0x5c(%rbp),%eax
  je     L271
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x8c(%rbp),%eax
  sub    $0x4,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L237
L271:
  nop
L237:
  addl   $0x1,-0x8c(%rbp)
L278:
  mov    -0x8c(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jle    L472
  mov    -0x7c(%rbp),%eax
  mov    %eax,-0x84(%rbp)
  jmp    L321
L569:
  movl   $0x0,-0x68(%rbp)
  mov    -0x84(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  je     L115
  movl   $0x1,-0x68(%rbp)
L115:
  mov    -0x84(%rbp),%eax
  cmp    -0x7c(%rbp),%eax
  je     L43
  mov    -0x84(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  mov    -0x84(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  mov    %eax,%edx
  mov    -0x84(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  cmpl   $0x0,-0x68(%rbp)
  je     L139
  mov    -0x84(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x48(%rbp)
L139:
  movsd  -0x58(%rbp),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm0
  movq   D136(%rip),%xmm2
  andpd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x28(%rbp),%xmm0
  jp     L119
  pxor   %xmm0,%xmm0
  ucomisd -0x28(%rbp),%xmm0
  je     L137
L119:
  movsd  -0x58(%rbp),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x50(%rbp),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x48(%rbp),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  jmp    L50
L43:
  nop
L50:
  movsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  movsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x84(%rbp),%eax
  cmp    -0x7c(%rbp),%eax
  je     L558
  movsd  -0x40(%rbp),%xmm0
  movq   D102(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0x84(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  mulsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L154
L558:
  mov    -0x80(%rbp),%eax
  cmp    -0x7c(%rbp),%eax
  je     L154
  mov    -0x84(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x84(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x2,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movq   D102(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L154:
  movsd  -0x58(%rbp),%xmm0
  addsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x58(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x50(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x48(%rbp),%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x50(%rbp),%xmm0
  divsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x48(%rbp),%xmm0
  divsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x84(%rbp),%eax
  mov    %eax,-0x88(%rbp)
  jmp    L19
L337:
  mov    -0x84(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  mov    %eax,%edx
  mov    -0x84(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x50(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  cmpl   $0x0,-0x68(%rbp)
  je     L216
  mov    -0x84(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x48(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  mov    -0x84(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  mulsd  -0x8(%rbp),%xmm1
  mov    -0x84(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L216:
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  mov    %eax,%edx
  mov    -0x84(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  mulsd  -0x20(%rbp),%xmm1
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  mov    %eax,%edx
  mov    -0x84(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x84(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  mulsd  -0x28(%rbp),%xmm1
  mov    -0x84(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x88(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x88(%rbp)
L19:
  mov    -0x88(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jle    L337
  mov    -0x84(%rbp),%eax
  lea    0x3(%rax),%edx
  mov    -0x78(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  mov    $0x0,%eax
  call   <imin@plt>
  mov    %eax,-0x88(%rbp)
  mov    -0x80(%rbp),%eax
  mov    %eax,-0x8c(%rbp)
  jmp    L335
L132:
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x28(%rbp),%xmm1
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  cmpl   $0x0,-0x68(%rbp)
  je     L11
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  add    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  add    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  mulsd  -0x48(%rbp),%xmm1
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  add    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L11:
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  mulsd  -0x50(%rbp),%xmm1
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x84(%rbp),%eax
  sub    $0x1,%eax
  imul   -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  subsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x8c(%rbp)
L335:
  mov    -0x8c(%rbp),%eax
  cmp    -0x88(%rbp),%eax
  jle    L132
  jmp    L125
L137:
  nop
L125:
  addl   $0x1,-0x84(%rbp)
L321:
  mov    -0x84(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jle    L569
  jmp    L319
L563:
  nop
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  add    $0x1,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x64(%rbp),%eax
  mov    %eax,-0x78(%rbp)
  jmp    L388
L340:
  nop
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x28(%rbp),%xmm0
  movsd  D239(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x58(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0x30(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movq   D136(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x8(%rbp)
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x50(%rbp),%xmm0
  ja     L264
  movsd  -0x58(%rbp),%xmm0
  mov    -0x8(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rax,%rdx
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0xb0(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  pxor   %xmm0,%xmm0
  ucomisd -0x8(%rbp),%xmm0
  jp     L242
  pxor   %xmm0,%xmm0
  ucomisd -0x8(%rbp),%xmm0
  je     L290
L242:
  movsd  -0x30(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x8(%rbp),%xmm1
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L290:
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  add    $0x1,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  add    $0x1,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L408
L264:
  nop
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x64(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  add    $0x1,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x78(%rbp),%eax
  sub    $0x1,%eax
  add    %eax,%eax
  add    $0x1,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x8(%rbp),%xmm0
  movq   D102(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L408:
  mov    -0x60(%rbp),%eax
  mov    %eax,-0x78(%rbp)
  jmp    L388
L407:
  nop
  mov    -0xb8(%rbp),%rax
  mov    -0x78(%rbp),%edx
  mov    %edx,(%rax)
  jmp    L33
L120:
  nop
L33:
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D96
type_guessed: byte[8]
name: extract_symbolic_data
arguments: label: D136
type_guessed: qword
name: extract_symbolic_data
arguments: label: D32
type_guessed: f64
name: extract_symbolic_data
arguments: label: D243
type_guessed: f32[2]
name: extract_symbolic_data
arguments: label: D102
type_guessed: f32[2]
name: extract_symbolic_data
arguments: label: D239
type_guessed: f64

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

tool: D136: 9223372036854775807

tool: D32: 0.75

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

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

tool: D239: 2

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;
}
```