user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<pearson>:
  endbr64
  sub    $0x48,%rsp
  test   %rdx,%rdx
  je     L54
  cmp    $0x1,%rdx
  je     L35
  mov    %rdx,%rcx
  pxor   %xmm8,%xmm8
  xor    %eax,%eax
  shr    %rcx
  movapd %xmm8,%xmm5
  movapd %xmm8,%xmm7
  shl    $0x4,%rcx
  movapd %xmm8,%xmm6
  movapd %xmm8,%xmm3
  movapd %xmm8,%xmm4
  nopl   (%rax)
L10:
  movupd (%rdi,%rax,1),%xmm2
  movupd (%rsi,%rax,1),%xmm1
  add    $0x10,%rax
  movapd %xmm2,%xmm0
  addsd  %xmm4,%xmm0
  movapd %xmm2,%xmm4
  unpckhpd %xmm4,%xmm4
  addsd  %xmm0,%xmm4
  movapd %xmm1,%xmm0
  addsd  %xmm3,%xmm0
  movapd %xmm1,%xmm3
  unpckhpd %xmm3,%xmm3
  addsd  %xmm0,%xmm3
  movapd %xmm2,%xmm0
  mulpd  %xmm1,%xmm0
  mulpd  %xmm2,%xmm2
  mulpd  %xmm1,%xmm1
  addsd  %xmm0,%xmm6
  unpckhpd %xmm0,%xmm0
  addsd  %xmm0,%xmm6
  movapd %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm7,%xmm0
  movapd %xmm2,%xmm7
  addsd  %xmm0,%xmm7
  movapd %xmm1,%xmm0
  unpckhpd %xmm1,%xmm1
  addsd  %xmm5,%xmm0
  movapd %xmm1,%xmm5
  addsd  %xmm0,%xmm5
  cmp    %rcx,%rax
  jne    L10
  mov    %rdx,%rax
  and    $0xfffffffffffffffe,%rax
  test   $0x1,%dl
  je     L22
L14:
  movsd  (%rsi,%rax,8),%xmm1
  movsd  (%rdi,%rax,8),%xmm2
  movapd %xmm1,%xmm0
  addsd  %xmm2,%xmm4
  addsd  %xmm1,%xmm3
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm1
  mulsd  %xmm2,%xmm2
  addsd  %xmm0,%xmm6
  addsd  %xmm1,%xmm5
  addsd  %xmm2,%xmm7
L22:
  movapd %xmm3,%xmm2
  mulsd  %xmm4,%xmm2
  mulsd  %xmm4,%xmm4
  movsd  %xmm2,0x38(%rsp)
L74:
  pxor   %xmm2,%xmm2
  cvtsi2sd %rdx,%xmm2
  mulsd  %xmm2,%xmm7
  movsd  %xmm2,0x8(%rsp)
  subsd  %xmm4,%xmm7
  comisd D29(%rip),%xmm7
  jae    L98
  cvttsd2si %xmm7,%rdi
L73:
  movsd  %xmm8,0x30(%rsp)
  movsd  %xmm3,0x28(%rsp)
  movsd  %xmm6,0x20(%rsp)
  movsd  %xmm5,0x18(%rsp)
  call   <sqrt@plt>
  movsd  0x28(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm5
  mulsd  0x8(%rsp),%xmm5
  movsd  0x20(%rsp),%xmm6
  movsd  %xmm0,0x10(%rsp)
  mulsd  %xmm3,%xmm3
  movsd  0x30(%rsp),%xmm8
  subsd  %xmm3,%xmm5
  comisd D29(%rip),%xmm5
  jae    L52
  cvttsd2si %xmm5,%rdi
L0:
  movsd  %xmm8,0x18(%rsp)
  movsd  %xmm6,0x20(%rsp)
  call   <sqrt@plt>
  movsd  0x10(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x10(%rsp)
  call   <isnan@plt>
  movsd  0x18(%rsp),%xmm8
  movapd %xmm8,%xmm0
  test   %rax,%rax
  jne    L16
  movsd  0x8(%rsp),%xmm0
  movsd  0x20(%rsp),%xmm6
  movsd  0x10(%rsp),%xmm1
  mulsd  %xmm6,%xmm0
  subsd  0x38(%rsp),%xmm0
  divsd  %xmm1,%xmm0
L16:
  add    $0x48,%rsp
  ret
  nopl   0x0(%rax,%rax,1)
L52:
  subsd  D29(%rip),%xmm5
  cvttsd2si %xmm5,%rdi
  btc    $0x3f,%rdi
  jmp    L0
  nopl   0x0(%rax)
L98:
  subsd  D29(%rip),%xmm7
  cvttsd2si %xmm7,%rdi
  btc    $0x3f,%rdi
  jmp    L73
  nop
L54:
  pxor   %xmm8,%xmm8
  movapd %xmm8,%xmm4
  movapd %xmm8,%xmm5
  movapd %xmm8,%xmm7
  movsd  %xmm8,0x38(%rsp)
  movapd %xmm8,%xmm6
  movapd %xmm8,%xmm3
  jmp    L74
L35:
  pxor   %xmm8,%xmm8
  xor    %eax,%eax
  movapd %xmm8,%xmm5
  movapd %xmm8,%xmm7
  movapd %xmm8,%xmm6
  movapd %xmm8,%xmm3
  movapd %xmm8,%xmm4
  jmp    L14
```

tool_calls: name: disassemble_data
arguments: label: D29
type: f32

tool: D29: 0

assistant: ```c
double pearson(double *x, double *y, uint32_t len) {
  double sumx;
  double sumy;
  double sumxy;
  double sumx2;
  double sumy2;
  double numer;
  double denom;
  uint64_t i;
  sumx = 0;
  sumy = 0;
  sumxy = 0;
  sumx2 = 0;
  sumy2 = 0;
  for (i = 0; i < len; i++) {
    sumx += x[i];
    sumy += y[i];
    sumxy += x[i] * y[i];
    sumx2 += x[i] * x[i];
    sumy2 += y[i] * y[i];
  }
  numer = (len * sumxy) - (sumx * sumy);
  denom = sqrt((len * sumx2) - sumx * sumx) * sqrt((len * sumy2) - sumy * sumy);
  if (isnan(denom))
    return 0.0;
  return numer / denom;
}
```