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

```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    L117
  pxor   %xmm7,%xmm7
  cvtsi2ssl 0x18(%rsp),%xmm7
  movss  %xmm7,0x1c(%rsp)
L110:
  mov    0x18(%rsp),%edx
  lea    -0x1(%rdx),%eax
  mov    %eax,0x2c(%rsp)
  cmp    $0x2,%eax
  jbe    L50
  shr    $0x2,%edx
  pxor   %xmm6,%xmm6
  lea    0x4(%r14),%rax
  mov    %edx,%edi
  movss  %xmm6,0x28(%rsp)
  movaps %xmm6,%xmm1
  sub    $0x1,%edi
  shl    $0x4,%rdi
  lea    0x14(%r14,%rdi,1),%rdi
  nopw   0x0(%rax,%rax,1)
L18:
  movups (%rax),%xmm0
  add    $0x10,%rax
  movaps %xmm0,%xmm2
  addss  %xmm1,%xmm2
  movaps %xmm0,%xmm1
  shufps $0x55,%xmm0,%xmm1
  addss  %xmm1,%xmm2
  movaps %xmm0,%xmm1
  unpckhps %xmm0,%xmm1
  shufps $0xff,%xmm0,%xmm0
  addss  %xmm2,%xmm1
  addss  %xmm0,%xmm1
  cmp    %rdi,%rax
  jne    L18
  mov    0x18(%rsp),%ecx
  mov    %ecx,%edi
  and    $0xfffffffc,%edi
  lea    0x1(%rdi),%eax
  cmp    %ecx,%edi
  je     L93
L109:
  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     L93
  add    $0x2,%eax
  addss  0x4(%r14,%r8,1),%xmm1
  cmp    %eax,%ecx
  jl     L93
  addss  0x8(%r14,%r8,1),%xmm1
L93:
  mov    0x10(%rsp),%rax
  divss  0x1c(%rsp),%xmm1
  mov    $0x1,%ebx
  pxor   %xmm2,%xmm2
  movss  %xmm1,(%rax)
  movl   $0x0,0x0(%rbp)
  movl   $0x0,0x0(%r13)
  movl   $0x0,(%r15)
  movl   $0x0,(%r12)
  nop
L90:
  mov    0x10(%rsp),%rax
  movss  (%r14,%rbx,4),%xmm0
  add    $0x1,%rbx
  movss  %xmm2,0xc(%rsp)
  subss  (%rax),%xmm0
  movss  %xmm0,0x8(%rsp)
  call   <fabs@plt>
  pxor   %xmm1,%xmm1
  movss  0x8(%rsp),%xmm0
  movss  0xc(%rsp),%xmm2
  cvtsi2ss %rax,%xmm1
  addss  (%r12),%xmm1
  addss  %xmm0,%xmm2
  movss  %xmm1,(%r12)
  movaps %xmm0,%xmm1
  movss  (%r15),%xmm3
  mulss  %xmm0,%xmm1
  addss  %xmm1,%xmm3
  mulss  %xmm0,%xmm1
  movss  %xmm3,(%r15)
  movss  0x0(%r13),%xmm3
  addss  %xmm1,%xmm3
  mulss  %xmm1,%xmm0
  movss  %xmm3,0x0(%r13)
  addss  0x0(%rbp),%xmm0
  movss  %xmm0,0x0(%rbp)
  cmp    %ebx,0x18(%rsp)
  jge    L90
  mulss  %xmm2,%xmm2
  movss  0x1c(%rsp),%xmm5
  movss  (%r12),%xmm0
  divss  %xmm5,%xmm0
  divss  %xmm5,%xmm2
L73:
  movss  %xmm0,(%r12)
  pxor   %xmm1,%xmm1
  movss  (%r15),%xmm0
  cvtsi2ssl 0x2c(%rsp),%xmm1
  subss  %xmm2,%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,(%r15)
  call   <sqrt@plt>
  mov    0x20(%rsp),%rax
  movss  %xmm0,(%rax)
  movss  (%r15),%xmm1
  ucomiss 0x28(%rsp),%xmm1
  jp     L11
  je     L83
L11:
  movss  0x1c(%rsp),%xmm4
  mulss  %xmm4,%xmm1
  mulss  %xmm0,%xmm1
  movss  0x0(%r13),%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,0x0(%r13)
  movss  (%r15),%xmm0
  mulss  %xmm0,%xmm4
  movaps %xmm4,%xmm1
  mulss  %xmm0,%xmm1
  movss  0x0(%rbp),%xmm0
  divss  %xmm1,%xmm0
  subss  D17(%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,%rax,1)
L83:
  add    $0x38,%rsp
  lea    D2(%rip),%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1090 <nrerror@plt>
  nopw   0x0(%rax,%rax,1)
L117:
  lea    D14(%rip),%rdi
  call   <nrerror@plt>
  mov    0x18(%rsp),%eax
  pxor   %xmm6,%xmm6
  cvtsi2ss %eax,%xmm6
  movss  %xmm6,0x1c(%rsp)
  test   %eax,%eax
  jg     L110
  mov    0x10(%rsp),%rax
  pxor   %xmm6,%xmm6
  movaps %xmm6,%xmm0
  divss  0x1c(%rsp),%xmm0
  movss  %xmm6,0x28(%rsp)
  movaps %xmm0,%xmm2
  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    L73
L50:
  pxor   %xmm6,%xmm6
  mov    $0x1,%eax
  movss  %xmm6,0x28(%rsp)
  movaps %xmm6,%xmm1
  jmp    L109
```

tool_calls: name: unpack_data
arguments: label: D17
type_guessed: f32
name: unpack_data
arguments: label: D2
type_guessed: string
name: unpack_data
arguments: label: D14
type_guessed: string

tool: D17: 3

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

tool: D14: 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));
    ep += s;
    *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)");
}
```