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

```asm
<svd_lsq>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %esi,-0x4c(%rbp)
  mov    %edx,-0x50(%rbp)
  mov    %rcx,-0x58(%rbp)
  mov    %r8,-0x60(%rbp)
  mov    %r9,-0x68(%rbp)
  mov    -0x4c(%rbp),%edx
  mov    -0x50(%rbp),%eax
  lea    (%rdx,%rax,1),%ecx
  mov    -0x4c(%rbp),%eax
  mov    $0x8,%edx
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <alloc2d@plt>
  mov    %rax,-0x28(%rbp)
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x20(%rbp)
  mov    -0x50(%rbp),%edx
  mov    -0x4c(%rbp),%eax
  cmp    %eax,%edx
  cmovle %edx,%eax
  mov    %eax,-0x2c(%rbp)
  cmpq   $0x0,-0x60(%rbp)
  je     L40
  movl   $0x0,-0x34(%rbp)
  jmp    L174
L110:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x0,-0x38(%rbp)
  jmp    L218
L84:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x38(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x38(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x38(%rbp)
L218:
  mov    -0x38(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jl     L84
  addl   $0x1,-0x34(%rbp)
L174:
  mov    -0x34(%rbp),%eax
  cmp    -0x50(%rbp),%eax
  jl     L110
L40:
  mov    -0x28(%rbp),%rdi
  mov    -0x68(%rbp),%rcx
  mov    -0x50(%rbp),%edx
  mov    -0x4c(%rbp),%esi
  mov    -0x48(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <svd@plt>
  movl   $0x0,-0x34(%rbp)
  jmp    L97
L68:
  movl   $0x0,-0x38(%rbp)
  jmp    L117
L154:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  jp     L159
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  je     L247
L159:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x38(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x38(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L206
L247:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x38(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
L206:
  addl   $0x1,-0x38(%rbp)
L117:
  mov    -0x38(%rbp),%eax
  cmp    -0x2c(%rbp),%eax
  jl     L154
  addl   $0x1,-0x34(%rbp)
L97:
  mov    -0x34(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jl     L68
  movl   $0x0,-0x38(%rbp)
  jmp    L82
L111:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    %rax,-0x18(%rbp)
  movl   $0x0,-0x34(%rbp)
  jmp    L86
L126:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x38(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  movl   $0x0,-0x30(%rbp)
  jmp    L60
L127:
  mov    -0x8(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x8(%rbp),%rax
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x30(%rbp)
L60:
  mov    -0x30(%rbp),%eax
  cmp    -0x2c(%rbp),%eax
  jl     L127
  addl   $0x1,-0x34(%rbp)
L86:
  mov    -0x34(%rbp),%eax
  cmp    -0x50(%rbp),%eax
  jl     L126
  addl   $0x1,-0x38(%rbp)
L82:
  mov    -0x38(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jl     L111
  movl   $0x0,-0x38(%rbp)
  jmp    L52
L92:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x0,-0x34(%rbp)
  jmp    L171
L61:
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x38(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x34(%rbp)
L171:
  mov    -0x34(%rbp),%eax
  cmp    -0x50(%rbp),%eax
  jl     L61
  addl   $0x1,-0x38(%rbp)
L52:
  mov    -0x38(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jl     L92
  cmpq   $0x0,-0x60(%rbp)
  je     L207
  movl   $0x0,-0x34(%rbp)
  jmp    L31
L227:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x34(%rbp)
L31:
  mov    -0x34(%rbp),%eax
  cmp    -0x50(%rbp),%eax
  jl     L227
L207:
  mov    -0x28(%rbp),%rax
  mov    %rax,%rdi
  call   <free2d@plt>
  nop
  leave
  ret
```

assistant: ```c
void svd_lsq(double **A, int ni, int nj, double *z, double *std, double *w, double *sol) {
  double **V = alloc2d(ni, ni + nj, sizeof(double));
  double **B = &V[ni];
  int nijmin = (ni < nj) ? ni : nj;
  int i, j, ii;
  if (std != NULL)
    for (j = 0; j < nj; ++j) {
      z[j] /= std[j];
      for (i = 0; i < ni; ++i)
        A[j][i] /= std[j];
    }
  svd(A, ni, nj, w, V);
  for (j = 0; j < ni; ++j)
    for (i = 0; i < nijmin; ++i)
      if (w[i] != 0.0)
        V[j][i] /= w[i];
      else
        V[j][i] = 0.0;
  for (i = 0; i < ni; ++i) {
    double *v = V[i];
    for (j = 0; j < nj; ++j) {
      double *a = A[j];
      double *b = &B[j][i];
      for (ii = 0; ii < nijmin; ++ii)
        *b += v[ii] * a[ii];
    }
  }
  for (i = 0; i < ni; ++i) {
    sol[i] = 0.0;
    for (j = 0; j < nj; ++j)
      sol[i] += B[j][i] * z[j];
  }
  if (std != NULL)
    for (j = 0; j < nj; ++j)
      z[j] *= std[j];
  free2d(V);
}
```