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

```asm
<relax_svd>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xe8,%rsp
  mov    %rdi,0x58(%rsp)
  mov    0x120(%rsp),%r14
  mov    %ecx,0x20(%rsp)
  mov    0x128(%rsp),%r10
  mov    %r8,0x68(%rsp)
  mov    %r9d,0x60(%rsp)
  cmp    %edx,%esi
  jl     L132
  lea    -0x1(%rdx),%eax
  mov    %esi,%ebp
  mov    %edx,%ebx
  mov    %eax,0x80(%rsp)
  test   %edx,%edx
  jle    L68
  lea    -0x2(%rdx),%eax
  pxor   %xmm7,%xmm7
  mov    %rdi,%r11
  mov    $0x1,%r15d
  mov    %eax,0x38(%rsp)
  lea    0x1(%rcx),%eax
  movapd %xmm7,%xmm3
  movapd %xmm7,%xmm2
  mov    %eax,0x78(%rsp)
  movslq %ecx,%rax
  mov    %r14,%r9
  mov    %r15,%rsi
  lea    0x0(,%rax,8),%rdi
  lea    0x8(,%rax,8),%rcx
  mov    %r11,%rax
  movsd  %xmm7,0x70(%rsp)
  movsd  %xmm7,0x18(%rsp)
  movq   D43(%rip),%xmm7
  add    $0x10,%rax
  mov    %rdi,%r12
  mov    %r11,0x48(%rsp)
  xor    %r11d,%r11d
  mov    %ebp,%edi
  mov    %r10,%rbp
  mov    %rcx,0x88(%rsp)
  mov    %r11,%r14
  movl   $0x1,0x40(%rsp)
  movl   $0x0,0x28(%rsp)
  mov    %rax,0xa0(%rsp)
  movaps %xmm7,(%rsp)
  jmp    L745
L860:
  movq   $0x0,(%r9)
  pxor   %xmm2,%xmm2
L958:
  movapd (%rsp),%xmm7
  movsd  0x0(%rbp,%r14,8),%xmm0
  pxor   %xmm3,%xmm3
  andpd  %xmm7,%xmm0
  andpd  %xmm7,%xmm2
  movsd  0x18(%rsp),%xmm7
  addsd  %xmm2,%xmm0
  comisd %xmm0,%xmm7
  ja     L367
L841:
  movsd  %xmm0,0x18(%rsp)
L367:
  cmp    %esi,%ebx
  jle    L593
L924:
  mov    0x28(%rsp),%eax
  mov    0x20(%rsp),%ecx
  pxor   %xmm2,%xmm2
  lea    (%rax,%rcx,1),%r10d
L194:
  mov    0x78(%rsp),%ecx
  add    %ecx,0x40(%rsp)
  add    $0x1,%r14
  add    $0x1,%rsi
  mov    0x88(%rsp),%rcx
  subl   $0x1,0x38(%rsp)
  add    $0x8,%r9
  add    %rcx,0x48(%rsp)
  mov    %r10d,0x28(%rsp)
L745:
  mulsd  %xmm3,%xmm2
  mov    %esi,0x50(%rsp)
  mov    %r14d,%r8d
  mov    %r14d,0x30(%rsp)
  movsd  %xmm2,0x0(%rbp,%r14,8)
  cmp    %r14d,%edi
  jle    L860
  movslq 0x28(%rsp),%rax
  mov    0x58(%rsp),%rcx
  pxor   %xmm2,%xmm2
  mov    %rax,0x90(%rsp)
  add    %r14,%rax
  lea    (%rcx,%rax,8),%r10
  mov    %r14d,%eax
  mov    %r10,%rdx
L376:
  movsd  (%rdx),%xmm0
  andpd  (%rsp),%xmm0
  mov    %eax,%r13d
  add    $0x1,%eax
  add    %r12,%rdx
  addsd  %xmm0,%xmm2
  cmp    %eax,%edi
  jne    L376
  movsd  0x70(%rsp),%xmm7
  ucomisd %xmm7,%xmm2
  jp     L170
  jne    L170
  mulsd  %xmm7,%xmm2
  movsd  %xmm2,(%r9)
  cmp    0x30(%rsp),%edi
  jle    L958
L425:
  mov    0x30(%rsp),%ecx
  cmp    %ecx,0x80(%rsp)
  jne    L855
  movsd  (%r9),%xmm2
  jmp    L958
L855:
  cmp    %esi,%ebx
  jle    L396
  mov    0x90(%rsp),%rcx
  mov    0x58(%rsp),%rdx
  pxor   %xmm3,%xmm3
  lea    (%rcx,%rsi,1),%rax
  lea    (%rdx,%rax,8),%r15
  mov    0x38(%rsp),%eax
  add    %rcx,%rax
  mov    0xa0(%rsp),%rcx
  add    %r14,%rax
  lea    (%rcx,%rax,8),%r13
  mov    %r15,%rax
L967:
  movsd  (%rax),%xmm0
  andpd  (%rsp),%xmm0
  add    $0x8,%rax
  addsd  %xmm0,%xmm3
  cmp    %r13,%rax
  jne    L967
  ucomisd 0x70(%rsp),%xmm3
  jp     L17
  jne    L17
  movapd (%rsp),%xmm7
  movsd  (%r9),%xmm0
  movsd  0x0(%rbp,%r14,8),%xmm1
  andpd  %xmm7,%xmm0
  andpd  %xmm7,%xmm1
  movsd  0x18(%rsp),%xmm7
  addsd  %xmm1,%xmm0
  comisd %xmm0,%xmm7
  ja     L924
  jmp    L841
L132:
  mov    $0xffffffff,%eax
L543:
  add    $0xe8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L593:
  mov    0x120(%rsp),%r14
  mov    %rbp,%r10
  mov    %edi,%ebp
L54:
  mov    0x80(%rsp),%r15d
  mov    %r15d,%ecx
  test   %r15d,%r15d
  js     L257
  movslq 0x60(%rsp),%r8
  mov    0x68(%rsp),%rsi
  movslq %r15d,%r13
  xor    %edi,%edi
  lea    -0x2(%rbx),%r11d
  mov    %ebp,0xb0(%rsp)
  mov    %edi,%ebp
  mov    %r8d,%eax
  mov    %r10,0x128(%rsp)
  lea    0x0(,%r8,8),%r12
  not    %eax
  mov    %r13,0x78(%rsp)
  cltq
  mov    %r11d,0x88(%rsp)
  shl    $0x3,%rax
  mov    %ebx,0xa8(%rsp)
  mov    %rax,0x40(%rsp)
  mov    %r8d,%eax
  imul   %r15d,%eax
  mov    %r14,0x120(%rsp)
  movslq %eax,%rdx
  mov    %eax,%r10d
  add    %r13,%rdx
  lea    (%rsi,%rdx,8),%r9
  mov    %r11d,%edx
  mov    0x20(%rsp),%esi
  imul   %r8d,%edx
  mov    %edx,0x38(%rsp)
  mov    %esi,%edx
  not    %esi
  imul   %r11d,%edx
  mov    %esi,0x28(%rsp)
  mov    0x68(%rsp),%rsi
  add    %edx,%ecx
  mov    %edx,0x30(%rsp)
  mov    %ecx,(%rsp)
  lea    0x8(%rsi),%rcx
  mov    %rcx,0x98(%rsp)
  mov    0x58(%rsp),%rcx
  add    $0x8,%rcx
  mov    %rcx,0xa0(%rsp)
  jmp    L646
L152:
  mov    0x60(%rsp),%ecx
  sub    %ecx,0x38(%rsp)
  add    $0x1,%ebp
  sub    $0x1,%r13
  add    0x40(%rsp),%r9
  mov    %edi,%r15d
  sub    %ecx,%r10d
  mov    0x20(%rsp),%ecx
  sub    %ecx,0x30(%rsp)
  mov    0x28(%rsp),%ecx
  add    %ecx,(%rsp)
L646:
  mov    D60(%rip),%rax
  lea    -0x1(%r15),%edi
  mov    %rax,(%r9)
  test   %r15d,%r15d
  je     L938
  cmp    %r13d,%ebx
  jle    L152
  mov    0x128(%rsp),%rax
  movsd  (%rax,%r13,8),%xmm1
  ucomisd 0x70(%rsp),%xmm1
  jp     L900
  movslq %edi,%rax
  movslq %r10d,%r11
  mov    %ebp,%r8d
  mov    %rax,0x90(%rsp)
  jne    L900
L379:
  mov    0x90(%rsp),%rax
  mov    0x68(%rsp),%rsi
  movslq 0x38(%rsp),%rcx
  add    %r11,%rax
  lea    (%rsi,%rax,8),%rdx
  lea    0x0(%r13,%rcx,1),%rax
  add    %r8,%rcx
  lea    (%rsi,%rax,8),%rax
  mov    0x98(%rsp),%rsi
  add    %r13,%rcx
  lea    (%rsi,%rcx,8),%rcx
  nopl   0x0(%rax)
L383:
  movq   $0x0,(%rdx)
  add    $0x8,%rax
  add    %r12,%rdx
  movq   $0x0,-0x8(%rax)
  cmp    %rax,%rcx
  jne    L383
  jmp    L152
L938:
  movslq 0x28(%rsp),%rax
  mov    0x20(%rsp),%edi
  mov    %ebx,0x38(%rsp)
  mov    0x80(%rsp),%r9d
  mov    0x78(%rsp),%rsi
  mov    %r12,0x80(%rsp)
  shl    $0x3,%rax
  mov    0x58(%rsp),%r15
  mov    0xa8(%rsp),%r11d
  mov    %rax,0x60(%rsp)
  mov    %edi,%eax
  mov    0xb0(%rsp),%ebp
  mov    %rsi,%r12
  imul   %r9d,%eax
  mov    %r11d,%r13d
  movslq %eax,%rdx
  mov    %eax,(%rsp)
  lea    0x1(%rdi),%eax
  imul   %r9d,%eax
  lea    (%rsi,%rdx,1),%rcx
  lea    (%r15,%rcx,8),%r8
  mov    %r8,%r14
  mov    %eax,0x30(%rsp)
  mov    %edi,%eax
  neg    %eax
  cltq
  lea    -0x8(,%rax,8),%rax
  mov    %rax,0x50(%rsp)
  lea    0x1(%rdx,%rsi,1),%rax
  lea    (%r15,%rax,8),%rdx
  movslq %edi,%rax
  mov    %rax,0xa8(%rsp)
  lea    0x0(,%rax,8),%r15
  mov    %rdx,%rcx
L658:
  mov    0x120(%rsp),%rax
  mov    %r12d,%ebx
  movsd  (%rax,%r12,8),%xmm0
  cmp    %r13d,0x38(%rsp)
  jle    L923
  mov    0x88(%rsp),%r8d
  xor    %esi,%esi
  mov    %rcx,%rdi
  movsd  %xmm0,0x48(%rsp)
  sub    %r12d,%r8d
  lea    0x8(,%r8,8),%rdx
  mov    %r8,0x40(%rsp)
  call   <memset@plt>
  movsd  0x48(%rsp),%xmm0
  mov    0x40(%rsp),%r8
  ucomisd 0x70(%rsp),%xmm0
  mov    %rax,%rcx
  jp     L545
  je     L689
L545:
  movsd  D60(%rip),%xmm2
  mov    0x58(%rsp),%rdi
  movslq %r13d,%rdx
  mov    $0xfffffffffffffffe,%r9
  movslq 0x30(%rsp),%rax
  sub    %r8,%r9
  divsd  %xmm0,%xmm2
  lea    (%rdi,%rax,8),%r10
  movslq (%rsp),%rax
  mov    %rax,%rsi
  add    %rdx,%rax
  mov    %rdi,%rdx
  lea    (%rdi,%rax,8),%rdi
  mov    0x20(%rsp),%eax
  add    %esi,%eax
  mov    $0x1,%esi
  cltq
  add    %r12,%rax
  lea    (%rdx,%rax,8),%r8
  mov    $0xffffffffffffffff,%rdx
  nopl   0x0(%rax)
L333:
  mov    %r8,%rax
  mov    %r13d,%r11d
  pxor   %xmm1,%xmm1
  cmp    %r13d,%ebp
  jle    L994
  nop
L823:
  movsd  (%rax),%xmm0
  mulsd  (%rax,%rsi,8),%xmm0
  add    $0x1,%r11d
  add    %r15,%rax
  addsd  %xmm0,%xmm1
  cmp    %r11d,%ebp
  jne    L823
L994:
  divsd  (%r10),%xmm1
  mulsd  %xmm2,%xmm1
  mov    %rdi,%rax
  mov    %ebx,%r11d
  cmp    %r12d,%ebp
  jle    L482
  nopl   (%rax)
L654:
  movsd  (%rax,%rdx,8),%xmm0
  add    $0x1,%r11d
  mulsd  %xmm1,%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    %r15,%rax
  cmp    %r11d,%ebp
  jne    L654
L482:
  sub    $0x1,%rdx
  add    $0x8,%rdi
  add    $0x1,%rsi
  cmp    %rdx,%r9
  jne    L333
L1008:
  cmp    %r12d,%ebp
  jle    L369
  movslq (%rsp),%rax
  mov    0x58(%rsp),%rdi
  add    %r12,%rax
  lea    (%rdi,%rax,8),%rax
  nopl   0x0(%rax,%rax,1)
L1020:
  movsd  (%rax),%xmm0
  add    $0x1,%ebx
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  add    %r15,%rax
  cmp    %ebx,%ebp
  jne    L1020
L369:
  movsd  D60(%rip),%xmm0
  sub    $0x1,%r12
  mov    0x20(%rsp),%edi
  sub    $0x1,%r13d
  addsd  (%r14),%xmm0
  sub    %edi,(%rsp)
  add    0x50(%rsp),%rcx
  mov    0x28(%rsp),%edi
  add    %edi,0x30(%rsp)
  movsd  %xmm0,(%r14)
  add    0x60(%rsp),%r14
  test   %r12d,%r12d
  jns    L658
  mov    0x80(%rsp),%r12
  mov    0x38(%rsp),%ebx
  mov    0x128(%rsp),%r10
  movq   D43(%rip),%xmm7
  mov    %r12,%r14
  mov    0x120(%rsp),%r12
  mov    %ebx,%r13d
  movaps %xmm7,(%rsp)
  mov    %r10,%rbx
L456:
  mov    0x78(%rsp),%rcx
  mov    %rbx,0x128(%rsp)
  movl   $0x1e,0x90(%rsp)
  mov    %ecx,%eax
  mov    %ecx,0x80(%rsp)
  sub    $0x1,%eax
  mov    %eax,0x60(%rsp)
  cltq
  shl    $0x3,%rax
  lea    (%r12,%rax,1),%rdi
  add    %rbx,%rax
  mov    %r13d,%ebx
  mov    %r12,%r13
  mov    %rdi,0x98(%rsp)
  mov    0x128(%rsp),%r12
  mov    %rax,0xa0(%rsp)
  xchg   %ax,%ax
L143:
  mov    0x78(%rsp),%rsi
  nopl   (%rax)
L1039:
  movsd  0x18(%rsp),%xmm7
  movsd  (%r12,%rsi,8),%xmm0
  mov    %esi,%r9d
  lea    -0x1(%rsi),%eax
  andpd  (%rsp),%xmm0
  addsd  %xmm7,%xmm0
  ucomisd %xmm7,%xmm0
  jp     L486
  je     L966
L486:
  movsd  0x18(%rsp),%xmm7
  movsd  -0x8(%r13,%rsi,8),%xmm0
  andpd  (%rsp),%xmm0
  addsd  %xmm7,%xmm0
  ucomisd %xmm7,%xmm0
  jp     L15
  je     L850
L15:
  sub    $0x1,%rsi
  mov    %esi,%r9d
  cmp    $0xffffffff,%esi
  jne    L1039
L850:
  cmp    0x78(%rsp),%r9d
  jge    L966
  movslq %r9d,%r10
  cltq
  mov    0x58(%rsp),%rcx
  mov    %r9d,0x38(%rsp)
  mov    %r10,%rsi
  mov    %r14,0x40(%rsp)
  pxor   %xmm1,%xmm1
  mov    %ebp,%r14d
  sub    %rax,%rsi
  add    0xa8(%rsp),%rax
  mov    %ebx,0x48(%rsp)
  mov    %r13,%rbp
  lea    (%rcx,%rax,8),%rax
  mov    %r12,%rbx
  mov    %r10,%r13
  mov    %r15,%r12
  mov    %rax,0x30(%rsp)
  mov    D60(%rip),%rax
  mov    %rsi,%r15
  movq   %rax,%xmm4
  nopl   (%rax)
L690:
  movsd  (%rbx,%r13,8),%xmm0
  movapd (%rsp),%xmm7
  mulsd  %xmm0,%xmm4
  mulsd  %xmm1,%xmm0
  andpd  %xmm4,%xmm7
  movsd  %xmm0,(%rbx,%r13,8)
  movapd %xmm7,%xmm0
  movsd  0x18(%rsp),%xmm7
  addsd  %xmm7,%xmm0
  ucomisd %xmm7,%xmm0
  jp     L24
  je     L741
L24:
  movsd  0x0(%rbp,%r13,8),%xmm1
  movapd %xmm4,%xmm0
  mov    $0x2,%eax
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm1,0x28(%rsp)
  call   <pythag@plt>
  pxor   %xmm2,%xmm2
  movsd  D60(%rip),%xmm0
  movsd  0x28(%rsp),%xmm1
  cvtsi2sd %eax,%xmm2
  movsd  0x20(%rsp),%xmm4
  mov    0x30(%rsp),%rax
  mov    $0x1,%r8d
  xorpd  D50(%rip),%xmm4
  divsd  %xmm2,%xmm0
  movsd  %xmm2,0x0(%rbp,%r13,8)
  mulsd  %xmm0,%xmm1
  mulsd  %xmm0,%xmm4
  cmp    $0x1,%r14d
  jle    L79
  nopl   0x0(%rax,%rax,1)
L344:
  movsd  (%rax),%xmm2
  movapd %xmm1,%xmm3
  movapd %xmm4,%xmm5
  add    $0x1,%r8d
  movsd  (%rax,%r15,8),%xmm0
  mulsd  %xmm2,%xmm3
  mulsd  %xmm0,%xmm5
  mulsd  %xmm4,%xmm2
  mulsd  %xmm1,%xmm0
  addsd  %xmm5,%xmm3
  subsd  %xmm2,%xmm0
  movsd  %xmm3,(%rax)
  movsd  %xmm0,(%rax,%r15,8)
  add    %r12,%rax
  cmp    %r8d,%r14d
  jne    L344
L79:
  add    $0x1,%r13
  add    $0x1,%r15
  cmp    %r13d,0x80(%rsp)
  jg     L690
L741:
  mov    0x78(%rsp),%rax
  mov    %rbp,%r13
  mov    %r12,%r15
  mov    %r14d,%ebp
  mov    %rbx,%r12
  mov    0x38(%rsp),%r9d
  mov    0x40(%rsp),%r14
  mov    0x48(%rsp),%ebx
  movsd  0x0(%r13,%rax,8),%xmm3
L660:
  subl   $0x1,0x90(%rsp)
  je     L132
  movslq %r9d,%rsi
  movapd %xmm3,%xmm5
  mov    %r9d,0x88(%rsp)
  lea    0x0(,%rsi,8),%rax
  movsd  0x0(%r13,%rsi,8),%xmm8
  mov    %rsi,0x50(%rsp)
  mov    %rax,0x20(%rsp)
  mov    0x98(%rsp),%rax
  movsd  %xmm8,0x48(%rsp)
  movsd  (%rax),%xmm2
  mov    0xa0(%rsp),%rax
  movsd  %xmm3,0x40(%rsp)
  movapd %xmm2,%xmm0
  addsd  %xmm2,%xmm5
  movsd  (%rax),%xmm1
  mov    0x78(%rsp),%rax
  subsd  %xmm3,%xmm0
  movsd  %xmm2,0x30(%rsp)
  movsd  (%r12,%rax,8),%xmm6
  mov    D60(%rip),%rax
  mulsd  %xmm5,%xmm0
  movapd %xmm1,%xmm5
  addsd  %xmm6,%xmm1
  movsd  %xmm6,0x38(%rsp)
  subsd  %xmm6,%xmm5
  mulsd  %xmm5,%xmm1
  addsd  %xmm1,%xmm0
  movapd %xmm6,%xmm1
  addsd  %xmm6,%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movq   %rax,%xmm1
  mov    $0x2,%eax
  movsd  %xmm0,0x28(%rsp)
  call   <pythag@plt>
  movsd  0x48(%rsp),%xmm8
  movsd  0x40(%rsp),%xmm3
  pxor   %xmm1,%xmm1
  movsd  0x28(%rsp),%xmm0
  comisd 0x70(%rsp),%xmm0
  cvtsi2sd %eax,%xmm1
  movapd %xmm8,%xmm5
  andpd  (%rsp),%xmm1
  movsd  0x30(%rsp),%xmm2
  subsd  %xmm3,%xmm5
  addsd  %xmm8,%xmm3
  movsd  0x38(%rsp),%xmm6
  mov    0x50(%rsp),%rsi
  mov    0x88(%rsp),%r9d
  mulsd  %xmm5,%xmm3
  jae    L639
  xorpd  D50(%rip),%xmm1
L639:
  addsd  %xmm0,%xmm1
  mov    0x20(%rsp),%rax
  add    %r12,%rax
  divsd  %xmm1,%xmm2
  mov    %rax,0x88(%rsp)
  subsd  %xmm6,%xmm2
  mulsd  %xmm2,%xmm6
  addsd  %xmm3,%xmm6
  divsd  %xmm8,%xmm6
  cmp    0x60(%rsp),%r9d
  jg     L461
  mov    D60(%rip),%rax
  mov    %r13,0x20(%rsp)
  mov    %ebx,%r13d
  mov    %r14,%rbx
  mov    %rsi,%r14
  movq   %rax,%xmm2
  movq   %rax,%xmm4
  mov    %ebp,%eax
  mov    %r15,%rbp
  mov    %eax,%r15d
  nopl   0x0(%rax)
L669:
  movsd  0x8(%r12,%r14,8),%xmm7
  mov    0x20(%rsp),%rax
  movapd %xmm6,%xmm0
  movsd  %xmm8,0x48(%rsp)
  movsd  %xmm6,0x28(%rsp)
  mulsd  %xmm7,%xmm2
  movsd  0x8(%rax,%r14,8),%xmm3
  mov    $0x2,%eax
  mulsd  %xmm4,%xmm7
  movsd  %xmm3,0x40(%rsp)
  movapd %xmm2,%xmm1
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm7,0x38(%rsp)
  call   <pythag@plt>
  pxor   %xmm0,%xmm0
  xor    %r8d,%r8d
  movsd  0x30(%rsp),%xmm1
  cvtsi2sd %eax,%xmm0
  movsd  0x28(%rsp),%xmm6
  movsd  0x38(%rsp),%xmm7
  movapd %xmm1,%xmm2
  movsd  0x40(%rsp),%xmm3
  mov    0x68(%rsp),%rax
  movsd  0x48(%rsp),%xmm8
  movapd %xmm3,%xmm1
  lea    (%rax,%r14,8),%rax
  movsd  %xmm0,(%r12,%r14,8)
  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   %r13d,%r13d
  jle    L1000
  nopl   0x0(%rax,%rax,1)
L176:
  movsd  (%rax),%xmm3
  movapd %xmm4,%xmm5
  movapd %xmm2,%xmm6
  add    $0x1,%r8d
  movsd  0x8(%rax),%xmm0
  mulsd  %xmm3,%xmm5
  mulsd  %xmm0,%xmm6
  mulsd  %xmm2,%xmm3
  mulsd  %xmm4,%xmm0
  addsd  %xmm6,%xmm5
  subsd  %xmm3,%xmm0
  movsd  %xmm5,(%rax)
  movsd  %xmm0,0x8(%rax)
  add    %rbx,%rax
  cmp    %r8d,%r13d
  jne    L176
L1000:
  movapd %xmm9,%xmm0
  mov    $0x2,%eax
  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>
  pxor   %xmm0,%xmm0
  mov    0x20(%rsp),%rdi
  movsd  0x28(%rsp),%xmm9
  cvtsi2sd %eax,%xmm0
  test   %eax,%eax
  movsd  0x30(%rsp),%xmm1
  movsd  0x38(%rsp),%xmm7
  movsd  0x40(%rsp),%xmm8
  movsd  %xmm0,(%rdi,%r14,8)
  jne    L635
  movsd  0x48(%rsp),%xmm4
  movsd  0x50(%rsp),%xmm2
L81:
  movapd %xmm4,%xmm6
  movapd %xmm2,%xmm0
  mov    0x58(%rsp),%rax
  xor    %r8d,%r8d
  mulsd  %xmm7,%xmm6
  mulsd  %xmm8,%xmm0
  lea    (%rax,%r14,8),%rax
  mulsd  %xmm4,%xmm8
  mulsd  %xmm2,%xmm7
  addsd  %xmm0,%xmm6
  subsd  %xmm7,%xmm8
  test   %r15d,%r15d
  jle    L892
  nopl   0x0(%rax,%rax,1)
L1029:
  movsd  (%rax),%xmm1
  movapd %xmm4,%xmm3
  movapd %xmm2,%xmm5
  add    $0x1,%r8d
  movsd  0x8(%rax),%xmm0
  mulsd  %xmm1,%xmm3
  mulsd  %xmm0,%xmm5
  mulsd  %xmm2,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm5,%xmm3
  subsd  %xmm1,%xmm0
  movsd  %xmm3,(%rax)
  movsd  %xmm0,0x8(%rax)
  add    %rbp,%rax
  cmp    %r8d,%r15d
  jne    L1029
L892:
  add    $0x1,%r14
  cmp    %r14d,0x60(%rsp)
  jge    L669
  mov    %r15d,%eax
  mov    %rbx,%r14
  mov    %r13d,%ebx
  mov    0x20(%rsp),%r13
  mov    %rbp,%r15
  mov    %eax,%ebp
L461:
  mov    0x88(%rsp),%rax
  movq   $0x0,(%rax)
  mov    0x78(%rsp),%rax
  movsd  %xmm6,(%r12,%rax,8)
  movsd  %xmm8,0x0(%r13,%rax,8)
  jmp    L143
  cs nopw 0x0(%rax,%rax,1)
L635:
  movsd  D60(%rip),%xmm2
  movapd %xmm9,%xmm4
  divsd  %xmm0,%xmm2
  mulsd  %xmm2,%xmm4
  mulsd  %xmm1,%xmm2
  jmp    L81
L966:
  mov    0x78(%rsp),%rax
  movsd  0x0(%r13,%rax,8),%xmm3
  cmp    %eax,%r9d
  jne    L660
  movsd  0x70(%rsp),%xmm7
  mov    %r12,%rax
  mov    %r13,%r12
  mov    %ebx,%r13d
  mov    %rax,%rbx
  comisd %xmm3,%xmm7
  ja     L26
L23:
  subq   $0x1,0x78(%rsp)
  mov    0x78(%rsp),%rax
  test   %eax,%eax
  jns    L456
L257:
  xor    %eax,%eax
  jmp    L543
L923:
  ucomisd 0x70(%rsp),%xmm0
  jp     L153
  jne    L153
L689:
  cmp    %r12d,%ebp
  jle    L369
  movslq (%rsp),%rax
  mov    0x58(%rsp),%rdi
  add    %r12,%rax
  lea    (%rdi,%rax,8),%rax
  nopl   0x0(%rax)
L394:
  add    $0x1,%ebx
  movq   $0x0,(%rax)
  add    %r15,%rax
  cmp    %ebx,%ebp
  jne    L394
  jmp    L369
L900:
  mov    0x58(%rsp),%rcx
  movslq (%rsp),%rax
  movslq %r10d,%r11
  mov    %ebp,%r8d
  lea    (%rcx,%rax,8),%rsi
  movslq 0x30(%rsp),%rax
  lea    0x0(%r13,%rax,1),%rdx
  add    %r8,%rax
  lea    (%rcx,%rdx,8),%r14
  movslq %edi,%rcx
  add    %r13,%rax
  lea    (%rcx,%r11,1),%rdx
  mov    %rcx,0x90(%rsp)
  mov    0x68(%rsp),%rcx
  mov    %r14,0x48(%rsp)
  lea    (%rcx,%rdx,8),%rdx
  mov    0xa0(%rsp),%rcx
  lea    (%rcx,%rax,8),%rcx
  mov    %r14,%rax
L16:
  movsd  (%rax),%xmm0
  divsd  (%rsi),%xmm0
  divsd  %xmm1,%xmm0
  add    $0x8,%rax
  movsd  %xmm0,(%rdx)
  add    %r12,%rdx
  cmp    %rax,%rcx
  jne    L16
  mov    0x68(%rsp),%rdx
  lea    0x0(%r13,%r11,1),%rax
  lea    (%rdx,%rax,8),%rsi
  mov    %r8,%rax
  xor    %edx,%edx
  not    %rax
  mov    %rax,0x50(%rsp)
L78:
  mov    0x48(%rsp),%rax
  mov    %rsi,%r14
  pxor   %xmm1,%xmm1
  nopl   0x0(%rax)
L1069:
  movsd  (%rax),%xmm0
  mulsd  (%r14),%xmm0
  add    $0x8,%rax
  add    %r12,%r14
  addsd  %xmm0,%xmm1
  cmp    %rax,%rcx
  jne    L1069
  mov    %rsi,%rax
  mov    %r15d,%r14d
  nopw   0x0(%rax,%rax,1)
L1057:
  movsd  -0x8(%rax,%rdx,8),%xmm0
  add    $0x1,%r14d
  mulsd  %xmm1,%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    %r12,%rax
  cmp    %r14d,%ebx
  jne    L1057
  add    $0x8,%rsi
  sub    $0x1,%rdx
  cmp    %rdx,0x50(%rsp)
  jne    L78
  jmp    L379
L170:
  mov    %r10,%rax
  mov    %r8d,%edx
  pxor   %xmm1,%xmm1
L509:
  movsd  (%rax),%xmm0
  mov    %edx,%ecx
  add    $0x1,%edx
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  mulsd  %xmm0,%xmm0
  add    %r12,%rax
  addsd  %xmm0,%xmm1
  cmp    %r13d,%ecx
  jne    L509
  mov    0x40(%rsp),%eax
  mov    0x58(%rsp),%rcx
  movsd  0x70(%rsp),%xmm4
  sub    $0x1,%eax
  cltq
  lea    (%rcx,%rax,8),%r15
  movsd  (%r15),%xmm7
  movsd  %xmm7,0x98(%rsp)
  comisd %xmm4,%xmm7
  jb     L740
  ucomisd %xmm1,%xmm4
  ja     L1065
  movapd %xmm1,%xmm3
  sqrtsd %xmm3,%xmm3
L534:
  andpd  (%rsp),%xmm3
  xorpd  D50(%rip),%xmm3
L1033:
  movsd  0x98(%rsp),%xmm7
  movapd %xmm7,%xmm4
  subsd  %xmm3,%xmm7
  mulsd  %xmm3,%xmm4
  movsd  %xmm7,(%r15)
  subsd  %xmm1,%xmm4
  cmp    %esi,%ebx
  jle    L829
  mov    0x90(%rsp),%rax
  mov    0x58(%rsp),%rcx
  mov    $0xfffffffffffffffe,%rdx
  mov    %rsi,0xa8(%rsp)
  mov    $0x1,%r11d
  add    %rsi,%rax
  lea    (%rcx,%rax,8),%r15
  mov    0x38(%rsp),%eax
  mov    $0xffffffffffffffff,%rcx
  mov    %r15,%rsi
  sub    %rax,%rdx
  mov    %rdx,0x98(%rsp)
L89:
  mov    0x48(%rsp),%rax
  mov    %r8d,%edx
  pxor   %xmm1,%xmm1
L788:
  movsd  (%rax),%xmm0
  mulsd  (%rax,%r11,8),%xmm0
  mov    %edx,%r15d
  add    %r12,%rax
  add    $0x1,%edx
  addsd  %xmm0,%xmm1
  cmp    %r13d,%r15d
  jne    L788
  divsd  %xmm4,%xmm1
  mov    %rsi,%rax
  mov    %r8d,%edx
  cs nopw 0x0(%rax,%rax,1)
L853:
  movsd  (%rax,%rcx,8),%xmm0
  mov    %edx,%r15d
  add    $0x1,%edx
  mulsd  %xmm1,%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    %r12,%rax
  cmp    %r13d,%r15d
  jne    L853
  add    $0x8,%rsi
  sub    $0x1,%rcx
  add    $0x1,%r11
  cmp    0x98(%rsp),%rcx
  jne    L89
  mov    0xa8(%rsp),%rsi
L829:
  movsd  (%r10),%xmm0
  mov    %r8d,%eax
  add    $0x1,%r8d
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,(%r10)
  add    %r12,%r10
  cmp    %r13d,%eax
  jne    L829
  mulsd  %xmm3,%xmm2
  movsd  %xmm2,(%r9)
  jmp    L425
L153:
  movsd  D60(%rip),%xmm2
  divsd  %xmm0,%xmm2
  jmp    L1008
L740:
  movsd  0x70(%rsp),%xmm7
  ucomisd %xmm1,%xmm7
  ja     L148
  movapd %xmm1,%xmm3
  sqrtsd %xmm3,%xmm3
L679:
  andpd  (%rsp),%xmm3
  jmp    L1033
L26:
  xorpd  D50(%rip),%xmm3
  mov    0x78(%rsp),%rax
  movsd  %xmm3,(%r12,%rax,8)
  test   %r13d,%r13d
  jle    L23
  mov    0x68(%rsp),%rcx
  xor    %edi,%edi
  lea    (%rcx,%rax,8),%rax
L1035:
  movsd  (%rax),%xmm0
  xorpd  D50(%rip),%xmm0
  add    $0x1,%edi
  movsd  %xmm0,(%rax)
  add    %r14,%rax
  cmp    %edi,%r13d
  jne    L1035
  jmp    L23
L17:
  mov    %r15,%rax
  pxor   %xmm1,%xmm1
L318:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  divsd  %xmm3,%xmm0
  movsd  %xmm0,-0x8(%rax)
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  cmp    %r13,%rax
  jne    L318
  movslq 0x40(%rsp),%rax
  mov    0x58(%rsp),%rcx
  movsd  0x70(%rsp),%xmm4
  lea    (%rcx,%rax,8),%rax
  movsd  (%rax),%xmm7
  mov    %rax,0x30(%rsp)
  movsd  %xmm7,0x98(%rsp)
  comisd %xmm4,%xmm7
  jb     L372
  ucomisd %xmm1,%xmm4
  ja     L133
  movapd %xmm1,%xmm2
  sqrtsd %xmm2,%xmm2
L462:
  andpd  (%rsp),%xmm2
  xorpd  D50(%rip),%xmm2
L677:
  movsd  0x98(%rsp),%xmm7
  mov    0x30(%rsp),%rax
  mov    0x90(%rsp),%rcx
  movapd %xmm7,%xmm4
  subsd  %xmm2,%xmm7
  mulsd  %xmm2,%xmm4
  movsd  %xmm7,(%rax)
  mov    0x58(%rsp),%rax
  movapd %xmm7,%xmm0
  lea    (%rax,%rcx,8),%r8
  mov    %rsi,%rax
  subsd  %xmm1,%xmm4
  jmp    L561
L207:
  movsd  (%r8,%rax,8),%xmm0
L561:
  divsd  %xmm4,%xmm0
  movsd  %xmm0,0x0(%rbp,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%ebx
  jg     L207
  mov    0x28(%rsp),%r10d
  mov    0x58(%rsp),%rcx
  add    0x20(%rsp),%r10d
  mov    0x50(%rsp),%edx
  movslq %r10d,%rax
  lea    (%rcx,%rax,8),%rcx
L123:
  mov    %rsi,%rax
  pxor   %xmm1,%xmm1
  nopl   0x0(%rax,%rax,1)
L831:
  movsd  (%rcx,%rax,8),%xmm0
  mulsd  (%r8,%rax,8),%xmm0
  add    $0x1,%rax
  addsd  %xmm0,%xmm1
  cmp    %eax,%ebx
  jg     L831
  mov    %rsi,%rax
  nopw   0x0(%rax,%rax,1)
L186:
  movsd  0x0(%rbp,%rax,8),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  (%rcx,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%ebx
  jg     L186
  add    $0x1,%edx
  add    %r12,%rcx
  cmp    %edx,%edi
  jne    L123
L437:
  movsd  (%r15),%xmm0
  add    $0x8,%r15
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,-0x8(%r15)
  cmp    %r13,%r15
  jne    L437
  movapd (%rsp),%xmm7
  movsd  (%r9),%xmm0
  movsd  0x0(%rbp,%r14,8),%xmm1
  andpd  %xmm7,%xmm0
  andpd  %xmm7,%xmm1
  movsd  0x18(%rsp),%xmm7
  addsd  %xmm1,%xmm0
  maxsd  %xmm0,%xmm7
  movsd  %xmm7,0x18(%rsp)
  jmp    L194
L372:
  movsd  0x70(%rsp),%xmm7
  ucomisd %xmm1,%xmm7
  ja     L812
  movapd %xmm1,%xmm2
  sqrtsd %xmm2,%xmm2
L572:
  andpd  (%rsp),%xmm2
  jmp    L677
L68:
  pxor   %xmm7,%xmm7
  movsd  %xmm7,0x70(%rsp)
  movsd  %xmm7,0x18(%rsp)
  jmp    L54
L396:
  movapd (%rsp),%xmm7
  movsd  (%r9),%xmm0
  movsd  0x0(%rbp,%r14,8),%xmm1
  andpd  %xmm7,%xmm0
  andpd  %xmm7,%xmm1
  movsd  0x18(%rsp),%xmm7
  addsd  %xmm1,%xmm0
  comisd %xmm0,%xmm7
  ja     L593
  pxor   %xmm3,%xmm3
  jmp    L841
L812:
  movapd %xmm1,%xmm0
  mov    %edi,0xc8(%rsp)
  mov    %r9,0xc0(%rsp)
  mov    %rsi,0xb8(%rsp)
  movsd  %xmm3,0xb0(%rsp)
  movsd  %xmm1,0xa8(%rsp)
  call   <sqrt@plt>
  mov    0xc8(%rsp),%edi
  mov    0xc0(%rsp),%r9
  mov    0xb8(%rsp),%rsi
  movapd %xmm0,%xmm2
  movsd  0xb0(%rsp),%xmm3
  movsd  0xa8(%rsp),%xmm1
  jmp    L572
L1065:
  movapd %xmm1,%xmm0
  mov    %edi,0xdc(%rsp)
  mov    %r9,0xd0(%rsp)
  mov    %rsi,0xc8(%rsp)
  mov    %r10,0xc0(%rsp)
  mov    %r8d,0xb0(%rsp)
  movsd  %xmm2,0xb8(%rsp)
  movsd  %xmm1,0xa8(%rsp)
  call   <sqrt@plt>
  mov    0xdc(%rsp),%edi
  mov    0xd0(%rsp),%r9
  mov    0xc8(%rsp),%rsi
  mov    0xc0(%rsp),%r10
  movapd %xmm0,%xmm3
  movsd  0xb8(%rsp),%xmm2
  mov    0xb0(%rsp),%r8d
  movsd  0xa8(%rsp),%xmm1
  jmp    L534
L148:
  movapd %xmm1,%xmm0
  mov    %edi,0xdc(%rsp)
  mov    %r9,0xd0(%rsp)
  mov    %rsi,0xc8(%rsp)
  mov    %r10,0xc0(%rsp)
  mov    %r8d,0xb0(%rsp)
  movsd  %xmm2,0xb8(%rsp)
  movsd  %xmm1,0xa8(%rsp)
  call   <sqrt@plt>
  mov    0xdc(%rsp),%edi
  mov    0xd0(%rsp),%r9
  mov    0xc8(%rsp),%rsi
  mov    0xc0(%rsp),%r10
  movapd %xmm0,%xmm3
  movsd  0xb8(%rsp),%xmm2
  mov    0xb0(%rsp),%r8d
  movsd  0xa8(%rsp),%xmm1
  jmp    L679
L133:
  movapd %xmm1,%xmm0
  mov    %edi,0xc8(%rsp)
  mov    %r9,0xc0(%rsp)
  mov    %rsi,0xb8(%rsp)
  movsd  %xmm3,0xb0(%rsp)
  movsd  %xmm1,0xa8(%rsp)
  call   <sqrt@plt>
  mov    0xc8(%rsp),%edi
  mov    0xc0(%rsp),%r9
  mov    0xb8(%rsp),%rsi
  movapd %xmm0,%xmm2
  movsd  0xb0(%rsp),%xmm3
  movsd  0xa8(%rsp),%xmm1
  jmp    L462
```

tool_calls: name: extract_symbolic_data
arguments: label: D43
type_infered: qword
name: extract_symbolic_data
arguments: label: D60
type_infered: f64
name: extract_symbolic_data
arguments: label: D50
type_infered: f32[4]

tool: D43: 9223372036854775807

tool: D60: 1

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

assistant: ```c
int relax_svd(double *U, int m, int n, int ustride, double *V, int vstride, double *w, double *rv1) {
  int flag, i, its, j, o, k, l, nm;
  double anorm, c, f, g, h, s, scale, x, y, z;
  if (m < n) {
    return -1;
  }
  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(U[k * ustride + i]);
      if (scale != 0.0) {
        for (k = i; k < m; k++) {
          U[k * ustride + i] /= scale;
          s += U[k * ustride + i] * U[k * ustride + i];
        }
        f = U[i * ustride + i];
        g = -((f) >= 0.0 ? fabs(sqrt(s)) : -fabs(sqrt(s)));
        h = f * g - s;
        U[i * ustride + i] = f - g;
        for (j = l; j < n; j++) {
          for (s = 0.0, k = i; k < m; k++)
            s += U[k * ustride + i] * U[k * ustride + j];
          f = s / h;
          for (k = i; k < m; k++)
            U[k * ustride + j] += f * U[k * ustride + i];
        }
        for (k = i; k < m; k++)
          U[k * ustride + 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(U[i * ustride + k]);
      if (scale) {
        for (k = l; k < n; k++) {
          U[i * ustride + k] /= scale;
          s += U[i * ustride + k] * U[i * ustride + k];
        }
        f = U[i * ustride + l];
        g = -((f) >= 0.0 ? fabs(sqrt(s)) : -fabs(sqrt(s)));
        h = f * g - s;
        U[i * ustride + l] = f - g;
        for (k = l; k < n; k++)
          rv1[k] = U[i * ustride + k] / h;
        for (j = l; j < m; j++) {
          for (s = 0.0, k = l; k < n; k++)
            s += U[j * ustride + k] * U[i * ustride + k];
          for (k = l; k < n; k++)
            U[j * ustride + k] += s * rv1[k];
        }
        for (k = l; k < n; k++)
          U[i * ustride + k] *= scale;
      }
    }
    double tmp = fabs(w[i]) + fabs(rv1[i]);
    anorm = anorm > tmp ? anorm : tmp;
  }
  for (i = n - 1; i >= 0; i--) {
    if (i < n - 1) {
      if (g != 0.0) {
        for (j = l; j < n; j++)
          V[j * vstride + i] = (U[i * ustride + j] / U[i * ustride + l]) / g;
        for (j = l; j < n; j++) {
          for (s = 0.0, k = l; k < n; k++)
            s += U[i * ustride + k] * V[k * vstride + j];
          for (k = l; k < n; k++)
            V[k * vstride + j] += s * V[k * vstride + i];
        }
      }
      for (j = l; j < n; j++)
        V[i * vstride + j] = V[j * vstride + i] = 0.0;
    }
    V[i * vstride + i] = 1.0;
    g = rv1[i];
    l = i;
  }
  for (i = n - 1; i >= 0; i--) {
    l = i + 1;
    g = w[i];
    for (j = l; j < n; j++)
      U[i * ustride + j] = 0.0;
    if (g != 0.0) {
      g = 1.0 / g;
      for (j = l; j < n; j++) {
        for (s = 0.0, k = l; k < m; k++)
          s += U[k * ustride + i] * U[k * ustride + j];
        f = (s / U[i * ustride + i]) * g;
        for (k = i; k < m; k++)
          U[k * ustride + j] += f * U[k * ustride + i];
      }
      for (j = i; j < m; j++)
        U[j * ustride + i] *= g;
    } else {
      for (j = i; j < m; j++)
        U[j * ustride + i] = 0.0;
    }
    U[i * ustride + i] += 1.0;
  }
  for (k = n - 1; k >= 0; k--) {
    for (its = 1; its <= 30; its++) {
      flag = 1;
      for (l = k; l >= 0; l--) {
        nm = l - 1;
        if (fabs(rv1[l]) + anorm == anorm) {
          flag = 0;
          break;
        }
        if (fabs(w[nm]) + anorm == anorm)
          break;
      }
      if (flag != 0) {
        c = 0.0;
        s = 1.0;
        for (i = l; i < k; i++) {
          f = s * rv1[i];
          rv1[i] = c * rv1[i];
          if (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 = 1; j < m; j++) {
            y = U[j * ustride + nm];
            z = U[j * ustride + i];
            U[j * ustride + nm] = y * c + z * s;
            U[j * ustride + 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 * vstride + k] = -V[j * vstride + k];
        }
        break;
      }
      if (its == 30)
        return -1;
      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.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 (o = 0; o < n; o++) {
          x = V[o * vstride + j];
          z = V[o * vstride + i];
          V[o * vstride + j] = x * c + z * s;
          V[o * vstride + i] = z * c - x * s;
        }
        z = pythag(f, h);
        w[j] = z;
        if (z != 0.0) {
          z = 1.0 / z;
          c = f * z;
          s = h * z;
        }
        f = c * g + s * y;
        x = c * y - s * g;
        for (o = 0; o < m; o++) {
          y = U[o * ustride + j];
          z = U[o * ustride + i];
          U[o * ustride + j] = y * c + z * s;
          U[o * ustride + i] = z * c - y * s;
        }
      }
      rv1[l] = 0.0;
      rv1[k] = f;
      w[k] = x;
    }
  }
  return 0;
}
```