#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 */
double log(double);
void transLogOdds(int *mdim, double *empTR, double *empTLO) {
  int i, j, k, pos;
  for (i = 0; i < mdim[2]; i++) {
    for (j = 0; j < *mdim; j++) {
      for (k = 0; k < *mdim; k++) {
        pos = *mdim * k + j + *mdim * *mdim * i;
        empTLO[pos] = log(empTR[pos] / (1.0 - empTR[pos]));
      }
    }
  }
}


x86 assembly (clang, -O0):
<transLogOdds>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x8(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    11f6 <transLogOdds+0xe6>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x8(%rbp),%rcx
cmp    (%rcx),%eax
jge    11e3 <transLogOdds+0xd3>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x8(%rbp),%rcx
cmp    (%rcx),%eax
jge    11d0 <transLogOdds+0xc0>
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
imul   -0x24(%rbp),%eax
add    -0x20(%rbp),%eax
mov    -0x8(%rbp),%rcx
mov    (%rcx),%ecx
mov    -0x8(%rbp),%rdx
imul   (%rdx),%ecx
imul   -0x1c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  0xe59(%rip),%xmm1
subsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
call   1030 <log@plt>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1158 <transLogOdds+0x48>
jmp    11d5 <transLogOdds+0xc5>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1142 <transLogOdds+0x32>
jmp    11e8 <transLogOdds+0xd8>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112b <transLogOdds+0x1b>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<transLogOdds>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    0x8(%rdi),%ecx
test   %ecx,%ecx
jle    11c4 <transLogOdds+0xb4>
mov    %rdx,%r15
mov    %rsi,%r12
mov    %rdi,%r13
mov    (%rdi),%r8d
xor    %ebx,%ebx
jmp    114b <transLogOdds+0x3b>
cs nopw 0x0(%rax,%rax,1)
mov    0x8(%r13),%ecx
add    $0x1,%ebx
cmp    %ecx,%ebx
jge    11c4 <transLogOdds+0xb4>
test   %r8d,%r8d
jle    1144 <transLogOdds+0x34>
xor    %eax,%eax
jmp    1168 <transLogOdds+0x58>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%eax
cmp    %r8d,%eax
jge    1140 <transLogOdds+0x30>
test   %r8d,%r8d
jle    1160 <transLogOdds+0x50>
xor    %ebp,%ebp
mov    %eax,0x4(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %ebx,%ecx
imul   %r8d,%ecx
add    %ebp,%ecx
imul   %r8d,%ecx
add    %eax,%ecx
movslq %ecx,%r14
movsd  (%r12,%r14,8),%xmm0
movsd  0xe61(%rip),%xmm1
subsd  %xmm0,%xmm1
divsd  %xmm1,%xmm0
call   1030 <log@plt>
mov    0x4(%rsp),%eax
movsd  %xmm0,(%r15,%r14,8)
add    $0x1,%ebp
mov    0x0(%r13),%r8d
cmp    %r8d,%ebp
jl     1180 <transLogOdds+0x70>
jmp    1160 <transLogOdds+0x50>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @transLogOdds(i32* %0, double* %1, double* %2) {
%4 = alloca i32*
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32* %0, i32** %4
store double* %1, double** %5
store double* %2, double** %6
store i32 0, i32* %7
br label %11
11: ; preds = %69, %3
%12 = load i32, i32* %7
%13 = load i32*, i32** %4
%14 = getelementptr inbounds i32, i32* %13, i64 2
%15 = load i32, i32* %14
%16 = icmp slt i32 %12, %15
br i1 %16, label %17, label %72
17: ; preds = %11
store i32 0, i32* %8
br label %18
18: ; preds = %65, %17
%19 = load i32, i32* %8
%20 = load i32*, i32** %4
%21 = load i32, i32* %20
%22 = icmp slt i32 %19, %21
br i1 %22, label %23, label %68
23: ; preds = %18
store i32 0, i32* %9
br label %24
24: ; preds = %61, %23
%25 = load i32, i32* %9
%26 = load i32*, i32** %4
%27 = load i32, i32* %26
%28 = icmp slt i32 %25, %27
br i1 %28, label %29, label %64
29: ; preds = %24
%30 = load i32*, i32** %4
%31 = load i32, i32* %30
%32 = load i32, i32* %9
%33 = mul nsw i32 %31, %32
%34 = load i32, i32* %8
%35 = add nsw i32 %33, %34
%36 = load i32*, i32** %4
%37 = load i32, i32* %36
%38 = load i32*, i32** %4
%39 = load i32, i32* %38
%40 = mul nsw i32 %37, %39
%41 = load i32, i32* %7
%42 = mul nsw i32 %40, %41
%43 = add nsw i32 %35, %42
store i32 %43, i32* %10
%44 = load double*, double** %5
%45 = load i32, i32* %10
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %44, i64 %46
%48 = load double, double* %47
%49 = load double*, double** %5
%50 = load i32, i32* %10
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
%53 = load double, double* %52
%54 = fsub double 1.000000e+00, %53
%55 = fdiv double %48, %54
%56 = call double @log(double %55)
%57 = load double*, double** %6
%58 = load i32, i32* %10
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %57, i64 %59
store double %56, double* %60
br label %61
61: ; preds = %29
%62 = load i32, i32* %9
%63 = add nsw i32 %62, 1
store i32 %63, i32* %9
br label %24
64: ; preds = %24
br label %65
65: ; preds = %64
%66 = load i32, i32* %8
%67 = add nsw i32 %66, 1
store i32 %67, i32* %8
br label %18
68: ; preds = %18
br label %69
69: ; preds = %68
%70 = load i32, i32* %7
%71 = add nsw i32 %70, 1
store i32 %71, i32* %7
br label %11
72: ; preds = %11
ret void
}

LLVM IR (-O2):
define void @transLogOdds(i32* nocapture %0, double* nocapture %1, double* nocapture writeonly %2) {
%4 = getelementptr inbounds i32, i32* %0, i64 2
%5 = load i32, i32* %4
%6 = icmp sgt i32 %5, 0
br i1 %6, label %7, label %49
7: ; preds = %3
%8 = load i32, i32* %0
br label %9
9: ; preds = %7, %43
%10 = phi i32 [ %44, %43 ], [ %5, %7 ]
%11 = phi i32 [ %45, %43 ], [ %8, %7 ]
%12 = phi i32 [ %46, %43 ], [ %8, %7 ]
%13 = phi i32 [ %47, %43 ], [ 0, %7 ]
%14 = icmp sgt i32 %12, 0
br i1 %14, label %15, label %43
15: ; preds = %9, %37
%16 = phi i32 [ %38, %37 ], [ %11, %9 ]
%17 = phi i32 [ %38, %37 ], [ %12, %9 ]
%18 = phi i32 [ %39, %37 ], [ 0, %9 ]
%19 = icmp sgt i32 %17, 0
br i1 %19, label %20, label %37
20: ; preds = %15, %20
%21 = phi i32 [ %35, %20 ], [ %17, %15 ]
%22 = phi i32 [ %34, %20 ], [ 0, %15 ]
%23 = mul i32 %21, %13
%24 = add i32 %23, %22
%25 = mul i32 %24, %21
%26 = add i32 %25, %18
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %1, i64 %27
%29 = load double, double* %28
%30 = fsub double 1.000000e+00, %29
%31 = fdiv double %29, %30
%32 = tail call double @log(double %31)
%33 = getelementptr inbounds double, double* %2, i64 %27
store double %32, double* %33
%34 = add nuw nsw i32 %22, 1
%35 = load i32, i32* %0
%36 = icmp slt i32 %34, %35
br i1 %36, label %20, label %37
37: ; preds = %20, %15
%38 = phi i32 [ %16, %15 ], [ %35, %20 ]
%39 = add nuw nsw i32 %18, 1
%40 = icmp slt i32 %39, %38
br i1 %40, label %15, label %41
41: ; preds = %37
%42 = load i32, i32* %4
br label %43
43: ; preds = %41, %9
%44 = phi i32 [ %42, %41 ], [ %10, %9 ]
%45 = phi i32 [ %38, %41 ], [ %11, %9 ]
%46 = phi i32 [ %38, %41 ], [ %12, %9 ]
%47 = add nuw nsw i32 %13, 1
%48 = icmp slt i32 %47, %44
br i1 %48, label %9, label %49
49: ; preds = %43, %3
ret void
}