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


x86 assembly (clang, -O0):
<ancestralGamma>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%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,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movslq -0x20(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1030 <rax_malloc@plt>
mov    %rax,-0x58(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11ae <ancestralGamma+0x8e>
mov    -0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1184 <ancestralGamma+0x64>
movl   $0x0,-0x2c(%rbp)
cmpl   $0x4,-0x2c(%rbp)
jge    1299 <ancestralGamma+0x179>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1286 <ancestralGamma+0x166>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1252 <ancestralGamma+0x132>
mov    -0x40(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
mov    -0x34(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11fa <ancestralGamma+0xda>
movsd  -0x60(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11df <ancestralGamma+0xbf>
jmp    128b <ancestralGamma+0x16b>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b5 <ancestralGamma+0x95>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    12d9 <ancestralGamma+0x1b9>
mov    -0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x50(%rbp),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12a0 <ancestralGamma+0x180>
mov    -0x58(%rbp),%rdi
call   1040 <rax_free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ancestralGamma>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r9d,%r15d
mov    %r8d,%ebp
mov    %ecx,%r14d
mov    %rdx,%rbx
mov    %rsi,0x8(%rsp)
mov    %rdi,%r13
lea    0x0(,%rbp,8),%edi
call   1030 <rax_malloc@plt>
mov    %rax,%r12
test   %ebp,%ebp
jle    1640 <ancestralGamma+0x510>
imul   %r14d,%r15d
movslq %r15d,%rcx
mov    0x8(%rsp),%rax
mov    %rcx,0x10(%rsp)
lea    (%rax,%rcx,8),%rax
mov    %rax,0x18(%rsp)
mov    %ebp,%r15d
lea    0x0(,%r15,8),%rdx
xor    %r14d,%r14d
mov    %r12,%rdi
xor    %esi,%esi
mov    %rbp,0x20(%rsp)
call   1050 <memset@plt>
mov    0x20(%rsp),%r10
lea    -0x1(%r15),%r11
mov    %r15d,%ecx
and    $0xfffffffe,%ecx
xorpd  %xmm0,%xmm0
mov    $0x1,%eax
xor    %edx,%edx
jmp    11e3 <ancestralGamma+0xb3>
nopl   0x0(%rax)
addsd  %xmm1,%xmm0
addsd  (%r12,%rdx,8),%xmm1
movsd  %xmm1,(%r12,%rdx,8)
add    $0x1,%rdx
add    %r15,%rax
add    %r15,%r14
cmp    %r15,%rdx
je     1270 <ancestralGamma+0x140>
test   %r11,%r11
je     1240 <ancestralGamma+0x110>
xorpd  %xmm2,%xmm2
xor    %esi,%esi
xchg   %ax,%ax
lea    (%r14,%rsi,1),%edi
movslq %edi,%rdi
movsd  0x0(%r13,%rsi,8),%xmm3
movsd  0x8(%r13,%rsi,8),%xmm1
mulsd  (%rbx,%rdi,8),%xmm3
lea    (%rax,%rsi,1),%edi
movslq %edi,%rdi
mulsd  (%rbx,%rdi,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rsi
movapd %xmm1,%xmm2
cmp    %rsi,%rcx
jne    11f0 <ancestralGamma+0xc0>
test   $0x1,%r15b
je     11c0 <ancestralGamma+0x90>
jmp    1250 <ancestralGamma+0x120>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xor    %esi,%esi
xorpd  %xmm1,%xmm1
test   $0x1,%r15b
je     11c0 <ancestralGamma+0x90>
mov    %edx,%edi
imul   %r10d,%edi
movsd  0x0(%r13,%rsi,8),%xmm2
add    %edi,%esi
movslq %esi,%rsi
mulsd  (%rbx,%rsi,8),%xmm2
addsd  %xmm2,%xmm1
jmp    11c0 <ancestralGamma+0x90>
test   %r10d,%r10d
jle    1640 <ancestralGamma+0x510>
lea    0x0(,%r15,8),%rdx
add    %r13,%rdx
mov    %r15d,%esi
and    $0xfffffffe,%esi
mov    %r10d,%r8d
imul   %r8d,%r8d
lea    0x1(%r8),%ecx
xor    %r9d,%r9d
mov    %r8,%rax
jmp    12bf <ancestralGamma+0x18f>
nopl   (%rax)
addsd  %xmm1,%xmm0
addsd  (%r12,%r9,8),%xmm1
movsd  %xmm1,(%r12,%r9,8)
add    $0x1,%r9
add    %r15,%rcx
add    %r15,%rax
cmp    %r15,%r9
je     133c <ancestralGamma+0x20c>
test   %r11,%r11
je     1310 <ancestralGamma+0x1e0>
xorpd  %xmm2,%xmm2
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
lea    (%rax,%rbp,1),%edi
movslq %edi,%rdi
movsd  (%rdx,%rbp,8),%xmm3
movsd  0x8(%rdx,%rbp,8),%xmm1
mulsd  (%rbx,%rdi,8),%xmm3
lea    (%rcx,%rbp,1),%edi
movslq %edi,%rdi
mulsd  (%rbx,%rdi,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rbp
movapd %xmm1,%xmm2
cmp    %rbp,%rsi
jne    12d0 <ancestralGamma+0x1a0>
test   $0x1,%r15b
je     12a0 <ancestralGamma+0x170>
jmp    131c <ancestralGamma+0x1ec>
xchg   %ax,%ax
xorpd  %xmm1,%xmm1
xor    %ebp,%ebp
test   $0x1,%r15b
je     12a0 <ancestralGamma+0x170>
lea    (%r9,%r10,1),%edi
imul   %r10d,%edi
movsd  (%rdx,%rbp,8),%xmm2
add    %edi,%ebp
movslq %ebp,%rdi
mulsd  (%rbx,%rdi,8),%xmm2
addsd  %xmm2,%xmm1
jmp    12a0 <ancestralGamma+0x170>
test   %r10d,%r10d
jle    1640 <ancestralGamma+0x510>
lea    (%r15,%r15,1),%rax
mov    %r10,%rcx
mov    %r15,%r10
shl    $0x4,%r10
add    %r13,%r10
lea    (%rcx,%rcx,1),%r9d
mov    %r15d,%edi
and    $0xfffffffe,%edi
lea    0x8(,%rax,8),%rcx
add    %r13,%rcx
lea    (%r8,%r8,1),%esi
lea    0x1(%rsi),%rdx
xor    %r14d,%r14d
jmp    139f <ancestralGamma+0x26f>
nopl   0x0(%rax,%rax,1)
addsd  %xmm1,%xmm0
addsd  (%r12,%r14,8),%xmm1
movsd  %xmm1,(%r12,%r14,8)
add    $0x1,%r14
add    %r15,%rdx
add    %r15,%rsi
cmp    %r15,%r14
je     141e <ancestralGamma+0x2ee>
test   %r11,%r11
je     13f0 <ancestralGamma+0x2c0>
xorpd  %xmm2,%xmm2
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
lea    (%rsi,%rbp,1),%eax
cltq
movsd  -0x8(%rcx,%rbp,8),%xmm3
movsd  (%rcx,%rbp,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm3
lea    (%rdx,%rbp,1),%eax
cltq
mulsd  (%rbx,%rax,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rbp
movapd %xmm1,%xmm2
cmp    %rbp,%rdi
jne    13b0 <ancestralGamma+0x280>
test   $0x1,%r15b
je     1380 <ancestralGamma+0x250>
jmp    13fc <ancestralGamma+0x2cc>
nopl   0x0(%rax)
xorpd  %xmm1,%xmm1
xor    %ebp,%ebp
test   $0x1,%r15b
je     1380 <ancestralGamma+0x250>
lea    (%r9,%r14,1),%eax
imul   0x20(%rsp),%eax
movsd  (%r10,%rbp,8),%xmm2
add    %eax,%ebp
movslq %ebp,%rax
mulsd  (%rbx,%rax,8),%xmm2
addsd  %xmm2,%xmm1
jmp    1380 <ancestralGamma+0x250>
mov    0x20(%rsp),%r10
test   %r10d,%r10d
jle    1640 <ancestralGamma+0x510>
lea    (%r15,%r15,2),%rax
lea    0x0(,%rax,8),%rdx
add    %r13,%rdx
lea    (%r10,%r10,2),%r9d
mov    %r15d,%edi
and    $0xfffffffe,%edi
lea    (%r8,%r8,2),%ecx
lea    (%r8,%r8,2),%esi
add    $0x1,%esi
xor    %r8d,%r8d
jmp    147f <ancestralGamma+0x34f>
cs nopw 0x0(%rax,%rax,1)
nop
addsd  %xmm1,%xmm0
addsd  (%r12,%r8,8),%xmm1
movsd  %xmm1,(%r12,%r8,8)
add    $0x1,%r8
add    %r15,%rsi
add    %r15,%rcx
cmp    %r15,%r8
je     14fc <ancestralGamma+0x3cc>
test   %r11,%r11
je     14d0 <ancestralGamma+0x3a0>
xorpd  %xmm2,%xmm2
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
lea    (%rcx,%rbp,1),%eax
cltq
movsd  (%rdx,%rbp,8),%xmm3
movsd  0x8(%rdx,%rbp,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm3
lea    (%rsi,%rbp,1),%eax
cltq
mulsd  (%rbx,%rax,8),%xmm1
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0x2,%rbp
movapd %xmm1,%xmm2
cmp    %rbp,%rdi
jne    1490 <ancestralGamma+0x360>
test   $0x1,%r15b
je     1460 <ancestralGamma+0x330>
jmp    14dc <ancestralGamma+0x3ac>
nopl   0x0(%rax)
xorpd  %xmm1,%xmm1
xor    %ebp,%ebp
test   $0x1,%r15b
je     1460 <ancestralGamma+0x330>
lea    (%r9,%r8,1),%eax
imul   %r10d,%eax
movsd  (%rdx,%rbp,8),%xmm2
add    %eax,%ebp
movslq %ebp,%rax
mulsd  (%rbx,%rax,8),%xmm2
addsd  %xmm2,%xmm1
jmp    1460 <ancestralGamma+0x330>
test   %r10d,%r10d
jle    1640 <ancestralGamma+0x510>
cmp    $0x1,%r10d
je     1530 <ancestralGamma+0x400>
lea    (%r12,%r15,8),%rax
cmp    %rax,0x18(%rsp)
jae    15ae <ancestralGamma+0x47e>
mov    0x10(%rsp),%rax
add    %r15,%rax
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
cmp    %rax,%r12
jae    15ae <ancestralGamma+0x47e>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
test   $0x1,%r15b
je     1556 <ancestralGamma+0x426>
movsd  (%r12,%rax,8),%xmm1
divsd  %xmm0,%xmm1
mov    0x18(%rsp),%rdx
movsd  %xmm1,(%rdx,%rax,8)
or     $0x1,%rax
add    %r15,%rcx
je     1640 <ancestralGamma+0x510>
mov    0x8(%rsp),%rcx
mov    0x10(%rsp),%rdx
lea    (%rcx,%rdx,8),%rcx
add    $0x8,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%r12,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%rcx,%rax,8)
movsd  0x8(%r12,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx,%rax,8)
add    $0x2,%rax
cmp    %rax,%r15
jne    1580 <ancestralGamma+0x450>
jmp    1640 <ancestralGamma+0x510>
mov    %r15d,%eax
and    $0xfffffffe,%eax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rax),%rdx
mov    %rdx,%rcx
shr    %rcx
add    $0x1,%rcx
test   %rdx,%rdx
je     161c <ancestralGamma+0x4ec>
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
mov    0x8(%rsp),%rdx
mov    0x10(%rsp),%rdi
lea    (%rdx,%rdi,8),%rdi
add    $0x10,%rdi
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movupd (%r12,%rdx,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rdi,%rdx,8)
movupd 0x10(%r12,%rdx,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdi,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    15f0 <ancestralGamma+0x4c0>
jmp    161e <ancestralGamma+0x4ee>
xor    %edx,%edx
test   $0x1,%cl
je     1637 <ancestralGamma+0x507>
movupd (%r12,%rdx,8),%xmm2
divpd  %xmm1,%xmm2
mov    0x18(%rsp),%rcx
movupd %xmm2,(%rcx,%rdx,8)
cmp    %r15,%rax
jne    1532 <ancestralGamma+0x402>
mov    %r12,%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <rax_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 i32
%16 = alloca i32
%17 = alloca double*
%18 = alloca double*
%19 = alloca double
%20 = alloca double*
%21 = 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
%22 = load i32, i32* %11
%23 = load i32, i32* %11
%24 = mul nsw i32 %22, %23
store i32 %24, i32* %13
%25 = load double*, double** %8
%26 = load i32, i32* %12
%27 = load i32, i32* %10
%28 = mul nsw i32 %26, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %25, i64 %29
store double* %30, double** %18
store double 0.000000e+00, double* %19
%31 = load i32, i32* %11
%32 = sext i32 %31 to i64
%33 = mul i64 8, %32
%34 = trunc i64 %33 to i32
%35 = call i64 @rax_malloc(i32 %34)
%36 = inttoptr i64 %35 to double*
store double* %36, double** %20
store i32 0, i32* %16
br label %37
37: ; preds = %46, %6
%38 = load i32, i32* %16
%39 = load i32, i32* %11
%40 = icmp slt i32 %38, %39
br i1 %40, label %41, label %49
41: ; preds = %37
%42 = load double*, double** %20
%43 = load i32, i32* %16
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %42, i64 %44
store double 0.000000e+00, double* %45
br label %46
46: ; preds = %41
%47 = load i32, i32* %16
%48 = add nsw i32 %47, 1
store i32 %48, i32* %16
br label %37
49: ; preds = %37
store i32 0, i32* %14
br label %50
50: ; preds = %108, %49
%51 = load i32, i32* %14
%52 = icmp slt i32 %51, 4
br i1 %52, label %53, label %111
53: ; preds = %50
%54 = load double*, double** %7
%55 = load i32, i32* %11
%56 = load i32, i32* %14
%57 = mul nsw i32 %55, %56
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %54, i64 %58
store double* %59, double** %17
store i32 0, i32* %16
br label %60
60: ; preds = %104, %53
%61 = load i32, i32* %16
%62 = load i32, i32* %11
%63 = icmp slt i32 %61, %62
br i1 %63, label %64, label %107
64: ; preds = %60
store double 0.000000e+00, double* %21
store i32 0, i32* %15
br label %65
65: ; preds = %90, %64
%66 = load i32, i32* %15
%67 = load i32, i32* %11
%68 = icmp slt i32 %66, %67
br i1 %68, label %69, label %93
69: ; preds = %65
%70 = load double*, double** %17
%71 = load i32, i32* %15
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %70, i64 %72
%74 = load double, double* %73
%75 = load double*, double** %9
%76 = load i32, i32* %14
%77 = load i32, i32* %13
%78 = mul nsw i32 %76, %77
%79 = load i32, i32* %16
%80 = load i32, i32* %11
%81 = mul nsw i32 %79, %80
%82 = add nsw i32 %78, %81
%83 = load i32, i32* %15
%84 = add nsw i32 %82, %83
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %75, i64 %85
%87 = load double, double* %86
%88 = load double, double* %21
%89 = call double @llvm.fmuladd.f64(double %74, double %87, double %88)
store double %89, double* %21
br label %90
90: ; preds = %69
%91 = load i32, i32* %15
%92 = add nsw i32 %91, 1
store i32 %92, i32* %15
br label %65
93: ; preds = %65
%94 = load double, double* %21
%95 = load double*, double** %20
%96 = load i32, i32* %16
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %95, i64 %97
%99 = load double, double* %98
%100 = fadd double %99, %94
store double %100, double* %98
%101 = load double, double* %21
%102 = load double, double* %19
%103 = fadd double %102, %101
store double %103, double* %19
br label %104
104: ; preds = %93
%105 = load i32, i32* %16
%106 = add nsw i32 %105, 1
store i32 %106, i32* %16
br label %60
107: ; preds = %60
br label %108
108: ; preds = %107
%109 = load i32, i32* %14
%110 = add nsw i32 %109, 1
store i32 %110, i32* %14
br label %50
111: ; preds = %50
store i32 0, i32* %16
br label %112
112: ; preds = %128, %111
%113 = load i32, i32* %16
%114 = load i32, i32* %11
%115 = icmp slt i32 %113, %114
br i1 %115, label %116, label %131
116: ; preds = %112
%117 = load double*, double** %20
%118 = load i32, i32* %16
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %117, i64 %119
%121 = load double, double* %120
%122 = load double, double* %19
%123 = fdiv double %121, %122
%124 = load double*, double** %18
%125 = load i32, i32* %16
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds double, double* %124, i64 %126
store double %123, double* %127
br label %128
128: ; preds = %116
%129 = load i32, i32* %16
%130 = add nsw i32 %129, 1
store i32 %130, i32* %16
br label %112
131: ; preds = %112
%132 = load double*, double** %20
%133 = call i32 @rax_free(double* %132)
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 = mul nsw i32 %5, %3
%8 = sext i32 %7 to i64
%9 = getelementptr double, double* %1, i64 %8
%10 = shl i32 %4, 3
%11 = tail call i64 @rax_malloc(i32 %10)
%12 = inttoptr i64 %11 to double*
%13 = inttoptr i64 %11 to i8*
%14 = icmp sgt i32 %4, 0
br i1 %14, label %83, label %339
15: ; preds = %325
%16 = zext i32 %4 to i64
%17 = icmp eq i32 %4, 1
br i1 %17, label %68, label %18
18: ; preds = %15
%19 = add nsw i64 %8, %88
%20 = getelementptr double, double* %1, i64 %19
%21 = bitcast double* %20 to i8*
%22 = getelementptr double, double* %12, i64 %88
%23 = icmp ult double* %9, %22
%24 = icmp ult i8* %13, %21
%25 = and i1 %23, %24
br i1 %25, label %68, label %26
26: ; preds = %18
%27 = and i64 %88, 4294967294
%28 = insertelement <2 x double> poison, double %322, i64 0
%29 = shufflevector <2 x double> %28, <2 x double> poison, <2 x i32> zeroinitializer
%30 = add nsw i64 %27, -2
%31 = lshr exact i64 %30, 1
%32 = add nuw i64 %31, 1
%33 = and i64 %32, 1
%34 = icmp eq i64 %30, 0
br i1 %34, label %56, label %35
35: ; preds = %26
%36 = and i64 %32, -2
br label %37
37: ; preds = %37, %35
%38 = phi i64 [ 0, %35 ], [ %53, %37 ]
%39 = phi i64 [ 0, %35 ], [ %54, %37 ]
%40 = getelementptr inbounds double, double* %12, i64 %38
%41 = bitcast double* %40 to <2 x double>*
%42 = load <2 x double>, <2 x double>* %41
%43 = fdiv <2 x double> %42, %29
%44 = getelementptr inbounds double, double* %9, i64 %38
%45 = bitcast double* %44 to <2 x double>*
store <2 x double> %43, <2 x double>* %45
%46 = or i64 %38, 2
%47 = getelementptr inbounds double, double* %12, i64 %46
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = fdiv <2 x double> %49, %29
%51 = getelementptr inbounds double, double* %9, i64 %46
%52 = bitcast double* %51 to <2 x double>*
store <2 x double> %50, <2 x double>* %52
%53 = add nuw i64 %38, 4
%54 = add nuw i64 %39, 2
%55 = icmp eq i64 %54, %36
br i1 %55, label %56, label %37
56: ; preds = %37, %26
%57 = phi i64 [ 0, %26 ], [ %53, %37 ]
%58 = icmp eq i64 %33, 0
br i1 %58, label %66, label %59
59: ; preds = %56
%60 = getelementptr inbounds double, double* %12, i64 %57
%61 = bitcast double* %60 to <2 x double>*
%62 = load <2 x double>, <2 x double>* %61
%63 = fdiv <2 x double> %62, %29
%64 = getelementptr inbounds double, double* %9, i64 %57
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> %63, <2 x double>* %65
br label %66
66: ; preds = %56, %59
%67 = icmp eq i64 %27, %88
br i1 %67, label %339, label %68
68: ; preds = %18, %15, %66
%69 = phi i64 [ 0, %18 ], [ 0, %15 ], [ %27, %66 ]
%70 = xor i64 %69, -1
%71 = and i64 %88, 1
%72 = icmp eq i64 %71, 0
br i1 %72, label %79, label %73
73: ; preds = %68
%74 = getelementptr inbounds double, double* %12, i64 %69
%75 = load double, double* %74
%76 = fdiv double %75, %322
%77 = getelementptr inbounds double, double* %9, i64 %69
store double %76, double* %77
%78 = or i64 %69, 1
br label %79
79: ; preds = %73, %68
%80 = phi i64 [ %69, %68 ], [ %78, %73 ]
%81 = sub nsw i64 0, %88
%82 = icmp eq i64 %70, %81
br i1 %82, label %339, label %326
83: ; preds = %6
%84 = inttoptr i64 %11 to i8*
%85 = zext i32 %4 to i64
%86 = shl nuw nsw i64 %85, 3
%87 = zext i32 %4 to i64
%88 = zext i32 %4 to i64
%89 = add nsw i64 %88, -1
%90 = and i64 %88, 1
%91 = icmp eq i64 %89, 0
%92 = and i64 %88, 4294967294
%93 = icmp eq i64 %90, 0
br label %94
94: ; preds = %136, %83
%95 = phi i64 [ 0, %83 ], [ %142, %136 ]
%96 = phi double [ 0.000000e+00, %83 ], [ %141, %136 ]
%97 = trunc i64 %95 to i32
%98 = mul i32 %97, %4
br i1 %91, label %123, label %99
99: ; preds = %94, %99
%100 = phi i64 [ %120, %99 ], [ 0, %94 ]
%101 = phi double [ %119, %99 ], [ 0.000000e+00, %94 ]
%102 = phi i64 [ %121, %99 ], [ 0, %94 ]
%103 = getelementptr inbounds double, double* %0, i64 %100
%104 = load double, double* %103
%105 = trunc i64 %100 to i32
%106 = add i32 %98, %105
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds double, double* %2, i64 %107
%109 = load double, double* %108
%110 = tail call double @llvm.fmuladd.f64(double %104, double %109, double %101)
%111 = or i64 %100, 1
%112 = getelementptr inbounds double, double* %0, i64 %111
%113 = load double, double* %112
%114 = trunc i64 %111 to i32
%115 = add i32 %98, %114
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %2, i64 %116
%118 = load double, double* %117
%119 = tail call double @llvm.fmuladd.f64(double %113, double %118, double %110)
%120 = add nuw nsw i64 %100, 2
%121 = add i64 %102, 2
%122 = icmp eq i64 %121, %92
br i1 %122, label %123, label %99
123: ; preds = %99, %94
%124 = phi double [ undef, %94 ], [ %119, %99 ]
%125 = phi i64 [ 0, %94 ], [ %120, %99 ]
%126 = phi double [ 0.000000e+00, %94 ], [ %119, %99 ]
br i1 %93, label %136, label %127
127: ; preds = %123
%128 = getelementptr inbounds double, double* %0, i64 %125
%129 = load double, double* %128
%130 = trunc i64 %125 to i32
%131 = add i32 %98, %130
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %2, i64 %132
%134 = load double, double* %133
%135 = tail call double @llvm.fmuladd.f64(double %129, double %134, double %126)
br label %136
136: ; preds = %123, %127
%137 = phi double [ %124, %123 ], [ %135, %127 ]
%138 = getelementptr inbounds double, double* %12, i64 %95
%139 = load double, double* %138
%140 = fadd double %137, %139
store double %140, double* %138
%141 = fadd double %96, %137
%142 = add nuw nsw i64 %95, 1
%143 = icmp eq i64 %142, %88
br i1 %143, label %144, label %94
144: ; preds = %136
%145 = getelementptr inbounds double, double* %0, i64 %87
br i1 %14, label %146, label %339
146: ; preds = %144
%147 = zext i32 %4 to i64
%148 = and i64 %88, 1
%149 = icmp eq i64 %89, 0
%150 = and i64 %88, 4294967294
%151 = icmp eq i64 %148, 0
br label %152
152: ; preds = %146, %195
%153 = phi i64 [ 0, %146 ], [ %201, %195 ]
%154 = phi double [ %141, %146 ], [ %200, %195 ]
%155 = trunc i64 %153 to i32
%156 = add i32 %155, %4
%157 = mul i32 %156, %4
br i1 %149, label %182, label %158
158: ; preds = %152, %158
%159 = phi i64 [ %179, %158 ], [ 0, %152 ]
%160 = phi double [ %178, %158 ], [ 0.000000e+00, %152 ]
%161 = phi i64 [ %180, %158 ], [ 0, %152 ]
%162 = getelementptr inbounds double, double* %145, i64 %159
%163 = load double, double* %162
%164 = trunc i64 %159 to i32
%165 = add i32 %157, %164
%166 = sext i32 %165 to i64
%167 = getelementptr inbounds double, double* %2, i64 %166
%168 = load double, double* %167
%169 = tail call double @llvm.fmuladd.f64(double %163, double %168, double %160)
%170 = or i64 %159, 1
%171 = getelementptr inbounds double, double* %145, i64 %170
%172 = load double, double* %171
%173 = trunc i64 %170 to i32
%174 = add i32 %157, %173
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds double, double* %2, i64 %175
%177 = load double, double* %176
%178 = tail call double @llvm.fmuladd.f64(double %172, double %177, double %169)
%179 = add nuw nsw i64 %159, 2
%180 = add i64 %161, 2
%181 = icmp eq i64 %180, %150
br i1 %181, label %182, label %158
182: ; preds = %158, %152
%183 = phi double [ undef, %152 ], [ %178, %158 ]
%184 = phi i64 [ 0, %152 ], [ %179, %158 ]
%185 = phi double [ 0.000000e+00, %152 ], [ %178, %158 ]
br i1 %151, label %195, label %186
186: ; preds = %182
%187 = getelementptr inbounds double, double* %145, i64 %184
%188 = load double, double* %187
%189 = trunc i64 %184 to i32
%190 = add i32 %157, %189
%191 = sext i32 %190 to i64
%192 = getelementptr inbounds double, double* %2, i64 %191
%193 = load double, double* %192
%194 = tail call double @llvm.fmuladd.f64(double %188, double %193, double %185)
br label %195
195: ; preds = %182, %186
%196 = phi double [ %183, %182 ], [ %194, %186 ]
%197 = getelementptr inbounds double, double* %12, i64 %153
%198 = load double, double* %197
%199 = fadd double %196, %198
store double %199, double* %197
%200 = fadd double %154, %196
%201 = add nuw nsw i64 %153, 1
%202 = icmp eq i64 %201, %147
br i1 %202, label %203, label %152
203: ; preds = %195
%204 = shl nuw nsw i64 %87, 1
%205 = getelementptr inbounds double, double* %0, i64 %204
br i1 %14, label %206, label %339
206: ; preds = %203
%207 = shl nuw i32 %4, 1
%208 = zext i32 %4 to i64
%209 = and i64 %88, 1
%210 = icmp eq i64 %89, 0
%211 = and i64 %88, 4294967294
%212 = icmp eq i64 %209, 0
br label %213
213: ; preds = %206, %256
%214 = phi i64 [ 0, %206 ], [ %262, %256 ]
%215 = phi double [ %200, %206 ], [ %261, %256 ]
%216 = trunc i64 %214 to i32
%217 = add i32 %207, %216
%218 = mul i32 %217, %4
br i1 %210, label %243, label %219
219: ; preds = %213, %219
%220 = phi i64 [ %240, %219 ], [ 0, %213 ]
%221 = phi double [ %239, %219 ], [ 0.000000e+00, %213 ]
%222 = phi i64 [ %241, %219 ], [ 0, %213 ]
%223 = getelementptr inbounds double, double* %205, i64 %220
%224 = load double, double* %223
%225 = trunc i64 %220 to i32
%226 = add i32 %218, %225
%227 = sext i32 %226 to i64
%228 = getelementptr inbounds double, double* %2, i64 %227
%229 = load double, double* %228
%230 = tail call double @llvm.fmuladd.f64(double %224, double %229, double %221)
%231 = or i64 %220, 1
%232 = getelementptr inbounds double, double* %205, i64 %231
%233 = load double, double* %232
%234 = trunc i64 %231 to i32
%235 = add i32 %218, %234
%236 = sext i32 %235 to i64
%237 = getelementptr inbounds double, double* %2, i64 %236
%238 = load double, double* %237
%239 = tail call double @llvm.fmuladd.f64(double %233, double %238, double %230)
%240 = add nuw nsw i64 %220, 2
%241 = add i64 %222, 2
%242 = icmp eq i64 %241, %211
br i1 %242, label %243, label %219
243: ; preds = %219, %213
%244 = phi double [ undef, %213 ], [ %239, %219 ]
%245 = phi i64 [ 0, %213 ], [ %240, %219 ]
%246 = phi double [ 0.000000e+00, %213 ], [ %239, %219 ]
br i1 %212, label %256, label %247
247: ; preds = %243
%248 = getelementptr inbounds double, double* %205, i64 %245
%249 = load double, double* %248
%250 = trunc i64 %245 to i32
%251 = add i32 %218, %250
%252 = sext i32 %251 to i64
%253 = getelementptr inbounds double, double* %2, i64 %252
%254 = load double, double* %253
%255 = tail call double @llvm.fmuladd.f64(double %249, double %254, double %246)
br label %256
256: ; preds = %243, %247
%257 = phi double [ %244, %243 ], [ %255, %247 ]
%258 = getelementptr inbounds double, double* %12, i64 %214
%259 = load double, double* %258
%260 = fadd double %257, %259
store double %260, double* %258
%261 = fadd double %215, %257
%262 = add nuw nsw i64 %214, 1
%263 = icmp eq i64 %262, %208
br i1 %263, label %264, label %213
264: ; preds = %256
%265 = mul nuw nsw i64 %87, 3
%266 = getelementptr inbounds double, double* %0, i64 %265
br i1 %14, label %267, label %339
267: ; preds = %264
%268 = mul i32 %4, 3
%269 = zext i32 %4 to i64
%270 = and i64 %88, 1
%271 = icmp eq i64 %89, 0
%272 = and i64 %88, 4294967294
%273 = icmp eq i64 %270, 0
br label %274
274: ; preds = %267, %317
%275 = phi i64 [ 0, %267 ], [ %323, %317 ]
%276 = phi double [ %261, %267 ], [ %322, %317 ]
%277 = trunc i64 %275 to i32
%278 = add i32 %268, %277
%279 = mul i32 %278, %4
br i1 %271, label %304, label %280
280: ; preds = %274, %280
%281 = phi i64 [ %301, %280 ], [ 0, %274 ]
%282 = phi double [ %300, %280 ], [ 0.000000e+00, %274 ]
%283 = phi i64 [ %302, %280 ], [ 0, %274 ]
%284 = getelementptr inbounds double, double* %266, i64 %281
%285 = load double, double* %284
%286 = trunc i64 %281 to i32
%287 = add i32 %279, %286
%288 = sext i32 %287 to i64
%289 = getelementptr inbounds double, double* %2, i64 %288
%290 = load double, double* %289
%291 = tail call double @llvm.fmuladd.f64(double %285, double %290, double %282)
%292 = or i64 %281, 1
%293 = getelementptr inbounds double, double* %266, i64 %292
%294 = load double, double* %293
%295 = trunc i64 %292 to i32
%296 = add i32 %279, %295
%297 = sext i32 %296 to i64
%298 = getelementptr inbounds double, double* %2, i64 %297
%299 = load double, double* %298
%300 = tail call double @llvm.fmuladd.f64(double %294, double %299, double %291)
%301 = add nuw nsw i64 %281, 2
%302 = add i64 %283, 2
%303 = icmp eq i64 %302, %272
br i1 %303, label %304, label %280
304: ; preds = %280, %274
%305 = phi double [ undef, %274 ], [ %300, %280 ]
%306 = phi i64 [ 0, %274 ], [ %301, %280 ]
%307 = phi double [ 0.000000e+00, %274 ], [ %300, %280 ]
br i1 %273, label %317, label %308
308: ; preds = %304
%309 = getelementptr inbounds double, double* %266, i64 %306
%310 = load double, double* %309
%311 = trunc i64 %306 to i32
%312 = add i32 %279, %311
%313 = sext i32 %312 to i64
%314 = getelementptr inbounds double, double* %2, i64 %313
%315 = load double, double* %314
%316 = tail call double @llvm.fmuladd.f64(double %310, double %315, double %307)
br label %317
317: ; preds = %304, %308
%318 = phi double [ %305, %304 ], [ %316, %308 ]
%319 = getelementptr inbounds double, double* %12, i64 %275
%320 = load double, double* %319
%321 = fadd double %318, %320
store double %321, double* %319
%322 = fadd double %276, %318
%323 = add nuw nsw i64 %275, 1
%324 = icmp eq i64 %323, %269
br i1 %324, label %325, label %274
325: ; preds = %317
br i1 %14, label %15, label %339
326: ; preds = %79, %326
%327 = phi i64 [ %337, %326 ], [ %80, %79 ]
%328 = getelementptr inbounds double, double* %12, i64 %327
%329 = load double, double* %328
%330 = fdiv double %329, %322
%331 = getelementptr inbounds double, double* %9, i64 %327
store double %330, double* %331
%332 = add nuw nsw i64 %327, 1
%333 = getelementptr inbounds double, double* %12, i64 %332
%334 = load double, double* %333
%335 = fdiv double %334, %322
%336 = getelementptr inbounds double, double* %9, i64 %332
store double %335, double* %336
%337 = add nuw nsw i64 %327, 2
%338 = icmp eq i64 %337, %16
br i1 %338, label %339, label %326
339: ; preds = %79, %326, %66, %6, %144, %203, %264, %325
%340 = tail call i32 @rax_free(double* %12)
ret void
}