#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ free(double *);
double *malloc(int);
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;
}


x86 assembly (clang, -O0):
<SNPHWE_pValue>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    114b <SNPHWE_pValue+0x2b>
mov    -0xc(%rbp),%eax
mov    %eax,-0x54(%rbp)
jmp    1151 <SNPHWE_pValue+0x31>
mov    -0x8(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    116e <SNPHWE_pValue+0x4e>
mov    -0x8(%rbp),%eax
mov    %eax,-0x58(%rbp)
jmp    1174 <SNPHWE_pValue+0x54>
mov    -0xc(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
shl    $0x1,%eax
add    -0x4(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x4(%rbp),%eax
add    -0x24(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x30(%rbp),%ecx
shl    $0x1,%ecx
sub    -0x2c(%rbp),%ecx
imul   %ecx,%eax
mov    -0x30(%rbp),%ecx
shl    $0x1,%ecx
cltd
idiv   %ecx
mov    %eax,-0x3c(%rbp)
mov    -0x2c(%rbp),%eax
and    $0x1,%eax
mov    -0x3c(%rbp),%ecx
and    $0x1,%ecx
xor    %ecx,%eax
cmp    $0x0,%eax
je     11e4 <SNPHWE_pValue+0xc4>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
sub    -0x3c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x44(%rbp)
mov    -0x30(%rbp),%eax
sub    -0x40(%rbp),%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  0xde9(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x1,-0x40(%rbp)
jle    12ed <SNPHWE_pValue+0x1cd>
mov    -0x38(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x40(%rbp),%xmm1
movsd  0xd9f(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x44(%rbp),%xmm2
movsd  0xd8a(%rip),%xmm1
addsd  %xmm1,%xmm2
movsd  0xd8e(%rip),%xmm1
mulsd  %xmm2,%xmm1
cvtsi2sdl -0x48(%rbp),%xmm2
movsd  0xd6d(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    -0x40(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x40(%rbp)
jmp    1234 <SNPHWE_pValue+0x114>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
sub    -0x3c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x44(%rbp)
mov    -0x30(%rbp),%eax
sub    -0x40(%rbp),%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jg     13c6 <SNPHWE_pValue+0x2a6>
mov    -0x38(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  0xcd9(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x44(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x48(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x40(%rbp),%xmm1
movsd  0xcad(%rip),%xmm2
addsd  %xmm2,%xmm1
cvtsi2sdl -0x40(%rbp),%xmm2
movsd  0xc94(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x44(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x44(%rbp)
mov    -0x48(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x40(%rbp)
jmp    1316 <SNPHWE_pValue+0x1f6>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1402 <SNPHWE_pValue+0x2e2>
movsd  -0x50(%rbp),%xmm1
mov    -0x38(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    13cd <SNPHWE_pValue+0x2ad>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
cmpl   $0x1,-0x10(%rbp)
jne    1476 <SNPHWE_pValue+0x356>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1471 <SNPHWE_pValue+0x351>
mov    -0x38(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x4(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    144c <SNPHWE_pValue+0x32c>
jmp    1463 <SNPHWE_pValue+0x343>
mov    -0x38(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    141b <SNPHWE_pValue+0x2fb>
jmp    14c1 <SNPHWE_pValue+0x3a1>
cmpl   $0x2,-0x10(%rbp)
jne    14bc <SNPHWE_pValue+0x39c>
mov    -0x4(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     14b7 <SNPHWE_pValue+0x397>
mov    -0x38(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1486 <SNPHWE_pValue+0x366>
jmp    14bc <SNPHWE_pValue+0x39c>
jmp    14c1 <SNPHWE_pValue+0x3a1>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
movsd  -0x20(%rbp),%xmm0
movsd  0xb29(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    14f3 <SNPHWE_pValue+0x3d3>
movsd  0xb17(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    14fd <SNPHWE_pValue+0x3dd>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<SNPHWE_pValue>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %ecx,%r15d
mov    %edx,%r12d
mov    %edi,%r14d
cmp    %edx,%esi
mov    %esi,%ebx
cmovl  %edx,%ebx
cmovl  %esi,%r12d
lea    (%r14,%r12,2),%r13d
add    %edi,%ebx
add    %r12d,%ebx
lea    0x8(,%r13,8),%edi
call   1040 <malloc@plt>
mov    %rax,%r9
lea    (%rbx,%rbx,1),%ecx
mov    %ecx,%eax
sub    %r13d,%eax
imul   %r13d,%eax
cltd
idiv   %ecx
mov    %eax,%ecx
xor    %r13d,%ecx
and    $0x1,%ecx
lea    (%rcx,%rax,1),%r8d
mov    %r13d,%ebp
sub    %r8d,%ebp
mov    %ebp,%edx
shr    $0x1f,%edx
add    %ebp,%edx
sar    %edx
sub    %r8d,%ebx
sub    %edx,%ebx
movslq %r8d,%rbp
movabs $0x3ff0000000000000,%rsi
mov    %rsi,(%r9,%rbp,8)
cmp    $0x2,%ebp
jl     123b <SNPHWE_pValue+0x11b>
mov    %r8d,%esi
movsd  (%r9,%rsi,8),%xmm5
movd   %ebx,%xmm0
movd   %r8d,%xmm1
punpckldq %xmm0,%xmm1
movsd  0xe3e(%rip),%xmm2
movapd 0xe56(%rip),%xmm3
movdqa 0xe5e(%rip),%xmm4
movapd %xmm2,%xmm0
mov    %edx,%ebp
nopl   0x0(%rax,%rax,1)
xorps  %xmm6,%xmm6
cvtsi2sd %ebp,%xmm6
addsd  %xmm2,%xmm6
unpcklpd %xmm6,%xmm5
cvtdq2pd %xmm1,%xmm6
movapd %xmm6,%xmm7
movhpd 0xe09(%rip),%xmm7
mulpd  %xmm5,%xmm7
addpd  %xmm3,%xmm6
mulpd  %xmm7,%xmm6
movapd %xmm6,%xmm5
unpckhpd %xmm6,%xmm5
divsd  %xmm5,%xmm6
movd   %xmm1,%esi
paddd  %xmm4,%xmm1
movd   %xmm1,%edi
movsd  %xmm6,(%r9,%rdi,8)
addsd  %xmm6,%xmm0
add    $0x1,%ebp
movapd %xmm6,%xmm5
cmp    $0x3,%esi
ja     11e0 <SNPHWE_pValue+0xc0>
jmp    1243 <SNPHWE_pValue+0x123>
movsd  0xdbd(%rip),%xmm0
lea    -0x2(%r13),%ebp
cmp    %ebp,%r8d
jg     12cf <SNPHWE_pValue+0x1af>
cltq
add    %rax,%rcx
movslq %ebp,%rax
movsd  (%r9,%rcx,8),%xmm2
movsd  0xda2(%rip),%xmm1
movsd  0xda2(%rip),%xmm3
movsd  0xd8a(%rip),%xmm4
cs nopw 0x0(%rax,%rax,1)
mulsd  %xmm1,%xmm2
xorps  %xmm5,%xmm5
cvtsi2sd %edx,%xmm5
mulsd  %xmm2,%xmm5
xorps  %xmm2,%xmm2
cvtsi2sd %ebx,%xmm2
mulsd  %xmm5,%xmm2
xorps  %xmm5,%xmm5
cvtsi2sd %ecx,%xmm5
movapd %xmm5,%xmm6
addsd  %xmm3,%xmm6
addsd  %xmm4,%xmm5
mulsd  %xmm6,%xmm5
divsd  %xmm5,%xmm2
movsd  %xmm2,0x10(%r9,%rcx,8)
add    $0x2,%rcx
addsd  %xmm2,%xmm0
add    $0xffffffff,%edx
add    $0xffffffff,%ebx
cmp    %rax,%rcx
jle    1280 <SNPHWE_pValue+0x160>
test   %r13d,%r13d
js     1379 <SNPHWE_pValue+0x259>
lea    0x1(%r13),%eax
test   %r13d,%r13d
je     1359 <SNPHWE_pValue+0x239>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rcx),%rdx
mov    %rdx,%rbp
shr    %rbp
add    $0x1,%rbp
test   %rdx,%rdx
je     1453 <SNPHWE_pValue+0x333>
mov    %rbp,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %esi,%esi
xchg   %ax,%ax
movupd (%r9,%rsi,8),%xmm2
movupd 0x10(%r9,%rsi,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,(%r9,%rsi,8)
divpd  %xmm1,%xmm3
movupd %xmm3,0x10(%r9,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdx
jne    1310 <SNPHWE_pValue+0x1f0>
test   $0x1,%bpl
je     1352 <SNPHWE_pValue+0x232>
movupd (%r9,%rsi,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%r9,%rsi,8)
cmp    %rax,%rcx
jne    1360 <SNPHWE_pValue+0x240>
jmp    1379 <SNPHWE_pValue+0x259>
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  (%r9,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r9,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1360 <SNPHWE_pValue+0x240>
xorpd  %xmm6,%xmm6
cmp    $0x1,%r15d
je     1424 <SNPHWE_pValue+0x304>
cmp    $0x2,%r15d
jne    155b <SNPHWE_pValue+0x43b>
test   %r12d,%r12d
js     155b <SNPHWE_pValue+0x43b>
movslq %r14d,%rax
cmp    %r13d,%r14d
cmovg  %r14d,%r13d
lea    0x1(%r13),%ecx
mov    %ecx,%edx
sub    %r14d,%edx
sub    %r14d,%r13d
xorpd  %xmm6,%xmm6
and    $0x7,%edx
je     13cf <SNPHWE_pValue+0x2af>
nopl   0x0(%rax)
addsd  (%r9,%rax,8),%xmm6
add    $0x1,%rax
add    $0xffffffff,%edx
jne    13c0 <SNPHWE_pValue+0x2a0>
cmp    $0x7,%r13d
jb     155b <SNPHWE_pValue+0x43b>
mov    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
addsd  (%r9,%rax,8),%xmm6
addsd  0x8(%r9,%rax,8),%xmm6
addsd  0x10(%r9,%rax,8),%xmm6
addsd  0x18(%r9,%rax,8),%xmm6
addsd  0x20(%r9,%rax,8),%xmm6
addsd  0x28(%r9,%rax,8),%xmm6
addsd  0x30(%r9,%rax,8),%xmm6
addsd  0x38(%r9,%rax,8),%xmm6
add    $0x8,%rax
cmp    %eax,%ecx
jne    13e0 <SNPHWE_pValue+0x2c0>
jmp    155b <SNPHWE_pValue+0x43b>
test   %r13d,%r13d
js     155b <SNPHWE_pValue+0x43b>
movslq %r14d,%rax
movsd  (%r9,%rax,8),%xmm0
mov    %r13d,%ecx
add    $0x1,%rcx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0x3,%r13d
jae    1464 <SNPHWE_pValue+0x344>
xorpd  %xmm6,%xmm6
xor    %edx,%edx
jmp    151c <SNPHWE_pValue+0x3fc>
xor    %esi,%esi
test   $0x1,%bpl
jne    1342 <SNPHWE_pValue+0x222>
jmp    1352 <SNPHWE_pValue+0x232>
and    $0xfffffffc,%ecx
xorpd  %xmm6,%xmm6
xor    %edx,%edx
movsd  0xba3(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r9,%rdx,8),%xmm2
movsd  0x8(%r9,%rdx,8),%xmm3
movapd %xmm0,%xmm4
cmpltsd %xmm2,%xmm4
movapd %xmm4,%xmm5
andnpd %xmm2,%xmm5
andpd  %xmm1,%xmm4
orpd   %xmm5,%xmm4
addsd  %xmm6,%xmm4
movapd %xmm0,%xmm2
cmpltsd %xmm3,%xmm2
movapd %xmm2,%xmm5
andnpd %xmm3,%xmm5
andpd  %xmm1,%xmm2
orpd   %xmm5,%xmm2
addsd  %xmm4,%xmm2
movsd  0x10(%r9,%rdx,8),%xmm3
movapd %xmm0,%xmm4
cmpltsd %xmm3,%xmm4
movapd %xmm4,%xmm5
andnpd %xmm3,%xmm5
andpd  %xmm1,%xmm4
orpd   %xmm5,%xmm4
addsd  %xmm2,%xmm4
movsd  0x18(%r9,%rdx,8),%xmm2
movapd %xmm0,%xmm6
cmpltsd %xmm2,%xmm6
movapd %xmm6,%xmm3
andnpd %xmm2,%xmm3
andpd  %xmm1,%xmm6
orpd   %xmm3,%xmm6
addsd  %xmm4,%xmm6
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1480 <SNPHWE_pValue+0x360>
test   %rax,%rax
je     155b <SNPHWE_pValue+0x43b>
lea    (%r9,%rdx,8),%rcx
xor    %edx,%edx
movsd  0xae9(%rip),%xmm1
nop
movsd  (%rcx,%rdx,8),%xmm2
movapd %xmm0,%xmm3
cmpltsd %xmm2,%xmm3
movapd %xmm3,%xmm4
andnpd %xmm2,%xmm4
andpd  %xmm1,%xmm3
orpd   %xmm4,%xmm3
addsd  %xmm3,%xmm6
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1530 <SNPHWE_pValue+0x410>
movapd %xmm6,(%rsp)
mov    %r9,%rdi
call   1030 <free@plt>
movsd  0xa90(%rip),%xmm0
minsd  (%rsp),%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @SNPHWE_pValue(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double*
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca double
store i32 %0, i32* %5
store i32 %1, i32* %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%21 = load i32, i32* %6
%22 = load i32, i32* %7
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %26
24: ; preds = %4
%25 = load i32, i32* %7
br label %28
26: ; preds = %4
%27 = load i32, i32* %6
br label %28
28: ; preds = %26, %24
%29 = phi i32 [ %25, %24 ], [ %27, %26 ]
store i32 %29, i32* %11
%30 = load i32, i32* %6
%31 = load i32, i32* %7
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %35
33: ; preds = %28
%34 = load i32, i32* %6
br label %37
35: ; preds = %28
%36 = load i32, i32* %7
br label %37
37: ; preds = %35, %33
%38 = phi i32 [ %34, %33 ], [ %36, %35 ]
store i32 %38, i32* %12
%39 = load i32, i32* %12
%40 = mul nsw i32 2, %39
%41 = load i32, i32* %5
%42 = add nsw i32 %40, %41
store i32 %42, i32* %13
%43 = load i32, i32* %5
%44 = load i32, i32* %11
%45 = add nsw i32 %43, %44
%46 = load i32, i32* %12
%47 = add nsw i32 %45, %46
store i32 %47, i32* %14
%48 = load i32, i32* %13
%49 = add nsw i32 %48, 1
%50 = sext i32 %49 to i64
%51 = mul i64 %50, 8
%52 = trunc i64 %51 to i32
%53 = call double* @malloc(i32 %52)
store double* %53, double** %15
%54 = load i32, i32* %13
%55 = load i32, i32* %14
%56 = mul nsw i32 2, %55
%57 = load i32, i32* %13
%58 = sub nsw i32 %56, %57
%59 = mul nsw i32 %54, %58
%60 = load i32, i32* %14
%61 = mul nsw i32 2, %60
%62 = sdiv i32 %59, %61
store i32 %62, i32* %16
%63 = load i32, i32* %13
%64 = and i32 %63, 1
%65 = load i32, i32* %16
%66 = and i32 %65, 1
%67 = xor i32 %64, %66
%68 = icmp ne i32 %67, 0
br i1 %68, label %69, label %72
69: ; preds = %37
%70 = load i32, i32* %16
%71 = add nsw i32 %70, 1
store i32 %71, i32* %16
br label %72
72: ; preds = %69, %37
%73 = load i32, i32* %16
store i32 %73, i32* %17
%74 = load i32, i32* %13
%75 = load i32, i32* %16
%76 = sub nsw i32 %74, %75
%77 = sdiv i32 %76, 2
store i32 %77, i32* %18
%78 = load i32, i32* %14
%79 = load i32, i32* %17
%80 = sub nsw i32 %78, %79
%81 = load i32, i32* %18
%82 = sub nsw i32 %80, %81
store i32 %82, i32* %19
%83 = load double*, double** %15
%84 = load i32, i32* %16
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
store double 1.000000e+00, double* %86
%87 = load double*, double** %15
%88 = load i32, i32* %16
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
store double %91, double* %20
%92 = load i32, i32* %16
store i32 %92, i32* %17
br label %93
93: ; preds = %135, %72
%94 = load i32, i32* %17
%95 = icmp sgt i32 %94, 1
br i1 %95, label %96, label %138
96: ; preds = %93
%97 = load double*, double** %15
%98 = load i32, i32* %17
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %97, i64 %99
%101 = load double, double* %100
%102 = load i32, i32* %17
%103 = sitofp i32 %102 to double
%104 = fmul double %101, %103
%105 = load i32, i32* %17
%106 = sitofp i32 %105 to double
%107 = fsub double %106, 1.000000e+00
%108 = fmul double %104, %107
%109 = load i32, i32* %18
%110 = sitofp i32 %109 to double
%111 = fadd double %110, 1.000000e+00
%112 = fmul double 4.000000e+00, %111
%113 = load i32, i32* %19
%114 = sitofp i32 %113 to double
%115 = fadd double %114, 1.000000e+00
%116 = fmul double %112, %115
%117 = fdiv double %108, %116
%118 = load double*, double** %15
%119 = load i32, i32* %17
%120 = sub nsw i32 %119, 2
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %118, i64 %121
store double %117, double* %122
%123 = load double*, double** %15
%124 = load i32, i32* %17
%125 = sub nsw i32 %124, 2
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds double, double* %123, i64 %126
%128 = load double, double* %127
%129 = load double, double* %20
%130 = fadd double %129, %128
store double %130, double* %20
%131 = load i32, i32* %18
%132 = add nsw i32 %131, 1
store i32 %132, i32* %18
%133 = load i32, i32* %19
%134 = add nsw i32 %133, 1
store i32 %134, i32* %19
br label %135
135: ; preds = %96
%136 = load i32, i32* %17
%137 = sub nsw i32 %136, 2
store i32 %137, i32* %17
br label %93
138: ; preds = %93
%139 = load i32, i32* %16
store i32 %139, i32* %17
%140 = load i32, i32* %13
%141 = load i32, i32* %16
%142 = sub nsw i32 %140, %141
%143 = sdiv i32 %142, 2
store i32 %143, i32* %18
%144 = load i32, i32* %14
%145 = load i32, i32* %17
%146 = sub nsw i32 %144, %145
%147 = load i32, i32* %18
%148 = sub nsw i32 %146, %147
store i32 %148, i32* %19
%149 = load i32, i32* %16
store i32 %149, i32* %17
br label %150
150: ; preds = %193, %138
%151 = load i32, i32* %17
%152 = load i32, i32* %13
%153 = sub nsw i32 %152, 2
%154 = icmp sle i32 %151, %153
br i1 %154, label %155, label %196
155: ; preds = %150
%156 = load double*, double** %15
%157 = load i32, i32* %17
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %156, i64 %158
%160 = load double, double* %159
%161 = fmul double %160, 4.000000e+00
%162 = load i32, i32* %18
%163 = sitofp i32 %162 to double
%164 = fmul double %161, %163
%165 = load i32, i32* %19
%166 = sitofp i32 %165 to double
%167 = fmul double %164, %166
%168 = load i32, i32* %17
%169 = sitofp i32 %168 to double
%170 = fadd double %169, 2.000000e+00
%171 = load i32, i32* %17
%172 = sitofp i32 %171 to double
%173 = fadd double %172, 1.000000e+00
%174 = fmul double %170, %173
%175 = fdiv double %167, %174
%176 = load double*, double** %15
%177 = load i32, i32* %17
%178 = add nsw i32 %177, 2
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds double, double* %176, i64 %179
store double %175, double* %180
%181 = load double*, double** %15
%182 = load i32, i32* %17
%183 = add nsw i32 %182, 2
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds double, double* %181, i64 %184
%186 = load double, double* %185
%187 = load double, double* %20
%188 = fadd double %187, %186
store double %188, double* %20
%189 = load i32, i32* %18
%190 = add nsw i32 %189, -1
store i32 %190, i32* %18
%191 = load i32, i32* %19
%192 = add nsw i32 %191, -1
store i32 %192, i32* %19
br label %193
193: ; preds = %155
%194 = load i32, i32* %17
%195 = add nsw i32 %194, 2
store i32 %195, i32* %17
br label %150
196: ; preds = %150
store i32 0, i32* %9
br label %197
197: ; preds = %209, %196
%198 = load i32, i32* %9
%199 = load i32, i32* %13
%200 = icmp sle i32 %198, %199
br i1 %200, label %201, label %212
201: ; preds = %197
%202 = load double, double* %20
%203 = load double*, double** %15
%204 = load i32, i32* %9
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds double, double* %203, i64 %205
%207 = load double, double* %206
%208 = fdiv double %207, %202
store double %208, double* %206
br label %209
209: ; preds = %201
%210 = load i32, i32* %9
%211 = add nsw i32 %210, 1
store i32 %211, i32* %9
br label %197
212: ; preds = %197
store double 0.000000e+00, double* %10
%213 = load i32, i32* %8
%214 = icmp eq i32 %213, 1
br i1 %214, label %215, label %245
215: ; preds = %212
store i32 0, i32* %9
br label %216
216: ; preds = %241, %215
%217 = load i32, i32* %9
%218 = load i32, i32* %13
%219 = icmp sle i32 %217, %218
br i1 %219, label %220, label %244
220: ; preds = %216
%221 = load double*, double** %15
%222 = load i32, i32* %9
%223 = sext i32 %222 to i64
%224 = getelementptr inbounds double, double* %221, i64 %223
%225 = load double, double* %224
%226 = load double*, double** %15
%227 = load i32, i32* %5
%228 = sext i32 %227 to i64
%229 = getelementptr inbounds double, double* %226, i64 %228
%230 = load double, double* %229
%231 = fcmp ogt double %225, %230
br i1 %231, label %232, label %233
232: ; preds = %220
br label %241
233: ; preds = %220
%234 = load double*, double** %15
%235 = load i32, i32* %9
%236 = sext i32 %235 to i64
%237 = getelementptr inbounds double, double* %234, i64 %236
%238 = load double, double* %237
%239 = load double, double* %10
%240 = fadd double %239, %238
store double %240, double* %10
br label %241
241: ; preds = %233, %232
%242 = load i32, i32* %9
%243 = add nsw i32 %242, 1
store i32 %243, i32* %9
br label %216
244: ; preds = %216
br label %267
245: ; preds = %212
%246 = load i32, i32* %8
%247 = icmp eq i32 %246, 2
br i1 %247, label %248, label %266
248: ; preds = %245
%249 = load i32, i32* %5
store i32 %249, i32* %9
br label %250
250: ; preds = %262, %248
%251 = load i32, i32* %9
%252 = load i32, i32* %13
%253 = icmp sle i32 %251, %252
br i1 %253, label %254, label %265
254: ; preds = %250
%255 = load double*, double** %15
%256 = load i32, i32* %9
%257 = sext i32 %256 to i64
%258 = getelementptr inbounds double, double* %255, i64 %257
%259 = load double, double* %258
%260 = load double, double* %10
%261 = fadd double %260, %259
store double %261, double* %10
br label %262
262: ; preds = %254
%263 = load i32, i32* %9
%264 = add nsw i32 %263, 1
store i32 %264, i32* %9
br label %250
265: ; preds = %250
br label %266
266: ; preds = %265, %245
br label %267
267: ; preds = %266, %244
%268 = load double*, double** %15
%269 = call i32 @free(double* %268)
%270 = load double, double* %10
%271 = fcmp ogt double %270, 1.000000e+00
br i1 %271, label %272, label %273
272: ; preds = %267
br label %275
273: ; preds = %267
%274 = load double, double* %10
br label %275
275: ; preds = %273, %272
%276 = phi double [ 1.000000e+00, %272 ], [ %274, %273 ]
ret double %276
}

LLVM IR (-O2):
define double @SNPHWE_pValue(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = icmp slt i32 %1, %2
%6 = select i1 %5, i32 %2, i32 %1
%7 = select i1 %5, i32 %1, i32 %2
%8 = shl nsw i32 %7, 1
%9 = add i32 %8, %0
%10 = add nsw i32 %6, %0
%11 = add nsw i32 %10, %7
%12 = shl i32 %9, 3
%13 = add i32 %12, 8
%14 = tail call double* @malloc(i32 %13)
%15 = shl nsw i32 %11, 1
%16 = sub nsw i32 %15, %9
%17 = mul nsw i32 %16, %9
%18 = sdiv i32 %17, %15
%19 = xor i32 %18, %9
%20 = and i32 %19, 1
%21 = add nsw i32 %20, %18
%22 = sub nsw i32 %9, %21
%23 = sdiv i32 %22, 2
%24 = add i32 %21, %23
%25 = sub i32 %11, %24
%26 = sext i32 %21 to i64
%27 = getelementptr inbounds double, double* %14, i64 %26
store double 1.000000e+00, double* %27
%28 = icmp sgt i32 %21, 1
br i1 %28, label %29, label %35
29: ; preds = %4
%30 = zext i32 %21 to i64
%31 = getelementptr inbounds double, double* %14, i64 %30
%32 = load double, double* %31
%33 = insertelement <2 x i32> poison, i32 %21, i64 0
%34 = insertelement <2 x i32> %33, i32 %25, i64 1
br label %47
35: ; preds = %47, %4
%36 = phi double [ 1.000000e+00, %4 ], [ %68, %47 ]
%37 = add nsw i32 %9, -2
%38 = icmp sgt i32 %21, %37
br i1 %38, label %72, label %39
39: ; preds = %35
%40 = sext i32 %18 to i64
%41 = and i32 %19, 1
%42 = zext i32 %41 to i64
%43 = add nsw i64 %40, %42
%44 = sext i32 %37 to i64
%45 = getelementptr inbounds double, double* %14, i64 %43
%46 = load double, double* %45
br label %120
47: ; preds = %29, %47
%48 = phi double [ %63, %47 ], [ %32, %29 ]
%49 = phi double [ %68, %47 ], [ 1.000000e+00, %29 ]
%50 = phi i32 [ %69, %47 ], [ %23, %29 ]
%51 = phi <2 x i32> [ %64, %47 ], [ %34, %29 ]
%52 = sitofp i32 %50 to double
%53 = fadd double %52, 1.000000e+00
%54 = sitofp <2 x i32> %51 to <2 x double>
%55 = insertelement <2 x double> poison, double %48, i64 0
%56 = insertelement <2 x double> %55, double %53, i64 1
%57 = insertelement <2 x double> %54, double 4.000000e+00, i64 1
%58 = fmul <2 x double> %56, %57
%59 = fadd <2 x double> %54, <double -1.000000e+00, double 1.000000e+00>
%60 = fmul <2 x double> %59, %58
%61 = extractelement <2 x double> %60, i64 0
%62 = extractelement <2 x double> %60, i64 1
%63 = fdiv double %61, %62
%64 = add nsw <2 x i32> %51, <i32 -2, i32 1>
%65 = extractelement <2 x i32> %64, i64 0
%66 = zext i32 %65 to i64
%67 = getelementptr inbounds double, double* %14, i64 %66
store double %63, double* %67
%68 = fadd double %49, %63
%69 = add nsw i32 %50, 1
%70 = extractelement <2 x i32> %51, i64 0
%71 = icmp ugt i32 %70, 3
br i1 %71, label %47, label %35
72: ; preds = %120, %35
%73 = phi double [ %36, %35 ], [ %139, %120 ]
%74 = icmp slt i32 %9, 0
br i1 %74, label %150, label %75
75: ; preds = %72
%76 = add nuw i32 %9, 1
%77 = zext i32 %76 to i64
%78 = icmp eq i32 %9, 0
br i1 %78, label %118, label %79
79: ; preds = %75
%80 = and i64 %77, 4294967294
%81 = insertelement <2 x double> poison, double %73, i64 0
%82 = shufflevector <2 x double> %81, <2 x double> poison, <2 x i32> zeroinitializer
%83 = add nsw i64 %80, -2
%84 = lshr exact i64 %83, 1
%85 = add nuw i64 %84, 1
%86 = and i64 %85, 1
%87 = icmp eq i64 %83, 0
br i1 %87, label %107, label %88
88: ; preds = %79
%89 = and i64 %85, -2
br label %90
90: ; preds = %90, %88
%91 = phi i64 [ 0, %88 ], [ %104, %90 ]
%92 = phi i64 [ 0, %88 ], [ %105, %90 ]
%93 = getelementptr inbounds double, double* %14, i64 %91
%94 = bitcast double* %93 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = fdiv <2 x double> %95, %82
%97 = bitcast double* %93 to <2 x double>*
store <2 x double> %96, <2 x double>* %97
%98 = or i64 %91, 2
%99 = getelementptr inbounds double, double* %14, i64 %98
%100 = bitcast double* %99 to <2 x double>*
%101 = load <2 x double>, <2 x double>* %100
%102 = fdiv <2 x double> %101, %82
%103 = bitcast double* %99 to <2 x double>*
store <2 x double> %102, <2 x double>* %103
%104 = add nuw i64 %91, 4
%105 = add nuw i64 %92, 2
%106 = icmp eq i64 %105, %89
br i1 %106, label %107, label %90
107: ; preds = %90, %79
%108 = phi i64 [ 0, %79 ], [ %104, %90 ]
%109 = icmp eq i64 %86, 0
br i1 %109, label %116, label %110
110: ; preds = %107
%111 = getelementptr inbounds double, double* %14, i64 %108
%112 = bitcast double* %111 to <2 x double>*
%113 = load <2 x double>, <2 x double>* %112
%114 = fdiv <2 x double> %113, %82
%115 = bitcast double* %111 to <2 x double>*
store <2 x double> %114, <2 x double>* %115
br label %116
116: ; preds = %107, %110
%117 = icmp eq i64 %80, %77
br i1 %117, label %150, label %118
118: ; preds = %75, %116
%119 = phi i64 [ 0, %75 ], [ %80, %116 ]
br label %143
120: ; preds = %39, %120
%121 = phi double [ %46, %39 ], [ %136, %120 ]
%122 = phi i64 [ %43, %39 ], [ %137, %120 ]
%123 = phi double [ %36, %39 ], [ %139, %120 ]
%124 = phi i32 [ %25, %39 ], [ %141, %120 ]
%125 = phi i32 [ %23, %39 ], [ %140, %120 ]
%126 = fmul double %121, 4.000000e+00
%127 = sitofp i32 %125 to double
%128 = fmul double %126, %127
%129 = sitofp i32 %124 to double
%130 = fmul double %128, %129
%131 = trunc i64 %122 to i32
%132 = sitofp i32 %131 to double
%133 = fadd double %132, 2.000000e+00
%134 = fadd double %132, 1.000000e+00
%135 = fmul double %133, %134
%136 = fdiv double %130, %135
%137 = add nsw i64 %122, 2
%138 = getelementptr inbounds double, double* %14, i64 %137
store double %136, double* %138
%139 = fadd double %123, %136
%140 = add nsw i32 %125, -1
%141 = add nsw i32 %124, -1
%142 = icmp sgt i64 %137, %44
br i1 %142, label %72, label %120
143: ; preds = %118, %143
%144 = phi i64 [ %148, %143 ], [ %119, %118 ]
%145 = getelementptr inbounds double, double* %14, i64 %144
%146 = load double, double* %145
%147 = fdiv double %146, %73
store double %147, double* %145
%148 = add nuw nsw i64 %144, 1
%149 = icmp eq i64 %148, %77
br i1 %149, label %150, label %143
150: ; preds = %143, %116, %72
switch i32 %3, label %272 [
i32 1, label %176
i32 2, label %151
]
151: ; preds = %150
%152 = icmp slt i32 %7, 0
br i1 %152, label %272, label %153
153: ; preds = %151
%154 = sext i32 %0 to i64
%155 = call i32 @llvm.smax.i32(i32 %0, i32 %9)
%156 = add i32 %155, 1
%157 = sub i32 %156, %0
%158 = sub i32 %155, %0
%159 = and i32 %157, 7
%160 = icmp eq i32 %159, 0
br i1 %160, label %171, label %161
161: ; preds = %153, %161
%162 = phi i64 [ %168, %161 ], [ %154, %153 ]
%163 = phi double [ %167, %161 ], [ 0.000000e+00, %153 ]
%164 = phi i32 [ %169, %161 ], [ 0, %153 ]
%165 = getelementptr inbounds double, double* %14, i64 %162
%166 = load double, double* %165
%167 = fadd double %163, %166
%168 = add nsw i64 %162, 1
%169 = add i32 %164, 1
%170 = icmp eq i32 %169, %159
br i1 %170, label %171, label %161
171: ; preds = %161, %153
%172 = phi double [ undef, %153 ], [ %167, %161 ]
%173 = phi i64 [ %154, %153 ], [ %168, %161 ]
%174 = phi double [ 0.000000e+00, %153 ], [ %167, %161 ]
%175 = icmp ult i32 %158, 7
br i1 %175, label %272, label %218
176: ; preds = %150
br i1 %74, label %272, label %177
177: ; preds = %176
%178 = sext i32 %0 to i64
%179 = getelementptr inbounds double, double* %14, i64 %178
%180 = load double, double* %179
%181 = add nuw i32 %9, 1
%182 = zext i32 %181 to i64
%183 = add nsw i64 %182, -1
%184 = and i64 %182, 3
%185 = icmp ult i64 %183, 3
br i1 %185, label %255, label %186
186: ; preds = %177
%187 = and i64 %182, 4294967292
br label %188
188: ; preds = %188, %186
%189 = phi i64 [ 0, %186 ], [ %215, %188 ]
%190 = phi double [ 0.000000e+00, %186 ], [ %214, %188 ]
%191 = phi i64 [ 0, %186 ], [ %216, %188 ]
%192 = getelementptr inbounds double, double* %14, i64 %189
%193 = load double, double* %192
%194 = fcmp ogt double %193, %180
%195 = select i1 %194, double -0.000000e+00, double %193
%196 = fadd double %190, %195
%197 = or i64 %189, 1
%198 = getelementptr inbounds double, double* %14, i64 %197
%199 = load double, double* %198
%200 = fcmp ogt double %199, %180
%201 = select i1 %200, double -0.000000e+00, double %199
%202 = fadd double %196, %201
%203 = or i64 %189, 2
%204 = getelementptr inbounds double, double* %14, i64 %203
%205 = load double, double* %204
%206 = fcmp ogt double %205, %180
%207 = select i1 %206, double -0.000000e+00, double %205
%208 = fadd double %202, %207
%209 = or i64 %189, 3
%210 = getelementptr inbounds double, double* %14, i64 %209
%211 = load double, double* %210
%212 = fcmp ogt double %211, %180
%213 = select i1 %212, double -0.000000e+00, double %211
%214 = fadd double %208, %213
%215 = add nuw nsw i64 %189, 4
%216 = add i64 %191, 4
%217 = icmp eq i64 %216, %187
br i1 %217, label %255, label %188
218: ; preds = %171, %218
%219 = phi i64 [ %252, %218 ], [ %173, %171 ]
%220 = phi double [ %251, %218 ], [ %174, %171 ]
%221 = getelementptr inbounds double, double* %14, i64 %219
%222 = load double, double* %221
%223 = fadd double %220, %222
%224 = add nsw i64 %219, 1
%225 = getelementptr inbounds double, double* %14, i64 %224
%226 = load double, double* %225
%227 = fadd double %223, %226
%228 = add nsw i64 %219, 2
%229 = getelementptr inbounds double, double* %14, i64 %228
%230 = load double, double* %229
%231 = fadd double %227, %230
%232 = add nsw i64 %219, 3
%233 = getelementptr inbounds double, double* %14, i64 %232
%234 = load double, double* %233
%235 = fadd double %231, %234
%236 = add nsw i64 %219, 4
%237 = getelementptr inbounds double, double* %14, i64 %236
%238 = load double, double* %237
%239 = fadd double %235, %238
%240 = add nsw i64 %219, 5
%241 = getelementptr inbounds double, double* %14, i64 %240
%242 = load double, double* %241
%243 = fadd double %239, %242
%244 = add nsw i64 %219, 6
%245 = getelementptr inbounds double, double* %14, i64 %244
%246 = load double, double* %245
%247 = fadd double %243, %246
%248 = add nsw i64 %219, 7
%249 = getelementptr inbounds double, double* %14, i64 %248
%250 = load double, double* %249
%251 = fadd double %247, %250
%252 = add nsw i64 %219, 8
%253 = trunc i64 %252 to i32
%254 = icmp eq i32 %156, %253
br i1 %254, label %272, label %218
255: ; preds = %188, %177
%256 = phi double [ undef, %177 ], [ %214, %188 ]
%257 = phi i64 [ 0, %177 ], [ %215, %188 ]
%258 = phi double [ 0.000000e+00, %177 ], [ %214, %188 ]
%259 = icmp eq i64 %184, 0
br i1 %259, label %272, label %260
260: ; preds = %255, %260
%261 = phi i64 [ %269, %260 ], [ %257, %255 ]
%262 = phi double [ %268, %260 ], [ %258, %255 ]
%263 = phi i64 [ %270, %260 ], [ 0, %255 ]
%264 = getelementptr inbounds double, double* %14, i64 %261
%265 = load double, double* %264
%266 = fcmp ogt double %265, %180
%267 = select i1 %266, double -0.000000e+00, double %265
%268 = fadd double %262, %267
%269 = add nuw nsw i64 %261, 1
%270 = add i64 %263, 1
%271 = icmp eq i64 %270, %184
br i1 %271, label %272, label %260
272: ; preds = %171, %218, %255, %260, %151, %176, %150
%273 = phi double [ 0.000000e+00, %150 ], [ 0.000000e+00, %176 ], [ 0.000000e+00, %151 ], [ %256, %255 ], [ %268, %260 ], [ %172, %171 ], [ %251, %218 ]
%274 = tail call i32 @free(double* %14)
%275 = fcmp ogt double %273, 1.000000e+00
%276 = select i1 %275, double 1.000000e+00, double %273
ret double %276
}