void transProbs(int *mpoints, int *nk, double *rwsum, double *empTR) {
  int i, j, k;
  for (i = 0; i < *mpoints; i++) {
    for (j = 0; j < *nk; j++) {
      for (k = 0; k < *nk; k++) {
        empTR[*nk * k + j + *nk * *nk * i] /= rwsum[*nk * i + j];
      }
    }
  }
}


x86 assembly (clang, -O0):
<transProbs>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x8(%rbp),%rcx
cmp    (%rcx),%eax
jge    11d8 <transProbs+0xd8>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    11c5 <transProbs+0xc5>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    11b2 <transProbs+0xb2>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    (%rcx),%ecx
imul   -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    (%rcx),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    -0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    -0x10(%rbp),%rsi
imul   (%rsi),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1147 <transProbs+0x47>
jmp    11b7 <transProbs+0xb7>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1131 <transProbs+0x31>
jmp    11ca <transProbs+0xca>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111b <transProbs+0x1b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<transProbs>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    (%rdi),%eax
mov    %rax,-0x8(%rsp)
test   %eax,%eax
jle    120d <transProbs+0x10d>
mov    (%rsi),%esi
movslq %esi,%rax
mov    %eax,%edi
and    $0xfffffffe,%edi
mov    %eax,%ebp
imul   %eax,%ebp
mov    %ebp,-0x14(%rsp)
mov    %rax,-0x10(%rsp)
add    %eax,%eax
xor    %ebp,%ebp
mov    %rbp,-0x20(%rsp)
xor    %r14d,%r14d
jmp    115d <transProbs+0x5d>
nop
add    $0x1,%r14
mov    -0x20(%rsp),%rbp
add    -0x14(%rsp),%ebp
mov    %rbp,-0x20(%rsp)
cmp    -0x8(%rsp),%r14
je     120d <transProbs+0x10d>
test   %esi,%esi
jle    1140 <transProbs+0x40>
mov    %r14,%r12
imul   -0x10(%rsp),%r12
mov    %esi,%r15d
imul   %r14d,%r15d
mov    -0x20(%rsp),%rbp
mov    %ebp,%r13d
xor    %r11d,%r11d
jmp    118d <transProbs+0x8d>
xchg   %ax,%ax
add    $0x1,%r11
add    $0x1,%r13d
cmp    %rsi,%r11
je     1140 <transProbs+0x40>
lea    (%r12,%r11,1),%rbp
cmp    $0x1,%esi
jne    11a0 <transProbs+0xa0>
xor    %r8d,%r8d
jmp    11e6 <transProbs+0xe6>
nopl   0x0(%rax,%rax,1)
mov    %r13d,%r9d
xor    %r8d,%r8d
cs nopw 0x0(%rax,%rax,1)
movslq %r9d,%r9
movsd  (%rcx,%r9,8),%xmm0
divsd  (%rdx,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%r9,8)
lea    (%rsi,%r9,1),%r10d
movslq %r10d,%rbx
movsd  (%rcx,%rbx,8),%xmm0
divsd  (%rdx,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%rbx,8)
add    $0x2,%r8
add    %eax,%r9d
cmp    %r8,%rdi
jne    11b0 <transProbs+0xb0>
test   $0x1,%sil
je     1180 <transProbs+0x80>
add    %r15d,%r8d
imul   %esi,%r8d
add    %r11d,%r8d
movslq %r8d,%rbx
movsd  (%rcx,%rbx,8),%xmm0
divsd  (%rdx,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%rbx,8)
jmp    1180 <transProbs+0x80>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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