user: What is the c source code of the assembly code below?

```asm
<svdsolve>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %rdx,-0x38(%rbp)
  mov    %ecx,-0x3c(%rbp)
  mov    %r8d,-0x40(%rbp)
  mov    %r9,-0x48(%rbp)
  movq   $0x0,-0x8(%rbp)
  mov    -0x40(%rbp),%eax
  mov    $0x8,%esi
  mov    %eax,%edi
  call   <calloc@plt>
  mov    %rax,-0x8(%rbp)
  cmpq   $0x0,-0x8(%rbp)
  jne    L29
  lea    D7(%rip),%rax
  mov    %rax,%rdi
  call   <perror@plt>
  mov    $0xffffffff,%eax
  jmp    L93
L29:
  movl   $0x0,-0x18(%rbp)
  jmp    L77
L75:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  jp     L56
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  je     L39
L56:
  movl   $0x0,-0x14(%rbp)
  jmp    L53
L78:
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x10(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  addl   $0x1,-0x14(%rbp)
L53:
  mov    -0x14(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jl     L78
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x10(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
L39:
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x18(%rbp)
L77:
  mov    -0x18(%rbp),%eax
  cmp    -0x40(%rbp),%eax
  jl     L75
  movl   $0x0,-0x18(%rbp)
  jmp    L102
L90:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movl   $0x0,-0x1c(%rbp)
  jmp    L62
L87:
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x1c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x10(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  addl   $0x1,-0x1c(%rbp)
L62:
  mov    -0x1c(%rbp),%eax
  cmp    -0x40(%rbp),%eax
  jl     L87
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x18(%rbp)
L102:
  mov    -0x18(%rbp),%eax
  cmp    -0x40(%rbp),%eax
  jl     L90
  mov    -0x8(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    $0x0,%eax
L93:
  leave
  ret
```

tool_calls: name: parse_data
arguments: label: D7
type: string

tool: D7: svdsolve

assistant: ```c
int svdsolve(double **u, double *w, double **v, int m, int n, double *b, double *x) {
  int jj, j, i;
  double s, *tmp = NULL;
  tmp = (double *)calloc(n, sizeof(double));
  if (NULL == tmp) {
    perror("svdsolve");
    return -1;
  }
  for (j = 0; j < n; j++) {
    s = 0.0;
    if (w[j]) {
      for (i = 0; i < m; i++) {
        s += u[i][j] * b[i];
      }
      s /= w[j];
    }
    tmp[j] = s;
  }
  for (j = 0; j < n; j++) {
    s = 0.0;
    for (jj = 0; jj < n; jj++) {
      s += v[j][jj] * tmp[jj];
    }
    x[j] = s;
  }
  free(tmp);
  return 0;
}
```