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

```asm
<SNPHWE_t>:
  endbr64
  mov    %rsi,%rcx
  cmp    %rdx,%rsi
  jl     L336
  mov    %rdx,%rax
  mov    %rsi,%rdx
  mov    %rax,%rcx
L336:
  lea    (%rdx,%rdi,1),%rax
  add    %rcx,%rax
  je     L62
  lea    (%rdi,%rcx,2),%rsi
  add    %rax,%rax
  pxor   %xmm1,%xmm1
  cvtsi2sd %rdi,%xmm1
  pxor   %xmm2,%xmm2
  cvtsi2sd %rcx,%xmm2
  pxor   %xmm3,%xmm3
  cvtsi2sd %rdx,%xmm3
  movsd  D63(%rip),%xmm8
  subsd  %xmm0,%xmm8
  divsd  %xmm0,%xmm8
  mov    %rdi,%r8
  imul   %rax,%r8
  sub    %rsi,%rax
  imul   %rsi,%rax
  cmp    %rax,%r8
  jle    L270
  mov    $0x0,%eax
  cmp    $0x1,%rdi
  jle    L62
  pxor   %xmm10,%xmm10
  cvtsi2sd %rsi,%xmm10
  mulsd  %xmm8,%xmm10
  mulsd  D224(%rip),%xmm10
  movapd %xmm3,%xmm5
  movapd %xmm2,%xmm4
  movsd  D21(%rip),%xmm0
  pxor   %xmm6,%xmm6
  movsd  D63(%rip),%xmm9
  movsd  D181(%rip),%xmm13
  movsd  D187(%rip),%xmm12
  movsd  D224(%rip),%xmm11
  movsd  D156(%rip),%xmm14
L46:
  addsd  %xmm9,%xmm4
  addsd  %xmm9,%xmm5
  movapd %xmm1,%xmm7
  subsd  %xmm9,%xmm7
  mulsd  %xmm1,%xmm7
  movapd %xmm4,%xmm15
  mulsd  %xmm13,%xmm15
  mulsd  %xmm5,%xmm15
  divsd  %xmm15,%xmm7
  mulsd  %xmm7,%xmm0
  subsd  %xmm12,%xmm1
  comisd %xmm0,%xmm11
  ja     L317
  addsd  %xmm0,%xmm6
  comisd %xmm10,%xmm6
  ja     L342
  comisd %xmm14,%xmm1
  ja     L46
  pxor   %xmm13,%xmm13
  jmp    L98
L317:
  movapd %xmm0,%xmm13
L98:
  divsd  %xmm8,%xmm6
  mov    $0x0,%eax
  movapd %xmm13,%xmm7
  addsd  D21(%rip),%xmm7
  comisd %xmm6,%xmm7
  jae    L62
  movsd  D63(%rip),%xmm7
  movapd %xmm1,%xmm9
  subsd  %xmm7,%xmm9
  mulsd  %xmm1,%xmm9
  movapd %xmm4,%xmm10
  addsd  %xmm7,%xmm10
  movsd  D181(%rip),%xmm8
  mulsd  %xmm8,%xmm10
  movapd %xmm5,%xmm11
  addsd  %xmm7,%xmm11
  mulsd  %xmm11,%xmm10
  divsd  %xmm10,%xmm9
  movapd %xmm7,%xmm10
  subsd  %xmm9,%xmm10
  movapd %xmm13,%xmm14
  divsd  %xmm10,%xmm14
  movsd  %xmm14,-0x8(%rsp)
  add    $0x2,%rdi
  pxor   %xmm10,%xmm10
  cvtsi2sd %rdi,%xmm10
  mulsd  %xmm2,%xmm8
  mulsd  %xmm3,%xmm8
  movapd %xmm10,%xmm9
  subsd  %xmm7,%xmm9
  mulsd  %xmm10,%xmm9
  divsd  %xmm9,%xmm8
  subsd  %xmm8,%xmm7
  movsd  D21(%rip),%xmm9
  divsd  %xmm7,%xmm9
  movapd %xmm9,%xmm7
  addsd  %xmm14,%xmm7
  mov    $0x1,%eax
  comisd %xmm7,%xmm6
  ja     L62
  movsd  D21(%rip),%xmm9
  mulsd  %xmm9,%xmm8
  addsd  %xmm8,%xmm9
  cmp    $0x1,%rcx
  jle    L197
  movapd %xmm6,%xmm15
  subsd  %xmm13,%xmm15
  movsd  D187(%rip),%xmm14
  movsd  D63(%rip),%xmm12
L282:
  addsd  %xmm14,%xmm10
  subsd  %xmm12,%xmm2
  subsd  %xmm12,%xmm3
  movsd  D181(%rip),%xmm7
  mulsd  %xmm2,%xmm7
  mulsd  %xmm3,%xmm7
  movapd %xmm10,%xmm11
  subsd  %xmm12,%xmm11
  mulsd  %xmm10,%xmm11
  divsd  %xmm11,%xmm7
  mulsd  %xmm7,%xmm8
  movapd %xmm9,%xmm7
  addsd  %xmm8,%xmm9
  comisd %xmm15,%xmm9
  ja     L89
  comisd %xmm9,%xmm7
  jae    L197
  comisd D156(%rip),%xmm2
  ja     L282
L197:
  movsd  -0x8(%rsp),%xmm2
  addsd  %xmm9,%xmm2
  mov    $0x1,%eax
  comisd %xmm2,%xmm6
  ja     L62
  subsd  %xmm9,%xmm6
  comisd D63(%rip),%xmm1
  jbe    L133
  movsd  D63(%rip),%xmm7
  movsd  D181(%rip),%xmm8
  movsd  D187(%rip),%xmm10
  movapd %xmm7,%xmm9
L340:
  addsd  %xmm7,%xmm4
  addsd  %xmm7,%xmm5
  movapd %xmm1,%xmm2
  subsd  %xmm7,%xmm2
  mulsd  %xmm1,%xmm2
  movapd %xmm4,%xmm3
  mulsd  %xmm8,%xmm3
  mulsd  %xmm5,%xmm3
  divsd  %xmm3,%xmm2
  mulsd  %xmm2,%xmm0
  movapd %xmm13,%xmm2
  addsd  %xmm0,%xmm13
  comisd %xmm6,%xmm13
  jae    L57
  comisd %xmm13,%xmm2
  jae    L118
  subsd  %xmm10,%xmm1
  comisd %xmm9,%xmm1
  ja     L340
  mov    $0x1,%eax
  ret
L270:
  mov    %rcx,%rax
  test   %rcx,%rcx
  je     L62
  pxor   %xmm10,%xmm10
  cvtsi2sd %rsi,%xmm10
  mulsd  %xmm8,%xmm10
  mulsd  D224(%rip),%xmm10
  movapd %xmm1,%xmm0
  movsd  D21(%rip),%xmm4
  pxor   %xmm5,%xmm5
  movsd  D187(%rip),%xmm13
  movsd  D181(%rip),%xmm12
  movsd  D63(%rip),%xmm9
  movsd  D224(%rip),%xmm11
  movsd  D302(%rip),%xmm14
L296:
  addsd  %xmm13,%xmm0
  movapd %xmm2,%xmm6
  mulsd  %xmm12,%xmm6
  mulsd  %xmm3,%xmm6
  movapd %xmm0,%xmm7
  subsd  %xmm9,%xmm7
  mulsd  %xmm0,%xmm7
  divsd  %xmm7,%xmm6
  mulsd  %xmm6,%xmm4
  subsd  %xmm9,%xmm2
  subsd  %xmm9,%xmm3
  comisd %xmm4,%xmm11
  ja     L168
  addsd  %xmm4,%xmm5
  comisd %xmm10,%xmm5
  ja     L93
  comisd %xmm14,%xmm2
  ja     L296
  pxor   %xmm14,%xmm14
  jmp    L125
L168:
  movapd %xmm4,%xmm14
L125:
  divsd  %xmm8,%xmm5
  mov    $0x0,%eax
  movapd %xmm14,%xmm6
  addsd  D21(%rip),%xmm6
  comisd %xmm5,%xmm6
  jae    L62
  movsd  D181(%rip),%xmm7
  movapd %xmm2,%xmm8
  mulsd  %xmm7,%xmm8
  mulsd  %xmm3,%xmm8
  movsd  D63(%rip),%xmm6
  movapd %xmm0,%xmm10
  addsd  %xmm6,%xmm10
  movapd %xmm0,%xmm9
  addsd  D187(%rip),%xmm9
  mulsd  %xmm10,%xmm9
  divsd  %xmm9,%xmm8
  movapd %xmm6,%xmm9
  subsd  %xmm8,%xmm9
  movapd %xmm14,%xmm13
  divsd  %xmm9,%xmm13
  add    $0x1,%rcx
  pxor   %xmm10,%xmm10
  cvtsi2sd %rcx,%xmm10
  add    $0x1,%rdx
  pxor   %xmm11,%xmm11
  cvtsi2sd %rdx,%xmm11
  movapd %xmm1,%xmm8
  subsd  %xmm6,%xmm8
  mulsd  %xmm1,%xmm8
  mulsd  %xmm10,%xmm7
  mulsd  %xmm11,%xmm7
  divsd  %xmm7,%xmm8
  subsd  %xmm8,%xmm6
  movsd  D21(%rip),%xmm7
  divsd  %xmm6,%xmm7
  movapd %xmm7,%xmm6
  addsd  %xmm13,%xmm6
  mov    $0x1,%eax
  comisd %xmm6,%xmm5
  ja     L62
  movsd  D21(%rip),%xmm6
  mulsd  %xmm6,%xmm8
  addsd  %xmm8,%xmm6
  cmp    $0x3,%rdi
  jle    L36
  movapd %xmm5,%xmm15
  subsd  %xmm14,%xmm15
  movsd  D63(%rip),%xmm12
L202:
  subsd  D187(%rip),%xmm1
  addsd  %xmm12,%xmm10
  addsd  %xmm12,%xmm11
  movapd %xmm1,%xmm7
  subsd  %xmm12,%xmm7
  mulsd  %xmm1,%xmm7
  movsd  D181(%rip),%xmm9
  mulsd  %xmm10,%xmm9
  mulsd  %xmm11,%xmm9
  divsd  %xmm9,%xmm7
  mulsd  %xmm7,%xmm8
  movapd %xmm6,%xmm7
  addsd  %xmm8,%xmm6
  comisd %xmm15,%xmm6
  ja     L91
  comisd %xmm6,%xmm7
  jae    L36
  comisd D286(%rip),%xmm1
  ja     L202
L36:
  addsd  %xmm6,%xmm13
  mov    $0x1,%eax
  comisd %xmm13,%xmm5
  ja     L62
  subsd  %xmm6,%xmm5
  movsd  D302(%rip),%xmm8
  movsd  D187(%rip),%xmm10
  movsd  D181(%rip),%xmm9
  movsd  D63(%rip),%xmm7
L144:
  comisd %xmm8,%xmm2
  jbe    L328
  addsd  %xmm10,%xmm0
  movapd %xmm2,%xmm1
  mulsd  %xmm9,%xmm1
  mulsd  %xmm3,%xmm1
  movapd %xmm0,%xmm6
  subsd  %xmm7,%xmm6
  mulsd  %xmm0,%xmm6
  divsd  %xmm6,%xmm1
  mulsd  %xmm1,%xmm4
  subsd  %xmm7,%xmm2
  subsd  %xmm7,%xmm3
  movapd %xmm14,%xmm1
  addsd  %xmm4,%xmm1
  comisd %xmm5,%xmm1
  jae    L86
  comisd %xmm1,%xmm14
  jae    L169
  movapd %xmm1,%xmm14
  jmp    L144
L328:
  mov    $0x1,%eax
  ret
L342:
  mov    $0x1,%eax
  ret
L89:
  mov    $0x0,%eax
  ret
L133:
  mov    $0x1,%eax
  ret
L57:
  mov    $0x0,%eax
  ret
L118:
  mov    $0x1,%eax
  ret
L93:
  mov    $0x1,%eax
  ret
L91:
  mov    $0x0,%eax
  ret
L86:
  mov    $0x0,%eax
  ret
L169:
  mov    $0x1,%eax
L62:
  ret
```

tool_calls: name: parse_data
arguments: data_tag: D63
datatype: f64[2]
name: parse_data
arguments: data_tag: D224
datatype: f64
name: parse_data
arguments: data_tag: D21
datatype: byte[8]
name: parse_data
arguments: data_tag: D181
datatype: f32[2]
name: parse_data
arguments: data_tag: D187
datatype: f32[2]
name: parse_data
arguments: data_tag: D156
datatype: f64[2]
name: parse_data
arguments: data_tag: D302
datatype: f64
name: parse_data
arguments: data_tag: D286
datatype: f64

tool: {"D63": [1.0,1.0339757656912846e-25]}

tool: D224: 1.0339757656912846e-25

tool: {"D21": [0x00,0xfe,0xff,0xff,0xff,0xff,0xbf,0x3a]}

tool: {"D181": [0.0,2.25]}

tool: {"D187": [0.0,2.0]}

tool: {"D156": [1.5,0.5]}

tool: D302: 0.5

tool: D286: 3.5

assistant: ```c
int32_t SNPHWE_t(int32_t obs_hets, int32_t obs_hom1, int32_t obs_hom2, double thresh) {
  intptr_t obs_homc;
  intptr_t obs_homr;
  if (obs_hom1 < obs_hom2) {
    obs_homc = obs_hom2;
    obs_homr = obs_hom1;
  } else {
    obs_homc = obs_hom1;
    obs_homr = obs_hom2;
  }
  int64_t rare_copies = 2LL * obs_homr + obs_hets;
  int64_t genotypes2 = (obs_hets + obs_homc + obs_homr) * 2LL;
  double curr_hets_t2 = obs_hets;
  double curr_homr_t2 = obs_homr;
  double curr_homc_t2 = obs_homc;
  double tailp1 = (1 - 0.00000000000005684341886080801486968994140625) * 0.00000000000000000000000010339757656912845935892608650874535669572651386260986328125;
  double centerp = 0;
  double lastp2 = tailp1;
  double tailp2 = 0;
  double tail1_ceil;
  double tail2_ceil;
  double lastp1;
  double curr_hets_t1;
  double curr_homr_t1;
  double curr_homc_t1;
  double exit_thresh;
  double exit_threshx;
  double ratio;
  double preaddp;
  if (!genotypes2) {
    return 0;
  }
  thresh = (1 - thresh) / thresh;
  if (obs_hets * genotypes2 > rare_copies * (genotypes2 - rare_copies)) {
    if (obs_hets < 2) {
      return 0;
    }
    exit_thresh = rare_copies * thresh * 0.00000000000000000000000010339757656912845935892608650874535669572651386260986328125;
    do {
      curr_homr_t2 += 1;
      curr_homc_t2 += 1;
      lastp2 *= (curr_hets_t2 * (curr_hets_t2 - 1)) / (4 * curr_homr_t2 * curr_homc_t2);
      curr_hets_t2 -= 2;
      if (lastp2 < 0.00000000000000000000000010339757656912845935892608650874535669572651386260986328125) {
        tailp2 = lastp2;
        break;
      }
      centerp += lastp2;
      if (centerp > exit_thresh) {
        return 1;
      }
    } while (curr_hets_t2 > 1.5);
    exit_thresh = centerp / thresh;
    if (tailp1 + tailp2 >= exit_thresh) {
      return 0;
    }
    ratio = (curr_hets_t2 * (curr_hets_t2 - 1)) / (4 * (curr_homr_t2 + 1) * (curr_homc_t2 + 1));
    tail2_ceil = tailp2 / (1 - ratio);
    curr_hets_t1 = obs_hets + 2;
    curr_homr_t1 = obs_homr;
    curr_homc_t1 = obs_homc;
    lastp1 = (4 * curr_homr_t1 * curr_homc_t1) / (curr_hets_t1 * (curr_hets_t1 - 1));
    tail1_ceil = tailp1 / (1 - lastp1);
    if (tail1_ceil + tail2_ceil < exit_thresh) {
      return 1;
    }
    lastp1 *= tailp1;
    tailp1 += lastp1;
    if (obs_homr > 1) {
      exit_threshx = exit_thresh - tailp2;
      do {
        curr_hets_t1 += 2;
        curr_homr_t1 -= 1;
        curr_homc_t1 -= 1;
        lastp1 *= (4 * curr_homr_t1 * curr_homc_t1) / (curr_hets_t1 * (curr_hets_t1 - 1));
        preaddp = tailp1;
        tailp1 += lastp1;
        if (tailp1 > exit_threshx) {
          return 0;
        }
        if (tailp1 <= preaddp) {
          break;
        }
      } while (curr_homr_t1 > 1.5);
    }
    if (tailp1 + tail2_ceil < exit_thresh) {
      return 1;
    }
    exit_threshx = exit_thresh - tailp1;
    while (curr_hets_t2 > 1) {
      curr_homr_t2 += 1;
      curr_homc_t2 += 1;
      lastp2 *= (curr_hets_t2 * (curr_hets_t2 - 1)) / (4 * curr_homr_t2 * curr_homc_t2);
      preaddp = tailp2;
      tailp2 += lastp2;
      if (tailp2 >= exit_threshx) {
        return 0;
      }
      if (tailp2 <= preaddp) {
        return 1;
      }
      curr_hets_t2 -= 2;
    }
    return 1;
  } else {
    if (!obs_homr) {
      return 0;
    }
    exit_thresh = rare_copies * thresh * 0.00000000000000000000000010339757656912845935892608650874535669572651386260986328125;
    do {
      curr_hets_t2 += 2;
      lastp2 *= (4 * curr_homr_t2 * curr_homc_t2) / (curr_hets_t2 * (curr_hets_t2 - 1));
      curr_homr_t2 -= 1;
      curr_homc_t2 -= 1;
      if (lastp2 < 0.00000000000000000000000010339757656912845935892608650874535669572651386260986328125) {
        tailp2 = lastp2;
        break;
      }
      centerp += lastp2;
      if (centerp > exit_thresh) {
        return 1;
      }
    } while (curr_homr_t2 > 0.5);
    exit_thresh = centerp / thresh;
    if (tailp1 + tailp2 >= exit_thresh) {
      return 0;
    }
    ratio = (4 * curr_homr_t2 * curr_homc_t2) / ((curr_hets_t2 + 2) * (curr_hets_t2 + 1));
    tail2_ceil = tailp2 / (1 - ratio);
    curr_hets_t1 = obs_hets;
    curr_homr_t1 = obs_homr + 1;
    curr_homc_t1 = obs_homc + 1;
    lastp1 = (curr_hets_t1 * (curr_hets_t1 - 1)) / (4 * curr_homr_t1 * curr_homc_t1);
    tail1_ceil = tailp1 / (1 - lastp1);
    lastp1 *= tailp1;
    tailp1 += lastp1;
    if (tail1_ceil + tail2_ceil < exit_thresh) {
      return 1;
    }
    if (obs_hets >= 4) {
      exit_threshx = exit_thresh - tailp2;
      do {
        curr_hets_t1 -= 2;
        curr_homr_t1 += 1;
        curr_homc_t1 += 1;
        lastp1 *= (curr_hets_t1 * (curr_hets_t1 - 1)) / (4 * curr_homr_t1 * curr_homc_t1);
        preaddp = tailp1;
        tailp1 += lastp1;
        if (tailp1 > exit_threshx) {
          return 0;
        }
        if (tailp1 <= preaddp) {
          break;
        }
      } while (curr_hets_t1 > 3.5);
    }
    if (tailp1 + tail2_ceil < exit_thresh) {
      return 1;
    }
    exit_threshx = exit_thresh - tailp1;
    while (curr_homr_t2 > 0.5) {
      curr_hets_t2 += 2;
      lastp2 *= (4 * curr_homr_t2 * curr_homc_t2) / (curr_hets_t2 * (curr_hets_t2 - 1));
      curr_homr_t2 -= 1;
      curr_homc_t2 -= 1;
      preaddp = tailp2;
      tailp2 += lastp2;
      if (tailp2 >= exit_threshx) {
        return 0;
      }
      if (tailp2 <= preaddp) {
        return 1;
      }
    }
    return 1;
  }
}
```