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

```asm
<add_point>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %ecx,-0x4c(%rbp)
  mov    %r8,-0x58(%rbp)
  mov    %r9d,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x28(%rbp)
  cmpl   $0x0,-0x4c(%rbp)
  je     L10
  movl   $0x1,-0x24(%rbp)
  jmp    L24
L15:
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  ucomiss %xmm1,%xmm0
  jp     L28
  pxor   %xmm1,%xmm1
  ucomiss %xmm1,%xmm0
  je     L18
L28:
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x18(%rbp)
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x14(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2ssl -0x4c(%rbp),%xmm0
  mulss  -0x18(%rbp),%xmm0
  mulss  -0x18(%rbp),%xmm0
  movss  -0x14(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x10(%rbp)
  movss  -0x14(%rbp),%xmm0
  divss  -0x10(%rbp),%xmm0
  movss  %xmm0,-0xc(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2ssl -0x4c(%rbp),%xmm0
  mulss  -0x18(%rbp),%xmm0
  divss  -0x10(%rbp),%xmm0
  movss  %xmm0,-0x8(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2ssl -0x4c(%rbp),%xmm0
  mulss  -0xc(%rbp),%xmm0
  cvttss2si %xmm0,%eax
  mov    %eax,-0x4c(%rbp)
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movss  -0x10(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x50(%rbp),%eax
  add    %eax,%eax
  sub    -0x24(%rbp),%eax
  imul   %edx,%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  add    $0x1,%eax
  mov    %eax,-0x1c(%rbp)
  mov    -0x24(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x20(%rbp)
  jmp    L40
L21:
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x4(%rbp)
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x18(%rbp),%xmm1
  mov    -0x20(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  -0x4(%rbp),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0xc(%rbp),%xmm1
  movss  -0x8(%rbp),%xmm0
  mulss  -0x4(%rbp),%xmm0
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x1c(%rbp)
  addl   $0x1,-0x20(%rbp)
L40:
  mov    -0x20(%rbp),%eax
  cmp    -0x50(%rbp),%eax
  jle    L21
  mov    -0x50(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x4(%rbp)
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x18(%rbp),%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  -0x4(%rbp),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0xc(%rbp),%xmm1
  movss  -0x8(%rbp),%xmm0
  mulss  -0x4(%rbp),%xmm0
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
L18:
  addl   $0x1,-0x24(%rbp)
L24:
  mov    -0x24(%rbp),%eax
  cmp    -0x50(%rbp),%eax
  jle    L15
  pxor   %xmm1,%xmm1
  cvtsi2ssl -0x4c(%rbp),%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x28(%rbp)
L10:
  movss  -0x28(%rbp),%xmm0
  pop    %rbp
  ret
```

assistant: ```c
float add_point(r_bar, diag, theta_bar, weight, drow, ivars)
float *r_bar, *diag, *theta_bar, *drow;
int ivars;
{
  float ss_error = 0.0, c_bar, s_bar, xi, xk, di, d_prime_i;
  int i, k, next_r;
  if (weight != 0.0) {
    for (i = 1; i <= ivars; i++) {
      if (drow[i] != 0) {
        xi = drow[i];
        di = diag[i];
        d_prime_i = di + weight * xi * xi;
        c_bar = di / d_prime_i;
        s_bar = weight * xi / d_prime_i;
        weight = c_bar * weight;
        diag[i] = d_prime_i;
        next_r = (i - 1) * (2 * ivars - i) / 2 + 1;
        for (k = i + 1; k <= ivars; k++) {
          xk = drow[k];
          drow[k] = xk - xi * r_bar[next_r];
          r_bar[next_r] = c_bar * r_bar[next_r] + s_bar * xk;
          next_r++;
        }
        xk = drow[ivars + 1];
        drow[ivars + 1] = xk - xi * theta_bar[i];
        theta_bar[i] = c_bar * theta_bar[i] + s_bar * xk;
      }
    }
    ss_error = weight * drow[ivars + 1] * drow[ivars + 1];
  }
  return (ss_error);
}
```