user: Please show me the C implementation that matches the assembly code listed below:

```asm
<dsvdcmp0>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rcx,%rbp
  push   %rbx
  sub    $0xd8,%rsp
  mov    %rdi,0x60(%rsp)
  lea    0x0(,%rdx,8),%edi
  mov    %esi,0x70(%rsp)
  mov    %edx,0x74(%rsp)
  mov    %r8,0x80(%rsp)
  call   <malloc@plt>
  mov    %rax,0x88(%rsp)
  test   %rax,%rax
  je     L301
  mov    0x74(%rsp),%eax
  lea    -0x1(%rax),%edi
  mov    %edi,0x20(%rsp)
  test   %eax,%eax
  jle    L321
L74:
  sub    $0x2,%eax
  pxor   %xmm7,%xmm7
  mov    $0x1,%r12d
  xor    %r15d,%r15d
  mov    %eax,0x38(%rsp)
  mov    0x70(%rsp),%eax
  movapd %xmm7,%xmm2
  movapd %xmm7,%xmm0
  movsd  %xmm7,0x58(%rsp)
  sub    $0x1,%eax
  movsd  %xmm7,0x18(%rsp)
  movq   D67(%rip),%xmm7
  mov    %eax,0x30(%rsp)
  mov    %eax,0x50(%rsp)
  mov    0x60(%rsp),%rax
  mov    %r13d,0x98(%rsp)
  mov    %r12,%r13
  mov    %rbp,%r12
  mov    %r15,%rbp
  mov    %rax,%r10
  add    $0x8,%rax
  movaps %xmm7,(%rsp)
  movq   $0x8,0x28(%rsp)
  mov    %r10,%r14
  mov    %rax,0x90(%rsp)
  jmp    L755
L201:
  movsd  0x58(%rsp),%xmm1
  movq   $0x0,(%r12,%rbp,8)
  movapd %xmm1,%xmm2
L1233:
  mov    0x48(%rsp),%rax
  movsd  0x18(%rsp),%xmm7
  movsd  (%rax),%xmm0
  andpd  (%rsp),%xmm0
  addsd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm7
  ja     L888
L171:
  movsd  %xmm1,0x18(%rsp)
L888:
  cmp    %r13d,0x74(%rsp)
  jle    L684
L331:
  addq   $0x8,0x28(%rsp)
  add    $0x1,%rbp
  add    $0x1,%r13
  add    $0x8,%r14
  subl   $0x1,0x38(%rsp)
  subl   $0x1,0x50(%rsp)
L755:
  mulsd  %xmm2,%xmm0
  mov    0x88(%rsp),%rax
  lea    0x0(,%rbp,8),%r15
  mov    %ebp,0x68(%rsp)
  add    %r15,%rax
  mov    %rax,0x48(%rsp)
  movsd  %xmm0,(%rax)
  cmp    %ebp,0x70(%rsp)
  jle    L201
  mov    0x50(%rsp),%ebx
  mov    0x90(%rsp),%rax
  pxor   %xmm4,%xmm4
  add    %rbp,%rbx
  shl    $0x3,%rbx
  lea    (%rax,%rbx,1),%rcx
  mov    %r14,%rax
L484:
  mov    (%rax),%rdx
  add    $0x8,%rax
  movsd  (%rdx,%r15,1),%xmm0
  andpd  (%rsp),%xmm0
  addsd  %xmm0,%xmm4
  cmp    %rax,%rcx
  jne    L484
  movsd  0x58(%rsp),%xmm7
  ucomisd %xmm7,%xmm4
  jp     L859
  movapd %xmm7,%xmm1
  jne    L859
L822:
  mulsd  %xmm4,%xmm1
  mov    0x20(%rsp),%edi
  movsd  %xmm1,(%r12,%rbp,8)
  cmp    %edi,0x68(%rsp)
  jne    L821
L843:
  andpd  (%rsp),%xmm1
  pxor   %xmm2,%xmm2
  jmp    L1233
L821:
  mov    %r13d,0x68(%rsp)
  cmp    0x74(%rsp),%r13d
  jge    L843
  mov    (%r14),%rbx
  mov    %r13,%rax
  pxor   %xmm2,%xmm2
L1094:
  movsd  (%rbx,%rax,8),%xmm0
  andpd  (%rsp),%xmm0
  add    $0x1,%rax
  addsd  %xmm0,%xmm2
  cmp    %eax,0x74(%rsp)
  jg     L1094
  ucomisd 0x58(%rsp),%xmm2
  jp     L761
  jne    L761
  andpd  (%rsp),%xmm1
  jmp    L1233
L859:
  mov    %r14,%rax
  pxor   %xmm2,%xmm2
L813:
  mov    (%rax),%rdx
  add    $0x8,%rax
  add    %r15,%rdx
  movsd  (%rdx),%xmm0
  divsd  %xmm4,%xmm0
  movsd  %xmm0,(%rdx)
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm2
  cmp    %rax,%rcx
  jne    L813
  mov    (%r14),%r9
  movsd  0x58(%rsp),%xmm6
  lea    (%r9,%r15,1),%r11
  movsd  (%r11),%xmm7
  comisd %xmm6,%xmm7
  movsd  %xmm7,0x40(%rsp)
  jbe    L1069
  ucomisd %xmm2,%xmm6
  ja     L183
  movapd %xmm2,%xmm1
  sqrtsd %xmm1,%xmm1
L174:
  andpd  (%rsp),%xmm1
  xorpd  D59(%rip),%xmm1
L237:
  movsd  0x40(%rsp),%xmm7
  movapd %xmm7,%xmm5
  subsd  %xmm1,%xmm7
  mulsd  %xmm1,%xmm5
  movapd %xmm7,%xmm0
  movsd  %xmm7,(%r11)
  subsd  %xmm2,%xmm5
  cmp    %r13d,0x74(%rsp)
  jle    L480
  mov    0x38(%rsp),%eax
  mov    0x28(%rsp),%r8
  add    0x60(%rsp),%rbx
  lea    0x2(%rbp,%rax,1),%rax
  lea    0x0(,%rax,8),%rdx
L184:
  mov    %r14,%rcx
  movapd %xmm0,%xmm2
  pxor   %xmm3,%xmm3
  mov    %r9,%rax
  jmp    L1057
L657:
  mov    0x8(%rcx),%rax
  add    $0x8,%rcx
  movsd  (%rax,%r15,1),%xmm2
L1057:
  mulsd  (%rax,%r8,1),%xmm2
  addsd  %xmm2,%xmm3
  cmp    %rcx,%rbx
  jne    L657
  divsd  %xmm5,%xmm3
  mov    %r14,%rsi
  mov    %r9,%rax
  jmp    L573
L288:
  mov    0x8(%rsi),%rax
  add    $0x8,%rsi
  movsd  (%rax,%r15,1),%xmm0
L573:
  mulsd  %xmm3,%xmm0
  add    %r8,%rax
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  cmp    %rsi,%rcx
  jne    L288
  add    $0x8,%r8
  movsd  (%r11),%xmm0
  cmp    %r8,%rdx
  jne    L184
L176:
  mov    %r14,%rax
  jmp    L409
L1178:
  mov    0x8(%rax),%r9
  add    $0x8,%rax
  movsd  (%r9,%r15,1),%xmm0
L409:
  mulsd  %xmm4,%xmm0
  movsd  %xmm0,(%r9,%r15,1)
  cmp    %rax,%rbx
  jne    L1178
  jmp    L822
L684:
  mov    0x98(%rsp),%r13d
  mov    %r12,%rbp
L789:
  movslq 0x20(%rsp),%rdi
  movsd  D25(%rip),%xmm7
  mov    0x80(%rsp),%r9
  mov    0x74(%rsp),%r10d
  mov    %rdi,%rax
  movsd  %xmm7,0x78(%rsp)
  shl    $0x3,%rdi
  lea    0x8(%r9),%r12
  jmp    L1012
L837:
  mov    %r14d,%eax
L1012:
  mov    (%r9,%rdi,1),%rdx
  movsd  0x78(%rsp),%xmm7
  lea    -0x1(%rax),%r14d
  mov    0x88(%rsp),%rsi
  movsd  %xmm7,(%rdx,%rdi,1)
  movsd  (%rsi,%rdi,1),%xmm1
  sub    $0x8,%rdi
  test   %eax,%eax
  je     L604
  cmp    %r10d,%eax
  jge    L837
  ucomisd 0x58(%rsp),%xmm1
  jp     L952
  movslq %eax,%rdx
  jne    L952
L382:
  mov    (%r9,%rdi,1),%rcx
L834:
  mov    (%r9,%rdx,8),%rax
  movq   $0x0,(%rax,%rdi,1)
  movq   $0x0,(%rcx,%rdx,8)
  add    $0x1,%rdx
  cmp    %edx,%r10d
  jg     L834
  jmp    L837
L604:
  mov    0x70(%rsp),%edi
  cmp    %edi,0x74(%rsp)
  jg     L488
  mov    0x20(%rsp),%eax
L261:
  mov    0x74(%rsp),%esi
  mov    0x60(%rsp),%rdi
  movslq %eax,%r11
  lea    0x1(%rax),%r9d
  mov    0x70(%rsp),%r12d
  mov    %r13d,0x38(%rsp)
  mov    %r11,%r15
  mov    %rbp,%rcx
  lea    -0x2(%rsi),%edx
  mov    0x30(%rsp),%esi
  lea    (%rdi,%r11,8),%r10
  mov    %rdi,%r13
  sub    %eax,%edx
  sub    %eax,%esi
  mov    %edx,0x28(%rsp)
  mov    %esi,(%rsp)
L627:
  mov    (%r10),%r14
  lea    0x0(,%r15,8),%rbx
  movsd  (%rcx,%r15,8),%xmm0
  lea    (%r14,%rbx,1),%rbp
  cmp    %r9d,0x74(%rsp)
  jle    L794
  mov    0x28(%rsp),%r8d
  lea    0x8(%r14,%rbx,1),%rax
  xor    %esi,%esi
  mov    %rcx,0x90(%rsp)
  mov    %rax,%rdi
  mov    %r10,0x68(%rsp)
  lea    0x8(,%r8,8),%rdx
  mov    %r9d,0x50(%rsp)
  mov    %r8,0x40(%rsp)
  movsd  %xmm0,0x48(%rsp)
  call   <memset@plt>
  movsd  0x48(%rsp),%xmm0
  ucomisd 0x58(%rsp),%xmm0
  mov    0x40(%rsp),%r8
  mov    0x50(%rsp),%r9d
  mov    0x68(%rsp),%r10
  mov    0x90(%rsp),%rcx
  jp     L1203
  je     L78
L1203:
  movsd  0x78(%rsp),%xmm3
  mov    (%rsp),%eax
  movslq %r9d,%r11
  lea    0x2(%r8,%r15,1),%r8
  lea    0x0(,%r11,8),%rsi
  shl    $0x3,%r8
  divsd  %xmm0,%xmm3
  add    %r15,%rax
  movsd  0x0(%rbp),%xmm0
  lea    0x0(%r13,%rax,8),%rdi
L166:
  mov    %r11,%rax
  pxor   %xmm2,%xmm2
  cmp    %r9d,%r12d
  jle    L384
  nopw   0x0(%rax,%rax,1)
L629:
  mov    0x0(%r13,%rax,8),%rdx
  add    $0x1,%rax
  movsd  (%rdx,%rbx,1),%xmm1
  mulsd  (%rdx,%rsi,1),%xmm1
  addsd  %xmm1,%xmm2
  cmp    %eax,%r12d
  jg     L629
L384:
  divsd  %xmm0,%xmm2
  mulsd  %xmm3,%xmm2
  cmp    %r15d,%r12d
  jle    L811
  mov    %r10,%rdx
  mov    %r14,%rax
  jmp    L709
  nopl   0x0(%rax)
L766:
  mov    0x8(%rdx),%rax
  add    $0x8,%rdx
  movsd  (%rax,%rbx,1),%xmm0
L709:
  mulsd  %xmm2,%xmm0
  add    %rsi,%rax
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  cmp    %rdx,%rdi
  jne    L766
  movsd  0x0(%rbp),%xmm0
L811:
  add    $0x8,%rsi
  cmp    %rsi,%r8
  jne    L166
L368:
  cmp    %r15d,%r12d
  jle    L1085
  mov    (%rsp),%eax
  add    %r15,%rax
  lea    0x0(%r13,%rax,8),%rdx
  mov    %r10,%rax
  jmp    L572
L497:
  mov    0x8(%rax),%r14
  add    $0x8,%rax
  movsd  (%r14,%rbx,1),%xmm0
L572:
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,(%r14,%rbx,1)
  cmp    %rax,%rdx
  jne    L497
L586:
  movsd  0x0(%rbp),%xmm0
L1085:
  sub    $0x1,%r15
  addl   $0x1,(%rsp)
  sub    $0x8,%r10
  sub    $0x1,%r9d
  addsd  0x78(%rsp),%xmm0
  addl   $0x1,0x28(%rsp)
  movsd  %xmm0,0x0(%rbp)
  test   %r15d,%r15d
  jns    L627
  mov    0x38(%rsp),%r13d
  mov    %rcx,%rbp
L675:
  mov    0x20(%rsp),%edi
  test   %edi,%edi
  js     L227
  movslq %edi,%rax
  mov    0x60(%rsp),%rsi
  movq   D67(%rip),%xmm7
  mov    %rbp,%r14
  mov    %rax,0x90(%rsp)
  mov    0x30(%rsp),%eax
  movaps %xmm7,(%rsp)
  lea    0x8(%rsi,%rax,8),%r15
  mov    %edi,%eax
  mov    0x80(%rsp),%rdi
  lea    0x8(%rdi,%rax,8),%r12
  mov    %r12,%rbp
  mov    %r13d,%r12d
L135:
  mov    0x90(%rsp),%rax
  mov    %r12d,%r9d
  mov    %r14,%r12
  movl   $0x1e,0xb8(%rsp)
  lea    0x0(,%rax,8),%rsi
  mov    %eax,0xb0(%rsp)
  sub    $0x1,%eax
  mov    %eax,0x68(%rsp)
  cltq
  shl    $0x3,%rax
  mov    %rsi,0xa8(%rsp)
  lea    (%r14,%rax,1),%rdi
  mov    %rdi,0xc0(%rsp)
  mov    0x88(%rsp),%rdi
  add    %rdi,%rax
  add    %rsi,%rdi
  mov    %rax,0xc8(%rsp)
  mov    %rdi,0xa0(%rsp)
  nopw   0x0(%rax,%rax,1)
L64:
  mov    0xa8(%rsp),%rax
  mov    0xb0(%rsp),%ebx
  mov    0x88(%rsp),%rsi
  jmp    L831
  nopl   0x0(%rax)
L972:
  lea    (%rsi,%rax,1),%rdx
  movsd  0x18(%rsp),%xmm7
  lea    -0x1(%rbx),%r9d
  movsd  (%rdx),%xmm0
  andpd  (%rsp),%xmm0
  addsd  %xmm7,%xmm0
  ucomisd %xmm7,%xmm0
  jp     L1072
  je     L84
L1072:
  movsd  0x18(%rsp),%xmm7
  movsd  -0x8(%r12,%rax,1),%xmm0
  lea    -0x8(%rax),%rcx
  andpd  (%rsp),%xmm0
  addsd  %xmm7,%xmm0
  ucomisd %xmm7,%xmm0
  jp     L527
  je     L556
L527:
  mov    %rcx,%rax
  mov    %r9d,%ebx
L831:
  test   %ebx,%ebx
  jne    L972
  mov    0x88(%rsp),%rax
  mov    0x90(%rsp),%rdi
  mov    %rax,0x98(%rsp)
  xor    %eax,%eax
  mov    %edi,0x30(%rsp)
  cmp    %edi,%ebx
  jg     L1224
L964:
  movslq %ebx,%r10
  movslq %r9d,%r8
  mov    %ebx,0x38(%rsp)
  movsd  0x78(%rsp),%xmm4
  mov    %r9d,0x40(%rsp)
  shl    $0x3,%r8
  mov    %r12,%r14
  pxor   %xmm1,%xmm1
  mov    0x88(%rsp),%rax
  mov    %rbp,0x48(%rsp)
  mov    %r8,%rbp
  lea    (%rax,%r10,8),%r13
  mov    %r13,%r12
  mov    %r15,%r13
  mov    %r10,%r15
  nopw   0x0(%rax,%rax,1)
L465:
  movsd  (%r12),%xmm0
  movapd (%rsp),%xmm7
  lea    0x0(,%r15,8),%rbx
  mulsd  %xmm0,%xmm4
  mulsd  %xmm1,%xmm0
  andpd  %xmm4,%xmm7
  movsd  %xmm0,(%r12)
  movapd %xmm7,%xmm0
  movsd  0x18(%rsp),%xmm7
  addsd  %xmm7,%xmm0
  ucomisd %xmm7,%xmm0
  jp     L795
  je     L214
L795:
  movsd  (%r14,%r15,8),%xmm1
  movapd %xmm4,%xmm0
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm1,0x28(%rsp)
  call   <pythag@plt>
  movsd  0x78(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm4
  mov    0x70(%rsp),%edi
  movsd  %xmm0,(%r14,%r15,8)
  divsd  %xmm0,%xmm2
  xorpd  D59(%rip),%xmm4
  mov    0x60(%rsp),%rdx
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm4
  test   %edi,%edi
  jle    L1154
  nopl   0x0(%rax,%rax,1)
L1202:
  mov    (%rdx),%rax
  movapd %xmm1,%xmm3
  movapd %xmm4,%xmm5
  add    $0x8,%rdx
  lea    (%rax,%rbp,1),%rcx
  add    %rbx,%rax
  movsd  (%rcx),%xmm2
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm3
  mulsd  %xmm0,%xmm5
  mulsd  %xmm4,%xmm2
  mulsd  %xmm1,%xmm0
  addsd  %xmm5,%xmm3
  subsd  %xmm2,%xmm0
  movsd  %xmm3,(%rcx)
  movsd  %xmm0,(%rax)
  cmp    %rdx,%r13
  jne    L1202
L1154:
  add    $0x1,%r15
  add    $0x8,%r12
  cmp    %r15d,0x30(%rsp)
  jge    L465
L214:
  mov    0x38(%rsp),%ebx
  mov    0x40(%rsp),%r9d
  mov    %r13,%r15
  mov    %r14,%r12
  mov    0x48(%rsp),%rbp
L84:
  mov    0x90(%rsp),%rax
  movsd  (%r12,%rax,8),%xmm4
  cmp    %eax,%ebx
  je     L1176
  mov    0x88(%rsp),%rdi
  movslq %ebx,%rax
  shl    $0x3,%rax
  add    %rax,%rdi
  mov    %rdi,0x98(%rsp)
L1052:
  movsd  (%r12,%rax,1),%xmm8
  movapd %xmm4,%xmm1
  mov    0xc0(%rsp),%rax
  movsd  %xmm4,0x38(%rsp)
  movsd  (%rax),%xmm2
  mov    0xc8(%rsp),%rax
  movsd  %xmm8,0x40(%rsp)
  movapd %xmm2,%xmm3
  addsd  %xmm2,%xmm1
  movsd  (%rax),%xmm0
  mov    0xa0(%rsp),%rax
  subsd  %xmm4,%xmm3
  movsd  %xmm2,0x28(%rsp)
  movsd  (%rax),%xmm6
  mulsd  %xmm1,%xmm3
  movapd %xmm0,%xmm1
  addsd  %xmm6,%xmm0
  movsd  %xmm6,0x30(%rsp)
  subsd  %xmm6,%xmm1
  mulsd  %xmm1,%xmm0
  movsd  0x78(%rsp),%xmm1
  addsd  %xmm0,%xmm3
  movapd %xmm6,%xmm0
  addsd  %xmm6,%xmm0
  mulsd  %xmm2,%xmm0
  divsd  %xmm0,%xmm3
  movapd %xmm3,%xmm0
  movsd  %xmm3,0x20(%rsp)
  call   <pythag@plt>
  movsd  0x38(%rsp),%xmm4
  andpd  (%rsp),%xmm0
  movsd  0x40(%rsp),%xmm8
  movsd  0x20(%rsp),%xmm3
  comisd 0x58(%rsp),%xmm3
  movsd  0x28(%rsp),%xmm2
  movapd %xmm8,%xmm1
  movsd  0x30(%rsp),%xmm6
  subsd  %xmm4,%xmm1
  addsd  %xmm8,%xmm4
  mulsd  %xmm1,%xmm4
  ja     L632
  xorpd  D59(%rip),%xmm0
L632:
  addsd  %xmm3,%xmm0
  divsd  %xmm0,%xmm2
  subsd  %xmm6,%xmm2
  mulsd  %xmm2,%xmm6
  addsd  %xmm4,%xmm6
  divsd  %xmm8,%xmm6
  cmp    0x68(%rsp),%ebx
  jg     L312
  mov    0x88(%rsp),%rax
  movsd  0x78(%rsp),%xmm4
  mov    %r12,0x20(%rsp)
  movslq %ebx,%r14
  mov    %r15,%rbx
  mov    %r14,%r15
  lea    (%rax,%r14,8),%r13
  movapd %xmm4,%xmm2
  nop
L1048:
  mov    %r13,0x48(%rsp)
  movapd %xmm2,%xmm1
  movapd %xmm6,%xmm0
  add    $0x8,%r13
  movsd  0x0(%r13),%xmm7
  mov    0x20(%rsp),%rax
  movsd  %xmm8,0x50(%rsp)
  lea    0x0(,%r15,8),%r12
  movsd  %xmm6,0x28(%rsp)
  lea    0x8(%r12),%r14
  mulsd  %xmm7,%xmm1
  movsd  0x8(%rax,%r15,8),%xmm3
  mulsd  %xmm4,%xmm7
  movsd  %xmm3,0x40(%rsp)
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm7,0x38(%rsp)
  call   <pythag@plt>
  movsd  0x30(%rsp),%xmm1
  movsd  0x28(%rsp),%xmm6
  movsd  0x38(%rsp),%xmm7
  mov    0x48(%rsp),%rax
  movapd %xmm1,%xmm2
  movsd  0x50(%rsp),%xmm8
  mov    0x74(%rsp),%ecx
  movsd  %xmm0,(%rax)
  movsd  0x40(%rsp),%xmm3
  mov    0x80(%rsp),%rdx
  movapd %xmm3,%xmm1
  divsd  %xmm0,%xmm6
  divsd  %xmm0,%xmm2
  movapd %xmm7,%xmm0
  mulsd  %xmm6,%xmm3
  movapd %xmm6,%xmm9
  movapd %xmm6,%xmm4
  mulsd  %xmm8,%xmm9
  mulsd  %xmm6,%xmm7
  mulsd  %xmm2,%xmm0
  mulsd  %xmm2,%xmm1
  addsd  %xmm0,%xmm9
  movapd %xmm8,%xmm0
  movapd %xmm3,%xmm8
  mulsd  %xmm2,%xmm0
  subsd  %xmm0,%xmm7
  test   %ecx,%ecx
  jle    L564
  nopl   (%rax)
L392:
  mov    (%rdx),%rax
  movapd %xmm4,%xmm5
  movapd %xmm2,%xmm6
  add    $0x8,%rdx
  lea    (%rax,%r12,1),%rcx
  add    %r14,%rax
  movsd  (%rcx),%xmm3
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm5
  mulsd  %xmm0,%xmm6
  mulsd  %xmm2,%xmm3
  mulsd  %xmm4,%xmm0
  addsd  %xmm6,%xmm5
  subsd  %xmm3,%xmm0
  movsd  %xmm5,(%rcx)
  movsd  %xmm0,(%rax)
  cmp    %rdx,%rbp
  jne    L392
L564:
  movapd %xmm9,%xmm0
  movsd  %xmm8,0x40(%rsp)
  movsd  %xmm7,0x38(%rsp)
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm9,0x28(%rsp)
  movsd  %xmm2,0x50(%rsp)
  movsd  %xmm4,0x48(%rsp)
  call   <pythag@plt>
  ucomisd 0x58(%rsp),%xmm0
  mov    0x20(%rsp),%rax
  movsd  0x28(%rsp),%xmm9
  movsd  0x30(%rsp),%xmm1
  movsd  0x38(%rsp),%xmm7
  movsd  0x40(%rsp),%xmm8
  movsd  %xmm0,(%rax,%r15,8)
  jp     L1039
  movsd  0x48(%rsp),%xmm4
  movsd  0x50(%rsp),%xmm2
  je     L993
L1039:
  movsd  0x78(%rsp),%xmm2
  movapd %xmm9,%xmm4
  divsd  %xmm0,%xmm2
  mulsd  %xmm2,%xmm4
  mulsd  %xmm1,%xmm2
L993:
  movapd %xmm4,%xmm6
  movapd %xmm2,%xmm0
  mov    0x70(%rsp),%eax
  mov    0x60(%rsp),%rdx
  mulsd  %xmm7,%xmm6
  mulsd  %xmm8,%xmm0
  mulsd  %xmm2,%xmm7
  mulsd  %xmm4,%xmm8
  addsd  %xmm0,%xmm6
  subsd  %xmm7,%xmm8
  test   %eax,%eax
  jle    L552
  nopl   (%rax)
L615:
  mov    (%rdx),%rax
  movapd %xmm4,%xmm3
  movapd %xmm2,%xmm5
  add    $0x8,%rdx
  lea    (%rax,%r12,1),%rcx
  add    %r14,%rax
  movsd  (%rcx),%xmm1
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm3
  mulsd  %xmm0,%xmm5
  mulsd  %xmm2,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm5,%xmm3
  subsd  %xmm1,%xmm0
  movsd  %xmm3,(%rcx)
  movsd  %xmm0,(%rax)
  cmp    %rdx,%rbx
  jne    L615
L552:
  add    $0x1,%r15
  cmp    %r15d,0x68(%rsp)
  jge    L1048
  mov    0x20(%rsp),%r12
  mov    %rbx,%r15
L312:
  mov    0x98(%rsp),%rax
  subl   $0x1,0xb8(%rsp)
  movq   $0x0,(%rax)
  mov    0xa0(%rsp),%rax
  movsd  %xmm6,(%rax)
  mov    0x90(%rsp),%rax
  movsd  %xmm8,(%r12,%rax,8)
  je     L1002
  mov    0x68(%rsp),%r9d
  jmp    L64
L952:
  cmp    %r10d,%eax
  jge    L837
  mov    0x60(%rsp),%rsi
  movslq %eax,%rdx
  lea    0x0(,%rdx,8),%r8
  mov    %rdx,%rcx
  mov    (%rsi,%rdi,1),%r11
  lea    (%r11,%r8,1),%rsi
L911:
  movsd  (%r11,%rcx,8),%xmm0
  divsd  (%rsi),%xmm0
  divsd  %xmm1,%xmm0
  mov    (%r9,%rcx,8),%rbx
  add    $0x1,%rcx
  movsd  %xmm0,(%rbx,%rdi,1)
  cmp    %ecx,%r10d
  jg     L911
  mov    0x20(%rsp),%ecx
  lea    (%r9,%r8,1),%r15
  sub    %eax,%ecx
  lea    0x1(%rcx,%rdx,1),%rax
  add    %rdx,%rcx
  shl    $0x3,%rax
  lea    (%r12,%rcx,8),%rbx
  mov    %rax,(%rsp)
L711:
  mov    %rdx,%rax
  pxor   %xmm1,%xmm1
L443:
  mov    (%r9,%rax,8),%rcx
  movsd  (%rcx,%r8,1),%xmm0
  mulsd  (%r11,%rax,8),%xmm0
  add    $0x1,%rax
  addsd  %xmm0,%xmm1
  cmp    %eax,%r10d
  jg     L443
  mov    %r15,%rax
L1232:
  mov    (%rax),%rsi
  add    $0x8,%rax
  movsd  (%rsi,%rdi,1),%xmm0
  lea    (%rsi,%r8,1),%rcx
  mulsd  %xmm1,%xmm0
  addsd  (%rcx),%xmm0
  movsd  %xmm0,(%rcx)
  cmp    %rax,%rbx
  jne    L1232
  add    $0x8,%r8
  cmp    %r8,(%rsp)
  jne    L711
  jmp    L382
L488:
  mov    0x30(%rsp),%edi
  mov    %edi,%eax
  test   %edi,%edi
  jns    L261
  jmp    L675
  nop
L556:
  mov    0x90(%rsp),%rdi
  mov    %rdx,0x98(%rsp)
  mov    %edi,0x30(%rsp)
  cmp    %edi,%ebx
  jle    L964
L1224:
  movsd  (%r12,%rdi,8),%xmm4
  jmp    L1052
L1176:
  movsd  0x58(%rsp),%xmm7
  mov    %r12,%r14
  mov    %r9d,%r12d
  comisd %xmm4,%xmm7
  ja     L869
L274:
  subq   $0x1,0x90(%rsp)
  mov    0x90(%rsp),%rax
  test   %eax,%eax
  jns    L135
L227:
  mov    0x88(%rsp),%rdi
  add    $0xd8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10a0 <free@plt>
L78:
  cmp    %r15d,%r12d
  jle    L586
  mov    (%rsp),%eax
  add    %r15,%rax
  lea    0x0(%r13,%rax,8),%rdx
  mov    %r10,%rax
  jmp    L45
L939:
  mov    0x8(%rax),%r14
  add    $0x8,%rax
L45:
  movq   $0x0,(%r14,%rbx,1)
  cmp    %rax,%rdx
  jne    L939
  jmp    L586
L1002:
  mov    %r12,%r14
  mov    0x68(%rsp),%r12d
  jmp    L274
L869:
  xorpd  D59(%rip),%xmm4
  mov    0x74(%rsp),%esi
  movsd  %xmm4,(%r14,%rax,8)
  test   %esi,%esi
  jle    L274
  mov    0x80(%rsp),%rax
L30:
  mov    0xa8(%rsp),%rdx
  add    (%rax),%rdx
  add    $0x8,%rax
  movsd  (%rdx),%xmm0
  xorpd  D59(%rip),%xmm0
  movsd  %xmm0,(%rdx)
  cmp    %rbp,%rax
  jne    L30
  jmp    L274
L761:
  mov    0x38(%rsp),%r15d
  mov    0x28(%rsp),%rax
  pxor   %xmm1,%xmm1
  add    %rbp,%r15
  add    %rbx,%rax
  shl    $0x3,%r15
  mov    %rax,%r8
  lea    0x10(%rbx,%r15,1),%rdi
  mov    %rdi,0x40(%rsp)
  mov    %rax,%rdi
L356:
  movsd  (%rdi),%xmm0
  add    $0x8,%rdi
  divsd  %xmm2,%xmm0
  movsd  %xmm0,-0x8(%rdi)
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  cmp    %rdi,0x40(%rsp)
  jne    L356
  movsd  (%rax),%xmm7
  movsd  0x58(%rsp),%xmm6
  comisd %xmm6,%xmm7
  movsd  %xmm7,0x78(%rsp)
  jbe    L207
  ucomisd %xmm1,%xmm6
  ja     L1073
  movapd %xmm1,%xmm0
  sqrtsd %xmm0,%xmm0
L944:
  andpd  (%rsp),%xmm0
  xorpd  D59(%rip),%xmm0
L304:
  movsd  0x78(%rsp),%xmm7
  mov    0x88(%rsp),%rsi
  mov    %rbx,%rdi
  movapd %xmm7,%xmm3
  subsd  %xmm0,%xmm7
  lea    0x10(%rsi,%r15,1),%r9
  sub    %rsi,%rdi
  mulsd  %xmm0,%xmm3
  movsd  %xmm7,(%rax)
  mov    0x28(%rsp),%rax
  add    %rsi,%rax
  subsd  %xmm1,%xmm3
  movapd %xmm7,%xmm1
  jmp    L4
L1191:
  movsd  (%rdi,%rax,1),%xmm1
L4:
  divsd  %xmm3,%xmm1
  add    $0x8,%rax
  movsd  %xmm1,-0x8(%rax)
  cmp    %rax,%r9
  jne    L1191
  mov    0x70(%rsp),%edi
  mov    %r13,%r9
  lea    0x10(%r15),%rdx
  cmp    %edi,0x68(%rsp)
  jge    L996
  mov    0x74(%rsp),%edi
L324:
  mov    0x60(%rsp),%rax
  pxor   %xmm3,%xmm3
  mov    (%rax,%r9,8),%rcx
  mov    %r13,%rax
L222:
  movsd  (%rcx,%rax,8),%xmm1
  mulsd  (%rbx,%rax,8),%xmm1
  add    $0x1,%rax
  addsd  %xmm1,%xmm3
  cmp    %eax,%edi
  jg     L222
  mov    0x28(%rsp),%rax
  mov    0x88(%rsp),%rsi
  lea    (%rcx,%rdx,1),%r10
  add    %rcx,%rax
  sub    %rcx,%rsi
L15:
  movsd  (%rsi,%rax,1),%xmm1
  add    $0x8,%rax
  mulsd  %xmm3,%xmm1
  addsd  -0x8(%rax),%xmm1
  movsd  %xmm1,-0x8(%rax)
  cmp    %rax,%r10
  jne    L15
  add    $0x1,%r9
  cmp    %r9d,0x70(%rsp)
  jg     L324
L996:
  movsd  (%r8),%xmm1
  add    $0x8,%r8
  mulsd  %xmm2,%xmm1
  movsd  %xmm1,-0x8(%r8)
  cmp    %r8,0x40(%rsp)
  jne    L996
  mov    0x48(%rsp),%rax
  movapd (%rsp),%xmm7
  movsd  (%r12,%rbp,8),%xmm1
  movsd  (%rax),%xmm3
  andpd  %xmm7,%xmm1
  andpd  %xmm7,%xmm3
  movsd  0x18(%rsp),%xmm7
  addsd  %xmm3,%xmm1
  comisd %xmm1,%xmm7
  ja     L331
  jmp    L171
L794:
  ucomisd 0x58(%rsp),%xmm0
  jp     L1160
  je     L78
L1160:
  movsd  0x78(%rsp),%xmm3
  divsd  %xmm0,%xmm3
  movsd  0x0(%rbp),%xmm0
  jmp    L368
L301:
  call   <svd_errm@plt>
  mov    0x74(%rsp),%eax
  lea    -0x1(%rax),%edi
  mov    %edi,0x20(%rsp)
  test   %eax,%eax
  jg     L74
L321:
  mov    0x20(%rsp),%r8d
  test   %r8d,%r8d
  jns    L309
  mov    0x70(%rsp),%edi
  cmp    %edi,0x74(%rsp)
  jle    L227
  mov    0x70(%rsp),%eax
  sub    $0x1,%eax
  mov    %eax,0x30(%rsp)
  js     L227
  pxor   %xmm7,%xmm7
  movsd  %xmm7,0x58(%rsp)
  movsd  %xmm7,0x18(%rsp)
  movsd  D25(%rip),%xmm7
  movsd  %xmm7,0x78(%rsp)
  jmp    L261
L1069:
  movsd  0x58(%rsp),%xmm7
  ucomisd %xmm2,%xmm7
  ja     L8
  movapd %xmm2,%xmm1
  sqrtsd %xmm1,%xmm1
L896:
  andpd  (%rsp),%xmm1
  jmp    L237
L207:
  movsd  0x58(%rsp),%xmm7
  ucomisd %xmm1,%xmm7
  ja     L21
  movapd %xmm1,%xmm0
  sqrtsd %xmm0,%xmm0
L441:
  andpd  (%rsp),%xmm0
  jmp    L304
L480:
  mov    0x60(%rsp),%rax
  add    %rax,%rbx
  jmp    L176
L183:
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x78(%rsp)
  mov    %r9,0xb0(%rsp)
  mov    %r11,0xa0(%rsp)
  movsd  %xmm4,0xa8(%rsp)
  call   <sqrt@plt>
  mov    0xb0(%rsp),%r9
  movsd  0xa8(%rsp),%xmm4
  mov    0xa0(%rsp),%r11
  movsd  0x78(%rsp),%xmm2
  movapd %xmm0,%xmm1
  jmp    L174
L21:
  movapd %xmm1,%xmm0
  mov    %r8,0xb0(%rsp)
  mov    %rax,0xa8(%rsp)
  movsd  %xmm2,0xb8(%rsp)
  movsd  %xmm1,0xa0(%rsp)
  call   <sqrt@plt>
  mov    0xb0(%rsp),%r8
  movsd  0xb8(%rsp),%xmm2
  mov    0xa8(%rsp),%rax
  movsd  0xa0(%rsp),%xmm1
  jmp    L441
L1073:
  movapd %xmm1,%xmm0
  mov    %r8,0xb0(%rsp)
  mov    %rax,0xa8(%rsp)
  movsd  %xmm2,0xb8(%rsp)
  movsd  %xmm1,0xa0(%rsp)
  call   <sqrt@plt>
  mov    0xb0(%rsp),%r8
  movsd  0xb8(%rsp),%xmm2
  mov    0xa8(%rsp),%rax
  movsd  0xa0(%rsp),%xmm1
  jmp    L944
L8:
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x78(%rsp)
  mov    %r9,0xb0(%rsp)
  mov    %r11,0xa0(%rsp)
  movsd  %xmm4,0xa8(%rsp)
  call   <sqrt@plt>
  mov    0xb0(%rsp),%r9
  movsd  0xa8(%rsp),%xmm4
  mov    0xa0(%rsp),%r11
  movsd  0x78(%rsp),%xmm2
  movapd %xmm0,%xmm1
  jmp    L896
L309:
  mov    0x70(%rsp),%eax
  pxor   %xmm7,%xmm7
  movsd  %xmm7,0x58(%rsp)
  sub    $0x1,%eax
  movsd  %xmm7,0x18(%rsp)
  mov    %eax,0x30(%rsp)
  jmp    L789
```

tool_calls: name: unpack_data
arguments: tag: D67
type_infered: qword
name: unpack_data
arguments: tag: D59
type_infered: f32[4]
name: unpack_data
arguments: tag: D25
type_infered: f64

tool: D67: 9223372036854775807

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

tool: D25: 1

assistant: ```c
void dsvdcmp0(double **a, int m, int n, double *w, double **v) {
  int flag, i, its, j, jj, k, l, nm;
  double anorm, c, f, g, h, s, scale, x, y, z, *rv1;
  if (!(rv1 = (double *)malloc(n * sizeof(double))))
    svd_errm();
  g = scale = anorm = 0.0;
  for (i = 0; i < n; i++) {
    l = i + 1;
    rv1[i] = scale * g;
    g = s = scale = 0.0;
    if (i < m) {
      for (k = i; k < m; k++)
        scale += fabs(a[k][i]);
      if (scale) {
        for (k = i; k < m; k++) {
          a[k][i] /= scale;
          s += a[k][i] * a[k][i];
        }
        f = a[i][i];
        g = -((f) > 0 ? fabs(sqrt(s)) : -fabs(sqrt(s)));
        h = f * g - s;
        a[i][i] = f - g;
        for (j = l; j < n; j++) {
          for (s = 0.0, k = i; k < m; k++)
            s += a[k][i] * a[k][j];
          f = s / h;
          for (k = i; k < m; k++)
            a[k][j] += f * a[k][i];
        }
        for (k = i; k < m; k++)
          a[k][i] *= scale;
      }
    }
    w[i] = scale * g;
    g = s = scale = 0.0;
    if (i < m && i != n - 1) {
      for (k = l; k < n; k++)
        scale += fabs(a[i][k]);
      if (scale) {
        for (k = l; k < n; k++) {
          a[i][k] /= scale;
          s += a[i][k] * a[i][k];
        }
        f = a[i][l];
        g = -((f) > 0 ? fabs(sqrt(s)) : -fabs(sqrt(s)));
        h = f * g - s;
        a[i][l] = f - g;
        for (k = l; k < n; k++)
          rv1[k] = a[i][k] / h;
        for (j = l; j < m; j++) {
          for (s = 0.0, k = l; k < n; k++)
            s += a[j][k] * a[i][k];
          for (k = l; k < n; k++)
            a[j][k] += s * rv1[k];
        }
        for (k = l; k < n; k++)
          a[i][k] *= scale;
      }
    }
    anorm = ((anorm) > ((fabs(w[i]) + fabs(rv1[i]))) ? (anorm) : ((fabs(w[i]) + fabs(rv1[i]))));
  }
  for (i = n - 1; i >= 0; i--) {
    if (i < n - 1) {
      if (g) {
        for (j = l; j < n; j++)
          v[j][i] = (a[i][j] / a[i][l]) / g;
        for (j = l; j < n; j++) {
          for (s = 0.0, k = l; k < n; k++)
            s += a[i][k] * v[k][j];
          for (k = l; k < n; k++)
            v[k][j] += s * v[k][i];
        }
      }
      for (j = l; j < n; j++)
        v[i][j] = v[j][i] = 0.0;
    }
    v[i][i] = 1.0;
    g = rv1[i];
    l = i;
  }
  for (i = ((m - 1) < (n - 1) ? (m - 1) : (n - 1)); i >= 0; i--) {
    l = i + 1;
    g = w[i];
    for (j = l; j < n; j++)
      a[i][j] = 0.0;
    if (g) {
      g = 1.0 / g;
      for (j = l; j < n; j++) {
        for (s = 0.0, k = l; k < m; k++)
          s += a[k][i] * a[k][j];
        f = (s / a[i][i]) * g;
        for (k = i; k < m; k++)
          a[k][j] += f * a[k][i];
      }
      for (j = i; j < m; j++)
        a[j][i] *= g;
    } else
      for (j = i; j < m; j++)
        a[j][i] = 0.0;
    a[i][i] += 1.0;
  }
  for (k = n - 1; k >= 0; k--) {
    for (its = 0; its < 30; its++) {
      flag = 1;
      for (l = k; l > 0; l--) {
        nm = l - 1;
        if ((double)(fabs(rv1[l]) + anorm) == anorm) {
          flag = 0;
          break;
        }
        if ((double)(fabs(w[nm]) + anorm) == anorm)
          break;
      }
      if (flag) {
        c = 0.0;
        s = 1.0;
        for (i = l; i <= k; i++) {
          f = s * rv1[i];
          rv1[i] = c * rv1[i];
          if ((double)(fabs(f) + anorm) == anorm)
            break;
          g = w[i];
          h = pythag(f, g);
          w[i] = h;
          h = 1.0 / h;
          c = g * h;
          s = -f * h;
          for (j = 0; j < m; j++) {
            y = a[j][nm];
            z = a[j][i];
            a[j][nm] = y * c + z * s;
            a[j][i] = z * c - y * s;
          }
        }
      }
      z = w[k];
      if (l == k) {
        if (z < 0.0) {
          w[k] = -z;
          for (j = 0; j < n; j++)
            v[j][k] = -v[j][k];
        }
        break;
      }
      if (its >= 30) {
        printf("dsvdcmp0: no convergence in 30 iterations\n");
      }
      x = w[l];
      nm = k - 1;
      y = w[nm];
      g = rv1[nm];
      h = rv1[k];
      f = ((y - z) * (y + z) + (g - h) * (g + h)) / (2.0 * h * y);
      g = pythag(f, 1.0);
      f = ((x - z) * (x + z) + h * ((y / (f + ((f) > 0 ? fabs(g) : -fabs(g)))) - h)) / x;
      c = s = 1.0;
      for (j = l; j <= nm; j++) {
        i = j + 1;
        g = rv1[i];
        y = w[i];
        h = s * g;
        g = c * g;
        z = pythag(f, h);
        rv1[j] = z;
        c = f / z;
        s = h / z;
        f = x * c + g * s;
        g = g * c - x * s;
        h = y * s;
        y *= c;
        for (jj = 0; jj < n; jj++) {
          x = v[jj][j];
          z = v[jj][i];
          v[jj][j] = x * c + z * s;
          v[jj][i] = z * c - x * s;
        }
        z = pythag(f, h);
        w[j] = z;
        if (z) {
          z = 1.0 / z;
          c = f * z;
          s = h * z;
        }
        f = c * g + s * y;
        x = c * y - s * g;
        for (jj = 0; jj < m; jj++) {
          y = a[jj][j];
          z = a[jj][i];
          a[jj][j] = y * c + z * s;
          a[jj][i] = z * c - y * s;
        }
      }
      rv1[l] = 0.0;
      rv1[k] = f;
      w[k] = x;
    }
  }
  free(rv1);
}
```