user: 请从给定的汇编指令推导出C源代码：

```asm
<add_point>:
  endbr64
  push   %r15
  pxor   %xmm0,%xmm0
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rsi,-0x30(%rsp)
  mov    %rdx,-0x28(%rsp)
  test   %ecx,%ecx
  je     L131
  lea    0x1(%r9),%eax
  mov    %ecx,%r13d
  cltq
  lea    (%r8,%rax,4),%rax
  mov    %rax,-0x38(%rsp)
  test   %r9d,%r9d
  jle    L64
  mov    %rdi,%rbx
  lea    0x1(%r9,%r9,1),%eax
  lea    0x8(%r8),%rdi
  mov    $0x1,%esi
  mov    %eax,-0x14(%rsp)
  lea    0x8(%rbx),%rax
  lea    -0x1(%r9),%r10d
  mov    $0x1,%edx
  mov    %rax,-0x10(%rsp)
  jmp    L162
  nopl   0x0(%rax,%rax,1)
L177:
  mov    %r11,%rsi
  add    $0x4,%rdi
  sub    $0x1,%r10d
  cmp    %ebp,%r9d
  je     L64
L162:
  movss  -0x4(%rdi),%xmm3
  mov    %edx,%ebp
  add    $0x1,%edx
  ucomiss %xmm0,%xmm3
  jp     L147
  lea    0x1(%rsi),%r11
  je     L177
L147:
  pxor   %xmm2,%xmm2
  movaps %xmm3,%xmm5
  movaps %xmm3,%xmm1
  mov    -0x30(%rsp),%rax
  cvtsi2ss %r13d,%xmm2
  lea    -0x1(%rbp),%ecx
  movss  (%rax,%rsi,4),%xmm4
  mulss  %xmm2,%xmm5
  mulss  %xmm5,%xmm1
  addss  %xmm4,%xmm1
  divss  %xmm1,%xmm4
  movss  %xmm1,(%rax,%rsi,4)
  mov    -0x14(%rsp),%eax
  sub    %edx,%eax
  imul   %ecx,%eax
  mov    %eax,%r11d
  shr    $0x1f,%r11d
  add    %eax,%r11d
  sar    %r11d
  lea    0x1(%r11),%ecx
  divss  %xmm1,%xmm5
  mulss  %xmm4,%xmm2
  cvttss2si %xmm2,%r13d
  cmp    %edx,%r9d
  jl     L81
  movslq %r11d,%r11
  lea    0x10(%rdi),%r14
  lea    0x14(,%r11,4),%r12
  lea    0x1(%r11),%rax
  mov    %rax,-0x20(%rsp)
  lea    -0x10(%rbx,%r12,1),%rax
  add    %rbx,%r12
  cmp    %r12,%rdi
  setae  %r15b
  cmp    %r14,%rax
  setae  %r12b
  or     %r15b,%r12b
  je     L92
  cmp    $0x1,%r10d
  je     L92
  mov    %r9d,%r11d
  sub    %edx,%r11d
  cmp    $0x2,%r11d
  jbe    L2
  movaps %xmm3,%xmm9
  movaps %xmm4,%xmm8
  movaps %xmm5,%xmm7
  mov    %r10d,%r11d
  shr    $0x2,%r11d
  shufps $0x0,%xmm9,%xmm9
  shufps $0x0,%xmm8,%xmm8
  shufps $0x0,%xmm7,%xmm7
  lea    -0x1(%r11),%r12d
  mov    %rdi,%r11
  shl    $0x4,%r12
  add    %r14,%r12
  nopl   0x0(%rax)
L0:
  movups (%rax),%xmm1
  movups (%r11),%xmm2
  add    $0x10,%r11
  add    $0x10,%rax
  mulps  %xmm9,%xmm1
  movaps %xmm2,%xmm6
  mulps  %xmm7,%xmm2
  subps  %xmm1,%xmm6
  movups %xmm6,-0x10(%r11)
  movups -0x10(%rax),%xmm1
  mulps  %xmm8,%xmm1
  addps  %xmm2,%xmm1
  movups %xmm1,-0x10(%rax)
  cmp    %r12,%r11
  jne    L0
  mov    %r10d,%r12d
  lea    0x1(%rsi),%r11
  and    $0xfffffffc,%r12d
  lea    (%rdx,%r12,1),%eax
  add    %r12d,%ecx
  cmp    %r10d,%r12d
  je     L48
  mov    %r10d,%r14d
  sub    %r12d,%r14d
  cmp    $0x1,%r14d
  je     L25
L171:
  lea    0x1(%rsi),%r11
  movaps %xmm3,%xmm2
  lea    (%r11,%r12,1),%r15
  add    -0x20(%rsp),%r12
  shufps $0xe0,%xmm2,%xmm2
  lea    (%rbx,%r12,4),%r12
  lea    (%r8,%r15,4),%r15
  movq   (%r12),%xmm1
  movq   (%r15),%xmm6
  mulps  %xmm2,%xmm1
  movaps %xmm6,%xmm2
  subps  %xmm1,%xmm2
  movlps %xmm2,(%r15)
  movq   (%r12),%xmm1
  movaps %xmm4,%xmm2
  shufps $0xe0,%xmm2,%xmm2
  mulps  %xmm2,%xmm1
  movaps %xmm5,%xmm2
  shufps $0xe0,%xmm2,%xmm2
  mulps  %xmm6,%xmm2
  addps  %xmm2,%xmm1
  movlps %xmm1,(%r12)
  mov    %r14d,%r12d
  and    $0xfffffffe,%r12d
  add    %r12d,%eax
  add    %r12d,%ecx
  cmp    %r12d,%r14d
  je     L48
L25:
  movslq %ecx,%rcx
  cltq
  lea    (%rbx,%rcx,4),%rcx
  lea    (%r8,%rax,4),%rax
  movss  (%rcx),%xmm6
  movss  (%rax),%xmm1
  mulss  %xmm3,%xmm6
  movaps %xmm1,%xmm2
  mulss  %xmm5,%xmm1
  subss  %xmm6,%xmm2
  movss  %xmm2,(%rax)
  movss  (%rcx),%xmm2
  mulss  %xmm4,%xmm2
  addss  %xmm2,%xmm1
  movss  %xmm1,(%rcx)
L48:
  mov    -0x28(%rsp),%rax
  mov    -0x38(%rsp),%rcx
  add    $0x4,%rdi
  sub    $0x1,%r10d
  mulss  (%rax,%rsi,4),%xmm3
  movss  (%rcx),%xmm1
  movaps %xmm1,%xmm2
  mulss  %xmm5,%xmm1
  subss  %xmm3,%xmm2
  movss  %xmm2,(%rcx)
  mulss  (%rax,%rsi,4),%xmm4
  addss  %xmm1,%xmm4
  movss  %xmm4,(%rax,%rsi,4)
  mov    %r11,%rsi
  cmp    %ebp,%r9d
  jne    L162
L64:
  pxor   %xmm0,%xmm0
  mov    -0x38(%rsp),%rax
  cvtsi2ss %r13d,%xmm0
  movss  (%rax),%xmm1
  mulss  %xmm1,%xmm0
  mulss  %xmm1,%xmm0
L131:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L92:
  mov    %r9d,%ecx
  mov    -0x10(%rsp),%r15
  sub    %edx,%ecx
  add    %r11,%rcx
  lea    (%r15,%rcx,4),%r11
  mov    %rdi,%rcx
  nopl   0x0(%rax)
L168:
  movss  (%rax),%xmm6
  movss  (%rcx),%xmm1
  add    $0x4,%rax
  add    $0x4,%rcx
  mulss  %xmm3,%xmm6
  movaps %xmm1,%xmm2
  mulss  %xmm5,%xmm1
  subss  %xmm6,%xmm2
  movss  %xmm2,-0x4(%rcx)
  movss  -0x4(%rax),%xmm2
  mulss  %xmm4,%xmm2
  addss  %xmm2,%xmm1
  movss  %xmm1,-0x4(%rax)
  cmp    %r11,%rax
  jne    L168
L81:
  lea    0x1(%rsi),%r11
  jmp    L48
L2:
  mov    %r10d,%r14d
  mov    %edx,%eax
  xor    %r12d,%r12d
  jmp    L171
```

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);
}
```