void normalMM(float c[1024][1024], float a[1024][1024], float b[1024][1024]) {
  int i, j, k;
  for (i = 0; i < 1024; i++)
    for (j = 0; j < 1024; j++)
      for (k = 0; k < 1024; k++)
        c[j][i] = c[j][i] + a[k][i] * b[j][k];
}


x86 assembly (clang, -O0):
<normalMM>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x400,-0x1c(%rbp)
jge    11e8 <normalMM+0xe8>
movl   $0x0,-0x20(%rbp)
cmpl   $0x400,-0x20(%rbp)
jge    11d5 <normalMM+0xd5>
movl   $0x0,-0x24(%rbp)
cmpl   $0x400,-0x24(%rbp)
jge    11c2 <normalMM+0xc2>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0xc,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0xc,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0xc,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
shl    $0xc,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    113f <normalMM+0x3f>
jmp    11c7 <normalMM+0xc7>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112b <normalMM+0x2b>
jmp    11da <normalMM+0xda>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1117 <normalMM+0x17>
pop    %rbp
ret

x86 assembly (clang, -O2):
<normalMM>:
add    $0x1000,%rsi
add    $0x4,%rdx
xor    %r8d,%r8d
xchg   %ax,%ax
mov    %rdx,%r10
xor    %r9d,%r9d
cs nopw 0x0(%rax,%rax,1)
mov    %r9,%rax
shl    $0xc,%rax
add    %rdi,%rax
lea    (%rax,%r8,4),%r11
movss  (%rax,%r8,4),%xmm0
mov    %rsi,%rcx
xor    %eax,%eax
nopl   0x0(%rax)
movss  -0x1000(%rcx),%xmm1
ff
mulss  -0x4(%r10,%rax,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r11)
movss  (%rcx),%xmm0
mulss  (%r10,%rax,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r11)
add    $0x2,%rax
add    $0x2000,%rcx
cmp    $0x400,%rax
jne    1140 <normalMM+0x40>
add    $0x1,%r9
add    $0x1000,%r10
cmp    $0x400,%r9
jne    1120 <normalMM+0x20>
add    $0x1,%r8
add    $0x4,%rsi
cmp    $0x400,%r8
jne    1110 <normalMM+0x10>
ret

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

LLVM IR (-O2):
define void @normalMM([1024 x float]* nocapture %0, [1024 x float]* nocapture %1, [1024 x float]* nocapture %2) {
br label %4
4: ; preds = %3, %29
%5 = phi i64 [ 0, %3 ], [ %30, %29 ]
br label %6
6: ; preds = %4, %26
%7 = phi i64 [ 0, %4 ], [ %27, %26 ]
%8 = getelementptr inbounds [1024 x float], [1024 x float]* %0, i64 %7, i64 %5
%9 = load float, float* %8
br label %10
10: ; preds = %10, %6
%11 = phi float [ %9, %6 ], [ %23, %10 ]
%12 = phi i64 [ 0, %6 ], [ %24, %10 ]
%13 = getelementptr inbounds [1024 x float], [1024 x float]* %1, i64 %12, i64 %5
%14 = load float, float* %13
%15 = getelementptr inbounds [1024 x float], [1024 x float]* %2, i64 %7, i64 %12
%16 = load float, float* %15
%17 = tail call float @llvm.fmuladd.f32(float %14, float %16, float %11)
store float %17, float* %8
%18 = or i64 %12, 1
%19 = getelementptr inbounds [1024 x float], [1024 x float]* %1, i64 %18, i64 %5
%20 = load float, float* %19
%21 = getelementptr inbounds [1024 x float], [1024 x float]* %2, i64 %7, i64 %18
%22 = load float, float* %21
%23 = tail call float @llvm.fmuladd.f32(float %20, float %22, float %17)
store float %23, float* %8
%24 = add nuw nsw i64 %12, 2
%25 = icmp eq i64 %24, 1024
br i1 %25, label %26, label %10
26: ; preds = %10
%27 = add nuw nsw i64 %7, 1
%28 = icmp eq i64 %27, 1024
br i1 %28, label %29, label %6
29: ; preds = %26
%30 = add nuw nsw i64 %5, 1
%31 = icmp eq i64 %30, 1024
br i1 %31, label %32, label %4
32: ; preds = %29
ret void
}