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

```asm
<SNPHWE_pValue>:
  endbr64
  push   %r14
  mov    %edx,%eax
  push   %r13
  push   %r12
  mov    %ecx,%r12d
  push   %rbp
  movslq %edi,%rbp
  push   %rbx
  sub    $0x10,%rsp
  cmp    %edx,%esi
  cmovle %esi,%eax
  cmovge %esi,%edx
  lea    (%rax,%rax,1),%r13d
  mov    %edx,%r14d
  lea    0x0(%r13,%rbp,1),%ebx
  add    %ebp,%r14d
  add    %eax,%r14d
  lea    0x8(,%rbx,8),%edi
  call   <malloc@plt>
  lea    (%r14,%r14,1),%ecx
  movsd  D44(%rip),%xmm4
  mov    %rax,%rdi
  mov    %ecx,%eax
  sub    %ebx,%eax
  imul   %ebx,%eax
  cltd
  idiv   %ecx
  mov    %ebx,%edx
  xor    %eax,%edx
  and    $0x1,%edx
  cmp    $0x1,%edx
  mov    %ebx,%edx
  sbb    $0xffffffff,%eax
  sub    %eax,%edx
  sub    %eax,%r14d
  mov    %edx,%esi
  mov    %r14d,%r9d
  shr    $0x1f,%esi
  add    %edx,%esi
  movslq %eax,%rdx
  sar    %esi
  movsd  %xmm4,(%rdi,%rdx,8)
  sub    %esi,%r9d
  cmp    $0x1,%eax
  jle    L19
  mov    %r9d,%r10d
  mov    %rdx,%rcx
  movapd %xmm4,%xmm0
  mov    %esi,%r8d
  movsd  D45(%rip),%xmm5
  movapd %xmm4,%xmm3
  sub    %esi,%r10d
  jmp    L76
  nopl   0x0(%rax)
L140:
  movsd  (%rdi,%rcx,8),%xmm0
L76:
  pxor   %xmm1,%xmm1
  pxor   %xmm2,%xmm2
  lea    (%r10,%r8,1),%r11d
  cvtsi2sd %ecx,%xmm1
  cvtsi2sd %r8d,%xmm2
  add    $0x1,%r8d
  mulsd  %xmm1,%xmm0
  subsd  %xmm4,%xmm1
  addsd  %xmm4,%xmm2
  mulsd  %xmm5,%xmm2
  mulsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %r11d,%xmm1
  addsd  %xmm4,%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rdi,%rcx,8)
  sub    $0x2,%rcx
  addsd  %xmm0,%xmm3
  cmp    $0x1,%ecx
  jg     L140
L106:
  lea    -0x1(%rbx),%ecx
  cmp    %ecx,%eax
  jge    L109
  movsd  D45(%rip),%xmm5
  movsd  D48(%rip),%xmm6
  sub    %esi,%r9d
  nopl   0x0(%rax)
L113:
  movsd  (%rdi,%rdx,8),%xmm0
  pxor   %xmm2,%xmm2
  lea    (%r9,%rsi,1),%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %esi,%xmm2
  sub    $0x1,%esi
  mulsd  %xmm5,%xmm0
  cvtsi2sd %edx,%xmm1
  mulsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  mulsd  %xmm2,%xmm0
  movapd %xmm1,%xmm2
  addsd  %xmm4,%xmm1
  addsd  %xmm6,%xmm2
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%rdi,%rdx,8)
  add    $0x2,%rdx
  addsd  %xmm0,%xmm3
  cmp    %edx,%ecx
  jg     L113
L109:
  test   %ebx,%ebx
  js     L165
  movslq %ebx,%rax
  mov    %rdi,%rdx
  lea    0x8(%rdi,%rax,8),%rcx
  mov    %rdi,%rax
  nopl   0x0(%rax)
L59:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  divsd  %xmm3,%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rcx
  jne    L59
  cmp    $0x1,%r12d
  je     L128
L134:
  cmp    $0x2,%r12d
  je     L68
L120:
  call   <free@plt>
  add    $0x10,%rsp
  pxor   %xmm0,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L68:
  cmp    %ebx,%ebp
  jg     L120
  lea    0x1(%rbp,%r13,1),%rdx
  lea    (%rdi,%rbp,8),%rax
  pxor   %xmm0,%xmm0
  lea    (%rdi,%rdx,8),%rdx
  nopl   (%rax)
L83:
  addsd  (%rax),%xmm0
  add    $0x8,%rax
  cmp    %rax,%rdx
  jne    L83
L169:
  movsd  %xmm0,0x8(%rsp)
  call   <free@plt>
  mov    D44(%rip),%rax
  movsd  0x8(%rsp),%xmm0
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  movq   %rax,%xmm4
  pop    %r12
  pop    %r13
  minsd  %xmm0,%xmm4
  pop    %r14
  movapd %xmm4,%xmm0
  ret
  nopw   0x0(%rax,%rax,1)
L128:
  movsd  (%rdi,%rbp,8),%xmm2
  pxor   %xmm0,%xmm0
  nopl   0x0(%rax)
L111:
  movsd  (%rdx),%xmm1
  comisd %xmm2,%xmm1
  ja     L143
  addsd  %xmm1,%xmm0
L143:
  add    $0x8,%rdx
  cmp    %rdx,%rcx
  jne    L111
  jmp    L169
  nopl   0x0(%rax)
L19:
  movapd %xmm4,%xmm3
  jmp    L106
  nopl   0x0(%rax)
L165:
  cmp    $0x1,%r12d
  jne    L134
  jmp    L120
```

tool_calls: name: retrieve_data
arguments: data_label: D44
type_infered: f64
name: retrieve_data
arguments: data_label: D45
type_infered: f64
name: retrieve_data
arguments: data_label: D48
type_infered: f64

tool: D44: 1

tool: D45: 4

tool: D48: 2

assistant: ```c
double SNPHWE_pValue(int obs_hets, int obs_hom1, int obs_hom2, int ctype) {
  int i;
  double p_hwe;
  int obs_homc = (obs_hom1 < obs_hom2) ? obs_hom2 : obs_hom1;
  int obs_homr = (obs_hom1 < obs_hom2) ? obs_hom1 : obs_hom2;
  int rare_copies = 2 * obs_homr + obs_hets;
  int genotypes = obs_hets + obs_homc + obs_homr;
  double *het_probs = malloc((rare_copies + 1) * sizeof(double));
  int mid = rare_copies * (2 * genotypes - rare_copies) / (2 * genotypes);
  if ((rare_copies & 1) ^ (mid & 1)) {
    mid++;
  }
  int curr_hets = mid;
  int curr_homr = (rare_copies - mid) / 2;
  int curr_homc = genotypes - curr_hets - curr_homr;
  het_probs[mid] = 1.0;
  double sum = het_probs[mid];
  for (curr_hets = mid; curr_hets > 1; curr_hets -= 2) {
    het_probs[curr_hets - 2] = het_probs[curr_hets] * curr_hets * (curr_hets - 1.0) / (4.0 * (curr_homr + 1.0) * (curr_homc + 1.0));
    sum += het_probs[curr_hets - 2];
    curr_homr++;
    curr_homc++;
  }
  curr_hets = mid;
  curr_homr = (rare_copies - mid) / 2;
  curr_homc = genotypes - curr_hets - curr_homr;
  for (curr_hets = mid; curr_hets <= rare_copies - 2; curr_hets += 2) {
    het_probs[curr_hets + 2] = het_probs[curr_hets] * 4.0 * curr_homr * curr_homc / ((curr_hets + 2.0) * (curr_hets + 1.0));
    sum += het_probs[curr_hets + 2];
    curr_homr--;
    curr_homc--;
  }
  for (i = 0; i <= rare_copies; i++) {
    het_probs[i] /= sum;
  }
  p_hwe = 0.0;
  if (ctype == 1) {
    for (i = 0; i <= rare_copies; i++) {
      if (het_probs[i] > het_probs[obs_hets])
        continue;
      p_hwe += het_probs[i];
    }
  } else if (ctype == 2) {
    for (i = obs_hets; i <= rare_copies; i++) {
      p_hwe += het_probs[i];
    }
  }
  free(het_probs);
  return (p_hwe > 1.0) ? 1.0 : p_hwe;
}
```