user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<c_corr_prune>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x50,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x48(%rsp)
  xor    %eax,%eax
  lea    D12(%rip),%rsi
  call   <fopen@plt>
  mov    %rax,%rbx
  test   %rax,%rax
  setne  %dil
  movzbl %dil,%edi
  call   <assert@plt>
  mov    $0x0,%r12d
  mov    D63(%rip),%rbp
  mov    %rbp,0x20(%rsp)
  mov    %rbp,0x8(%rsp)
  mov    %rbp,0x18(%rsp)
  mov    %rbp,0x10(%rsp)
  lea    0x38(%rsp),%r13
  lea    0x40(%rsp),%r14
L93:
  mov    %rbx,%rcx
  mov    $0x1,%edx
  mov    $0x8,%esi
  mov    %r13,%rdi
  call   <fread@plt>
  cmp    $0x1,%eax
  jne    L133
  mov    %rbx,%rcx
  mov    $0x1,%edx
  mov    $0x8,%esi
  mov    %r14,%rdi
  call   <fread@plt>
  cmp    $0x1,%eax
  jne    L133
  mov    %rbx,%rdi
  call   <feof@plt>
  test   %rax,%rax
  jne    L133
  movsd  0x38(%rsp),%xmm1
  comisd D33(%rip),%xmm1
  jae    L93
  movsd  0x40(%rsp),%xmm0
  comisd D33(%rip),%xmm0
  jae    L93
  movsd  0x10(%rsp),%xmm3
  addsd  %xmm1,%xmm3
  movsd  %xmm3,0x10(%rsp)
  movapd %xmm1,%xmm2
  mulsd  %xmm1,%xmm2
  addsd  0x18(%rsp),%xmm2
  movsd  %xmm2,0x18(%rsp)
  movsd  0x8(%rsp),%xmm4
  addsd  %xmm0,%xmm4
  movsd  %xmm4,0x8(%rsp)
  mulsd  %xmm0,%xmm1
  movq   %rbp,%xmm5
  addsd  %xmm1,%xmm5
  movq   %xmm5,%rbp
  mulsd  %xmm0,%xmm0
  addsd  0x20(%rsp),%xmm0
  movsd  %xmm0,0x20(%rsp)
  add    $0x1,%r12
  jmp    L93
L133:
  mov    %rbx,%rdi
  call   <fclose@plt>
  test   %r12,%r12
  js     L86
  pxor   %xmm4,%xmm4
  cvtsi2sd %r12,%xmm4
  movsd  %xmm4,0x28(%rsp)
L20:
  movq   %rbp,%xmm1
  movsd  0x28(%rsp),%xmm7
  mulsd  %xmm7,%xmm1
  movsd  0x10(%rsp),%xmm6
  movapd %xmm6,%xmm0
  mulsd  0x8(%rsp),%xmm0
  subsd  %xmm0,%xmm1
  movq   %xmm1,%rbp
  movsd  0x18(%rsp),%xmm0
  mulsd  %xmm7,%xmm0
  movapd %xmm6,%xmm1
  mulsd  %xmm6,%xmm1
  subsd  %xmm1,%xmm0
  comisd D6(%rip),%xmm0
  jae    L54
  cvttsd2si %xmm0,%rdi
L60:
  call   <sqrt@plt>
  test   %rax,%rax
  js     L16
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
L113:
  movq   %rbp,%xmm3
  divsd  %xmm0,%xmm3
  movq   %xmm3,%rbp
  movsd  0x28(%rsp),%xmm1
  mulsd  0x20(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  mulsd  %xmm0,%xmm0
  subsd  %xmm0,%xmm1
  comisd D6(%rip),%xmm1
  jae    L69
  cvttsd2si %xmm1,%rdi
L9:
  call   <sqrt@plt>
  test   %rax,%rax
  js     L81
  pxor   %xmm1,%xmm1
  cvtsi2sd %rax,%xmm1
L40:
  movq   %rbp,%xmm0
  divsd  %xmm1,%xmm0
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L102
  add    $0x50,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L86:
  mov    %r12,%rax
  shr    %rax
  and    $0x1,%r12d
  or     %r12,%rax
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,0x28(%rsp)
  jmp    L20
L54:
  subsd  D6(%rip),%xmm0
  cvttsd2si %xmm0,%rdi
  btc    $0x3f,%rdi
  jmp    L60
L16:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
  jmp    L113
L69:
  subsd  D6(%rip),%xmm1
  cvttsd2si %xmm1,%rdi
  btc    $0x3f,%rdi
  jmp    L9
L81:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm1,%xmm1
  cvtsi2sd %rdx,%xmm1
  addsd  %xmm1,%xmm1
  jmp    L40
L102:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: data_label: D12
type_guessed: string
name: disassemble_data
arguments: data_label: D63
type_guessed: f64
name: disassemble_data
arguments: data_label: D33
type_guessed: f64
name: disassemble_data
arguments: data_label: D6
type_guessed: byte[16]

tool: D12: rb

tool: D63: 0

tool: D33: 0.25

tool: {"D6": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x43,0x00,0x00,0x00,0x00,0x00,0x00,0xd0,0x3f]}

assistant: ```c
double c_corr_prune(const char f_name[]) {
  FILE *const pf = fopen(f_name, "rb");
  assert(pf != NULL);
  double sx = 0, sxx = 0, sy = 0, syy = 0, sxy = 0;
  size_t n = 0;
  while (1) {
    double x, y;
    if (fread(&x, sizeof(double), 1, pf) != 1 || fread(&y, sizeof(double), 1, pf) != 1 || feof(pf)) {
      fclose(pf);
      break;
    }
    if (x >= 0.25 || y >= 0.25)
      continue;
    sx += x;
    sxx += x * x;
    sy += y;
    sxy += x * y;
    syy += y * y;
    ++n;
  }
  return (n * sxy - sx * sy) / sqrt(n * sxx - sx * sx) / sqrt(n * syy - sy * sy);
}
```