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

```asm
<stat_pearson>:
  endbr64
  push   %r12
  xor    %eax,%eax
  mov    %edx,%r12d
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x20,%rsp
  call   <assert@plt>
  mov    %rbp,%rdi
  xor    %eax,%eax
  call   <assert@plt>
  xor    %edi,%edi
  test   %r12d,%r12d
  setg   %dil
  xor    %eax,%eax
  call   <assert@plt>
  mov    %r12d,%esi
  mov    %rbx,%rdi
  xor    %eax,%eax
  call   <stat_double_average@plt>
  pxor   %xmm5,%xmm5
  mov    %r12d,%esi
  mov    %rbx,%rdi
  cvtsi2sd %eax,%xmm5
  xor    %eax,%eax
  movsd  %xmm5,0x8(%rsp)
  call   <stat_double_variance@plt>
  pxor   %xmm3,%xmm3
  pxor   %xmm4,%xmm4
  movsd  0x8(%rsp),%xmm5
  cvtsi2sd %eax,%xmm3
  ucomisd %xmm3,%xmm4
  ja     L18
  sqrtsd %xmm3,%xmm3
L20:
  mov    %r12d,%esi
  mov    %rbp,%rdi
  xor    %eax,%eax
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm5,0x10(%rsp)
  call   <stat_double_average@plt>
  pxor   %xmm6,%xmm6
  mov    %r12d,%esi
  mov    %rbp,%rdi
  cvtsi2sd %eax,%xmm6
  xor    %eax,%eax
  movsd  %xmm6,0x8(%rsp)
  call   <stat_double_variance@plt>
  pxor   %xmm1,%xmm1
  pxor   %xmm4,%xmm4
  movsd  0x8(%rsp),%xmm6
  cvtsi2sd %eax,%xmm1
  movsd  0x10(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm3
  ucomisd %xmm1,%xmm4
  ja     L23
  sqrtsd %xmm1,%xmm1
L4:
  test   %r12d,%r12d
  jle    L46
  lea    -0x1(%r12),%eax
  pxor   %xmm4,%xmm4
  cvtsi2sd %eax,%xmm4
  cmp    $0x1,%r12d
  je     L10
  mov    %r12d,%edx
  movapd %xmm5,%xmm11
  movapd %xmm3,%xmm10
  xor    %eax,%eax
  movapd %xmm6,%xmm9
  movapd %xmm1,%xmm8
  movapd %xmm4,%xmm7
  shr    %edx
  unpcklpd %xmm11,%xmm11
  unpcklpd %xmm10,%xmm10
  unpcklpd %xmm9,%xmm9
  shl    $0x4,%rdx
  unpcklpd %xmm8,%xmm8
  unpcklpd %xmm7,%xmm7
  pxor   %xmm2,%xmm2
  nopl   (%rax)
L54:
  movupd 0x0(%rbp,%rax,1),%xmm0
  movupd (%rbx,%rax,1),%xmm12
  add    $0x10,%rax
  subpd  %xmm9,%xmm0
  subpd  %xmm11,%xmm12
  divpd  %xmm8,%xmm0
  divpd  %xmm10,%xmm12
  mulpd  %xmm12,%xmm0
  divpd  %xmm7,%xmm0
  movapd %xmm0,%xmm12
  unpckhpd %xmm0,%xmm0
  addsd  %xmm2,%xmm12
  movapd %xmm0,%xmm2
  addsd  %xmm12,%xmm2
  cmp    %rax,%rdx
  jne    L54
  mov    %r12d,%eax
  and    $0xfffffffe,%eax
  and    $0x1,%r12d
  je     L8
L59:
  cltq
  movsd  0x0(%rbp,%rax,8),%xmm0
  subsd  %xmm6,%xmm0
  divsd  %xmm1,%xmm0
  movsd  (%rbx,%rax,8),%xmm1
  subsd  %xmm5,%xmm1
  divsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm0
  divsd  %xmm4,%xmm0
  addsd  %xmm0,%xmm2
L8:
  add    $0x20,%rsp
  movapd %xmm2,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
  nopl   0x0(%rax,%rax,1)
L46:
  add    $0x20,%rsp
  pxor   %xmm2,%xmm2
  pop    %rbx
  movapd %xmm2,%xmm0
  pop    %rbp
  pop    %r12
  ret
L10:
  xor    %eax,%eax
  pxor   %xmm2,%xmm2
  jmp    L59
L18:
  movapd %xmm3,%xmm0
  movsd  %xmm5,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm5
  movapd %xmm0,%xmm3
  jmp    L20
L23:
  movapd %xmm1,%xmm0
  movsd  %xmm6,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm5,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm6
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm5
  movapd %xmm0,%xmm1
  jmp    L4
```

assistant: ```c
double stat_pearson(const double *x, const double *y, int n) {
  assert(x);
  assert(y);
  assert(n > 0);
  double x_mean = stat_double_average(x, n);
  double x_stddev = sqrt(stat_double_variance(x, n));
  double y_mean = stat_double_average(y, n);
  double y_stddev = sqrt(stat_double_variance(y, n));
  double r = 0.0;
  int i;
  for (i = 0; i < n; i++) {
    double x_score = (x[i] - x_mean) / x_stddev;
    double y_score = (y[i] - y_mean) / y_stddev;
    r += x_score * y_score / (n - 1);
  }
  return r;
}
```