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

```asm
<moment>:
  endbr64
  push   %r15
  mov    %r9,%r15
  push   %r14
  mov    %rdi,%r14
  push   %r13
  push   %r12
  mov    %rcx,%r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    0x70(%rsp),%r13
  mov    0x78(%rsp),%rbp
  mov    %esi,0x18(%rsp)
  mov    %rdx,0x10(%rsp)
  mov    %r8,0x20(%rsp)
  cmp    $0x1,%esi
  jle    L31
  pxor   %xmm7,%xmm7
  cvtsi2ssl 0x18(%rsp),%xmm7
  movss  %xmm7,0x1c(%rsp)
L8:
  mov    0x18(%rsp),%edx
  lea    -0x1(%rdx),%eax
  mov    %eax,0x2c(%rsp)
  cmp    $0x2,%eax
  jbe    L57
  shr    $0x2,%edx
  pxor   %xmm5,%xmm5
  lea    0x4(%r14),%rax
  mov    %edx,%edi
  movss  %xmm5,0x28(%rsp)
  movaps %xmm5,%xmm1
  sub    $0x1,%edi
  shl    $0x4,%rdi
  lea    0x14(%r14,%rdi,1),%rdi
  nopw   0x0(%rax,%rax,1)
L39:
  movups (%rax),%xmm0
  add    $0x10,%rax
  movaps %xmm0,%xmm2
  addss  %xmm1,%xmm2
  movaps %xmm0,%xmm1
  shufps $0x55,%xmm0,%xmm1
  addss  %xmm2,%xmm1
  movaps %xmm0,%xmm2
  unpckhps %xmm0,%xmm2
  shufps $0xff,%xmm0,%xmm0
  addss  %xmm2,%xmm1
  addss  %xmm0,%xmm1
  cmp    %rdi,%rax
  jne    L39
  mov    0x18(%rsp),%ecx
  mov    %ecx,%edi
  and    $0xfffffffc,%edi
  lea    0x1(%rdi),%eax
  cmp    %ecx,%edi
  je     L0
L7:
  movslq %eax,%rdi
  mov    0x18(%rsp),%ecx
  addss  (%r14,%rdi,4),%xmm1
  lea    0x0(,%rdi,4),%r8
  lea    0x1(%rax),%edi
  cmp    %edi,%ecx
  jl     L0
  add    $0x2,%eax
  addss  0x4(%r14,%r8,1),%xmm1
  cmp    %eax,%ecx
  jl     L0
  addss  0x8(%r14,%r8,1),%xmm1
L0:
  mov    0x10(%rsp),%rax
  divss  0x1c(%rsp),%xmm1
  mov    $0x1,%ebx
  movss  %xmm1,(%rax)
  movl   $0x0,0x0(%rbp)
  movl   $0x0,0x0(%r13)
  movl   $0x0,(%r15)
  movl   $0x0,(%r12)
  nopl   0x0(%rax,%rax,1)
L50:
  mov    0x10(%rsp),%rax
  movss  (%r14,%rbx,4),%xmm0
  add    $0x1,%rbx
  subss  (%rax),%xmm0
  movss  %xmm0,0xc(%rsp)
  call   <fabs@plt>
  pxor   %xmm1,%xmm1
  movss  0xc(%rsp),%xmm0
  cvtsi2ss %rax,%xmm1
  addss  (%r12),%xmm1
  movss  %xmm1,(%r12)
  movaps %xmm0,%xmm1
  movss  (%r15),%xmm2
  mulss  %xmm0,%xmm1
  addss  %xmm1,%xmm2
  mulss  %xmm0,%xmm1
  movss  %xmm2,(%r15)
  movss  0x0(%r13),%xmm2
  addss  %xmm1,%xmm2
  mulss  %xmm1,%xmm0
  movss  %xmm2,0x0(%r13)
  addss  0x0(%rbp),%xmm0
  movss  %xmm0,0x0(%rbp)
  cmp    %ebx,0x18(%rsp)
  jge    L50
  movss  0x1c(%rsp),%xmm4
  movss  (%r12),%xmm0
  movss  0x28(%rsp),%xmm1
  divss  %xmm4,%xmm0
  divss  %xmm4,%xmm1
L86:
  movss  %xmm0,(%r12)
  movss  (%r15),%xmm0
  subss  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2ssl 0x2c(%rsp),%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,(%r15)
  call   <sqrt@plt>
  mov    0x20(%rsp),%rax
  movss  %xmm0,(%rax)
  movss  (%r15),%xmm1
  ucomiss 0x28(%rsp),%xmm1
  jp     L113
  je     L45
L113:
  movss  0x1c(%rsp),%xmm3
  mulss  %xmm3,%xmm1
  mulss  %xmm0,%xmm1
  movss  0x0(%r13),%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,0x0(%r13)
  movss  (%r15),%xmm0
  mulss  %xmm0,%xmm3
  movaps %xmm3,%xmm1
  mulss  %xmm0,%xmm1
  movss  0x0(%rbp),%xmm0
  divss  %xmm1,%xmm0
  subss  D27(%rip),%xmm0
  movss  %xmm0,0x0(%rbp)
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L45:
  add    $0x38,%rsp
  lea    D11(%rip),%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1090 <nrerror@plt>
  nopw   0x0(%rax,%rax,1)
L31:
  lea    D23(%rip),%rdi
  call   <nrerror@plt>
  mov    0x18(%rsp),%eax
  pxor   %xmm6,%xmm6
  cvtsi2ss %eax,%xmm6
  movss  %xmm6,0x1c(%rsp)
  test   %eax,%eax
  jg     L8
  mov    0x10(%rsp),%rax
  pxor   %xmm6,%xmm6
  movaps %xmm6,%xmm0
  divss  0x1c(%rsp),%xmm0
  movss  %xmm6,0x28(%rsp)
  movaps %xmm0,%xmm1
  movss  %xmm0,(%rax)
  mov    0x18(%rsp),%eax
  movl   $0x0,0x0(%rbp)
  sub    $0x1,%eax
  movl   $0x0,0x0(%r13)
  mov    %eax,0x2c(%rsp)
  movl   $0x0,(%r15)
  jmp    L86
L57:
  pxor   %xmm6,%xmm6
  mov    $0x1,%eax
  movss  %xmm6,0x28(%rsp)
  movaps %xmm6,%xmm1
  jmp    L7
```

tool_calls: name: extract_symbolic_data
arguments: label: D27
type_guessed: f32
name: extract_symbolic_data
arguments: label: D11
type_guessed: string
name: extract_symbolic_data
arguments: label: D23
type_guessed: string

tool: D27: 3

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

tool: D23: n must be at least 2 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)");
}
```