double lr(table, rows, cols)
int **table;
int rows, cols;
{
  double sum;
  int *rowsums, *colsums;
  int total;
  int i, j;
  total = 0;
  rowsums = (int *)calloc(rows, sizeof(rowsums[0]));
  colsums = (int *)calloc(cols, sizeof(colsums[0]));
  for (j = 0; j < cols; j++) {
    colsums[j] = 0;
  }
  for (i = 0; i < rows; i++) {
    rowsums[i] = 0;
    for (j = 0; j < cols; j++) {
      rowsums[i] += table[i][j];
      colsums[j] += table[i][j];
    }
    total += rowsums[i];
  }
  sum = 0;
  for (i = 0; i < rows; i++) {
    sum += klogP(table[i], rowsums[i], table[i], cols) -
           klogP(colsums, total, table[i], cols);
  }
  free(rowsums);
  free(colsums);
  return 2 * sum;
}


x86 assembly (clang, -O0):
<lr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x0,-0x2c(%rbp)
movslq -0xc(%rbp),%rdi
mov    $0x4,%esi
call   1050 <calloc@plt>
mov    %rax,-0x20(%rbp)
movslq -0x10(%rbp),%rdi
mov    $0x4,%esi
call   1050 <calloc@plt>
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    119d <lr+0x6d>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1174 <lr+0x44>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1241 <lr+0x111>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1222 <lr+0xf2>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
add    (%rax,%rcx,4),%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
add    (%rax,%rcx,4),%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11c6 <lr+0x96>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11a4 <lr+0x74>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12cc <lr+0x19c>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%esi
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x10(%rbp),%ecx
mov    $0x0,%al
call   1040 <klogP@plt>
mov    %eax,-0x38(%rbp)
mov    -0x28(%rbp),%rdi
mov    -0x2c(%rbp),%esi
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x10(%rbp),%ecx
mov    $0x0,%al
call   1040 <klogP@plt>
mov    %eax,%ecx
mov    -0x38(%rbp),%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1250 <lr+0x120>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
movsd  0xd1a(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lr>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %edx,%r12d
mov    %esi,%r14d
mov    %rdi,0x10(%rsp)
movslq %esi,%rbp
mov    $0x4,%esi
mov    %rbp,%rdi
call   1050 <calloc@plt>
mov    %rax,%r15
movslq %r12d,%rdi
mov    $0x4,%esi
call   1050 <calloc@plt>
mov    %rax,%rbx
xorpd  %xmm8,%xmm8
test   %ebp,%ebp
jle    13f1 <lr+0x2c1>
mov    %r14d,0x1c(%rsp)
mov    %r14d,%ecx
mov    %r12d,%eax
lea    (%rbx,%rax,4),%rdx
mov    %rdx,0x8(%rsp)
mov    %eax,%r9d
and    $0xfffffff8,%r9d
lea    -0x8(%r9),%r11
mov    %r11,0x30(%rsp)
shr    $0x3,%r11
add    $0x1,%r11
mov    %r11,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x28(%rsp)
mov    %rax,%r14
neg    %r14
xor    %r10d,%r10d
xor    %r13d,%r13d
mov    0x10(%rsp),%r8
mov    %rcx,0x20(%rsp)
jmp    11e9 <lr+0xb9>
nopl   0x0(%rax)
mov    %edi,(%r15,%r10,4)
mov    0x20(%rsp),%rcx
add    %edi,%r13d
add    $0x1,%r10
cmp    %rcx,%r10
je     1372 <lr+0x242>
movl   $0x0,(%r15,%r10,4)
mov    $0x0,%edi
test   %r12d,%r12d
jle    11d9 <lr+0xa9>
mov    (%r8,%r10,8),%rcx
cmp    $0x8,%r12d
jb     1215 <lr+0xe5>
lea    (%rcx,%rax,4),%rdx
cmp    %rdx,%rbx
jae    1270 <lr+0x140>
cmp    %rcx,0x8(%rsp)
jbe    1270 <lr+0x140>
xor    %ebp,%ebp
xor    %edi,%edi
test   $0x1,%al
jne    1230 <lr+0x100>
mov    %rbp,%rdx
jmp    123f <lr+0x10f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%rcx,%rbp,4),%edx
add    %edx,%edi
add    %edx,(%rbx,%rbp,4)
mov    %rbp,%rdx
or     $0x1,%rdx
not    %rbp
cmp    %r14,%rbp
je     11d0 <lr+0xa0>
nopw   0x0(%rax,%rax,1)
mov    (%rcx,%rdx,4),%ebp
add    %ebp,(%rbx,%rdx,4)
mov    0x4(%rcx,%rdx,4),%esi
add    %esi,0x4(%rbx,%rdx,4)
add    %ebp,%edi
add    %esi,%edi
add    $0x2,%rdx
cmp    %rdx,%rax
jne    1250 <lr+0x120>
jmp    11d0 <lr+0xa0>
cmpq   $0x0,0x30(%rsp)
je     1301 <lr+0x1d1>
pxor   %xmm0,%xmm0
mov    0x28(%rsp),%r8
xor    %edi,%edi
pxor   %xmm1,%xmm1
nopl   0x0(%rax,%rax,1)
movdqu (%rcx,%rdi,4),%xmm2
movdqu 0x10(%rcx,%rdi,4),%xmm3
paddd  %xmm2,%xmm0
paddd  %xmm3,%xmm1
movdqu (%rbx,%rdi,4),%xmm4
movdqu 0x10(%rbx,%rdi,4),%xmm5
movdqu 0x20(%rbx,%rdi,4),%xmm6
movdqu 0x30(%rbx,%rdi,4),%xmm7
paddd  %xmm2,%xmm4
paddd  %xmm3,%xmm5
movdqu %xmm4,(%rbx,%rdi,4)
movdqu %xmm5,0x10(%rbx,%rdi,4)
movdqu 0x20(%rcx,%rdi,4),%xmm2
movdqu 0x30(%rcx,%rdi,4),%xmm3
paddd  %xmm2,%xmm0
paddd  %xmm3,%xmm1
paddd  %xmm2,%xmm6
paddd  %xmm3,%xmm7
movdqu %xmm6,0x20(%rbx,%rdi,4)
movdqu %xmm7,0x30(%rbx,%rdi,4)
add    $0x10,%rdi
add    $0xfffffffffffffffe,%r8
jne    1290 <lr+0x160>
jmp    130b <lr+0x1db>
pxor   %xmm0,%xmm0
xor    %edi,%edi
pxor   %xmm1,%xmm1
test   $0x1,%r11b
mov    0x10(%rsp),%r8
je     1347 <lr+0x217>
movdqu (%rcx,%rdi,4),%xmm2
movdqu 0x10(%rcx,%rdi,4),%xmm3
paddd  %xmm2,%xmm0
paddd  %xmm3,%xmm1
movdqu (%rbx,%rdi,4),%xmm4
movdqu 0x10(%rbx,%rdi,4),%xmm5
paddd  %xmm2,%xmm4
paddd  %xmm3,%xmm5
movdqu %xmm4,(%rbx,%rdi,4)
movdqu %xmm5,0x10(%rbx,%rdi,4)
paddd  %xmm1,%xmm0
pshufd $0xee,%xmm0,%xmm1
paddd  %xmm0,%xmm1
pshufd $0x55,%xmm1,%xmm0
paddd  %xmm1,%xmm0
movd   %xmm0,%edi
mov    %r9,%rbp
cmp    %rax,%r9
je     11d0 <lr+0xa0>
jmp    1219 <lr+0xe9>
cmpl   $0x0,0x1c(%rsp)
jle    13f1 <lr+0x2c1>
xorpd  %xmm8,%xmm8
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  %xmm8,0x8(%rsp)
mov    (%r8,%r14,8),%rdi
mov    (%r15,%r14,4),%esi
mov    %rdi,%rdx
mov    %r12d,%ecx
xor    %eax,%eax
call   1040 <klogP@plt>
mov    %eax,%ebp
mov    0x10(%rsp),%rax
mov    (%rax,%r14,8),%rdx
mov    %rbx,%rdi
mov    %r13d,%esi
mov    %r12d,%ecx
xor    %eax,%eax
call   1040 <klogP@plt>
movsd  0x8(%rsp),%xmm8
mov    0x10(%rsp),%r8
sub    %eax,%ebp
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
addsd  %xmm0,%xmm8
add    $0x1,%r14
cmp    %r14,0x20(%rsp)
jne    1390 <lr+0x260>
addsd  %xmm8,%xmm8
movsd  %xmm8,0x8(%rsp)
mov    %r15,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
movsd  0x8(%rsp),%xmm0
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @lr(i32** nocapture %0, i32 %1, i32 %2) {
%4 = sext i32 %1 to i64
%5 = tail call i8* @calloc(i64 %4, i64 4)
%6 = bitcast i8* %5 to i32*
%7 = sext i32 %2 to i64
%8 = tail call i8* @calloc(i64 %7, i64 4)
%9 = bitcast i8* %8 to i32*
%10 = icmp sgt i32 %2, 0
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %183
12: ; preds = %3
%13 = zext i32 %1 to i64
%14 = zext i32 %2 to i64
%15 = shl nuw nsw i64 %14, 2
%16 = getelementptr i8, i8* %8, i64 %15
%17 = and i64 %14, 4294967288
%18 = add nsw i64 %17, -8
%19 = lshr exact i64 %18, 3
%20 = add nuw nsw i64 %19, 1
%21 = icmp ult i32 %2, 8
%22 = and i64 %14, 4294967288
%23 = and i64 %20, 1
%24 = icmp eq i64 %18, 0
%25 = and i64 %20, 4611686018427387902
%26 = icmp eq i64 %23, 0
%27 = icmp eq i64 %22, %14
%28 = and i64 %14, 1
%29 = icmp eq i64 %28, 0
%30 = sub nsw i64 0, %14
br label %34
31: ; preds = %161
br i1 %11, label %32, label %183
32: ; preds = %31
%33 = zext i32 %1 to i64
br label %166
34: ; preds = %12, %161
%35 = phi i64 [ 0, %12 ], [ %164, %161 ]
%36 = phi i32 [ 0, %12 ], [ %163, %161 ]
%37 = getelementptr inbounds i32, i32* %6, i64 %35
store i32 0, i32* %37
br i1 %10, label %38, label %161
38: ; preds = %34
%39 = getelementptr inbounds i32*, i32** %0, i64 %35
%40 = load i32*, i32** %39
br i1 %21, label %124, label %41
41: ; preds = %38
%42 = bitcast i32* %40 to i8*
%43 = getelementptr i32, i32* %40, i64 %14
%44 = bitcast i32* %43 to i8*
%45 = icmp ult i8* %8, %44
%46 = icmp ugt i8* %16, %42
%47 = and i1 %45, %46
br i1 %47, label %124, label %48
48: ; preds = %41
br i1 %24, label %94, label %49
49: ; preds = %48, %49
%50 = phi i64 [ %91, %49 ], [ 0, %48 ]
%51 = phi <4 x i32> [ %79, %49 ], [ zeroinitializer, %48 ]
%52 = phi <4 x i32> [ %80, %49 ], [ zeroinitializer, %48 ]
%53 = phi i64 [ %92, %49 ], [ 0, %48 ]
%54 = getelementptr inbounds i32, i32* %40, i64 %50
%55 = bitcast i32* %54 to <4 x i32>*
%56 = load <4 x i32>, <4 x i32>* %55
%57 = getelementptr inbounds i32, i32* %54, i64 4
%58 = bitcast i32* %57 to <4 x i32>*
%59 = load <4 x i32>, <4 x i32>* %58
%60 = add <4 x i32> %51, %56
%61 = add <4 x i32> %52, %59
%62 = getelementptr inbounds i32, i32* %9, i64 %50
%63 = bitcast i32* %62 to <4 x i32>*
%64 = load <4 x i32>, <4 x i32>* %63
%65 = getelementptr inbounds i32, i32* %62, i64 4
%66 = bitcast i32* %65 to <4 x i32>*
%67 = load <4 x i32>, <4 x i32>* %66
%68 = add nsw <4 x i32> %64, %56
%69 = add nsw <4 x i32> %67, %59
%70 = bitcast i32* %62 to <4 x i32>*
store <4 x i32> %68, <4 x i32>* %70
%71 = bitcast i32* %65 to <4 x i32>*
store <4 x i32> %69, <4 x i32>* %71
%72 = or i64 %50, 8
%73 = getelementptr inbounds i32, i32* %40, i64 %72
%74 = bitcast i32* %73 to <4 x i32>*
%75 = load <4 x i32>, <4 x i32>* %74
%76 = getelementptr inbounds i32, i32* %73, i64 4
%77 = bitcast i32* %76 to <4 x i32>*
%78 = load <4 x i32>, <4 x i32>* %77
%79 = add <4 x i32> %60, %75
%80 = add <4 x i32> %61, %78
%81 = getelementptr inbounds i32, i32* %9, i64 %72
%82 = bitcast i32* %81 to <4 x i32>*
%83 = load <4 x i32>, <4 x i32>* %82
%84 = getelementptr inbounds i32, i32* %81, i64 4
%85 = bitcast i32* %84 to <4 x i32>*
%86 = load <4 x i32>, <4 x i32>* %85
%87 = add nsw <4 x i32> %83, %75
%88 = add nsw <4 x i32> %86, %78
%89 = bitcast i32* %81 to <4 x i32>*
store <4 x i32> %87, <4 x i32>* %89
%90 = bitcast i32* %84 to <4 x i32>*
store <4 x i32> %88, <4 x i32>* %90
%91 = add nuw i64 %50, 16
%92 = add i64 %53, 2
%93 = icmp eq i64 %92, %25
br i1 %93, label %94, label %49
94: ; preds = %49, %48
%95 = phi <4 x i32> [ undef, %48 ], [ %79, %49 ]
%96 = phi <4 x i32> [ undef, %48 ], [ %80, %49 ]
%97 = phi i64 [ 0, %48 ], [ %91, %49 ]
%98 = phi <4 x i32> [ zeroinitializer, %48 ], [ %79, %49 ]
%99 = phi <4 x i32> [ zeroinitializer, %48 ], [ %80, %49 ]
br i1 %26, label %119, label %100
100: ; preds = %94
%101 = getelementptr inbounds i32, i32* %40, i64 %97
%102 = bitcast i32* %101 to <4 x i32>*
%103 = load <4 x i32>, <4 x i32>* %102
%104 = getelementptr inbounds i32, i32* %101, i64 4
%105 = bitcast i32* %104 to <4 x i32>*
%106 = load <4 x i32>, <4 x i32>* %105
%107 = add <4 x i32> %98, %103
%108 = add <4 x i32> %99, %106
%109 = getelementptr inbounds i32, i32* %9, i64 %97
%110 = bitcast i32* %109 to <4 x i32>*
%111 = load <4 x i32>, <4 x i32>* %110
%112 = getelementptr inbounds i32, i32* %109, i64 4
%113 = bitcast i32* %112 to <4 x i32>*
%114 = load <4 x i32>, <4 x i32>* %113
%115 = add nsw <4 x i32> %111, %103
%116 = add nsw <4 x i32> %114, %106
%117 = bitcast i32* %109 to <4 x i32>*
store <4 x i32> %115, <4 x i32>* %117
%118 = bitcast i32* %112 to <4 x i32>*
store <4 x i32> %116, <4 x i32>* %118
br label %119
119: ; preds = %94, %100
%120 = phi <4 x i32> [ %95, %94 ], [ %107, %100 ]
%121 = phi <4 x i32> [ %96, %94 ], [ %108, %100 ]
%122 = add <4 x i32> %121, %120
%123 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %122)
br i1 %27, label %159, label %124
124: ; preds = %41, %38, %119
%125 = phi i64 [ 0, %41 ], [ 0, %38 ], [ %22, %119 ]
%126 = phi i32 [ 0, %41 ], [ 0, %38 ], [ %123, %119 ]
%127 = xor i64 %125, -1
br i1 %29, label %136, label %128
128: ; preds = %124
%129 = getelementptr inbounds i32, i32* %40, i64 %125
%130 = load i32, i32* %129
%131 = add nsw i32 %126, %130
%132 = getelementptr inbounds i32, i32* %9, i64 %125
%133 = load i32, i32* %132
%134 = add nsw i32 %133, %130
store i32 %134, i32* %132
%135 = or i64 %125, 1
br label %136
136: ; preds = %128, %124
%137 = phi i32 [ undef, %124 ], [ %131, %128 ]
%138 = phi i64 [ %125, %124 ], [ %135, %128 ]
%139 = phi i32 [ %126, %124 ], [ %131, %128 ]
%140 = icmp eq i64 %127, %30
br i1 %140, label %159, label %141
141: ; preds = %136, %141
%142 = phi i64 [ %157, %141 ], [ %138, %136 ]
%143 = phi i32 [ %153, %141 ], [ %139, %136 ]
%144 = getelementptr inbounds i32, i32* %40, i64 %142
%145 = load i32, i32* %144
%146 = add nsw i32 %143, %145
%147 = getelementptr inbounds i32, i32* %9, i64 %142
%148 = load i32, i32* %147
%149 = add nsw i32 %148, %145
store i32 %149, i32* %147
%150 = add nuw nsw i64 %142, 1
%151 = getelementptr inbounds i32, i32* %40, i64 %150
%152 = load i32, i32* %151
%153 = add nsw i32 %146, %152
%154 = getelementptr inbounds i32, i32* %9, i64 %150
%155 = load i32, i32* %154
%156 = add nsw i32 %155, %152
store i32 %156, i32* %154
%157 = add nuw nsw i64 %142, 2
%158 = icmp eq i64 %157, %14
br i1 %158, label %159, label %141
159: ; preds = %136, %141, %119
%160 = phi i32 [ %123, %119 ], [ %137, %136 ], [ %153, %141 ]
store i32 %160, i32* %37
br label %161
161: ; preds = %159, %34
%162 = phi i32 [ %160, %159 ], [ 0, %34 ]
%163 = add nsw i32 %162, %36
%164 = add nuw nsw i64 %35, 1
%165 = icmp eq i64 %164, %13
br i1 %165, label %31, label %34
166: ; preds = %32, %166
%167 = phi i64 [ 0, %32 ], [ %179, %166 ]
%168 = phi double [ 0.000000e+00, %32 ], [ %178, %166 ]
%169 = getelementptr inbounds i32*, i32** %0, i64 %167
%170 = load i32*, i32** %169
%171 = getelementptr inbounds i32, i32* %6, i64 %167
%172 = load i32, i32* %171
%173 = tail call i32 (i32*, i32, i32*, i32, ...) bitcast (i32 (...)* @klogP to i32 (i32*, i32, i32*, i32, ...)*)(i32* %170, i32 %172, i32* %170, i32 %2)
%174 = load i32*, i32** %169
%175 = tail call i32 (i32*, i32, i32*, i32, ...) bitcast (i32 (...)* @klogP to i32 (i32*, i32, i32*, i32, ...)*)(i32* %9, i32 %163, i32* %174, i32 %2)
%176 = sub nsw i32 %173, %175
%177 = sitofp i32 %176 to double
%178 = fadd double %168, %177
%179 = add nuw nsw i64 %167, 1
%180 = icmp eq i64 %179, %33
br i1 %180, label %181, label %166
181: ; preds = %166
%182 = fmul double %178, 2.000000e+00
br label %183
183: ; preds = %3, %181, %31
%184 = phi double [ 0.000000e+00, %31 ], [ %182, %181 ], [ 0.000000e+00, %3 ]
tail call void @free(i8* %5)
tail call void @free(i8* %8)
ret double %184
}