user: Converte this assembly snippet into its equivalent C code:

```asm
<qrfac>:
  endbr64
  push   %r15
  mov    %rdx,%r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x118,%rsp
  mov    %edi,0x24(%rsp)
  mov    0x160(%rsp),%rbx
  mov    %esi,0x58(%rsp)
  mov    %ecx,0x20(%rsp)
  mov    %r8d,0x40(%rsp)
  mov    %r9,0x88(%rsp)
  test   %esi,%esi
  jle    L545
  movslq 0x20(%rsp),%rax
  mov    %rdx,%r14
  lea    0x0(,%rax,8),%r12
  mov    0x40(%rsp),%eax
  test   %eax,%eax
  jne    L272
  xor    %r13d,%r13d
  mov    %rdx,0x8(%rsp)
  movslq 0x58(%rsp),%rbp
  mov    %rbx,0x160(%rsp)
  mov    0x24(%rsp),%r15d
  mov    %r13,%rbx
  mov    0x160(%rsp),%r13
  nopw   0x0(%rax,%rax,1)
L404:
  mov    %r14,%rsi
  mov    %r15d,%edi
  add    %r12,%r14
  call   <enorm@plt>
  mov    0x158(%rsp),%rax
  movsd  %xmm0,0x0(%r13,%rbx,8)
  movsd  %xmm0,(%rax,%rbx,8)
  mov    0x168(%rsp),%rax
  movsd  %xmm0,(%rax,%rbx,8)
  add    $0x1,%rbx
  cmp    %rbp,%rbx
  jne    L404
  mov    0x8(%rsp),%r15
L545:
  mov    0x58(%rsp),%ecx
  mov    0x24(%rsp),%eax
  cmp    %eax,%ecx
  cmovle %ecx,%eax
  mov    %eax,0x5c(%rsp)
  test   %eax,%eax
  jle    L462
  mov    0x24(%rsp),%eax
  test   %eax,%eax
  jle    L384
  mov    0x5c(%rsp),%ecx
  pxor   %xmm5,%xmm5
  mov    %r15,0x60(%rsp)
  mov    %r15,%r12
  movq   D41(%rip),%xmm4
  movsd  %xmm5,0x8(%rsp)
  movq   $0x0,0x68(%rsp)
  cmp    %ecx,%eax
  cmovle %eax,%ecx
  xor    %r10d,%r10d
  movaps %xmm4,0x30(%rsp)
  movq   $0x1,0x18(%rsp)
  mov    %r10,%rbp
  mov    %ecx,0x70(%rsp)
  movslq 0x20(%rsp),%rcx
  lea    0x8(,%rcx,8),%rsi
  mov    %rcx,0x28(%rsp)
  mov    %rcx,%r14
  shl    $0x3,%rcx
  mov    %rcx,0x50(%rsp)
  mov    %eax,%ecx
  shr    %eax
  sub    $0x1,%eax
  mov    %ecx,%r9d
  mov    %rsi,0x78(%rsp)
  add    $0x1,%rax
  and    $0xfffffffe,%r9d
  shl    $0x4,%rax
  mov    %r9d,0xd0(%rsp)
  mov    %ecx,%r9d
  mov    %rax,0xc0(%rsp)
  lea    -0x1(%rcx),%eax
  mov    %eax,0xbc(%rsp)
  jmp    L140
  nopl   0x0(%rax)
L554:
  mov    0x18(%rsp),%eax
  mov    %r9d,%edi
  mov    %r12,%rsi
  mov    %r9d,0x44(%rsp)
  mov    %eax,0x80(%rsp)
  lea    -0x1(%r9),%eax
  mov    %eax,0x48(%rsp)
  call   <enorm@plt>
  ucomisd 0x8(%rsp),%xmm0
  mov    0x44(%rsp),%r9d
  jp     L446
  jne    L446
  movapd %xmm0,%xmm3
  xorpd  0x30(%rsp),%xmm3
L587:
  mov    0x18(%rsp),%rax
  mov    0x158(%rsp),%rcx
  add    $0x1,%rbp
  mov    0x50(%rsp),%rsi
  add    0x78(%rsp),%r12
  movsd  %xmm3,-0x8(%rcx,%rax,8)
  add    $0x1,%rax
  mov    0x20(%rsp),%ecx
  mov    %rax,0x18(%rsp)
  mov    0x10(%rsp),%eax
  add    %rsi,0x60(%rsp)
  mov    0x28(%rsp),%rsi
  lea    (%rax,%rcx,2),%r14d
  add    %rsi,0x68(%rsp)
  cmp    %ebp,0x70(%rsp)
  jle    L202
  mov    0x48(%rsp),%r9d
L140:
  mov    %r14d,%eax
  sub    0x20(%rsp),%eax
  mov    %ebp,0x84(%rsp)
  mov    %ebp,%r13d
  mov    %eax,0x10(%rsp)
  mov    0x40(%rsp),%eax
  test   %eax,%eax
  je     L554
  mov    0x58(%rsp),%r11d
  mov    %ebp,%r10d
  cmp    %ebp,%r11d
  jle    L554
  mov    0x158(%rsp),%rax
  mov    0x18(%rsp),%rcx
  mov    %rbp,%rdi
  mov    0x158(%rsp),%rbx
  movsd  -0x8(%rax,%rcx,8),%xmm1
  lea    (%rax,%rbp,8),%r8
  mov    %ebp,%eax
  jmp    L341
L801:
  movslq %eax,%rsi
  add    $0x1,%rdi
  add    $0x8,%r8
  lea    0x0(,%rsi,8),%rdx
  lea    (%rbx,%rdx,1),%rcx
  cmp    %edi,%r11d
  jle    L240
L341:
  movsd  (%r8),%xmm0
  lea    0x0(,%rdi,8),%rdx
  mov    %r8,%rcx
  comisd %xmm1,%xmm0
  jbe    L801
  mov    %rdi,%rsi
  mov    %edi,%eax
  add    $0x1,%rdi
  add    $0x8,%r8
  movapd %xmm0,%xmm1
  cmp    %edi,%r11d
  jg     L341
L240:
  cmp    %r10d,%eax
  je     L554
  imul   0x20(%rsp),%eax
  mov    0x68(%rsp),%rdi
  lea    0x0(,%rdi,8),%r10
  movslq %eax,%r8
  lea    0x10(%r10),%r11
  lea    0x0(,%r8,8),%rdi
  lea    0x10(,%r8,8),%rbx
  cmp    %r11,%rdi
  setge  %r11b
  cmp    %r10,%rbx
  setle  %r10b
  or     %r11b,%r10b
  je     L845
  cmpl   $0x1,0xbc(%rsp)
  jbe    L845
  mov    0xc0(%rsp),%r11
  mov    0x60(%rsp),%r10
  add    %r15,%rdi
  xor    %r8d,%r8d
  nopw   0x0(%rax,%rax,1)
L819:
  movupd (%r10,%r8,1),%xmm0
  movupd (%rdi,%r8,1),%xmm6
  movups %xmm6,(%r10,%r8,1)
  movups %xmm0,(%rdi,%r8,1)
  add    $0x10,%r8
  cmp    %r11,%r8
  jne    L819
  mov    0xd0(%rsp),%ebx
  mov    0x24(%rsp),%edi
  cmp    %edi,%ebx
  je     L242
  mov    0x10(%rsp),%edi
  add    %ebx,%eax
  cltq
  add    %ebx,%edi
  lea    (%r15,%rax,8),%rax
  movslq %edi,%rdi
  movsd  (%rax),%xmm1
  lea    (%r15,%rdi,8),%rdi
  movsd  (%rdi),%xmm0
  movsd  %xmm1,(%rdi)
  movsd  %xmm0,(%rax)
L242:
  mov    0x158(%rsp),%rax
  mov    0x18(%rsp),%rbx
  movsd  -0x8(%rax,%rbx,8),%xmm0
  mov    0x168(%rsp),%rax
  mov    0x88(%rsp),%rbx
  movsd  %xmm0,(%rcx)
  movsd  (%rax,%rbp,8),%xmm0
  movsd  %xmm0,(%rax,%rdx,1)
  lea    (%rbx,%rsi,4),%rdx
  mov    (%rbx,%rbp,4),%eax
  mov    (%rdx),%ecx
  mov    %ecx,(%rbx,%rbp,4)
  mov    %eax,(%rdx)
  jmp    L554
  nop
L446:
  movsd  (%r12),%xmm1
  movsd  0x8(%rsp),%xmm7
  movapd 0x30(%rsp),%xmm6
  comisd %xmm1,%xmm7
  xorpd  %xmm0,%xmm6
  movapd %xmm6,%xmm3
  jbe    L646
  movq   %xmm0,%rax
  movapd %xmm6,%xmm0
  movq   %rax,%xmm3
L646:
  mov    %ebp,0xb8(%rsp)
  cmp    %ebp,0x24(%rsp)
  jle    L574
  cmp    $0x1,%r9d
  je     L715
  mov    %r9d,%eax
  movapd %xmm0,%xmm2
  shr    %eax
  unpcklpd %xmm2,%xmm2
  sub    $0x1,%eax
  shl    $0x4,%rax
  lea    0x10(%r12,%rax,1),%rdx
  mov    %r12,%rax
  nopl   0x0(%rax)
L663:
  movupd (%rax),%xmm1
  add    $0x10,%rax
  divpd  %xmm2,%xmm1
  movups %xmm1,-0x10(%rax)
  cmp    %rdx,%rax
  jne    L663
  mov    %r9d,%edx
  mov    %r13d,%eax
  and    $0xfffffffe,%edx
  add    %edx,%eax
  cmp    %r9d,%edx
  je     L467
L359:
  add    0x10(%rsp),%eax
  cltq
  lea    (%r15,%rax,8),%rax
  movsd  (%rax),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax)
L467:
  movsd  (%r12),%xmm1
L574:
  addsd  D62(%rip),%xmm1
  mov    0x18(%rsp),%ecx
  movsd  %xmm1,(%r12)
  cmp    %ecx,0x58(%rsp)
  jle    L587
  mov    0x18(%rsp),%rax
  mov    0x58(%rsp),%esi
  mov    %r13d,%r11d
  mov    %rbp,0xb0(%rsp)
  mov    0x158(%rsp),%rcx
  mov    0x84(%rsp),%edi
  shl    $0x3,%rax
  add    %rax,%rcx
  add    0x168(%rsp),%rax
  mov    %rax,0xc8(%rsp)
  movslq %r14d,%rax
  lea    (%rax,%rbp,1),%rdx
  lea    (%r15,%rax,8),%r8
  mov    %rcx,0xa8(%rsp)
  lea    -0x2(%rsi),%eax
  mov    0x158(%rsp),%rsi
  lea    (%r15,%rdx,8),%rbx
  sub    %edi,%eax
  add    %rbp,%rax
  lea    0x10(%rsi,%rax,8),%rdx
  mov    $0x1,%esi
  mov    %rdx,0x98(%rsp)
  mov    0x24(%rsp),%edx
  cmp    %r13d,%edx
  cmovg  %r9d,%esi
  cmp    0xb8(%rsp),%edx
  setle  %dl
  cmp    $0x1,%r9d
  mov    %esi,%ecx
  sete   %al
  mov    %esi,0x44(%rsp)
  and    $0xfffffffe,%esi
  shr    %ecx
  mov    %esi,%r10d
  mov    %esi,0x84(%rsp)
  mov    0x10(%rsp),%esi
  lea    -0x1(%rcx),%r13d
  add    %edi,%r10d
  add    $0x1,%r13
  shl    $0x4,%r13
  or     %al,%dl
  mov    %dl,0x90(%rsp)
  mov    %r11d,%edx
  cmove  %r10d,%edx
  lea    (%rdx,%rsi,1),%eax
  cltq
  lea    (%r15,%rax,8),%rcx
  lea    (%r10,%rsi,1),%eax
  cltq
  mov    %rcx,%rbp
  mov    0xa8(%rsp),%rcx
  mov    %r9d,0xa8(%rsp)
  lea    (%r15,%rax,8),%r11
  mov    %r12,%rax
  mov    0xc8(%rsp),%r9
  sub    %r15,%rax
  lea    0x8(%rax),%rsi
  mov    %rsi,0xa0(%rsp)
  lea    0x1(%rdi),%esi
  mov    %esi,0xd4(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L885:
  xor    %eax,%eax
  cmpb   $0x0,0x90(%rsp)
  pxor   %xmm0,%xmm0
  jne    L449
  nopl   0x0(%rax)
L807:
  movupd (%r12,%rax,1),%xmm2
  movupd (%rbx,%rax,1),%xmm5
  add    $0x10,%rax
  mulpd  %xmm5,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %rax,%r13
  jne    L807
  mov    0x84(%rsp),%esi
  cmp    %esi,0x44(%rsp)
  je     L244
L498:
  lea    (%r14,%rdx,1),%eax
  cltq
  movsd  (%r15,%rax,8),%xmm1
  mulsd  0x0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
L244:
  mov    %rbx,%rax
  divsd  (%r12),%xmm0
  sub    %r15,%rax
  cmp    %rax,0xa0(%rsp)
  je     L432
  mov    0xb8(%rsp),%esi
  cmp    %esi,0x24(%rsp)
  movapd %xmm0,%xmm4
  setg   %dil
  cmpl   $0x1,0xa8(%rsp)
  unpcklpd %xmm4,%xmm4
  setne  %sil
  xor    %eax,%eax
  test   %sil,%dil
  je     L432
  cs nopw 0x0(%rax,%rax,1)
L573:
  movupd (%r12,%rax,1),%xmm2
  movupd (%rbx,%rax,1),%xmm1
  mulpd  %xmm4,%xmm2
  subpd  %xmm2,%xmm1
  movups %xmm1,(%rbx,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%r13
  jne    L573
  mov    0x84(%rsp),%esi
  cmp    %esi,0x44(%rsp)
  je     L24
  mulsd  (%r11),%xmm0
  lea    (%r14,%r10,1),%eax
  cltq
  lea    (%r15,%rax,8),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax)
L24:
  mov    0x40(%rsp),%eax
  test   %eax,%eax
  je     L378
  movsd  (%rcx),%xmm1
  ucomisd 0x8(%rsp),%xmm1
  jp     L171
  je     L378
L171:
  movsd  (%rbx),%xmm0
  movsd  D62(%rip),%xmm2
  movsd  0x8(%rsp),%xmm7
  divsd  %xmm1,%xmm0
  mulsd  %xmm0,%xmm0
  subsd  %xmm0,%xmm2
  comisd %xmm2,%xmm7
  jae    L736
  ucomisd %xmm2,%xmm7
  ja     L891
L394:
  sqrtsd %xmm2,%xmm2
L760:
  mulsd  %xmm2,%xmm1
  movsd  D80(%rip),%xmm5
  movsd  %xmm1,(%rcx)
  divsd  (%r9),%xmm1
  mulsd  %xmm1,%xmm1
  mulsd  D72(%rip),%xmm1
  comisd %xmm1,%xmm5
  jae    L100
  nop
L378:
  mov    0x50(%rsp),%rax
  add    $0x8,%rcx
  add    $0x8,%r9
  add    0x20(%rsp),%r14d
  add    %rax,%rbx
  add    %rax,%r8
  cmp    0x98(%rsp),%rcx
  jne    L885
  mov    0xb0(%rsp),%rbp
  jmp    L587
  nopl   0x0(%rax,%rax,1)
L202:
  mov    0x5c(%rsp),%esi
  cmp    %esi,0x80(%rsp)
  jl     L206
L462:
  add    $0x118,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L384:
  movslq 0x20(%rsp),%rax
  movq   D41(%rip),%xmm5
  pxor   %xmm7,%xmm7
  movl   $0x0,0x80(%rsp)
  movsd  %xmm7,0x8(%rsp)
  mov    %rax,0x28(%rsp)
  movaps %xmm5,0x30(%rsp)
L206:
  movslq 0x80(%rsp),%rbx
  mov    0x20(%rsp),%esi
  mov    0x28(%rsp),%rax
  mov    0x24(%rsp),%edi
  mov    %ebx,%r12d
  lea    0x1(%rbx),%ebp
  lea    0x1(%rsi),%edx
  mov    %ebx,%r14d
  imul   %esi,%r12d
  lea    0x8(,%rax,8),%rax
  mov    %edx,0x50(%rsp)
  mov    0x158(%rsp),%r11
  mov    %rax,0x48(%rsp)
  imul   %ebp,%edx
  movslq %r12d,%rax
  mov    %rax,0x18(%rsp)
  add    %rbx,%rax
  lea    (%r15,%rax,8),%r13
  mov    %edi,%eax
  and    $0xfffffffe,%edi
  mov    %edx,0x44(%rsp)
  shr    %eax
  mov    %edi,0x60(%rsp)
  sub    $0x1,%eax
  add    $0x1,%rax
  shl    $0x4,%rax
  mov    %rax,0x78(%rsp)
  mov    %esi,%eax
  imul   %ebp,%eax
  add    %ebx,%eax
  sub    %edx,%eax
  mov    %eax,0x68(%rsp)
  jmp    L505
  nopl   (%rax)
L581:
  mov    0x24(%rsp),%eax
  mov    %r11,0x158(%rsp)
  mov    %eax,%esi
  sub    %ebx,%eax
  sub    %ebp,%esi
  mov    %eax,%edi
  mov    %eax,0x10(%rsp)
  mov    %esi,0x80(%rsp)
  mov    %r13,%rsi
  call   <enorm@plt>
  ucomisd 0x8(%rsp),%xmm0
  mov    0x10(%rsp),%edx
  mov    0x158(%rsp),%r11
  jp     L264
  jne    L264
  movapd %xmm0,%xmm2
  xorpd  0x30(%rsp),%xmm2
L70:
  mov    0x50(%rsp),%ecx
  mov    0x28(%rsp),%rsi
  movsd  %xmm2,(%r11,%rbx,8)
  add    $0x1,%rbx
  add    %ecx,0x44(%rsp)
  add    0x48(%rsp),%r13
  add    %rsi,0x18(%rsp)
  add    0x20(%rsp),%r12d
  cmp    %ebp,0x5c(%rsp)
  jle    L462
  mov    %ebp,%r14d
  add    $0x1,%ebp
L505:
  mov    0x40(%rsp),%r8d
  test   %r8d,%r8d
  je     L581
  mov    0x58(%rsp),%r9d
  cmp    %r14d,%r9d
  jle    L581
  movsd  (%r11,%rbx,8),%xmm2
  lea    (%r11,%rbx,8),%r8
  mov    %rbx,%rdi
  mov    %r14d,%eax
  movapd %xmm2,%xmm1
  jmp    L466
  nopw   0x0(%rax,%rax,1)
L734:
  movslq %eax,%rsi
  add    $0x1,%rdi
  add    $0x8,%r8
  lea    0x0(,%rsi,8),%rdx
  lea    (%r11,%rdx,1),%rcx
  cmp    %edi,%r9d
  jle    L596
L466:
  movsd  (%r8),%xmm0
  lea    0x0(,%rdi,8),%rdx
  mov    %r8,%rcx
  comisd %xmm1,%xmm0
  jbe    L734
  mov    %rdi,%rsi
  mov    %edi,%eax
  add    $0x1,%rdi
  add    $0x8,%r8
  movapd %xmm0,%xmm1
  cmp    %edi,%r9d
  jg     L466
L596:
  cmp    %r14d,%eax
  je     L581
  mov    0x24(%rsp),%edi
  test   %edi,%edi
  jle    L726
  imul   0x20(%rsp),%eax
  mov    0x24(%rsp),%edi
  mov    0x18(%rsp),%r8
  sub    $0x1,%edi
  movslq %eax,%r9
  mov    %edi,0x80(%rsp)
  lea    0x0(,%r8,8),%rdi
  lea    0x0(,%r9,8),%r8
  lea    0x10(,%r9,8),%r10
  mov    %r8,0x10(%rsp)
  lea    0x10(%rdi),%r8
  cmp    0x10(%rsp),%r8
  setle  0x70(%rsp)
  cmp    %r10,%rdi
  movzbl 0x70(%rsp),%r8d
  setge  %r10b
  add    %r15,%rdi
  or     %r8b,%r10b
  je     L127
  cmpl   $0x1,0x80(%rsp)
  jbe    L127
  mov    0x10(%rsp),%r8
  mov    0x78(%rsp),%r10
  xor    %r9d,%r9d
  add    %r15,%r8
  nopw   0x0(%rax,%rax,1)
L31:
  movupd (%rdi,%r9,1),%xmm0
  movupd (%r8,%r9,1),%xmm7
  movups %xmm7,(%rdi,%r9,1)
  movups %xmm0,(%r8,%r9,1)
  add    $0x10,%r9
  cmp    %r10,%r9
  jne    L31
  mov    0x60(%rsp),%r10d
  cmp    %r10d,0x24(%rsp)
  je     L156
  lea    (%r12,%r10,1),%edi
  add    %r10d,%eax
  movslq %edi,%rdi
  cltq
  lea    (%r15,%rdi,8),%rdi
  lea    (%r15,%rax,8),%rax
  movsd  (%rdi),%xmm0
  movsd  (%rax),%xmm1
  movsd  %xmm1,(%rdi)
  movsd  %xmm0,(%rax)
L156:
  movsd  (%r11,%rbx,8),%xmm2
L726:
  mov    0x168(%rsp),%rax
  movsd  %xmm2,(%rcx)
  mov    0x88(%rsp),%rdi
  movsd  (%rax,%rbx,8),%xmm0
  movsd  %xmm0,(%rax,%rdx,1)
  lea    (%rdi,%rsi,4),%rdx
  mov    (%rdi,%rbx,4),%eax
  mov    (%rdx),%ecx
  mov    %ecx,(%rdi,%rbx,4)
  mov    %eax,(%rdx)
  jmp    L581
  cs nopw 0x0(%rax,%rax,1)
L264:
  movsd  0x0(%r13),%xmm1
  movsd  0x8(%rsp),%xmm4
  movapd 0x30(%rsp),%xmm7
  comisd %xmm1,%xmm4
  xorpd  %xmm0,%xmm7
  movapd %xmm7,%xmm2
  jbe    L153
  movq   %xmm0,%rax
  movapd %xmm7,%xmm0
  movq   %rax,%xmm2
L153:
  cmp    %ebx,0x24(%rsp)
  jle    L62
  cmp    $0x1,%edx
  je     L872
  mov    %edx,%eax
  movapd %xmm0,%xmm3
  shr    %eax
  unpcklpd %xmm3,%xmm3
  sub    $0x1,%eax
  shl    $0x4,%rax
  lea    0x10(%r13,%rax,1),%rcx
  mov    %r13,%rax
  nopl   0x0(%rax)
L429:
  movupd (%rax),%xmm1
  add    $0x10,%rax
  divpd  %xmm3,%xmm1
  movups %xmm1,-0x10(%rax)
  cmp    %rax,%rcx
  jne    L429
  mov    %edx,%eax
  and    $0xfffffffe,%eax
  add    %eax,%r14d
  cmp    %eax,%edx
  je     L237
L872:
  add    %r12d,%r14d
  movslq %r14d,%r14
  lea    (%r15,%r14,8),%rax
  movsd  (%rax),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax)
L237:
  movsd  0x0(%r13),%xmm1
L62:
  addsd  D62(%rip),%xmm1
  movsd  %xmm1,0x0(%r13)
  cmp    %ebp,0x58(%rsp)
  jle    L70
  mov    0x58(%rsp),%ecx
  movslq %ebp,%rax
  mov    %r13,0x70(%rsp)
  shl    $0x3,%rax
  mov    0x20(%rsp),%r13d
  mov    %ebp,0xb8(%rsp)
  lea    -0x1(%rcx),%edx
  lea    (%r11,%rax,1),%r14
  add    0x168(%rsp),%rax
  mov    %rbx,0x10(%rsp)
  sub    %ebp,%edx
  mov    %rax,%rbp
  mov    %r15,%rax
  mov    0x40(%rsp),%r15d
  add    %rbx,%rdx
  mov    %r12d,0x84(%rsp)
  lea    0x10(%r11,%rdx,8),%rcx
  mov    0x44(%rsp),%edx
  mov    %rcx,%r12
  mov    %edx,%ebx
  nopl   0x0(%rax,%rax,1)
L619:
  test   %r15d,%r15d
  je     L854
  movsd  (%r14),%xmm3
  ucomisd 0x8(%rsp),%xmm3
  jp     L876
  je     L854
L876:
  mov    0x68(%rsp),%ecx
  movsd  0x8(%rsp),%xmm6
  lea    (%rcx,%rbx,1),%esi
  movslq %esi,%rsi
  movsd  (%rax,%rsi,8),%xmm0
  divsd  %xmm3,%xmm0
  mulsd  %xmm0,%xmm0
  movapd %xmm0,%xmm1
  movsd  D62(%rip),%xmm0
  subsd  %xmm1,%xmm0
  comisd %xmm0,%xmm6
  jae    L595
  ucomisd %xmm0,%xmm6
  ja     L49
L783:
  sqrtsd %xmm0,%xmm0
L211:
  mulsd  %xmm3,%xmm0
  movsd  D80(%rip),%xmm4
  movsd  %xmm0,(%r14)
  divsd  0x0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  D72(%rip),%xmm0
  comisd %xmm0,%xmm4
  jae    L395
  nopl   0x0(%rax)
L854:
  add    $0x8,%r14
  add    %r13d,%ebx
  add    $0x8,%rbp
  cmp    %r12,%r14
  jne    L619
  mov    0x10(%rsp),%rbx
  mov    0x70(%rsp),%r13
  mov    %rax,%r15
  mov    0x84(%rsp),%r12d
  mov    0xb8(%rsp),%ebp
  jmp    L70
  nopl   (%rax)
L432:
  mov    0xb0(%rsp),%rax
  mov    0x60(%rsp),%rsi
  mov    0x24(%rsp),%edi
  nopl   0x0(%rax)
L103:
  movsd  (%rsi,%rax,8),%xmm2
  movsd  (%r8,%rax,8),%xmm1
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  %xmm1,(%r8,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%edi
  jg     L103
  jmp    L24
  xchg   %ax,%ax
L272:
  mov    0x58(%rsp),%eax
  mov    $0x1,%r13d
  mov    %rdx,0x8(%rsp)
  mov    %rbx,0x160(%rsp)
  mov    %r13,%rbx
  mov    0x160(%rsp),%r13
  lea    0x1(%rax),%ebp
  mov    %rbp,%r15
  mov    0x24(%rsp),%ebp
  nopl   0x0(%rax)
L299:
  mov    %r14,%rsi
  mov    %ebp,%edi
  add    %r12,%r14
  call   <enorm@plt>
  mov    0x158(%rsp),%rax
  movsd  %xmm0,-0x8(%r13,%rbx,8)
  movsd  %xmm0,-0x8(%rax,%rbx,8)
  mov    0x168(%rsp),%rax
  movsd  %xmm0,-0x8(%rax,%rbx,8)
  mov    0x88(%rsp),%rax
  mov    %ebx,-0x4(%rax,%rbx,4)
  add    $0x1,%rbx
  cmp    %r15,%rbx
  jne    L299
  mov    0x8(%rsp),%r15
  jmp    L545
  nop
L595:
  movapd %xmm6,%xmm0
  jmp    L783
  nopl   0x0(%rax)
L736:
  movapd %xmm7,%xmm2
  jmp    L394
  nopl   0x0(%rax)
L449:
  pxor   %xmm0,%xmm0
  jmp    L498
L395:
  mov    0x80(%rsp),%edi
  movslq %ebx,%rsi
  mov    %r11,0x158(%rsp)
  lea    (%rax,%rsi,8),%rsi
  mov    %rax,0x90(%rsp)
  movsd  %xmm2,0x98(%rsp)
  call   <enorm@plt>
  mov    0x158(%rsp),%r11
  movsd  0x98(%rsp),%xmm2
  movsd  %xmm0,(%r14)
  mov    0x90(%rsp),%rax
  movsd  %xmm0,0x0(%rbp)
  jmp    L854
L100:
  mov    0xd4(%rsp),%eax
  mov    0x48(%rsp),%edi
  mov    %r8,0x100(%rsp)
  mov    %r9,0xf8(%rsp)
  add    %r14d,%eax
  mov    %edx,0xf0(%rsp)
  cltq
  mov    %r11,0xe8(%rsp)
  lea    (%r15,%rax,8),%rsi
  mov    %rcx,0xd8(%rsp)
  mov    %r10d,0xc8(%rsp)
  movsd  %xmm3,0xe0(%rsp)
  call   <enorm@plt>
  mov    0xd8(%rsp),%rcx
  mov    0xf8(%rsp),%r9
  mov    0x100(%rsp),%r8
  mov    0xf0(%rsp),%edx
  movsd  %xmm0,(%rcx)
  mov    0xe8(%rsp),%r11
  movsd  0xe0(%rsp),%xmm3
  mov    0xc8(%rsp),%r10d
  movsd  %xmm0,(%r9)
  jmp    L378
L715:
  mov    %r13d,%eax
  jmp    L359
L127:
  mov    0x18(%rsp),%r10
  mov    0x80(%rsp),%eax
  add    %r10,%rax
  sub    %r10,%r9
  lea    0x8(%r15,%rax,8),%rax
  nopl   0x0(%rax,%rax,1)
L636:
  movsd  (%rdi),%xmm0
  movsd  (%rdi,%r9,8),%xmm1
  movsd  %xmm1,(%rdi)
  movsd  %xmm0,(%rdi,%r9,8)
  add    $0x8,%rdi
  cmp    %rdi,%rax
  jne    L636
  jmp    L156
L845:
  mov    0x24(%rsp),%eax
  mov    0x60(%rsp),%rbx
  sub    0x68(%rsp),%r8
  lea    (%rbx,%rax,8),%rdi
  mov    %rbx,%rax
  nop
L889:
  movsd  (%rax),%xmm0
  movsd  (%rax,%r8,8),%xmm1
  movsd  %xmm1,(%rax)
  movsd  %xmm0,(%rax,%r8,8)
  add    $0x8,%rax
  cmp    %rdi,%rax
  jne    L889
  jmp    L242
L891:
  movapd %xmm2,%xmm0
  mov    %r8,0x108(%rsp)
  mov    %r9,0x100(%rsp)
  mov    %edx,0xf8(%rsp)
  mov    %r11,0xf0(%rsp)
  mov    %rcx,0xe0(%rsp)
  mov    %r10d,0xd8(%rsp)
  movsd  %xmm3,0xe8(%rsp)
  movsd  %xmm1,0xc8(%rsp)
  call   <sqrt@plt>
  mov    0x108(%rsp),%r8
  mov    0x100(%rsp),%r9
  mov    0xf8(%rsp),%edx
  mov    0xf0(%rsp),%r11
  movapd %xmm0,%xmm2
  mov    0xe0(%rsp),%rcx
  movsd  0xe8(%rsp),%xmm3
  mov    0xd8(%rsp),%r10d
  movsd  0xc8(%rsp),%xmm1
  jmp    L760
L49:
  mov    %r11,0x158(%rsp)
  mov    %rax,0xa0(%rsp)
  movsd  %xmm2,0x98(%rsp)
  movsd  %xmm3,0x90(%rsp)
  call   <sqrt@plt>
  mov    0xa0(%rsp),%rax
  movsd  0x90(%rsp),%xmm3
  movsd  0x98(%rsp),%xmm2
  mov    0x158(%rsp),%r11
  jmp    L211
```

tool_calls: name: read_data
arguments: tag: D41
type: f32[2]
name: read_data
arguments: tag: D62
type: f64
name: read_data
arguments: tag: D80
type: f32[2]
name: read_data
arguments: tag: D72
type: byte[8]

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

tool: D62: 1

tool: {"D80": [0.0,0.021484375]}

tool: {"D72": [0x9a,0x99,0x99,0x99,0x99,0x99,0xa9,0x3f]}

assistant: ```c
void qrfac(int m, int n, double *a, int lda, int pivot, int *ipvt, int lipvt, double *rdiag, double *acnorm, double *wa) {
  double d1;
  int i, j, k, jp1;
  double sum;
  double temp;
  int minmn;
  double epsmch;
  double ajnorm;
  (void)lipvt;
  epsmch = 2.2204460492503131e-16;
  for (j = 0; j < n; ++j) {
    acnorm[j] = enorm(m, &a[j * lda + 0]);
    rdiag[j] = acnorm[j];
    wa[j] = rdiag[j];
    if (pivot) {
      ipvt[j] = j + 1;
    }
  }
  minmn = ((m) <= (n) ? (m) : (n));
  for (j = 0; j < minmn; ++j) {
    if (pivot) {
      int kmax = j;
      for (k = j; k < n; ++k) {
        if (rdiag[k] > rdiag[kmax]) {
          kmax = k;
        }
      }
      if (kmax != j) {
        for (i = 0; i < m; ++i) {
          temp = a[i + j * lda];
          a[i + j * lda] = a[i + kmax * lda];
          a[i + kmax * lda] = temp;
        }
        rdiag[kmax] = rdiag[j];
        wa[kmax] = wa[j];
        k = ipvt[j];
        ipvt[j] = ipvt[kmax];
        ipvt[kmax] = k;
      }
    }
    ajnorm = enorm(m - (j + 1) + 1, &a[j + j * lda]);
    if (ajnorm != 0.) {
      if (a[j + j * lda] < 0.) {
        ajnorm = -ajnorm;
      }
      for (i = j; i < m; ++i) {
        a[i + j * lda] /= ajnorm;
      }
      a[j + j * lda] += 1.;
      jp1 = j + 1;
      if (n > jp1) {
        for (k = jp1; k < n; ++k) {
          sum = 0.;
          for (i = j; i < m; ++i) {
            sum += a[i + j * lda] * a[i + k * lda];
          }
          temp = sum / a[j + j * lda];
          for (i = j; i < m; ++i) {
            a[i + k * lda] -= temp * a[i + j * lda];
          }
          if (pivot && rdiag[k] != 0.) {
            temp = a[j + k * lda] / rdiag[k];
            d1 = 1. - temp * temp;
            rdiag[k] *= sqrt(((((double)0.) >= (d1) ? ((double)0.) : (d1))));
            d1 = rdiag[k] / wa[k];
            if (.05 * (d1 * d1) <= epsmch) {
              rdiag[k] = enorm(m - (j + 1), &a[jp1 + k * lda]);
              wa[k] = rdiag[k];
            }
          }
        }
      }
    }
    rdiag[j] = -ajnorm;
  }
}
```