void bmod(float row[32][32], float col[32][32], float inner[32][32]) {
  int i, j, k;
  for (i = 0; i < 32; i++)
    for (j = 0; j < 32; j++)
      for (k = 0; k < 32; k++)
        inner[i][j] -= row[i][k] * col[k][j];
}


x86 assembly (clang, -O0):
<bmod>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x20,-0x1c(%rbp)
jge    11da <bmod+0xda>
movl   $0x0,-0x20(%rbp)
cmpl   $0x20,-0x20(%rbp)
jge    11c7 <bmod+0xc7>
movl   $0x0,-0x24(%rbp)
cmpl   $0x20,-0x24(%rbp)
jge    11b4 <bmod+0xb4>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1139 <bmod+0x39>
jmp    11b9 <bmod+0xb9>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1128 <bmod+0x28>
jmp    11cc <bmod+0xcc>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1117 <bmod+0x17>
pop    %rbp
ret

x86 assembly (clang, -O2):
<bmod>:
add    $0x4,%rdi
sub    $0xffffffffffffff80,%rsi
xor    %r8d,%r8d
nopl   0x0(%rax,%rax,1)
mov    %rsi,%r9
xor    %r10d,%r10d
cs nopw 0x0(%rax,%rax,1)
mov    %r8,%rax
shl    $0x7,%rax
add    %rdx,%rax
lea    (%rax,%r10,4),%r11
movss  (%rax,%r10,4),%xmm0
mov    %r9,%rax
xor    %ecx,%ecx
nopl   0x0(%rax)
movss  -0x4(%rdi,%rcx,4),%xmm1
mulss  -0x80(%rax),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,(%r11)
movss  (%rdi,%rcx,4),%xmm1
mulss  (%rax),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,(%r11)
add    $0x2,%rcx
add    $0x100,%rax
cmp    $0x20,%rcx
jne    1140 <bmod+0x40>
add    $0x1,%r10
add    $0x4,%r9
cmp    $0x20,%r10
jne    1120 <bmod+0x20>
add    $0x1,%r8
sub    $0xffffffffffffff80,%rdi
cmp    $0x20,%r8
jne    1110 <bmod+0x10>
ret

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

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