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

```asm
<svd_pseudoinverse>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x30,%rsp
  mov    %rdi,-0x18(%rbp)
  mov    %rsi,-0x20(%rbp)
  mov    %edx,-0x24(%rbp)
  mov    %ecx,-0x28(%rbp)
  mov    -0x28(%rbp),%ecx
  mov    -0x24(%rbp),%edx
  mov    -0x20(%rbp),%rsi
  mov    -0x18(%rbp),%rax
  mov    %rax,%rdi
  call   <svd@plt>
  mov    %eax,-0x4(%rbp)
  movl   $0x0,-0xc(%rbp)
  jmp    L95
L31:
  movl   $0x0,-0x10(%rbp)
  jmp    L60
L92:
  mov    -0x24(%rbp),%edx
  mov    -0xc(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x10(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x24(%rbp),%edx
  mov    -0xc(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x10(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x10(%rbp)
L60:
  mov    -0x10(%rbp),%eax
  cmp    -0x4(%rbp),%eax
  jl     L92
  addl   $0x1,-0xc(%rbp)
L95:
  mov    -0xc(%rbp),%eax
  cmp    -0x28(%rbp),%eax
  jl     L31
  movl   $0x0,-0x10(%rbp)
  jmp    L25
L55:
  mov    -0x28(%rbp),%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edx
  mov    -0x20(%rbp),%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  movl   $0x0,-0xc(%rbp)
  jmp    L69
L18:
  movl   $0x0,-0x8(%rbp)
  jmp    L6
L74:
  mov    -0xc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x24(%rbp),%edx
  mov    -0xc(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0xc(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x8(%rbp)
L6:
  mov    -0x8(%rbp),%eax
  cmp    -0x4(%rbp),%eax
  jl     L74
  addl   $0x1,-0xc(%rbp)
L69:
  mov    -0xc(%rbp),%eax
  cmp    -0x28(%rbp),%eax
  jl     L18
  mov    -0x28(%rbp),%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edx
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x18(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rax
  mov    -0x20(%rbp),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
  addl   $0x1,-0x10(%rbp)
L25:
  mov    -0x10(%rbp),%eax
  cmp    -0x24(%rbp),%eax
  jl     L55
  mov    -0x4(%rbp),%eax
  leave
  ret
```

assistant: ```c
int svd_pseudoinverse(double **_w, double *_s, int _nrows, int _ncols) {
  int rank;
  int i;
  int j;
  int k;
  rank = svd(_w, _s, _nrows, _ncols);
  for (j = 0; j < _ncols; j++)
    for (i = 0; i < rank; i++)
      _w[_nrows + j][i] /= _s[i];
  for (i = 0; i < _nrows; i++) {
    memset(_s, 0, _ncols * sizeof(*_s));
    for (j = 0; j < _ncols; j++)
      for (k = 0; k < rank; k++)
        _s[j] += _w[i][k] * _w[_nrows + j][k];
    memcpy(_w[i], _s, _ncols * sizeof(*_s));
  }
  return rank;
}
```