user: 下面的汇编指令对应的C源代码是什么:

```asm
<ssvdinv>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %rdi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  movslq %edx,%rbx
  lea    0x0(,%rbx,4),%ebp
  mov    %rbx,%r14
  sub    $0xe8,%rsp
  mov    %rdi,0x58(%rsp)
  mov    %ebp,%edi
  mov    %rsi,0xd8(%rsp)
  mov    %rbx,0x78(%rsp)
  call   <malloc@plt>
  mov    %ebp,%edi
  mov    %rax,0x20(%rsp)
  mov    %rax,%r15
  call   <malloc@plt>
  mov    %ebx,%ecx
  imul   %ebx,%ecx
  mov    %rax,0xc0(%rsp)
  lea    0x0(,%rcx,4),%edi
  mov    %ecx,0xc8(%rsp)
  call   <malloc@plt>
  mov    %ebp,%edi
  mov    %rax,0x68(%rsp)
  call   <malloc@plt>
  lea    -0x1(%rbx),%edi
  mov    %rax,0x70(%rsp)
  mov    %edi,0xa0(%rsp)
  test   %ebx,%ebx
  jle    L751
  mov    %r15,%rcx
  lea    0x1(%rbx),%edx
  sub    $0x4,%rax
  mov    %r13,0x40(%rsp)
  mov    %edx,0x88(%rsp)
  sub    %rcx,%rax
  lea    -0x2(%rbx),%edi
  mov    %ebx,%r13d
  lea    0x4(,%rbx,4),%rdx
  pxor   %xmm7,%xmm7
  movl   $0x1,0x80(%rsp)
  lea    0x0(,%rbx,4),%r12
  mov    %rdx,0x90(%rsp)
  movaps %xmm7,%xmm2
  movaps %xmm7,%xmm3
  movq   $0x1,0x18(%rsp)
  movq   $0x0,0x10(%rsp)
  mov    %rax,0x98(%rsp)
  mov    %edi,0xb8(%rsp)
  lea    0x4(%r15),%rdi
  xor    %r15d,%r15d
  mov    %r15d,0x28(%rsp)
  mov    %r12,%r15
  mov    %rdi,%r12
  movss  %xmm7,0xc(%rsp)
  movss  %xmm7,0x38(%rsp)
L641:
  mulss  %xmm3,%xmm2
  mov    0x18(%rsp),%edi
  mov    0x10(%rsp),%rax
  pxor   %xmm3,%xmm3
  mov    0x40(%rsp),%r14
  mov    %edi,0x50(%rsp)
  lea    -0x4(%r12),%rdi
  mov    %eax,%ebx
  mov    %eax,%ebp
  mov    %rdi,0x48(%rsp)
  movss  %xmm2,-0x4(%r12)
  nopw   0x0(%rax,%rax,1)
L23:
  movss  (%r14),%xmm0
  movss  %xmm3,(%rsp)
  add    $0x1,%ebp
  add    %r15,%r14
  call   <fabs@plt>
  movss  (%rsp),%xmm3
  addss  %xmm0,%xmm3
  cmp    %ebp,%r13d
  jg     L23
  movss  0xc(%rsp),%xmm6
  ucomiss %xmm6,%xmm3
  jp     L290
  mov    0x28(%rsp),%eax
  movaps %xmm6,%xmm4
  lea    (%rax,%r13,1),%eax
  mov    %eax,0x30(%rsp)
  jne    L290
L120:
  mulss  %xmm3,%xmm4
  mov    0x98(%rsp),%rax
  mov    0x18(%rsp),%rcx
  add    %r12,%rax
  mov    %rax,0x60(%rsp)
  movaps %xmm4,%xmm0
  movss  %xmm4,(%rax)
  cmp    %ecx,%r13d
  je     L38
  jle    L38
  movslq 0x28(%rsp),%rdi
  mov    0xb8(%rsp),%edx
  pxor   %xmm2,%xmm2
  mov    %rdi,%rax
  mov    %rdi,0x28(%rsp)
  add    %rcx,%rax
  mov    0x58(%rsp),%rcx
  mov    %rdx,0xa8(%rsp)
  lea    (%rcx,%rax,4),%rbp
  mov    %rdi,%rax
  add    %rdx,%rax
  add    0x10(%rsp),%rax
  mov    %rbp,%r14
  lea    0x8(%rcx,%rax,4),%rbx
  nopl   (%rax)
L1143:
  movss  (%r14),%xmm0
  movss  %xmm2,(%rsp)
  add    $0x4,%r14
  call   <fabs@plt>
  movss  (%rsp),%xmm2
  addss  %xmm0,%xmm2
  cmp    %rbx,%r14
  jne    L1143
  movss  0xc(%rsp),%xmm3
  ucomiss %xmm3,%xmm2
  jp     L54
  jne    L54
L1083:
  mov    0x60(%rsp),%rax
  movss  (%rax),%xmm0
L248:
  movss  %xmm2,0x60(%rsp)
  add    $0x4,%r12
  movss  %xmm3,0x50(%rsp)
  call   <fabs@plt>
  mov    0x48(%rsp),%rax
  movss  %xmm0,(%rsp)
  movss  (%rax),%xmm0
  call   <fabs@plt>
  addss  (%rsp),%xmm0
  mov    0x30(%rsp),%edi
  movss  0x38(%rsp),%xmm6
  addq   $0x1,0x10(%rsp)
  mov    0x88(%rsp),%ecx
  mov    0x10(%rsp),%rax
  mov    %edi,0x28(%rsp)
  maxss  %xmm0,%xmm6
  add    %ecx,0x80(%rsp)
  mov    0x90(%rsp),%rcx
  addq   $0x1,0x18(%rsp)
  movss  0x50(%rsp),%xmm3
  subl   $0x1,0xb8(%rsp)
  movss  0x60(%rsp),%xmm2
  add    %rcx,0x40(%rsp)
  cmp    0x78(%rsp),%rax
  movss  %xmm6,0x38(%rsp)
  jne    L641
  movss  D181(%rip),%xmm6
  movss  0x38(%rsp),%xmm7
  mov    %r13d,%r14d
  mov    %r13d,%ebp
  movss  %xmm6,0xbc(%rsp)
  mulss  %xmm6,%xmm7
  movss  %xmm7,0x10(%rsp)
L1088:
  movslq 0xa0(%rsp),%rcx
  mov    %r14d,%eax
  mov    0x68(%rsp),%r15
  not    %eax
  mov    0x20(%rsp),%r12
  mov    %ecx,%edx
  mov    %eax,0x48(%rsp)
  cltq
  mov    %ecx,%edi
  imul   %r14d,%edx
  shl    $0x2,%rax
  lea    0x0(,%rcx,4),%rbx
  mov    %rcx,0xb0(%rsp)
  mov    %rax,0x18(%rsp)
  add    %rbx,%r12
  lea    0x4(%rbx),%rsi
  mov    %rbx,0x40(%rsp)
  movslq %edx,%rax
  mov    %edx,(%rsp)
  mov    %edx,%r13d
  mov    %rax,0x60(%rsp)
  add    %rcx,%rax
  shl    $0x2,%rax
  mov    %edx,0xb8(%rsp)
  mov    %rax,0x50(%rsp)
  add    %rax,%r15
  mov    0x78(%rsp),%rax
  mov    %rdi,0xd0(%rsp)
  sub    $0x2,%rax
  sub    %rdi,%rax
  mov    %rax,0xa8(%rsp)
  mov    0x58(%rsp),%rax
  add    $0x4,%rax
  mov    %rax,0x80(%rsp)
  jmp    L832
L243:
  movl   $0x3f800000,(%r15)
  sub    $0x1,%rcx
  movss  (%r12),%xmm3
  sub    %r14d,%r13d
  add    0x18(%rsp),%r15
  sub    $0x4,%r12
  cmp    %rcx,0xa8(%rsp)
  je     L571
L832:
  lea    0x1(%rcx),%eax
  mov    %ebp,%edi
  mov    %ecx,%ebp
  cmp    %r14d,%eax
  jge    L243
  ucomiss 0xc(%rsp),%xmm3
  jp     L540
  jne    L540
  movslq %edi,%rax
  mov    %rax,0x30(%rsp)
  movslq %r13d,%rax
  mov    %rax,0x38(%rsp)
  cmp    %edi,%r14d
  jle    L243
L686:
  mov    %r14d,%eax
  mov    0x68(%rsp),%rbx
  mov    0x38(%rsp),%rdx
  imul   %edi,%eax
  add    0x30(%rsp),%rdx
  lea    (%rbx,%rdx,4),%rdx
  cltq
  add    %rcx,%rax
  lea    (%rbx,%rax,4),%rax
  nopw   0x0(%rax,%rax,1)
L767:
  add    $0x1,%edi
  movl   $0x0,(%rax)
  add    $0x4,%rdx
  add    %rsi,%rax
  movl   $0x0,-0x4(%rdx)
  cmp    %edi,%r14d
  jg     L767
  jmp    L243
L38:
  movss  0xc(%rsp),%xmm3
  movaps %xmm3,%xmm2
  jmp    L248
L290:
  movslq 0x28(%rsp),%rcx
  mov    0x10(%rsp),%rax
  mov    %ebx,%edx
  pxor   %xmm1,%xmm1
  mov    0x58(%rsp),%rdi
  add    %rcx,%rax
  lea    (%rdi,%rax,4),%rbp
  mov    %rbp,%rax
L953:
  movss  (%rax),%xmm0
  add    $0x1,%edx
  divss  %xmm3,%xmm0
  movss  %xmm0,(%rax)
  mulss  %xmm0,%xmm0
  add    %r15,%rax
  addss  %xmm0,%xmm1
  cmp    %edx,%r13d
  jg     L953
  mov    0x80(%rsp),%eax
  mov    0x58(%rsp),%rdi
  pxor   %xmm0,%xmm0
  mov    %rcx,0xb0(%rsp)
  movss  %xmm3,0xa8(%rsp)
  cvtss2sd %xmm1,%xmm0
  sub    $0x1,%eax
  movss  %xmm1,0x60(%rsp)
  cltq
  lea    (%rdi,%rax,4),%rsi
  movss  (%rsi),%xmm6
  comiss 0xc(%rsp),%xmm6
  mov    %rsi,0x30(%rsp)
  movss  %xmm6,(%rsp)
  jb     L528
  call   <sqrt@plt>
  call   <fabs@plt>
  mov    0x30(%rsp),%rsi
  movss  0x60(%rsp),%xmm1
  movaps %xmm0,%xmm4
  mov    0xb0(%rsp),%rcx
  xorps  D83(%rip),%xmm4
  movss  0xa8(%rsp),%xmm3
L667:
  movss  (%rsp),%xmm6
  movaps %xmm6,%xmm5
  subss  %xmm4,%xmm6
  mulss  %xmm4,%xmm5
  movaps %xmm6,%xmm0
  movss  %xmm6,(%rsi)
  subss  %xmm1,%xmm5
  cmp    0x18(%rsp),%r13d
  jle    L887
  mov    0x18(%rsp),%rdi
  mov    0x58(%rsp),%rdx
  mov    0x10(%rsp),%r8
  mov    0x50(%rsp),%r11d
  lea    (%rcx,%rdi,1),%rax
  lea    (%rdx,%rax,4),%r10
  mov    0x28(%rsp),%eax
  sub    %rdi,%r8
  sub    %rcx,%r8
  add    %r13d,%eax
  mov    %eax,0x30(%rsp)
  cltq
  add    %rax,%r8
L341:
  mov    %r10,%rax
  movaps %xmm0,%xmm1
  pxor   %xmm2,%xmm2
  mov    %ebx,%edx
  jmp    L1130
L530:
  movss  (%rax,%r8,4),%xmm1
  mov    %ecx,%edx
  add    %r15,%rax
L1130:
  mulss  (%rax),%xmm1
  lea    0x1(%rdx),%ecx
  addss  %xmm1,%xmm2
  cmp    %ecx,%r13d
  jne    L530
  divss  %xmm5,%xmm2
  mov    %r10,%rax
  mov    %ebx,%ecx
  jmp    L187
  cs nopw 0x0(%rax,%rax,1)
L301:
  movss  (%rax,%r8,4),%xmm0
  mov    %edi,%ecx
  add    %r15,%rax
L187:
  mulss  %xmm2,%xmm0
  lea    0x1(%rcx),%edi
  addss  (%rax),%xmm0
  movss  %xmm0,(%rax)
  cmp    %ecx,%edx
  jne    L301
  movss  (%rsi),%xmm0
  lea    0x1(%r11),%eax
  add    $0x4,%r10
  sub    $0x1,%r8
  cmp    %edx,%r11d
  je     L814
  mov    %eax,%r11d
  jmp    L341
L833:
  movss  0x0(%rbp),%xmm0
L814:
  mulss  %xmm3,%xmm0
  add    $0x1,%ebx
  movss  %xmm0,0x0(%rbp)
  add    %r15,%rbp
  cmp    %ebx,%r13d
  jg     L833
  jmp    L120
L54:
  mov    %rbp,%rax
  pxor   %xmm1,%xmm1
L587:
  movss  (%rax),%xmm0
  add    $0x4,%rax
  divss  %xmm2,%xmm0
  movss  %xmm0,-0x4(%rax)
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  cmp    %rbx,%rax
  jne    L587
  movslq 0x80(%rsp),%rax
  mov    0x58(%rsp),%rdi
  pxor   %xmm0,%xmm0
  movss  %xmm2,0xbc(%rsp)
  movss  %xmm1,0xb0(%rsp)
  cvtss2sd %xmm1,%xmm0
  lea    (%rdi,%rax,4),%r14
  movss  (%r14),%xmm7
  comiss 0xc(%rsp),%xmm7
  movss  %xmm7,(%rsp)
  jb     L567
  call   <sqrt@plt>
  call   <fabs@plt>
  movss  0xb0(%rsp),%xmm1
  movss  0xbc(%rsp),%xmm2
  movaps %xmm0,%xmm3
  xorps  D83(%rip),%xmm3
L1239:
  movss  (%rsp),%xmm6
  mov    0x28(%rsp),%r9
  mov    0x20(%rsp),%rdi
  mov    0xa8(%rsp),%rax
  movaps %xmm6,%xmm4
  subss  %xmm3,%xmm6
  shl    $0x2,%r9
  add    0x10(%rsp),%rax
  mulss  %xmm3,%xmm4
  mov    %r9,%rdx
  lea    0x8(%rdi,%rax,4),%r8
  mov    %r12,%rax
  sub    %rdi,%rdx
  add    0x58(%rsp),%rdx
  movaps %xmm6,%xmm0
  movss  %xmm6,(%r14)
  subss  %xmm1,%xmm4
  jmp    L415
L850:
  movss  (%rdx,%rax,1),%xmm0
L415:
  divss  %xmm4,%xmm0
  add    $0x4,%rax
  movss  %xmm0,-0x4(%rax)
  cmp    %r8,%rax
  jne    L850
  movslq 0x30(%rsp),%rcx
  mov    0x58(%rsp),%rax
  mov    0x50(%rsp),%esi
  mov    0x18(%rsp),%rdi
  shl    $0x2,%rcx
  add    %rax,%r9
  mov    %rcx,%rdx
  sub    0x20(%rsp),%rdx
  add    %rax,%rcx
  add    %rax,%rdx
L1255:
  mov    %rdi,%rax
  pxor   %xmm1,%xmm1
  nopl   0x0(%rax,%rax,1)
L668:
  movss  (%rcx,%rax,4),%xmm0
  mulss  (%r9,%rax,4),%xmm0
  add    $0x1,%rax
  addss  %xmm0,%xmm1
  cmp    %eax,%r13d
  jg     L668
  mov    %r12,%rax
  nopl   0x0(%rax,%rax,1)
L1072:
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  addss  (%rdx,%rax,1),%xmm0
  movss  %xmm0,(%rdx,%rax,1)
  add    $0x4,%rax
  cmp    %r8,%rax
  jne    L1072
  add    $0x1,%esi
  add    %r15,%rdx
  add    %r15,%rcx
  cmp    %esi,%r13d
  jne    L1255
L799:
  movss  0x0(%rbp),%xmm0
  add    $0x4,%rbp
  mulss  %xmm2,%xmm0
  movss  %xmm0,-0x4(%rbp)
  cmp    %rbx,%rbp
  jne    L799
  jmp    L1083
L571:
  mov    0x40(%rsp),%rax
  add    0x70(%rsp),%rax
  mov    0x58(%rsp),%rdi
  mov    %rax,0x30(%rsp)
  mov    0x50(%rsp),%rax
  movslq 0xc8(%rsp),%rdx
  mov    0xb0(%rsp),%rcx
  mov    0x78(%rsp),%r15
  add    %rdi,%rax
  mov    %rax,0x28(%rsp)
  mov    %r14d,%eax
  shl    $0x2,%r15
  mov    %rcx,%rbp
  neg    %eax
  cltq
  lea    -0x4(,%rax,4),%rax
  mov    %rax,0x50(%rsp)
  lea    (%rcx,%rdx,1),%rax
  lea    (%rdi,%rax,4),%r8
  mov    0xb8(%rsp),%eax
  lea    -0x2(%r14),%edi
  add    0xa0(%rsp),%eax
  mov    %eax,0x38(%rsp)
  mov    0x60(%rsp),%rax
  mov    %r8,%r12
  mov    %edi,0x60(%rsp)
  sub    %rdx,%rax
  lea    0x4(,%rax,4),%rax
  mov    %rax,0x78(%rsp)
L469:
  mov    0x30(%rsp),%rax
  lea    0x1(%rbp),%r13d
  mov    %ebp,%ebx
  movss  (%rax),%xmm0
  cmp    %r13d,%r14d
  jle    L672
  mov    0x60(%rsp),%eax
  xor    %esi,%esi
  movss  %xmm0,0x40(%rsp)
  sub    %ebp,%eax
  lea    0x4(,%rax,4),%rdx
  mov    0x78(%rsp),%rax
  lea    (%rax,%r12,1),%rdi
  call   <memset@plt>
  movss  0x40(%rsp),%xmm0
  ucomiss 0xc(%rsp),%xmm0
  jp     L1228
  je     L1117
L1228:
  movslq 0x38(%rsp),%rax
  mov    0x58(%rsp),%rdi
  mov    %rbp,%rdx
  mov    %r13d,%r8d
  movss  D130(%rip),%xmm2
  mov    $0x1,%r9d
  lea    (%rdi,%rax,4),%r11
  movslq (%rsp),%rax
  divss  %xmm0,%xmm2
  mov    %rax,%rcx
  lea    0x1(%rbp,%rax,1),%rax
  add    %r14d,%ecx
  sub    %rax,%rdx
  lea    (%rdi,%rax,4),%r10
  movslq %ecx,%rcx
  add    %rdx,%rcx
  nopl   0x0(%rax)
L532:
  mov    %r12,%rax
  mov    %r13d,%esi
  pxor   %xmm1,%xmm1
  nopw   0x0(%rax,%rax,1)
L42:
  movss  (%rax),%xmm0
  mulss  (%rax,%r9,4),%xmm0
  mov    %esi,%edx
  lea    0x1(%rsi),%esi
  add    %r15,%rax
  addss  %xmm0,%xmm1
  cmp    %esi,%r14d
  jne    L42
  movss  (%r11),%xmm0
  divss  %xmm0,%xmm1
  mulss  %xmm2,%xmm1
  cmp    %ebp,%r14d
  jle    L48
  mov    %r10,%rax
  mov    %ebx,%esi
  jmp    L390
  nopl   0x0(%rax)
L125:
  movss  (%rax,%rcx,4),%xmm0
  mov    %edi,%esi
  add    %r15,%rax
L390:
  mulss  %xmm1,%xmm0
  lea    0x1(%rsi),%edi
  addss  (%rax),%xmm0
  movss  %xmm0,(%rax)
  cmp    %edx,%esi
  jne    L125
L48:
  lea    0x1(%r8),%eax
  add    $0x4,%r10
  sub    $0x1,%rcx
  add    $0x1,%r9
  cmp    %edx,%r8d
  je     L629
  mov    %eax,%r8d
  jmp    L532
L672:
  ucomiss 0xc(%rsp),%xmm0
  jp     L585
  jne    L585
L1117:
  cmp    %ebp,%r14d
  jle    L1212
  movslq (%rsp),%rax
  mov    0x58(%rsp),%rdi
  add    %rbp,%rax
  lea    (%rdi,%rax,4),%rax
  nopl   0x0(%rax)
L1203:
  add    $0x1,%ebx
  movl   $0x0,(%rax)
  add    %r15,%rax
  cmp    %ebx,%r14d
  jne    L1203
L1212:
  mov    0x28(%rsp),%rax
  mov    0x18(%rsp),%rdi
  sub    $0x1,%rbp
  movss  D130(%rip),%xmm0
  subq   $0x4,0x30(%rsp)
  addss  (%rax),%xmm0
  sub    %r14d,(%rsp)
  add    0x50(%rsp),%r12
  movss  %xmm0,(%rax)
  add    %rdi,%rax
  mov    0x48(%rsp),%edi
  mov    %rax,0x28(%rsp)
  add    %edi,0x38(%rsp)
  cmp    %rbp,0xa8(%rsp)
  jne    L469
L523:
  mov    0xb0(%rsp),%rax
  mov    0x70(%rsp),%rcx
  mov    %r15,%rbp
  mov    %r14d,%r15d
  movl   $0x1e,0x90(%rsp)
  mov    %eax,%edi
  mov    %rcx,%rdx
  add    $0x1,%edi
  mov    %edi,0xb8(%rsp)
  lea    0x0(,%rax,4),%rdi
  sub    $0x1,%eax
  mov    %eax,0x60(%rsp)
  cltq
  add    %rdi,%rdx
  shl    $0x2,%rax
  mov    %rdi,0xc8(%rsp)
  add    %rax,%rcx
  mov    %rdx,0x80(%rsp)
  mov    %rcx,0x98(%rsp)
  mov    0x20(%rsp),%rcx
  add    %rcx,%rax
  add    %rdi,%rcx
  mov    %rax,0xa0(%rsp)
  mov    %rcx,0x78(%rsp)
  nopl   0x0(%rax)
L519:
  mov    0x70(%rsp),%rax
  mov    0x78(%rsp),%r13
  mov    %rbp,(%rsp)
  mov    %r15d,0x18(%rsp)
  mov    0xb8(%rsp),%ebx
  sub    $0x4,%rax
  mov    0x20(%rsp),%r15
  mov    %rax,%rbp
  jmp    L694
L684:
  test   %ebx,%ebx
  je     L7
  mov    %rbp,%rax
  sub    %r15,%rax
  movss  (%rax,%r13,1),%xmm0
  sub    $0x4,%r13
  call   <fabs@plt>
  movss  0x10(%rsp),%xmm7
  comiss %xmm0,%xmm7
  jae    L626
L694:
  movss  0x0(%r13),%xmm0
  mov    %ebx,%r12d
  sub    $0x1,%ebx
  lea    -0x2(%r12),%r14d
  call   <fabs@plt>
  movss  0x10(%rsp),%xmm7
  comiss %xmm0,%xmm7
  jb     L684
  mov    (%rsp),%rbp
  mov    0x18(%rsp),%r15d
L660:
  mov    0x80(%rsp),%rax
  movss  (%rax),%xmm2
  cmp    0xb8(%rsp),%r12d
  je     L443
  subl   $0x1,0x90(%rsp)
  je     L320
  mov    0x70(%rsp),%rdi
  movslq %ebx,%rax
  movss  %xmm2,0x40(%rsp)
  lea    0x0(,%rax,4),%r12
  movss  (%rdi,%rax,4),%xmm7
  mov    0x98(%rsp),%rax
  movss  (%rax),%xmm6
  mov    0xa0(%rsp),%rax
  movss  %xmm7,0x18(%rsp)
  movaps %xmm6,%xmm0
  movaps %xmm6,%xmm4
  movss  (%rax),%xmm1
  mov    0x78(%rsp),%rax
  addss  %xmm2,%xmm4
  subss  %xmm2,%xmm0
  movss  %xmm6,0x30(%rsp)
  movss  (%rax),%xmm3
  mulss  %xmm4,%xmm0
  movaps %xmm1,%xmm4
  addss  %xmm3,%xmm1
  movss  %xmm3,0x38(%rsp)
  subss  %xmm3,%xmm4
  mulss  %xmm4,%xmm1
  pxor   %xmm4,%xmm4
  cvtss2sd %xmm6,%xmm4
  pxor   %xmm6,%xmm6
  addss  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd %xmm3,%xmm1
  addsd  %xmm1,%xmm1
  cvtss2sd %xmm0,%xmm0
  mulsd  %xmm4,%xmm1
  divsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm6
  movaps %xmm6,%xmm0
  movss  %xmm6,(%rsp)
  call   <fabs@plt>
  movss  %xmm0,0x28(%rsp)
  movss  D130(%rip),%xmm0
  call   <fabs@plt>
  movss  0x28(%rsp),%xmm4
  movss  0x38(%rsp),%xmm3
  movss  0x40(%rsp),%xmm2
  movaps %xmm0,%xmm1
  comiss %xmm0,%xmm4
  ja     L25
  movss  0xc(%rsp),%xmm7
  ucomiss %xmm7,%xmm0
  jp     L558
  movaps %xmm7,%xmm0
  jne    L558
L59:
  movss  0x18(%rsp),%xmm6
  movss  (%rsp),%xmm7
  movss  %xmm3,0x38(%rsp)
  comiss 0xc(%rsp),%xmm7
  movaps %xmm6,%xmm1
  subss  %xmm2,%xmm1
  addss  %xmm6,%xmm2
  mulss  %xmm2,%xmm1
  movss  %xmm1,0x28(%rsp)
  jb     L307
  call   <fabs@plt>
  movss  0x38(%rsp),%xmm3
L1175:
  addss  (%rsp),%xmm0
  movss  0x30(%rsp),%xmm6
  mov    0x20(%rsp),%rax
  divss  %xmm0,%xmm6
  add    %r12,%rax
  mov    %rax,0x88(%rsp)
  subss  %xmm3,%xmm6
  mulss  %xmm3,%xmm6
  addss  0x28(%rsp),%xmm6
  divss  0x18(%rsp),%xmm6
  cmp    0x60(%rsp),%ebx
  jg     L149
  movss  D130(%rip),%xmm1
  mov    %rax,%r13
  add    0x70(%rsp),%r12
  movaps %xmm1,%xmm3
  nopl   0x0(%rax,%rax,1)
L715:
  mov    %r12,0x38(%rsp)
  mov    %r13,%r14
  movss  0x4(%r12),%xmm7
  add    $0x1,%ebx
  movss  0x4(%r13),%xmm0
  movss  %xmm6,0x48(%rsp)
  add    $0x4,%r13
  add    $0x4,%r12
  movss  %xmm7,0x28(%rsp)
  mulss  %xmm0,%xmm1
  movaps %xmm0,%xmm7
  movaps %xmm6,%xmm0
  mulss  %xmm3,%xmm7
  movss  %xmm1,(%rsp)
  movss  %xmm7,0x30(%rsp)
  call   <fabs@plt>
  movss  %xmm0,0x40(%rsp)
  movss  (%rsp),%xmm0
  call   <fabs@plt>
  movss  0x40(%rsp),%xmm2
  movss  0x48(%rsp),%xmm6
  movaps %xmm0,%xmm1
  comiss %xmm0,%xmm2
  ja     L614
  movss  0xc(%rsp),%xmm7
  ucomiss %xmm7,%xmm0
  jp     L1023
  movaps %xmm7,%xmm0
  jne    L1023
L12:
  movss  (%rsp),%xmm1
  divss  %xmm0,%xmm6
  movss  %xmm0,(%r14)
  xor    %esi,%esi
  movss  0x18(%rsp),%xmm7
  movss  0x30(%rsp),%xmm5
  mov    0x68(%rsp),%rax
  sub    0x20(%rsp),%rax
  add    %r14,%rax
  divss  %xmm0,%xmm1
  movaps %xmm5,%xmm0
  movaps %xmm6,%xmm3
  movaps %xmm7,%xmm6
  mulss  %xmm3,%xmm5
  mulss  %xmm3,%xmm6
  mulss  %xmm1,%xmm0
  mulss  %xmm1,%xmm7
  addss  %xmm0,%xmm6
  movaps %xmm5,%xmm0
  subss  %xmm7,%xmm0
  movss  0x28(%rsp),%xmm7
  movaps %xmm7,%xmm5
  mulss  %xmm1,%xmm5
  movss  %xmm0,0x18(%rsp)
  mulss  %xmm3,%xmm7
  movss  %xmm5,(%rsp)
  movss  %xmm7,0x28(%rsp)
  test   %r15d,%r15d
  jle    L495
  nopw   0x0(%rax,%rax,1)
L351:
  movss  (%rax),%xmm2
  movss  0x4(%rax),%xmm0
  movaps %xmm3,%xmm4
  movaps %xmm1,%xmm5
  add    $0x1,%esi
  mulss  %xmm2,%xmm4
  mulss  %xmm0,%xmm5
  mulss  %xmm1,%xmm2
  mulss  %xmm3,%xmm0
  addss  %xmm5,%xmm4
  subss  %xmm2,%xmm0
  movss  %xmm4,(%rax)
  movss  %xmm0,0x4(%rax)
  add    %rbp,%rax
  cmp    %esi,%r15d
  jne    L351
L495:
  movaps %xmm6,%xmm0
  movss  %xmm1,0x50(%rsp)
  movss  %xmm3,0x48(%rsp)
  movss  %xmm6,0x40(%rsp)
  call   <fabs@plt>
  movss  %xmm0,0x30(%rsp)
  movss  (%rsp),%xmm0
  call   <fabs@plt>
  movss  0x30(%rsp),%xmm4
  movss  0x40(%rsp),%xmm6
  movss  0x48(%rsp),%xmm3
  movss  0x50(%rsp),%xmm1
  movaps %xmm0,%xmm2
  comiss %xmm0,%xmm4
  ja     L95
  movss  0xc(%rsp),%xmm7
  ucomiss %xmm7,%xmm0
  jp     L818
  movaps %xmm7,%xmm0
  jne    L818
L875:
  mov    0x38(%rsp),%rax
  movss  0x28(%rsp),%xmm5
  xor    %esi,%esi
  movss  0x18(%rsp),%xmm7
  movss  %xmm0,(%rax)
  movaps %xmm5,%xmm0
  mov    0x58(%rsp),%rax
  sub    0x20(%rsp),%rax
  mulss  %xmm3,%xmm5
  movaps %xmm7,%xmm6
  add    %r14,%rax
  mulss  %xmm1,%xmm7
  mulss  %xmm3,%xmm6
  mulss  %xmm1,%xmm0
  movaps %xmm5,%xmm2
  subss  %xmm7,%xmm2
  addss  %xmm0,%xmm6
  movss  %xmm2,0x18(%rsp)
  test   %r15d,%r15d
  jle    L191
  nopl   0x0(%rax)
L327:
  movss  (%rax),%xmm2
  movss  0x4(%rax),%xmm0
  movaps %xmm3,%xmm4
  movaps %xmm1,%xmm5
  add    $0x1,%esi
  mulss  %xmm2,%xmm4
  mulss  %xmm0,%xmm5
  mulss  %xmm1,%xmm2
  mulss  %xmm3,%xmm0
  addss  %xmm5,%xmm4
  subss  %xmm2,%xmm0
  movss  %xmm4,(%rax)
  movss  %xmm0,0x4(%rax)
  add    %rbp,%rax
  cmp    %esi,%r15d
  jne    L327
L191:
  cmp    0x60(%rsp),%ebx
  jle    L715
L149:
  mov    0x88(%rsp),%rax
  movl   $0x0,(%rax)
  mov    0x78(%rsp),%rax
  movss  %xmm6,(%rax)
  mov    0x80(%rsp),%rax
  movss  0x18(%rsp),%xmm6
  movss  %xmm6,(%rax)
  jmp    L519
  nopl   0x0(%rax,%rax,1)
L95:
  divss  %xmm4,%xmm2
  pxor   %xmm0,%xmm0
  movss  %xmm1,0x50(%rsp)
  movss  %xmm3,0x48(%rsp)
  movss  %xmm6,0x40(%rsp)
  movss  %xmm4,0x30(%rsp)
  mulss  %xmm2,%xmm2
  cvtss2sd %xmm2,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x30(%rsp),%xmm4
  movss  0x40(%rsp),%xmm6
  movss  0x48(%rsp),%xmm3
  movss  0x50(%rsp),%xmm1
  mulss  %xmm4,%xmm0
L941:
  ucomiss 0xc(%rsp),%xmm0
  jp     L644
  je     L875
L644:
  movss  D130(%rip),%xmm1
  divss  %xmm0,%xmm1
  mulss  %xmm1,%xmm6
  mulss  (%rsp),%xmm1
  movaps %xmm6,%xmm3
  jmp    L875
  nopl   0x0(%rax,%rax,1)
L614:
  divss  %xmm2,%xmm1
  pxor   %xmm0,%xmm0
  mulss  %xmm1,%xmm1
  cvtss2sd %xmm1,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x40(%rsp),%xmm2
  movss  0x48(%rsp),%xmm6
  mulss  %xmm2,%xmm0
  jmp    L12
  nopw   0x0(%rax,%rax,1)
L818:
  divss  %xmm2,%xmm4
  pxor   %xmm0,%xmm0
  movss  %xmm1,0x50(%rsp)
  movss  %xmm3,0x48(%rsp)
  movss  %xmm6,0x40(%rsp)
  movss  %xmm2,0x30(%rsp)
  mulss  %xmm4,%xmm4
  cvtss2sd %xmm4,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x30(%rsp),%xmm2
  movss  0x50(%rsp),%xmm1
  movss  0x48(%rsp),%xmm3
  movss  0x40(%rsp),%xmm6
  mulss  %xmm2,%xmm0
  jmp    L941
  cs nopw 0x0(%rax,%rax,1)
L1023:
  divss  %xmm1,%xmm2
  pxor   %xmm0,%xmm0
  movss  %xmm6,0x48(%rsp)
  movss  %xmm1,0x40(%rsp)
  mulss  %xmm2,%xmm2
  cvtss2sd %xmm2,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x40(%rsp),%xmm1
  movss  0x48(%rsp),%xmm6
  mulss  %xmm1,%xmm0
  jmp    L12
L307:
  call   <fabs@plt>
  movss  0x38(%rsp),%xmm3
  xorps  D83(%rip),%xmm0
  jmp    L1175
L7:
  mov    (%rsp),%rbp
  mov    0x18(%rsp),%r15d
  mov    $0x1,%r12d
  jmp    L660
L626:
  mov    0xb0(%rsp),%rax
  mov    (%rsp),%rbp
  mov    0x18(%rsp),%r15d
  mov    %eax,0x30(%rsp)
  cmp    %eax,%ebx
  jg     L660
  mov    0x20(%rsp),%rax
  movslq %ebx,%rdi
  movslq %r14d,%r14
  mov    %ebx,0x48(%rsp)
  movss  D130(%rip),%xmm4
  mov    %r12d,0x50(%rsp)
  pxor   %xmm3,%xmm3
  mov    %r15d,%r12d
  lea    (%rax,%rdi,4),%r8
  sub    %r14,%rdi
  mov    %rbp,%r15
  mov    %rdi,%r13
  mov    0x58(%rsp),%rdi
  mov    %r13,%rbp
  mov    %r8,%r13
  lea    (%rdi,%r14,4),%rdi
  mov    %rdi,0x38(%rsp)
  mov    0x70(%rsp),%rdi
  sub    %rax,%rdi
  mov    %rdi,0x40(%rsp)
  nopw   0x0(%rax,%rax,1)
L718:
  movss  0x0(%r13),%xmm0
  mulss  %xmm0,%xmm4
  mulss  %xmm3,%xmm0
  movss  %xmm4,0x18(%rsp)
  movss  %xmm0,0x0(%r13)
  movaps %xmm4,%xmm0
  call   <fabs@plt>
  movss  0x10(%rsp),%xmm7
  comiss %xmm0,%xmm7
  jae    L1241
  mov    0x40(%rsp),%rax
  movss  0x18(%rsp),%xmm4
  lea    (%rax,%r13,1),%r14
  movaps %xmm4,%xmm0
  movss  %xmm4,0x28(%rsp)
  movss  (%r14),%xmm7
  movss  %xmm7,(%rsp)
  call   <fabs@plt>
  movss  %xmm0,0x18(%rsp)
  movss  (%rsp),%xmm0
  call   <fabs@plt>
  movss  0x18(%rsp),%xmm2
  movss  0x28(%rsp),%xmm4
  movaps %xmm0,%xmm1
  comiss %xmm0,%xmm2
  ja     L1147
  movss  0xc(%rsp),%xmm7
  ucomiss %xmm7,%xmm0
  jp     L400
  movaps %xmm7,%xmm0
  jne    L400
L1069:
  movss  D130(%rip),%xmm1
  movss  (%rsp),%xmm3
  movss  %xmm0,(%r14)
  xor    %edx,%edx
  xorps  D83(%rip),%xmm4
  mov    0x38(%rsp),%rax
  divss  %xmm0,%xmm1
  mulss  %xmm1,%xmm3
  mulss  %xmm1,%xmm4
  test   %r12d,%r12d
  jle    L945
  nopl   0x0(%rax)
L1170:
  movss  (%rax),%xmm1
  movss  (%rax,%rbp,4),%xmm0
  movaps %xmm3,%xmm2
  movaps %xmm4,%xmm5
  add    $0x1,%edx
  mulss  %xmm1,%xmm2
  mulss  %xmm0,%xmm5
  mulss  %xmm4,%xmm1
  mulss  %xmm3,%xmm0
  addss  %xmm5,%xmm2
  subss  %xmm1,%xmm0
  movss  %xmm2,(%rax)
  movss  %xmm0,(%rax,%rbp,4)
  add    %r15,%rax
  cmp    %edx,%r12d
  jne    L1170
L945:
  add    $0x1,%ebx
  add    $0x4,%r13
  add    $0x1,%rbp
  cmp    0x30(%rsp),%ebx
  jle    L718
L1241:
  mov    %r15,%rbp
  mov    0x48(%rsp),%ebx
  mov    %r12d,%r15d
  mov    0x50(%rsp),%r12d
  jmp    L660
  nopl   0x0(%rax)
L1147:
  divss  %xmm2,%xmm1
  pxor   %xmm0,%xmm0
  mulss  %xmm1,%xmm1
  cvtss2sd %xmm1,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x18(%rsp),%xmm2
  movss  0x28(%rsp),%xmm4
  mulss  %xmm2,%xmm0
  jmp    L1069
L400:
  divss  %xmm1,%xmm2
  pxor   %xmm0,%xmm0
  movss  %xmm4,0x28(%rsp)
  movss  %xmm1,0x18(%rsp)
  mulss  %xmm2,%xmm2
  cvtss2sd %xmm2,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x18(%rsp),%xmm1
  movss  0x28(%rsp),%xmm4
  mulss  %xmm1,%xmm0
  jmp    L1069
L25:
  divss  %xmm4,%xmm1
  pxor   %xmm0,%xmm0
  mulss  %xmm1,%xmm1
  cvtss2sd %xmm1,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x28(%rsp),%xmm4
  movss  0x38(%rsp),%xmm3
  movss  0x40(%rsp),%xmm2
  mulss  %xmm4,%xmm0
  jmp    L59
L558:
  divss  %xmm1,%xmm4
  pxor   %xmm0,%xmm0
  movss  %xmm2,0x40(%rsp)
  movss  %xmm3,0x38(%rsp)
  movss  %xmm1,0x28(%rsp)
  mulss  %xmm4,%xmm4
  cvtss2sd %xmm4,%xmm0
  addsd  D76(%rip),%xmm0
  call   <sqrt@plt>
  movss  0x28(%rsp),%xmm1
  movss  0x40(%rsp),%xmm2
  movss  0x38(%rsp),%xmm3
  mulss  %xmm1,%xmm0
  jmp    L59
L585:
  movss  D130(%rip),%xmm2
  divss  %xmm0,%xmm2
L629:
  cmp    %ebp,%r14d
  jle    L1212
  movslq (%rsp),%rax
  mov    0x58(%rsp),%rdi
  add    %rbp,%rax
  lea    (%rdi,%rax,4),%rax
  nop
L658:
  movss  (%rax),%xmm0
  add    $0x1,%ebx
  mulss  %xmm2,%xmm0
  movss  %xmm0,(%rax)
  add    %r15,%rax
  cmp    %ebx,%r14d
  jne    L658
  jmp    L1212
L443:
  movss  0xc(%rsp),%xmm6
  mov    %r15d,%r14d
  mov    %rbp,%r15
  comiss %xmm2,%xmm6
  ja     L1084
L333:
  subq   $0x1,0xb0(%rsp)
  mov    0xb0(%rsp),%rax
  cmp    %rax,0xa8(%rsp)
  jne    L523
  mov    0x70(%rsp),%rax
  mov    0xd0(%rsp),%rdi
  pxor   %xmm0,%xmm0
  lea    0x4(%rax,%rdi,4),%rsi
L1137:
  movss  (%rax),%xmm1
  add    $0x4,%rax
  maxss  %xmm0,%xmm1
  movaps %xmm1,%xmm0
  cmp    %rsi,%rax
  jne    L1137
  mulss  0xbc(%rsp),%xmm0
  mov    0x70(%rsp),%rax
  xor    %r12d,%r12d
L61:
  movss  (%rax),%xmm2
  pxor   %xmm1,%xmm1
  comiss %xmm2,%xmm0
  jae    L939
  movss  D130(%rip),%xmm1
  add    $0x1,%r12d
  divss  %xmm2,%xmm1
L939:
  movss  %xmm1,(%rax)
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L61
  test   %r14d,%r14d
  jle    L250
  mov    0x68(%rsp),%rdi
  mov    0x70(%rsp),%r9
  xor    %r11d,%r11d
  mov    0xd8(%rsp),%rbx
  sub    %r9,%rdi
L1090:
  mov    0x58(%rsp),%rdx
  mov    %rbx,%rcx
  xor    %r8d,%r8d
  sub    %r9,%rdx
L1127:
  mov    %r9,%rax
  pxor   %xmm1,%xmm1
  nopl   0x0(%rax,%rax,1)
L266:
  movss  (%rax,%rdi,1),%xmm0
  mulss  (%rdx,%rax,1),%xmm0
  add    $0x4,%rax
  mulss  -0x4(%rax),%xmm0
  addss  %xmm0,%xmm1
  cmp    %rsi,%rax
  jne    L266
  lea    0x1(%r8),%eax
  movss  %xmm1,(%rcx)
  add    %r15,%rdx
  add    %r15,%rcx
  cmp    %eax,%r14d
  je     L2
  mov    %eax,%r8d
  jmp    L1127
L320:
  xor    %r12d,%r12d
L859:
  add    $0xe8,%rsp
  mov    %r12d,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L528:
  call   <sqrt@plt>
  call   <fabs@plt>
  movss  0x60(%rsp),%xmm1
  mov    0x30(%rsp),%rsi
  mov    0xb0(%rsp),%rcx
  movaps %xmm0,%xmm4
  movss  0xa8(%rsp),%xmm3
  jmp    L667
L540:
  cmp    %edi,%r14d
  jle    L243
  mov    0x58(%rsp),%rdx
  lea    (%rdi,%r13,1),%eax
  movslq %edi,%rbx
  mov    0x80(%rsp),%r11
  cltq
  mov    %rbx,0x30(%rsp)
  lea    (%rdx,%rax,4),%r8
  movslq %r13d,%rax
  mov    %rax,0x38(%rsp)
  add    %rbx,%rax
  lea    (%rdx,%rax,4),%rbx
  mov    %rax,%r9
  mov    %r14d,%eax
  mov    0x68(%rsp),%rdx
  imul   %edi,%eax
  mov    %rbx,0x28(%rsp)
  movslq %eax,%r10
  lea    (%rcx,%r10,1),%rax
  lea    (%rdx,%rax,4),%rdx
  mov    0xa0(%rsp),%eax
  sub    %edi,%eax
  add    %r9,%rax
  lea    (%r11,%rax,4),%r9
  mov    %rbx,%rax
L854:
  movss  (%rax),%xmm0
  divss  (%r8),%xmm0
  divss  %xmm3,%xmm0
  add    $0x4,%rax
  movss  %xmm0,(%rdx)
  add    %rsi,%rdx
  cmp    %rax,%r9
  jne    L854
  mov    0x30(%rsp),%rax
  mov    0x68(%rsp),%rdx
  mov    %rcx,%r8
  mov    %edi,%ebx
  add    %rax,%r10
  sub    %rax,%r8
  lea    (%rdx,%r10,4),%r11
  shl    $0x2,%r8
L92:
  mov    0x28(%rsp),%rax
  mov    %r11,%rdx
  pxor   %xmm1,%xmm1
  nopl   0x0(%rax)
L1196:
  movss  (%rax),%xmm0
  mulss  (%rdx),%xmm0
  add    $0x4,%rax
  add    %rsi,%rdx
  addss  %xmm0,%xmm1
  cmp    %rax,%r9
  jne    L1196
  mov    %r11,%rax
  mov    %edi,%edx
  nopl   (%rax)
L326:
  movss  (%r8,%rax,1),%xmm0
  mov    %edx,%r10d
  add    $0x1,%edx
  mulss  %xmm1,%xmm0
  addss  (%rax),%xmm0
  movss  %xmm0,(%rax)
  add    %rsi,%rax
  cmp    %edx,%r14d
  jne    L326
  lea    0x1(%rbx),%eax
  add    $0x4,%r11
  sub    $0x4,%r8
  cmp    %ebx,%r10d
  je     L686
  mov    %eax,%ebx
  jmp    L92
L567:
  call   <sqrt@plt>
  call   <fabs@plt>
  movss  0xbc(%rsp),%xmm2
  movss  0xb0(%rsp),%xmm1
  movaps %xmm0,%xmm3
  jmp    L1239
L1084:
  movss  D83(%rip),%xmm1
  xorps  %xmm1,%xmm2
  movss  %xmm2,(%rax)
  test   %r14d,%r14d
  jle    L333
  mov    0xc8(%rsp),%rax
  xor    %esi,%esi
  add    0x68(%rsp),%rax
L1215:
  movss  (%rax),%xmm0
  add    $0x1,%esi
  xorps  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  add    %r15,%rax
  cmp    %esi,%r14d
  jne    L1215
  jmp    L333
L887:
  mov    0x28(%rsp),%eax
  add    %r13d,%eax
  mov    %eax,0x30(%rsp)
  jmp    L814
L751:
  jne    L34
  xor    %r12d,%r12d
L250:
  mov    0xc0(%rsp),%rdi
  call   <free@plt>
  mov    0x20(%rsp),%rdi
  call   <free@plt>
  mov    0x68(%rsp),%rdi
  call   <free@plt>
  mov    0x70(%rsp),%rdi
  call   <free@plt>
  jmp    L859
L2:
  lea    0x1(%r11),%eax
  add    %r15,%rdi
  add    $0x4,%rbx
  cmp    %r11d,%r8d
  je     L250
  mov    %eax,%r11d
  jmp    L1090
L34:
  movss  D181(%rip),%xmm7
  pxor   %xmm6,%xmm6
  xor    %ebp,%ebp
  movss  %xmm6,0xc(%rsp)
  movaps %xmm6,%xmm3
  movss  %xmm7,0xbc(%rsp)
  movss  %xmm7,0x10(%rsp)
  jmp    L1088
```

tool_calls: name: disassemble_data
arguments: label: D181
type_infered: byte[4]
name: disassemble_data
arguments: label: D83
type_infered: byte[16]
name: disassemble_data
arguments: label: D130
type_infered: f32
name: disassemble_data
arguments: label: D76
type_infered: f32[2]

tool: {"D181": [0xbd,0x37,0x86,0x35]}

tool: {"D83": [0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: D130: 1

tool: {"D76": [0.0,1.875]}

assistant: ```c
int ssvdinv(float *a, float *b, int m) {
  int flag, i, its, j, jj, k, l, nm, rank;
  float *vmat, *wmat, *w, *rv1, *tmp, c, f, h, s, x, y, z, anorm, g, scale, at, bt, ct, maxarg1, maxarg2, thresh, wmax, tol, tanorm;
  anorm = g = scale = 0.0;
  tanorm = tol = sizeof(anorm) > 4 ? 1.0e-9 : 1.0e-6;
  rv1 = (float *)malloc(m * sizeof(float));
  tmp = (float *)malloc(m * sizeof(float));
  vmat = (float *)malloc(m * m * sizeof(float));
  wmat = (float *)malloc(m * sizeof(float));
  l = nm = 0;
  for (i = 0; i < m; i++) {
    l = i + 1;
    rv1[i] = scale * g;
    g = s = scale = 0.0;
    for (k = i; k < m; k++)
      scale += fabs(a[k * m + i]);
    if (scale) {
      for (k = i; k < m; k++) {
        a[k * m + i] /= scale;
        s += a[k * m + i] * a[k * m + i];
      }
      f = a[i * m + i];
      g = -((f) >= 0.0 ? fabs(sqrt(s)) : -fabs(sqrt(s)));
      h = f * g - s;
      a[i * m + i] = f - g;
      for (j = l; j < m; j++) {
        s = 0.0;
        for (k = i; k < m; k++)
          s += a[k * m + i] * a[k * m + j];
        f = s / h;
        for (k = i; k < m; k++)
          a[k * m + j] += f * a[k * m + i];
      }
      for (k = i; k < m; k++)
        a[k * m + i] *= scale;
    }
    wmat[i] = scale * g;
    g = s = scale = 0.0;
    if (l != m) {
      for (k = l; k < m; k++)
        scale += fabs(a[i * m + k]);
      if (scale) {
        for (k = l; k < m; k++) {
          a[i * m + k] /= scale;
          s += a[i * m + k] * a[i * m + k];
        }
        f = a[i * m + l];
        g = -((f) >= 0.0 ? fabs(sqrt(s)) : -fabs(sqrt(s)));
        h = f * g - s;
        a[i * m + l] = f - g;
        for (k = l; k < m; k++)
          rv1[k] = a[i * m + k] / h;
        for (j = l; j < m; j++) {
          s = 0.0;
          for (k = l; k < m; k++)
            s += a[j * m + k] * a[i * m + k];
          for (k = l; k < m; k++)
            a[j * m + k] += s * rv1[k];
        }
        for (k = l; k < m; k++)
          a[i * m + k] *= scale;
      }
    }
    anorm = (maxarg1 = (anorm), maxarg2 = ((fabs(wmat[i]) + fabs(rv1[i]))), (maxarg1) > (maxarg2) ? (maxarg1) : (maxarg2));
    tanorm = tol * anorm;
  }
  for (i = m; i--;) {
    if (i < m - 1) {
      if (g) {
        for (j = l; j < m; j++)
          vmat[j * m + i] = (a[i * m + j] / a[i * m + l]) / g;
        for (j = l; j < m; j++) {
          s = 0.0;
          for (k = l; k < m; k++)
            s += a[i * m + k] * vmat[k * m + j];
          for (k = l; k < m; k++)
            vmat[k * m + j] += s * vmat[k * m + i];
        }
      }
      for (j = l; j < m; j++)
        vmat[i * m + j] = vmat[j * m + i] = 0.0;
    }
    vmat[i * m + i] = 1.0;
    g = rv1[i];
    l = i;
  }
  for (i = m; i--;) {
    l = i + 1;
    g = wmat[i];
    for (j = l; j < m; j++)
      a[i * m + j] = 0.0;
    if (g) {
      g = 1.0 / g;
      for (j = l; j < m; j++) {
        s = 0.0;
        for (k = l; k < m; k++)
          s += a[k * m + i] * a[k * m + j];
        f = (s / a[i * m + i]) * g;
        for (k = i; k < m; k++)
          a[k * m + j] += f * a[k * m + i];
      }
      for (j = i; j < m; j++)
        a[j * m + i] *= g;
    } else
      for (j = i; j < m; j++)
        a[j * m + i] = 0.0;
    ++a[i * m + i];
  }
  for (k = m; k--;) {
    for (its = 0; its < 30; its++) {
      flag = 1;
      for (l = k + 1; l--;) {
        nm = l - 1;
        if (fabs(rv1[l]) <= tanorm || !l) {
          flag = 0;
          break;
        }
        if (fabs(wmat[nm]) <= tanorm)
          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 (fabs(f) <= tanorm)
            break;
          g = wmat[i];
          h = ((at = fabs(f)) > (bt = fabs(g)) ? (ct = bt / at, at * sqrt(1.0 + ct * ct)) : (bt ? (ct = at / bt, bt * sqrt(1.0 + ct * ct)) : 0.0));
          wmat[i] = h;
          h = 1.0 / h;
          c = g * h;
          s = -f * h;
          for (j = 0; j < m; j++) {
            y = a[j * m + nm];
            z = a[j * m + i];
            a[j * m + nm] = y * c + z * s;
            a[j * m + i] = z * c - y * s;
          }
        }
      }
      z = wmat[k];
      if (l == k) {
        if (z < 0.0) {
          wmat[k] = -z;
          for (j = 0; j < m; j++)
            vmat[j * m + k] = -vmat[j * m + k];
        }
        break;
      }
      if (its == 29)
        return 0;
      x = wmat[l];
      nm = k - 1;
      y = wmat[nm];
      g = rv1[nm];
      h = rv1[k];
      f = ((y - z) * (y + z) + (g - h) * (g + h)) / (2.0 * h * y);
      g = ((at = fabs(f)) > (bt = fabs(1.0)) ? (ct = bt / at, at * sqrt(1.0 + ct * ct)) : (bt ? (ct = at / bt, bt * sqrt(1.0 + ct * ct)) : 0.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 = wmat[i];
        h = s * g;
        g = c * g;
        z = ((at = fabs(f)) > (bt = fabs(h)) ? (ct = bt / at, at * sqrt(1.0 + ct * ct)) : (bt ? (ct = at / bt, bt * sqrt(1.0 + ct * ct)) : 0.0));
        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 < m; jj++) {
          x = vmat[jj * m + j];
          z = vmat[jj * m + i];
          vmat[jj * m + j] = x * c + z * s;
          vmat[jj * m + i] = z * c - x * s;
        }
        z = ((at = fabs(f)) > (bt = fabs(h)) ? (ct = bt / at, at * sqrt(1.0 + ct * ct)) : (bt ? (ct = at / bt, bt * sqrt(1.0 + ct * ct)) : 0.0));
        wmat[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 * m + j];
          z = a[jj * m + i];
          a[jj * m + j] = y * c + z * s;
          a[jj * m + i] = z * c - y * s;
        }
      }
      rv1[l] = 0.0;
      rv1[k] = f;
      wmat[k] = x;
    }
  }
  wmax = 0.0;
  w = wmat;
  for (j = m; j--; w++)
    if (*w > wmax)
      wmax = *w;
  thresh = tol * wmax;
  rank = 0;
  w = wmat;
  for (j = m; j--; w++)
    if (*w <= thresh)
      *w = 0.0;
    else {
      rank++;
      *w = 1.0 / *w;
    }
  for (j = 0; j < m; j++)
    for (i = 0; i < m; i++) {
      s = 0.0;
      for (k = 0; k < m; k++)
        s += vmat[k + j * m] * a[k + i * m] * wmat[k];
      b[j + i * m] = s;
    }
  free(tmp);
  free(rv1);
  free(vmat);
  free(wmat);
  return rank;
}
```