#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 *);
scalar_t__ malloc(int);
void ancestralGamma(double *x3, double *ancestralBuffer, double *diagptable,
                    const int n, const int numStates, const int gammaStates) {
  int i;
  const int statesSquare = numStates * numStates;
  double *term = (double *)malloc(sizeof(double) * numStates);
  for (i = 0; i < n; i++) {
    double sum = 0.0, *_v = &x3[gammaStates * i],
           *ancestral = &ancestralBuffer[numStates * i];
    int k, j, l;
    for (l = 0; l < numStates; l++)
      term[l] = 0.0;
    for (k = 0; k < 4; k++) {
      double *v = &(_v[numStates * k]);
      for (l = 0; l < numStates; l++) {
        double al = 0.0;
        for (j = 0; j < numStates; j++)
          al += v[j] * diagptable[k * statesSquare + l * numStates + j];
        term[l] += al;
        sum += al;
      }
    }
    for (l = 0; l < numStates; l++)
      ancestral[l] = term[l] / sum;
  }
  free(term);
}


x86 assembly (clang, -O0):
<ancestralGamma>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    -0x20(%rbp),%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movslq -0x20(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1318 <ancestralGamma+0x1f8>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11da <ancestralGamma+0xba>
mov    -0x38(%rbp),%rax
movslq -0x5c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    11b0 <ancestralGamma+0x90>
movl   $0x0,-0x54(%rbp)
cmpl   $0x4,-0x54(%rbp)
jge    12c5 <ancestralGamma+0x1a5>
mov    -0x48(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x54(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x68(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    12b2 <ancestralGamma+0x192>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    127e <ancestralGamma+0x15e>
mov    -0x68(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x54(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
mov    -0x5c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
add    -0x58(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1226 <ancestralGamma+0x106>
movsd  -0x70(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x5c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    120b <ancestralGamma+0xeb>
jmp    12b7 <ancestralGamma+0x197>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    11e1 <ancestralGamma+0xc1>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1305 <ancestralGamma+0x1e5>
mov    -0x38(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    12cc <ancestralGamma+0x1ac>
jmp    130a <ancestralGamma+0x1ea>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1163 <ancestralGamma+0x43>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ancestralGamma>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %r9d,%r14d
mov    %r8d,%ebp
mov    %ecx,%r12d
mov    %rdx,%rbx
mov    %rsi,0x28(%rsp)
mov    %rdi,0x20(%rsp)
lea    0x0(,%rbp,8),%edi
call   1050 <malloc@plt>
mov    %rax,%rdi
test   %r12d,%r12d
jle    183e <ancestralGamma+0x70e>
mov    %rbp,%r9
mov    %r9d,%r11d
lea    0x0(,%r11,8),%rax
mov    %rax,0xb8(%rsp)
movslq %r9d,%rdx
movslq %r14d,%r8
mov    %r12d,%eax
mov    %rax,0xe8(%rsp)
lea    (%rdx,%rdx,1),%rax
mov    %rax,0x98(%rsp)
lea    0x0(%rbp,%rbp,1),%eax
mov    %rax,0x100(%rsp)
lea    (%rdx,%rdx,2),%r10
mov    %r11d,%r15d
and    $0xfffffffe,%r15d
lea    -0x2(%r15),%rax
mov    %rax,0x58(%rsp)
mov    %rax,%rcx
shr    %rcx
add    $0x1,%rcx
mov    %rcx,%rax
mov    %rcx,0x50(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x48(%rsp)
mov    %r11,%rax
neg    %rax
mov    %rax,0x68(%rsp)
mov    %r9d,%esi
imul   %esi,%esi
mov    %rdx,%rax
shl    $0x4,%rax
mov    0x20(%rsp),%rcx
lea    (%rax,%rcx,1),%r12
add    $0x8,%r12
lea    0x0(%rbp,%rbp,2),%eax
mov    %rax,0xf8(%rsp)
lea    (%rdi,%r11,8),%rax
mov    %rax,0x60(%rsp)
lea    -0x1(%r11),%rax
mov    %rax,0x8(%rsp)
lea    0x8(%rcx),%r14
mov    %r8,0xb0(%rsp)
lea    0x0(,%r8,8),%rax
mov    %rax,0xe0(%rsp)
lea    0x8(%rcx,%rdx,8),%r13
lea    0x1(%rsi),%eax
mov    %rax,0xa0(%rsp)
lea    (%rsi,%rsi,1),%eax
mov    %rax,0x90(%rsp)
lea    0x1(%rax),%rax
mov    %rax,0x88(%rsp)
mov    %r10,0x80(%rsp)
lea    0x8(%rcx,%r10,8),%rcx
lea    (%rsi,%rsi,2),%eax
mov    %rax,0x78(%rsp)
mov    %rsi,0xa8(%rsp)
lea    0x1(%rsi,%rsi,2),%eax
mov    %rax,0x70(%rsp)
mov    0x28(%rsp),%rax
lea    0x10(%rax),%rsi
mov    %rsi,0x30(%rsp)
mov    %rdx,0x38(%rsp)
lea    0x0(,%rdx,8),%rdx
mov    %rdx,0xd8(%rsp)
lea    0x8(%rax),%rsi
xor    %eax,%eax
mov    %rax,0x18(%rsp)
mov    %rbp,0xc8(%rsp)
mov    %rdi,0xc0(%rsp)
jmp    1318 <ancestralGamma+0x1e8>
nopl   0x0(%rax)
mov    0x18(%rsp),%rdx
add    $0x1,%rdx
mov    0xe0(%rsp),%rax
add    %rax,%r14
add    %rax,%r13
add    %rax,%r12
mov    0x40(%rsp),%rcx
add    %rax,%rcx
mov    0xd8(%rsp),%rax
add    %rax,0x30(%rsp)
add    %rax,%rsi
mov    %rdx,%rax
mov    %rdx,0x18(%rsp)
cmp    0xe8(%rsp),%rdx
je     183e <ancestralGamma+0x70e>
mov    %rcx,0x40(%rsp)
test   %r9d,%r9d
jle    12d0 <ancestralGamma+0x1a0>
mov    %rsi,(%rsp)
mov    0x18(%rsp),%rcx
mov    %rcx,%rax
imul   0x38(%rsp),%rax
imul   0xb0(%rsp),%rcx
mov    0x28(%rsp),%rdx
lea    (%rdx,%rax,8),%rsi
mov    %rsi,0x10(%rsp)
add    %r11,%rax
lea    (%rdx,%rax,8),%rax
mov    %rax,0xd0(%rsp)
mov    0x20(%rsp),%rax
lea    (%rax,%rcx,8),%rax
mov    %rax,0xf0(%rsp)
xor    %esi,%esi
mov    0xb8(%rsp),%rdx
mov    %r11,%rbp
call   1040 <memset@plt>
mov    0xf0(%rsp),%r8
mov    %rbp,%r11
mov    0xc0(%rsp),%rdi
mov    0xc8(%rsp),%r9
xorpd  %xmm0,%xmm0
mov    $0x1,%eax
xor    %ecx,%ecx
xor    %edx,%edx
jmp    13d1 <ancestralGamma+0x2a1>
nopw   0x0(%rax,%rax,1)
addsd  %xmm1,%xmm0
addsd  (%rdi,%rdx,8),%xmm1
movsd  %xmm1,(%rdi,%rdx,8)
add    $0x1,%rdx
add    %r11,%rax
add    %r11,%rcx
cmp    %r11,%rdx
je     1450 <ancestralGamma+0x320>
cmpq   $0x0,0x8(%rsp)
je     1420 <ancestralGamma+0x2f0>
xorpd  %xmm2,%xmm2
xor    %esi,%esi
nop
lea    (%rcx,%rsi,1),%ebp
movslq %ebp,%rbp
movsd  -0x8(%r14,%rsi,8),%xmm3
movsd  (%r14,%rsi,8),%xmm1
mulsd  (%rbx,%rbp,8),%xmm3
lea    (%rax,%rsi,1),%ebp
movslq %ebp,%rbp
mulsd  (%rbx,%rbp,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rsi
movapd %xmm1,%xmm2
cmp    %rsi,%r15
jne    13e0 <ancestralGamma+0x2b0>
test   $0x1,%r11b
je     13b0 <ancestralGamma+0x280>
jmp    142c <ancestralGamma+0x2fc>
xor    %esi,%esi
xorpd  %xmm1,%xmm1
test   $0x1,%r11b
je     13b0 <ancestralGamma+0x280>
mov    %edx,%ebp
imul   %r9d,%ebp
movsd  (%r8,%rsi,8),%xmm2
add    %ebp,%esi
movslq %esi,%rsi
mulsd  (%rbx,%rsi,8),%xmm2
addsd  %xmm2,%xmm1
jmp    13b0 <ancestralGamma+0x280>
nopl   0x0(%rax,%rax,1)
test   %r9d,%r9d
mov    (%rsp),%rsi
jle    12d0 <ancestralGamma+0x1a0>
mov    0x38(%rsp),%rax
mov    %r8,%rbp
lea    (%r8,%rax,8),%r8
mov    0xa8(%rsp),%rcx
mov    0xa0(%rsp),%rdx
xor    %r10d,%r10d
jmp    14a3 <ancestralGamma+0x373>
xchg   %ax,%ax
addsd  %xmm1,%xmm0
addsd  (%rdi,%r10,8),%xmm1
movsd  %xmm1,(%rdi,%r10,8)
add    $0x1,%r10
add    %r11,%rdx
add    %r11,%rcx
cmp    %r11,%r10
je     1540 <ancestralGamma+0x410>
cmpq   $0x0,0x8(%rsp)
je     1500 <ancestralGamma+0x3d0>
xorpd  %xmm2,%xmm2
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%rcx,%rsi,1),%eax
cltq
movsd  -0x8(%r13,%rsi,8),%xmm3
movsd  0x0(%r13,%rsi,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm3
lea    (%rdx,%rsi,1),%eax
cltq
mulsd  (%rbx,%rax,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rsi
movapd %xmm1,%xmm2
cmp    %rsi,%r15
jne    14c0 <ancestralGamma+0x390>
test   $0x1,%r11b
je     1480 <ancestralGamma+0x350>
jmp    1510 <ancestralGamma+0x3e0>
nop
xorpd  %xmm1,%xmm1
xor    %esi,%esi
test   $0x1,%r11b
je     1480 <ancestralGamma+0x350>
lea    (%r10,%r9,1),%eax
imul   %r9d,%eax
movsd  (%r8,%rsi,8),%xmm2
add    %eax,%esi
movslq %esi,%rax
mulsd  (%rbx,%rax,8),%xmm2
addsd  %xmm2,%xmm1
jmp    1480 <ancestralGamma+0x350>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
test   %r9d,%r9d
mov    (%rsp),%rsi
jle    12d0 <ancestralGamma+0x1a0>
mov    %rbp,%rcx
mov    0x98(%rsp),%rax
lea    0x0(,%rax,8),%r8
add    %rbp,%r8
mov    0x90(%rsp),%rcx
mov    0x88(%rsp),%rdx
xor    %r10d,%r10d
jmp    15a3 <ancestralGamma+0x473>
nopl   0x0(%rax,%rax,1)
addsd  %xmm1,%xmm0
addsd  (%rdi,%r10,8),%xmm1
movsd  %xmm1,(%rdi,%r10,8)
add    $0x1,%r10
add    %r11,%rdx
add    %r11,%rcx
cmp    %r11,%r10
je     1640 <ancestralGamma+0x510>
cmpq   $0x0,0x8(%rsp)
je     1600 <ancestralGamma+0x4d0>
xorpd  %xmm2,%xmm2
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%rcx,%rsi,1),%eax
cltq
movsd  -0x8(%r12,%rsi,8),%xmm3
movsd  (%r12,%rsi,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm3
lea    (%rdx,%rsi,1),%eax
cltq
mulsd  (%rbx,%rax,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rsi
movapd %xmm1,%xmm2
cmp    %rsi,%r15
jne    15c0 <ancestralGamma+0x490>
test   $0x1,%r11b
je     1580 <ancestralGamma+0x450>
jmp    1610 <ancestralGamma+0x4e0>
xchg   %ax,%ax
xorpd  %xmm1,%xmm1
xor    %esi,%esi
test   $0x1,%r11b
je     1580 <ancestralGamma+0x450>
mov    0x100(%rsp),%rax
add    %r10d,%eax
imul   %r9d,%eax
movsd  (%r8,%rsi,8),%xmm2
add    %eax,%esi
movslq %esi,%rax
mulsd  (%rbx,%rax,8),%xmm2
addsd  %xmm2,%xmm1
jmp    1580 <ancestralGamma+0x450>
nopl   0x0(%rax,%rax,1)
test   %r9d,%r9d
mov    (%rsp),%rsi
jle    12d0 <ancestralGamma+0x1a0>
mov    0x80(%rsp),%rax
lea    0x0(,%rax,8),%r8
add    %rbp,%r8
mov    0x78(%rsp),%rcx
mov    0x70(%rsp),%rdx
xor    %r10d,%r10d
jmp    1693 <ancestralGamma+0x563>
nop
addsd  %xmm1,%xmm0
addsd  (%rdi,%r10,8),%xmm1
movsd  %xmm1,(%rdi,%r10,8)
add    $0x1,%r10
add    %r11,%rdx
add    %r11,%rcx
cmp    %r11,%r10
je     1728 <ancestralGamma+0x5f8>
cmpq   $0x0,0x8(%rsp)
je     16f0 <ancestralGamma+0x5c0>
xorpd  %xmm2,%xmm2
xor    %esi,%esi
mov    0x40(%rsp),%rbp
cs nopw 0x0(%rax,%rax,1)
lea    (%rcx,%rsi,1),%eax
cltq
movsd  -0x8(%rbp,%rsi,8),%xmm3
movsd  0x0(%rbp,%rsi,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm3
lea    (%rdx,%rsi,1),%eax
cltq
mulsd  (%rbx,%rax,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rsi
movapd %xmm1,%xmm2
cmp    %rsi,%r15
jne    16b0 <ancestralGamma+0x580>
test   $0x1,%r11b
je     1670 <ancestralGamma+0x540>
jmp    1700 <ancestralGamma+0x5d0>
nopl   (%rax)
xorpd  %xmm1,%xmm1
xor    %esi,%esi
test   $0x1,%r11b
je     1670 <ancestralGamma+0x540>
mov    0xf8(%rsp),%rax
add    %r10d,%eax
imul   %r9d,%eax
movsd  (%r8,%rsi,8),%xmm2
add    %eax,%esi
movslq %esi,%rax
mulsd  (%rbx,%rax,8),%xmm2
addsd  %xmm2,%xmm1
jmp    1670 <ancestralGamma+0x540>
test   %r9d,%r9d
mov    (%rsp),%rsi
jle    12d0 <ancestralGamma+0x1a0>
cmp    $0x1,%r9d
je     1751 <ancestralGamma+0x621>
mov    0x10(%rsp),%rax
cmp    0x60(%rsp),%rax
jae    17bc <ancestralGamma+0x68c>
cmp    0xd0(%rsp),%rdi
jae    17bc <ancestralGamma+0x68c>
xor    %ecx,%ecx
mov    %rcx,%rax
test   $0x1,%r11b
je     1776 <ancestralGamma+0x646>
movsd  (%rdi,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
mov    0x10(%rsp),%rax
movsd  %xmm1,(%rax,%rcx,8)
mov    %rcx,%rax
or     $0x1,%rax
not    %rcx
cmp    0x68(%rsp),%rcx
je     12d0 <ancestralGamma+0x1a0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rdi,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%rsi,%rax,8)
movsd  0x8(%rdi,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rax,8)
add    $0x2,%rax
cmp    %rax,%r11
jne    1790 <ancestralGamma+0x660>
jmp    12d0 <ancestralGamma+0x1a0>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
cmpq   $0x0,0x58(%rsp)
je     1833 <ancestralGamma+0x703>
mov    0x48(%rsp),%rcx
xor    %eax,%eax
mov    0x30(%rsp),%rdx
nopl   0x0(%rax,%rax,1)
movupd (%rdi,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rdx,%rax,8)
movupd 0x10(%rdi,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdx,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rcx
jne    17e0 <ancestralGamma+0x6b0>
testb  $0x1,0x50(%rsp)
je     1822 <ancestralGamma+0x6f2>
movupd (%rdi,%rax,8),%xmm2
divpd  %xmm1,%xmm2
mov    0x10(%rsp),%rcx
movupd %xmm2,(%rcx,%rax,8)
mov    %r15,%rcx
cmp    %r11,%r15
je     12d0 <ancestralGamma+0x1a0>
jmp    1753 <ancestralGamma+0x623>
xor    %eax,%eax
testb  $0x1,0x50(%rsp)
jne    180f <ancestralGamma+0x6df>
jmp    1822 <ancestralGamma+0x6f2>
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @ancestralGamma(double* %0, double* %1, double* %2, i32 %3, i32 %4, i32 %5) {
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double*
%16 = alloca double
%17 = alloca double*
%18 = alloca double*
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double*
%23 = alloca double
store double* %0, double** %7
store double* %1, double** %8
store double* %2, double** %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store i32 %5, i32* %12
%24 = load i32, i32* %11
%25 = load i32, i32* %11
%26 = mul nsw i32 %24, %25
store i32 %26, i32* %14
%27 = load i32, i32* %11
%28 = sext i32 %27 to i64
%29 = mul i64 8, %28
%30 = trunc i64 %29 to i32
%31 = call i64 @malloc(i32 %30)
%32 = inttoptr i64 %31 to double*
store double* %32, double** %15
store i32 0, i32* %13
br label %33
33: ; preds = %145, %6
%34 = load i32, i32* %13
%35 = load i32, i32* %10
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %148
37: ; preds = %33
store double 0.000000e+00, double* %16
%38 = load double*, double** %7
%39 = load i32, i32* %12
%40 = load i32, i32* %13
%41 = mul nsw i32 %39, %40
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %38, i64 %42
store double* %43, double** %17
%44 = load double*, double** %8
%45 = load i32, i32* %11
%46 = load i32, i32* %13
%47 = mul nsw i32 %45, %46
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %44, i64 %48
store double* %49, double** %18
store i32 0, i32* %21
br label %50
50: ; preds = %59, %37
%51 = load i32, i32* %21
%52 = load i32, i32* %11
%53 = icmp slt i32 %51, %52
br i1 %53, label %54, label %62
54: ; preds = %50
%55 = load double*, double** %15
%56 = load i32, i32* %21
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %55, i64 %57
store double 0.000000e+00, double* %58
br label %59
59: ; preds = %54
%60 = load i32, i32* %21
%61 = add nsw i32 %60, 1
store i32 %61, i32* %21
br label %50
62: ; preds = %50
store i32 0, i32* %19
br label %63
63: ; preds = %121, %62
%64 = load i32, i32* %19
%65 = icmp slt i32 %64, 4
br i1 %65, label %66, label %124
66: ; preds = %63
%67 = load double*, double** %17
%68 = load i32, i32* %11
%69 = load i32, i32* %19
%70 = mul nsw i32 %68, %69
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double, double* %67, i64 %71
store double* %72, double** %22
store i32 0, i32* %21
br label %73
73: ; preds = %117, %66
%74 = load i32, i32* %21
%75 = load i32, i32* %11
%76 = icmp slt i32 %74, %75
br i1 %76, label %77, label %120
77: ; preds = %73
store double 0.000000e+00, double* %23
store i32 0, i32* %20
br label %78
78: ; preds = %103, %77
%79 = load i32, i32* %20
%80 = load i32, i32* %11
%81 = icmp slt i32 %79, %80
br i1 %81, label %82, label %106
82: ; preds = %78
%83 = load double*, double** %22
%84 = load i32, i32* %20
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
%87 = load double, double* %86
%88 = load double*, double** %9
%89 = load i32, i32* %19
%90 = load i32, i32* %14
%91 = mul nsw i32 %89, %90
%92 = load i32, i32* %21
%93 = load i32, i32* %11
%94 = mul nsw i32 %92, %93
%95 = add nsw i32 %91, %94
%96 = load i32, i32* %20
%97 = add nsw i32 %95, %96
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %88, i64 %98
%100 = load double, double* %99
%101 = load double, double* %23
%102 = call double @llvm.fmuladd.f64(double %87, double %100, double %101)
store double %102, double* %23
br label %103
103: ; preds = %82
%104 = load i32, i32* %20
%105 = add nsw i32 %104, 1
store i32 %105, i32* %20
br label %78
106: ; preds = %78
%107 = load double, double* %23
%108 = load double*, double** %15
%109 = load i32, i32* %21
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %108, i64 %110
%112 = load double, double* %111
%113 = fadd double %112, %107
store double %113, double* %111
%114 = load double, double* %23
%115 = load double, double* %16
%116 = fadd double %115, %114
store double %116, double* %16
br label %117
117: ; preds = %106
%118 = load i32, i32* %21
%119 = add nsw i32 %118, 1
store i32 %119, i32* %21
br label %73
120: ; preds = %73
br label %121
121: ; preds = %120
%122 = load i32, i32* %19
%123 = add nsw i32 %122, 1
store i32 %123, i32* %19
br label %63
124: ; preds = %63
store i32 0, i32* %21
br label %125
125: ; preds = %141, %124
%126 = load i32, i32* %21
%127 = load i32, i32* %11
%128 = icmp slt i32 %126, %127
br i1 %128, label %129, label %144
129: ; preds = %125
%130 = load double*, double** %15
%131 = load i32, i32* %21
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %130, i64 %132
%134 = load double, double* %133
%135 = load double, double* %16
%136 = fdiv double %134, %135
%137 = load double*, double** %18
%138 = load i32, i32* %21
%139 = sext i32 %138 to i64
%140 = getelementptr inbounds double, double* %137, i64 %139
store double %136, double* %140
br label %141
141: ; preds = %129
%142 = load i32, i32* %21
%143 = add nsw i32 %142, 1
store i32 %143, i32* %21
br label %125
144: ; preds = %125
br label %145
145: ; preds = %144
%146 = load i32, i32* %13
%147 = add nsw i32 %146, 1
store i32 %147, i32* %13
br label %33
148: ; preds = %33
%149 = load double*, double** %15
%150 = call i32 @free(double* %149)
ret void
}

LLVM IR (-O2):
define void @ancestralGamma(double* nocapture %0, double* nocapture writeonly %1, double* nocapture %2, i32 %3, i32 %4, i32 %5) {
%7 = shl i32 %4, 3
%8 = tail call i64 @malloc(i32 %7)
%9 = inttoptr i64 %8 to double*
%10 = inttoptr i64 %8 to i8*
%11 = inttoptr i64 %8 to i8*
%12 = icmp sgt i32 %3, 0
br i1 %12, label %13, label %343
13: ; preds = %6
%14 = icmp sgt i32 %4, 0
%15 = zext i32 %4 to i64
%16 = shl nuw nsw i64 %15, 3
%17 = sext i32 %4 to i64
%18 = sext i32 %5 to i64
%19 = zext i32 %3 to i64
%20 = shl nsw i64 %17, 1
%21 = shl nuw i32 %4, 1
%22 = mul nsw i64 %17, 3
%23 = mul i32 %4, 3
%24 = getelementptr double, double* %9, i64 %15
%25 = add nsw i64 %15, -1
%26 = and i64 %15, 4294967294
%27 = add nsw i64 %26, -2
%28 = lshr exact i64 %27, 1
%29 = add nuw i64 %28, 1
%30 = and i64 %15, 1
%31 = icmp eq i64 %25, 0
%32 = and i64 %15, 4294967294
%33 = icmp eq i64 %30, 0
%34 = and i64 %15, 1
%35 = icmp eq i64 %25, 0
%36 = and i64 %15, 4294967294
%37 = icmp eq i64 %34, 0
%38 = and i64 %15, 1
%39 = icmp eq i64 %25, 0
%40 = and i64 %15, 4294967294
%41 = icmp eq i64 %38, 0
%42 = and i64 %15, 1
%43 = icmp eq i64 %25, 0
%44 = and i64 %15, 4294967294
%45 = icmp eq i64 %42, 0
%46 = icmp eq i32 %4, 1
%47 = and i64 %15, 4294967294
%48 = and i64 %29, 1
%49 = icmp eq i64 %27, 0
%50 = and i64 %29, -2
%51 = icmp eq i64 %48, 0
%52 = icmp eq i64 %47, %15
%53 = and i64 %15, 1
%54 = icmp eq i64 %53, 0
%55 = sub nsw i64 0, %15
br label %56
56: ; preds = %13, %340
%57 = phi i64 [ 0, %13 ], [ %341, %340 ]
%58 = mul i64 %57, %17
%59 = getelementptr double, double* %1, i64 %58
%60 = add i64 %58, %15
%61 = getelementptr double, double* %1, i64 %60
%62 = bitcast double* %61 to i8*
%63 = mul nsw i64 %57, %18
%64 = getelementptr inbounds double, double* %0, i64 %63
%65 = mul nsw i64 %57, %17
%66 = getelementptr inbounds double, double* %1, i64 %65
br i1 %14, label %67, label %340
67: ; preds = %56
br label %68
68: ; preds = %110, %67
%69 = phi i64 [ 0, %67 ], [ %116, %110 ]
%70 = phi double [ 0.000000e+00, %67 ], [ %115, %110 ]
%71 = trunc i64 %69 to i32
%72 = mul i32 %71, %4
br i1 %31, label %97, label %73
73: ; preds = %68, %73
%74 = phi i64 [ %94, %73 ], [ 0, %68 ]
%75 = phi double [ %93, %73 ], [ 0.000000e+00, %68 ]
%76 = phi i64 [ %95, %73 ], [ 0, %68 ]
%77 = getelementptr inbounds double, double* %64, i64 %74
%78 = load double, double* %77
%79 = trunc i64 %74 to i32
%80 = add i32 %72, %79
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %2, i64 %81
%83 = load double, double* %82
%84 = tail call double @llvm.fmuladd.f64(double %78, double %83, double %75)
%85 = or i64 %74, 1
%86 = getelementptr inbounds double, double* %64, i64 %85
%87 = load double, double* %86
%88 = trunc i64 %85 to i32
%89 = add i32 %72, %88
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %2, i64 %90
%92 = load double, double* %91
%93 = tail call double @llvm.fmuladd.f64(double %87, double %92, double %84)
%94 = add nuw nsw i64 %74, 2
%95 = add i64 %76, 2
%96 = icmp eq i64 %95, %32
br i1 %96, label %97, label %73
97: ; preds = %73, %68
%98 = phi double [ undef, %68 ], [ %93, %73 ]
%99 = phi i64 [ 0, %68 ], [ %94, %73 ]
%100 = phi double [ 0.000000e+00, %68 ], [ %93, %73 ]
br i1 %33, label %110, label %101
101: ; preds = %97
%102 = getelementptr inbounds double, double* %64, i64 %99
%103 = load double, double* %102
%104 = trunc i64 %99 to i32
%105 = add i32 %72, %104
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds double, double* %2, i64 %106
%108 = load double, double* %107
%109 = tail call double @llvm.fmuladd.f64(double %103, double %108, double %100)
br label %110
110: ; preds = %97, %101
%111 = phi double [ %98, %97 ], [ %109, %101 ]
%112 = getelementptr inbounds double, double* %9, i64 %69
%113 = load double, double* %112
%114 = fadd double %111, %113
store double %114, double* %112
%115 = fadd double %70, %111
%116 = add nuw nsw i64 %69, 1
%117 = icmp eq i64 %116, %15
br i1 %117, label %118, label %68
118: ; preds = %110
%119 = getelementptr inbounds double, double* %64, i64 %17
br i1 %14, label %120, label %340
120: ; preds = %118, %163
%121 = phi i64 [ %169, %163 ], [ 0, %118 ]
%122 = phi double [ %168, %163 ], [ %115, %118 ]
%123 = trunc i64 %121 to i32
%124 = add i32 %123, %4
%125 = mul i32 %124, %4
br i1 %35, label %150, label %126
126: ; preds = %120, %126
%127 = phi i64 [ %147, %126 ], [ 0, %120 ]
%128 = phi double [ %146, %126 ], [ 0.000000e+00, %120 ]
%129 = phi i64 [ %148, %126 ], [ 0, %120 ]
%130 = getelementptr inbounds double, double* %119, i64 %127
%131 = load double, double* %130
%132 = trunc i64 %127 to i32
%133 = add i32 %125, %132
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %2, i64 %134
%136 = load double, double* %135
%137 = tail call double @llvm.fmuladd.f64(double %131, double %136, double %128)
%138 = or i64 %127, 1
%139 = getelementptr inbounds double, double* %119, i64 %138
%140 = load double, double* %139
%141 = trunc i64 %138 to i32
%142 = add i32 %125, %141
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %2, i64 %143
%145 = load double, double* %144
%146 = tail call double @llvm.fmuladd.f64(double %140, double %145, double %137)
%147 = add nuw nsw i64 %127, 2
%148 = add i64 %129, 2
%149 = icmp eq i64 %148, %36
br i1 %149, label %150, label %126
150: ; preds = %126, %120
%151 = phi double [ undef, %120 ], [ %146, %126 ]
%152 = phi i64 [ 0, %120 ], [ %147, %126 ]
%153 = phi double [ 0.000000e+00, %120 ], [ %146, %126 ]
br i1 %37, label %163, label %154
154: ; preds = %150
%155 = getelementptr inbounds double, double* %119, i64 %152
%156 = load double, double* %155
%157 = trunc i64 %152 to i32
%158 = add i32 %125, %157
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds double, double* %2, i64 %159
%161 = load double, double* %160
%162 = tail call double @llvm.fmuladd.f64(double %156, double %161, double %153)
br label %163
163: ; preds = %150, %154
%164 = phi double [ %151, %150 ], [ %162, %154 ]
%165 = getelementptr inbounds double, double* %9, i64 %121
%166 = load double, double* %165
%167 = fadd double %164, %166
store double %167, double* %165
%168 = fadd double %122, %164
%169 = add nuw nsw i64 %121, 1
%170 = icmp eq i64 %169, %15
br i1 %170, label %171, label %120
171: ; preds = %163
%172 = getelementptr inbounds double, double* %64, i64 %20
br i1 %14, label %173, label %340
173: ; preds = %171, %216
%174 = phi i64 [ %222, %216 ], [ 0, %171 ]
%175 = phi double [ %221, %216 ], [ %168, %171 ]
%176 = trunc i64 %174 to i32
%177 = add i32 %21, %176
%178 = mul i32 %177, %4
br i1 %39, label %203, label %179
179: ; preds = %173, %179
%180 = phi i64 [ %200, %179 ], [ 0, %173 ]
%181 = phi double [ %199, %179 ], [ 0.000000e+00, %173 ]
%182 = phi i64 [ %201, %179 ], [ 0, %173 ]
%183 = getelementptr inbounds double, double* %172, i64 %180
%184 = load double, double* %183
%185 = trunc i64 %180 to i32
%186 = add i32 %178, %185
%187 = sext i32 %186 to i64
%188 = getelementptr inbounds double, double* %2, i64 %187
%189 = load double, double* %188
%190 = tail call double @llvm.fmuladd.f64(double %184, double %189, double %181)
%191 = or i64 %180, 1
%192 = getelementptr inbounds double, double* %172, i64 %191
%193 = load double, double* %192
%194 = trunc i64 %191 to i32
%195 = add i32 %178, %194
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds double, double* %2, i64 %196
%198 = load double, double* %197
%199 = tail call double @llvm.fmuladd.f64(double %193, double %198, double %190)
%200 = add nuw nsw i64 %180, 2
%201 = add i64 %182, 2
%202 = icmp eq i64 %201, %40
br i1 %202, label %203, label %179
203: ; preds = %179, %173
%204 = phi double [ undef, %173 ], [ %199, %179 ]
%205 = phi i64 [ 0, %173 ], [ %200, %179 ]
%206 = phi double [ 0.000000e+00, %173 ], [ %199, %179 ]
br i1 %41, label %216, label %207
207: ; preds = %203
%208 = getelementptr inbounds double, double* %172, i64 %205
%209 = load double, double* %208
%210 = trunc i64 %205 to i32
%211 = add i32 %178, %210
%212 = sext i32 %211 to i64
%213 = getelementptr inbounds double, double* %2, i64 %212
%214 = load double, double* %213
%215 = tail call double @llvm.fmuladd.f64(double %209, double %214, double %206)
br label %216
216: ; preds = %203, %207
%217 = phi double [ %204, %203 ], [ %215, %207 ]
%218 = getelementptr inbounds double, double* %9, i64 %174
%219 = load double, double* %218
%220 = fadd double %217, %219
store double %220, double* %218
%221 = fadd double %175, %217
%222 = add nuw nsw i64 %174, 1
%223 = icmp eq i64 %222, %15
br i1 %223, label %224, label %173
224: ; preds = %216
%225 = getelementptr inbounds double, double* %64, i64 %22
br i1 %14, label %226, label %340
226: ; preds = %224, %269
%227 = phi i64 [ %275, %269 ], [ 0, %224 ]
%228 = phi double [ %274, %269 ], [ %221, %224 ]
%229 = trunc i64 %227 to i32
%230 = add i32 %23, %229
%231 = mul i32 %230, %4
br i1 %43, label %256, label %232
232: ; preds = %226, %232
%233 = phi i64 [ %253, %232 ], [ 0, %226 ]
%234 = phi double [ %252, %232 ], [ 0.000000e+00, %226 ]
%235 = phi i64 [ %254, %232 ], [ 0, %226 ]
%236 = getelementptr inbounds double, double* %225, i64 %233
%237 = load double, double* %236
%238 = trunc i64 %233 to i32
%239 = add i32 %231, %238
%240 = sext i32 %239 to i64
%241 = getelementptr inbounds double, double* %2, i64 %240
%242 = load double, double* %241
%243 = tail call double @llvm.fmuladd.f64(double %237, double %242, double %234)
%244 = or i64 %233, 1
%245 = getelementptr inbounds double, double* %225, i64 %244
%246 = load double, double* %245
%247 = trunc i64 %244 to i32
%248 = add i32 %231, %247
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds double, double* %2, i64 %249
%251 = load double, double* %250
%252 = tail call double @llvm.fmuladd.f64(double %246, double %251, double %243)
%253 = add nuw nsw i64 %233, 2
%254 = add i64 %235, 2
%255 = icmp eq i64 %254, %44
br i1 %255, label %256, label %232
256: ; preds = %232, %226
%257 = phi double [ undef, %226 ], [ %252, %232 ]
%258 = phi i64 [ 0, %226 ], [ %253, %232 ]
%259 = phi double [ 0.000000e+00, %226 ], [ %252, %232 ]
br i1 %45, label %269, label %260
260: ; preds = %256
%261 = getelementptr inbounds double, double* %225, i64 %258
%262 = load double, double* %261
%263 = trunc i64 %258 to i32
%264 = add i32 %231, %263
%265 = sext i32 %264 to i64
%266 = getelementptr inbounds double, double* %2, i64 %265
%267 = load double, double* %266
%268 = tail call double @llvm.fmuladd.f64(double %262, double %267, double %259)
br label %269
269: ; preds = %256, %260
%270 = phi double [ %257, %256 ], [ %268, %260 ]
%271 = getelementptr inbounds double, double* %9, i64 %227
%272 = load double, double* %271
%273 = fadd double %270, %272
store double %273, double* %271
%274 = fadd double %228, %270
%275 = add nuw nsw i64 %227, 1
%276 = icmp eq i64 %275, %15
br i1 %276, label %277, label %226
277: ; preds = %269
br i1 %14, label %278, label %340
278: ; preds = %277
br i1 %46, label %315, label %279
279: ; preds = %278
%280 = icmp ult double* %59, %24
%281 = icmp ult i8* %10, %62
%282 = and i1 %280, %281
br i1 %282, label %315, label %283
283: ; preds = %279
%284 = insertelement <2 x double> poison, double %274, i64 0
%285 = shufflevector <2 x double> %284, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %49, label %305, label %286
286: ; preds = %283, %286
%287 = phi i64 [ %302, %286 ], [ 0, %283 ]
%288 = phi i64 [ %303, %286 ], [ 0, %283 ]
%289 = getelementptr inbounds double, double* %9, i64 %287
%290 = bitcast double* %289 to <2 x double>*
%291 = load <2 x double>, <2 x double>* %290
%292 = fdiv <2 x double> %291, %285
%293 = getelementptr inbounds double, double* %66, i64 %287
%294 = bitcast double* %293 to <2 x double>*
store <2 x double> %292, <2 x double>* %294
%295 = or i64 %287, 2
%296 = getelementptr inbounds double, double* %9, i64 %295
%297 = bitcast double* %296 to <2 x double>*
%298 = load <2 x double>, <2 x double>* %297
%299 = fdiv <2 x double> %298, %285
%300 = getelementptr inbounds double, double* %66, i64 %295
%301 = bitcast double* %300 to <2 x double>*
store <2 x double> %299, <2 x double>* %301
%302 = add nuw i64 %287, 4
%303 = add nuw i64 %288, 2
%304 = icmp eq i64 %303, %50
br i1 %304, label %305, label %286
305: ; preds = %286, %283
%306 = phi i64 [ 0, %283 ], [ %302, %286 ]
br i1 %51, label %314, label %307
307: ; preds = %305
%308 = getelementptr inbounds double, double* %9, i64 %306
%309 = bitcast double* %308 to <2 x double>*
%310 = load <2 x double>, <2 x double>* %309
%311 = fdiv <2 x double> %310, %285
%312 = getelementptr inbounds double, double* %66, i64 %306
%313 = bitcast double* %312 to <2 x double>*
store <2 x double> %311, <2 x double>* %313
br label %314
314: ; preds = %305, %307
br i1 %52, label %340, label %315
315: ; preds = %279, %278, %314
%316 = phi i64 [ 0, %279 ], [ 0, %278 ], [ %47, %314 ]
%317 = xor i64 %316, -1
br i1 %54, label %324, label %318
318: ; preds = %315
%319 = getelementptr inbounds double, double* %9, i64 %316
%320 = load double, double* %319
%321 = fdiv double %320, %274
%322 = getelementptr inbounds double, double* %66, i64 %316
store double %321, double* %322
%323 = or i64 %316, 1
br label %324
324: ; preds = %318, %315
%325 = phi i64 [ %316, %315 ], [ %323, %318 ]
%326 = icmp eq i64 %317, %55
br i1 %326, label %340, label %327
327: ; preds = %324, %327
%328 = phi i64 [ %338, %327 ], [ %325, %324 ]
%329 = getelementptr inbounds double, double* %9, i64 %328
%330 = load double, double* %329
%331 = fdiv double %330, %274
%332 = getelementptr inbounds double, double* %66, i64 %328
store double %331, double* %332
%333 = add nuw nsw i64 %328, 1
%334 = getelementptr inbounds double, double* %9, i64 %333
%335 = load double, double* %334
%336 = fdiv double %335, %274
%337 = getelementptr inbounds double, double* %66, i64 %333
store double %336, double* %337
%338 = add nuw nsw i64 %328, 2
%339 = icmp eq i64 %338, %15
br i1 %339, label %340, label %327
340: ; preds = %324, %327, %314, %56, %118, %171, %224, %277
%341 = add nuw nsw i64 %57, 1
%342 = icmp eq i64 %341, %19
br i1 %342, label %343, label %56
343: ; preds = %340, %6
%344 = tail call i32 @free(double* %9)
ret void
}