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

```asm
<moment>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rdi,%r15
  mov    %esi,0x10(%rsp)
  mov    %rdx,0x8(%rsp)
  mov    %rcx,%r12
  mov    %r8,0x18(%rsp)
  mov    %r9,%rbp
  mov    0x60(%rsp),%r13
  mov    0x68(%rsp),%r14
  cmp    $0x1,%esi
  jle    L82
L72:
  mov    0x10(%rsp),%ecx
  test   %ecx,%ecx
  jle    L25
  lea    0x4(%r15),%rax
  lea    -0x1(%rcx),%edx
  lea    0x8(%r15,%rdx,4),%rdx
  pxor   %xmm0,%xmm0
L85:
  addss  (%rax),%xmm0
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L85
  pxor   %xmm5,%xmm5
  cvtsi2ssl 0x10(%rsp),%xmm5
  movss  %xmm5,0x14(%rsp)
  divss  %xmm5,%xmm0
  mov    0x8(%rsp),%rax
  movss  %xmm0,(%rax)
  movl   $0x0,(%r14)
  movl   $0x0,0x0(%r13)
  movl   $0x0,0x0(%rbp)
  movl   $0x0,(%r12)
  mov    $0x1,%ebx
L73:
  movss  (%r15,%rbx,4),%xmm0
  mov    0x8(%rsp),%rax
  subss  (%rax),%xmm0
  movss  %xmm0,0x4(%rsp)
  call   <fabs@plt>
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  addss  (%r12),%xmm0
  movss  %xmm0,(%r12)
  movss  0x4(%rsp),%xmm2
  movaps %xmm2,%xmm0
  mulss  %xmm2,%xmm0
  movaps %xmm0,%xmm1
  addss  0x0(%rbp),%xmm1
  movss  %xmm1,0x0(%rbp)
  mulss  %xmm2,%xmm0
  movaps %xmm0,%xmm1
  addss  0x0(%r13),%xmm1
  movss  %xmm1,0x0(%r13)
  mulss  %xmm2,%xmm0
  addss  (%r14),%xmm0
  movss  %xmm0,(%r14)
  add    $0x1,%rbx
  cmp    %ebx,0x10(%rsp)
  jge    L73
L78:
  movss  (%r12),%xmm0
  movss  0x14(%rsp),%xmm3
  divss  %xmm3,%xmm0
  movss  %xmm0,(%r12)
  pxor   %xmm1,%xmm1
  divss  %xmm3,%xmm1
  movss  0x0(%rbp),%xmm0
  subss  %xmm1,%xmm0
  mov    0x10(%rsp),%eax
  sub    $0x1,%eax
  pxor   %xmm1,%xmm1
  cvtsi2ss %eax,%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,0x0(%rbp)
  call   <sqrt@plt>
  mov    0x18(%rsp),%rax
  movss  %xmm0,(%rax)
  movss  0x0(%rbp),%xmm1
  pxor   %xmm6,%xmm6
  ucomiss %xmm6,%xmm1
  jp     L4
  je     L21
L4:
  movss  0x14(%rsp),%xmm4
  mulss  %xmm4,%xmm1
  mulss  %xmm0,%xmm1
  movss  0x0(%r13),%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,0x0(%r13)
  movss  0x0(%rbp),%xmm0
  mulss  %xmm0,%xmm4
  movaps %xmm4,%xmm1
  mulss  %xmm0,%xmm1
  movss  (%r14),%xmm0
  divss  %xmm1,%xmm0
  subss  D18(%rip),%xmm0
  movss  %xmm0,(%r14)
L74:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L82:
  lea    D11(%rip),%rdi
  call   <nrerror@plt>
  jmp    L72
L21:
  lea    D25(%rip),%rdi
  call   <nrerror@plt>
  jmp    L74
L25:
  pxor   %xmm7,%xmm7
  cvtsi2ssl 0x10(%rsp),%xmm7
  movss  %xmm7,0x14(%rsp)
  pxor   %xmm0,%xmm0
  divss  %xmm7,%xmm0
  mov    0x8(%rsp),%rax
  movss  %xmm0,(%rax)
  movl   $0x0,(%r14)
  movl   $0x0,0x0(%r13)
  movl   $0x0,0x0(%rbp)
  movl   $0x0,(%r12)
  jmp    L78
```

tool_calls: name: disassemble_data
arguments: tag: D18
datatype: f32
name: disassemble_data
arguments: tag: D11
datatype: string
name: disassemble_data
arguments: tag: D25
datatype: string

tool: D18: 3

tool: D11: n must be at least 2 in moment

tool: D25: No skew/kurtosis when variance = 0 (in moment)

assistant: ```c
void moment(float data[], int n, float *ave, float *adev, float *sdev, float *var, float *skew, float *curt) {
  void nrerror(char error_text[]);
  int j;
  float ep = 0.0, s, p;
  if (n <= 1)
    nrerror("n must be at least 2 in moment");
  s = 0.0;
  for (j = 1; j <= n; j++)
    s += data[j];
  *ave = s / n;
  *adev = (*var) = (*skew) = (*curt) = 0.0;
  for (j = 1; j <= n; j++) {
    *adev += fabs(s = data[j] - (*ave));
    *var += (p = s * s);
    *skew += (p *= s);
    *curt += (p *= s);
  }
  *adev /= n;
  *var = (*var - ep * ep / n) / (n - 1);
  *sdev = sqrt(*var);
  if (*var) {
    *skew /= (n * (*var) * (*sdev));
    *curt = (*curt) / (n * (*var) * (*var)) - 3.0;
  } else
    nrerror("No skew/kurtosis when variance = 0 (in moment)");
}
```