double mutual(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 = -entropy_term(rowsums, total, cols, cols) -
        entropy_term(colsums, total, rows, rows);
  for (i = 0; i < rows; i++) {
    sum += entropy_term(table[i], total, cols, rows * cols);
  }
  free(rowsums);
  free(colsums);
  return sum;
}


x86 assembly (clang, -O0):
<mutual>:
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   1040 <calloc@plt>
mov    %rax,-0x20(%rbp)
movslq -0x10(%rbp),%rdi
mov    $0x4,%esi
call   1040 <calloc@plt>
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    119d <mutual+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 <mutual+0x44>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1241 <mutual+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 <mutual+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 <mutual+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 <mutual+0x74>
mov    -0x20(%rbp),%rdi
mov    -0x2c(%rbp),%esi
mov    -0x10(%rbp),%edx
mov    -0x10(%rbp),%ecx
mov    $0x0,%al
call   1050 <entropy_term@plt>
mov    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x38(%rbp)
mov    -0x28(%rbp),%rdi
mov    -0x2c(%rbp),%esi
mov    -0xc(%rbp),%edx
mov    -0xc(%rbp),%ecx
mov    $0x0,%al
call   1050 <entropy_term@plt>
mov    %eax,%ecx
mov    -0x38(%rbp),%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12d1 <mutual+0x1a1>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x2c(%rbp),%esi
mov    -0x10(%rbp),%edx
mov    -0xc(%rbp),%ecx
imul   -0x10(%rbp),%ecx
mov    $0x0,%al
call   1050 <entropy_term@plt>
cvtsi2sd %eax,%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1289 <mutual+0x159>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
movsd  -0x18(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mutual>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %edx,%r12d
mov    %esi,%r15d
mov    %rdi,%r14
movslq %esi,%r13
mov    $0x4,%esi
mov    %r13,%rdi
call   1040 <calloc@plt>
mov    %rax,%rbp
movslq %r12d,%rdi
mov    $0x4,%esi
call   1040 <calloc@plt>
mov    %rax,%rbx
mov    %r15d,%r10d
test   %r13d,%r13d
jle    11fb <mutual+0xcb>
mov    %r12d,%r8d
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
xor    %edx,%edx
xor    %r13d,%r13d
jmp    119c <mutual+0x6c>
nopw   0x0(%rax,%rax,1)
add    %esi,%r13d
add    $0x1,%rdx
cmp    %r10,%rdx
je     11fe <mutual+0xce>
movl   $0x0,0x0(%rbp,%rdx,4)
mov    $0x0,%esi
test   %r12d,%r12d
jle    1190 <mutual+0x60>
mov    (%r14,%rdx,8),%rdi
xor    %esi,%esi
xor    %eax,%eax
cmp    $0x1,%r12d
je     11e6 <mutual+0xb6>
nopl   0x0(%rax)
add    (%rdi,%rax,4),%esi
mov    %esi,0x0(%rbp,%rdx,4)
mov    (%rdi,%rax,4),%ecx
add    %ecx,(%rbx,%rax,4)
add    0x4(%rdi,%rax,4),%esi
mov    %esi,0x0(%rbp,%rdx,4)
mov    0x4(%rdi,%rax,4),%ecx
add    %ecx,0x4(%rbx,%rax,4)
add    $0x2,%rax
cmp    %rax,%r9
jne    11c0 <mutual+0x90>
test   $0x1,%r8b
je     1190 <mutual+0x60>
add    (%rdi,%rax,4),%esi
mov    %esi,0x0(%rbp,%rdx,4)
mov    (%rdi,%rax,4),%ecx
add    %ecx,(%rbx,%rax,4)
jmp    1190 <mutual+0x60>
xor    %r13d,%r13d
mov    %r10,0x18(%rsp)
mov    %rbp,%rdi
mov    %r13d,%esi
mov    %r12d,%edx
mov    %r12d,%ecx
xor    %eax,%eax
call   1050 <entropy_term@plt>
mov    %eax,0x8(%rsp)
mov    %rbx,%rdi
mov    %r13d,%esi
mov    %r15d,%edx
mov    %r15d,%ecx
xor    %eax,%eax
call   1050 <entropy_term@plt>
add    0x8(%rsp),%eax
neg    %eax
cvtsi2sd %eax,%xmm0
test   %r15d,%r15d
movsd  %xmm0,0x8(%rsp)
jle    129a <mutual+0x16a>
imul   %r12d,%r15d
mov    %r15d,0x14(%rsp)
mov    %r14,0x20(%rsp)
xor    %r15d,%r15d
mov    0x18(%rsp),%r14
nopl   0x0(%rax,%rax,1)
mov    0x20(%rsp),%rax
mov    (%rax,%r15,8),%rdi
mov    %r13d,%esi
mov    %r12d,%edx
mov    0x14(%rsp),%ecx
xor    %eax,%eax
call   1050 <entropy_term@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
add    $0x1,%r15
cmp    %r15,%r14
jne    1260 <mutual+0x130>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
movsd  0x8(%rsp),%xmm0
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @mutual(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
%93 = load i32*, i32** %8
%94 = load i32, i32* %10
%95 = load i32, i32* %6
%96 = load i32, i32* %6
%97 = call i32 (i32*, i32, i32, i32, ...) bitcast (i32 (...)* @entropy_term to i32 (i32*, i32, i32, i32, ...)*)(i32* %93, i32 %94, i32 %95, i32 %96)
%98 = sub nsw i32 0, %97
%99 = load i32*, i32** %9
%100 = load i32, i32* %10
%101 = load i32, i32* %5
%102 = load i32, i32* %5
%103 = call i32 (i32*, i32, i32, i32, ...) bitcast (i32 (...)* @entropy_term to i32 (i32*, i32, i32, i32, ...)*)(i32* %99, i32 %100, i32 %101, i32 %102)
%104 = sub nsw i32 %98, %103
%105 = sitofp i32 %104 to double
store double %105, double* %7
store i32 0, i32* %11
br label %106
106: ; preds = %125, %92
%107 = load i32, i32* %11
%108 = load i32, i32* %5
%109 = icmp slt i32 %107, %108
br i1 %109, label %110, label %128
110: ; preds = %106
%111 = load i32**, i32*** %4
%112 = load i32, i32* %11
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds i32*, i32** %111, i64 %113
%115 = load i32*, i32** %114
%116 = load i32, i32* %10
%117 = load i32, i32* %6
%118 = load i32, i32* %5
%119 = load i32, i32* %6
%120 = mul nsw i32 %118, %119
%121 = call i32 (i32*, i32, i32, i32, ...) bitcast (i32 (...)* @entropy_term to i32 (i32*, i32, i32, i32, ...)*)(i32* %115, i32 %116, i32 %117, i32 %120)
%122 = sitofp i32 %121 to double
%123 = load double, double* %7
%124 = fadd double %123, %122
store double %124, double* %7
br label %125
125: ; preds = %110
%126 = load i32, i32* %11
%127 = add nsw i32 %126, 1
store i32 %127, i32* %11
br label %106
128: ; preds = %106
%129 = load i32*, i32** %8
%130 = bitcast i32* %129 to i8*
%131 = load i32*, i32** %9
%132 = bitcast i32* %131 to i8*
%133 = load double, double* %7
ret double %133
}

LLVM IR (-O2):
define double @mutual(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 %65
12: ; preds = %3
%13 = zext i32 %1 to i64
%14 = zext i32 %2 to i64
%15 = and i64 %14, 1
%16 = icmp eq i32 %2, 1
%17 = and i64 %14, 4294967294
%18 = icmp eq i64 %15, 0
br label %19
19: ; preds = %12, %60
%20 = phi i64 [ 0, %12 ], [ %63, %60 ]
%21 = phi i32 [ 0, %12 ], [ %62, %60 ]
%22 = getelementptr inbounds i32, i32* %6, i64 %20
store i32 0, i32* %22
br i1 %10, label %23, label %60
23: ; preds = %19
%24 = getelementptr inbounds i32*, i32** %0, i64 %20
%25 = load i32*, i32** %24
br i1 %16, label %48, label %26
26: ; preds = %23, %26
%27 = phi i32 [ %40, %26 ], [ 0, %23 ]
%28 = phi i64 [ %45, %26 ], [ 0, %23 ]
%29 = phi i64 [ %46, %26 ], [ 0, %23 ]
%30 = getelementptr inbounds i32, i32* %25, i64 %28
%31 = load i32, i32* %30
%32 = add nsw i32 %27, %31
store i32 %32, i32* %22
%33 = load i32, i32* %30
%34 = getelementptr inbounds i32, i32* %9, i64 %28
%35 = load i32, i32* %34
%36 = add nsw i32 %35, %33
store i32 %36, i32* %34
%37 = or i64 %28, 1
%38 = getelementptr inbounds i32, i32* %25, i64 %37
%39 = load i32, i32* %38
%40 = add nsw i32 %32, %39
store i32 %40, i32* %22
%41 = load i32, i32* %38
%42 = getelementptr inbounds i32, i32* %9, i64 %37
%43 = load i32, i32* %42
%44 = add nsw i32 %43, %41
store i32 %44, i32* %42
%45 = add nuw nsw i64 %28, 2
%46 = add i64 %29, 2
%47 = icmp eq i64 %46, %17
br i1 %47, label %48, label %26
48: ; preds = %26, %23
%49 = phi i32 [ undef, %23 ], [ %40, %26 ]
%50 = phi i32 [ 0, %23 ], [ %40, %26 ]
%51 = phi i64 [ 0, %23 ], [ %45, %26 ]
br i1 %18, label %60, label %52
52: ; preds = %48
%53 = getelementptr inbounds i32, i32* %25, i64 %51
%54 = load i32, i32* %53
%55 = add nsw i32 %50, %54
store i32 %55, i32* %22
%56 = load i32, i32* %53
%57 = getelementptr inbounds i32, i32* %9, i64 %51
%58 = load i32, i32* %57
%59 = add nsw i32 %58, %56
store i32 %59, i32* %57
br label %60
60: ; preds = %52, %48, %19
%61 = phi i32 [ 0, %19 ], [ %49, %48 ], [ %55, %52 ]
%62 = add nsw i32 %61, %21
%63 = add nuw nsw i64 %20, 1
%64 = icmp eq i64 %63, %13
br i1 %64, label %65, label %19
65: ; preds = %60, %3
%66 = phi i32 [ 0, %3 ], [ %62, %60 ]
%67 = tail call i32 (i32*, i32, i32, i32, ...) bitcast (i32 (...)* @entropy_term to i32 (i32*, i32, i32, i32, ...)*)(i32* %6, i32 %66, i32 %2, i32 %2)
%68 = tail call i32 (i32*, i32, i32, i32, ...) bitcast (i32 (...)* @entropy_term to i32 (i32*, i32, i32, i32, ...)*)(i32* %9, i32 %66, i32 %1, i32 %1)
%69 = add i32 %67, %68
%70 = sub i32 0, %69
%71 = sitofp i32 %70 to double
br i1 %11, label %72, label %85
72: ; preds = %65
%73 = mul nsw i32 %2, %1
%74 = zext i32 %1 to i64
br label %75
75: ; preds = %72, %75
%76 = phi i64 [ 0, %72 ], [ %83, %75 ]
%77 = phi double [ %71, %72 ], [ %82, %75 ]
%78 = getelementptr inbounds i32*, i32** %0, i64 %76
%79 = load i32*, i32** %78
%80 = tail call i32 (i32*, i32, i32, i32, ...) bitcast (i32 (...)* @entropy_term to i32 (i32*, i32, i32, i32, ...)*)(i32* %79, i32 %66, i32 %2, i32 %73)
%81 = sitofp i32 %80 to double
%82 = fadd double %77, %81
%83 = add nuw nsw i64 %76, 1
%84 = icmp eq i64 %83, %74
br i1 %84, label %85, label %75
85: ; preds = %75, %65
%86 = phi double [ %71, %65 ], [ %82, %75 ]
tail call void @free(i8* %5)
tail call void @free(i8* %8)
ret double %86
}