char matrixMultiply(float **output, char outRows, char outCols, float **inputA,
                    char rowsA, char colsA, float **inputB, char rowsB,
                    char colsB) {
  char i = 0;
  char j = 0;
  char k = 0;
  float sum = 0;
  if (colsA != rowsB || outRows != rowsA || outCols != colsB) {
    return 0;
  } else {
    for (i = 0; i < outRows; i++) {
      for (j = 0; j < outCols; j++) {
        for (k = 0; k < colsA; k++) {
          sum = sum + inputA[i][k] * inputB[k][j];
        }
        output[i][j] = sum;
        sum = 0;
      }
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<matrixMultiply>:
push   %rbp
mov    %rsp,%rbp
mov    %r8d,-0x44(%rbp)
mov    %rcx,%rax
mov    -0x44(%rbp),%ecx
mov    %rax,-0x40(%rbp)
mov    %edx,%eax
mov    -0x40(%rbp),%rdx
mov    %eax,-0x34(%rbp)
mov    %esi,%eax
mov    -0x34(%rbp),%esi
mov    %eax,-0x30(%rbp)
mov    %rdi,%r8
mov    -0x30(%rbp),%edi
mov    %r9b,%al
mov    0x20(%rbp),%r9b
mov    0x18(%rbp),%r9b
mov    0x10(%rbp),%r9
mov    %r8,-0x10(%rbp)
mov    %dil,-0x11(%rbp)
mov    %sil,-0x12(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %cl,-0x21(%rbp)
mov    %al,-0x22(%rbp)
movb   $0x0,-0x23(%rbp)
movb   $0x0,-0x24(%rbp)
movb   $0x0,-0x25(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movsbl -0x22(%rbp),%eax
movsbl 0x18(%rbp),%ecx
cmp    %ecx,%eax
jne    1192 <matrixMultiply+0x92>
movsbl -0x11(%rbp),%eax
movsbl -0x21(%rbp),%ecx
cmp    %ecx,%eax
jne    1192 <matrixMultiply+0x92>
movsbl -0x12(%rbp),%eax
movsbl 0x20(%rbp),%ecx
cmp    %ecx,%eax
je     119b <matrixMultiply+0x9b>
movb   $0x0,-0x1(%rbp)
jmp    1270 <matrixMultiply+0x170>
movb   $0x0,-0x23(%rbp)
movsbl -0x23(%rbp),%eax
movsbl -0x11(%rbp),%ecx
cmp    %ecx,%eax
jge    1267 <matrixMultiply+0x167>
movb   $0x0,-0x24(%rbp)
movsbl -0x24(%rbp),%eax
movsbl -0x12(%rbp),%ecx
cmp    %ecx,%eax
jge    1255 <matrixMultiply+0x155>
movb   $0x0,-0x25(%rbp)
movsbl -0x25(%rbp),%eax
movsbl -0x22(%rbp),%ecx
cmp    %ecx,%eax
jge    1224 <matrixMultiply+0x124>
movss  -0x2c(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movsbq -0x23(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsbq -0x25(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x10(%rbp),%rax
movsbq -0x25(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsbq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x25(%rbp),%al
add    $0x1,%al
mov    %al,-0x25(%rbp)
jmp    11c7 <matrixMultiply+0xc7>
movss  -0x2c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsbq -0x23(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsbq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x24(%rbp),%al
add    $0x1,%al
mov    %al,-0x24(%rbp)
jmp    11b3 <matrixMultiply+0xb3>
jmp    125a <matrixMultiply+0x15a>
mov    -0x23(%rbp),%al
add    $0x1,%al
mov    %al,-0x23(%rbp)
jmp    119f <matrixMultiply+0x9f>
jmp    126c <matrixMultiply+0x16c>
movb   $0x1,-0x1(%rbp)
movsbl -0x1(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<matrixMultiply>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edx,-0x14(%rsp)
mov    %rdi,-0x10(%rsp)
xor    %eax,%eax
cmp    %r8b,%sil
jne    1263 <matrixMultiply+0x163>
cmp    0x40(%rsp),%r9b
jne    1263 <matrixMultiply+0x163>
mov    -0x14(%rsp),%edx
cmp    0x48(%rsp),%dl
jne    1263 <matrixMultiply+0x163>
mov    $0x1,%al
test   %sil,%sil
jle    1263 <matrixMultiply+0x163>
mov    0x38(%rsp),%r8
movzbl %sil,%eax
mov    %rax,-0x8(%rsp)
movzbl -0x14(%rsp),%r11d
movzbl %r9b,%r14d
lea    -0x1(%r14),%rdx
mov    %r14d,%r13d
and    $0x3,%r13d
and    $0xfffffffc,%r14d
xor    %r15d,%r15d
jmp    117f <matrixMultiply+0x7f>
xchg   %ax,%ax
add    $0x1,%r15
cmp    -0x8(%rsp),%r15
je     1261 <matrixMultiply+0x161>
cmpb   $0x0,-0x14(%rsp)
jle    1170 <matrixMultiply+0x70>
mov    -0x10(%rsp),%rax
mov    (%rax,%r15,8),%r12
xor    %ebx,%ebx
jmp    11b2 <matrixMultiply+0xb2>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorps  %xmm0,%xmm0
movss  %xmm0,(%r12,%rbx,4)
add    $0x1,%rbx
cmp    %r11,%rbx
je     1170 <matrixMultiply+0x70>
test   %r9b,%r9b
jle    11a0 <matrixMultiply+0xa0>
mov    (%rcx,%r15,8),%rbp
xorps  %xmm0,%xmm0
xor    %r10d,%r10d
cmp    $0x3,%rdx
jb     122c <matrixMultiply+0x12c>
nopw   0x0(%rax,%rax,1)
movss  0x0(%rbp,%r10,4),%xmm1
movss  0x4(%rbp,%r10,4),%xmm2
mov    (%r8,%r10,8),%rax
mov    0x8(%r8,%r10,8),%rdi
mulss  (%rax,%rbx,4),%xmm1
addss  %xmm0,%xmm1
mulss  (%rdi,%rbx,4),%xmm2
addss  %xmm1,%xmm2
movss  0x8(%rbp,%r10,4),%xmm1
mov    0x10(%r8,%r10,8),%rax
mulss  (%rax,%rbx,4),%xmm1
movss  0xc(%rbp,%r10,4),%xmm0
mov    0x18(%r8,%r10,8),%rax
mulss  (%rax,%rbx,4),%xmm0
addss  %xmm2,%xmm1
addss  %xmm1,%xmm0
add    $0x4,%r10
cmp    %r10,%r14
jne    11d0 <matrixMultiply+0xd0>
test   %r13,%r13
je     11a3 <matrixMultiply+0xa3>
lea    (%r8,%r10,8),%rax
lea    0x0(%rbp,%r10,4),%rbp
xor    %edi,%edi
movss  0x0(%rbp,%rdi,4),%xmm1
mov    (%rax,%rdi,8),%rsi
mulss  (%rsi,%rbx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rdi
cmp    %rdi,%r13
jne    1240 <matrixMultiply+0x140>
jmp    11a3 <matrixMultiply+0xa3>
mov    $0x1,%al
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define signext i8 @matrixMultiply(float** %0, i8 signext %1, i8 signext %2, float** %3, i8 signext %4, i8 signext %5, float** %6, i8 signext %7, i8 signext %8) {
%10 = alloca i8
%11 = alloca float**
%12 = alloca i8
%13 = alloca i8
%14 = alloca float**
%15 = alloca i8
%16 = alloca i8
%17 = alloca float**
%18 = alloca i8
%19 = alloca i8
%20 = alloca i8
%21 = alloca i8
%22 = alloca i8
%23 = alloca float
store float** %0, float*** %11
store i8 %1, i8* %12
store i8 %2, i8* %13
store float** %3, float*** %14
store i8 %4, i8* %15
store i8 %5, i8* %16
store float** %6, float*** %17
store i8 %7, i8* %18
store i8 %8, i8* %19
store i8 0, i8* %20
store i8 0, i8* %21
store i8 0, i8* %22
store float 0.000000e+00, float* %23
%24 = load i8, i8* %16
%25 = sext i8 %24 to i32
%26 = load i8, i8* %18
%27 = sext i8 %26 to i32
%28 = icmp ne i32 %25, %27
br i1 %28, label %41, label %29
29: ; preds = %9
%30 = load i8, i8* %12
%31 = sext i8 %30 to i32
%32 = load i8, i8* %15
%33 = sext i8 %32 to i32
%34 = icmp ne i32 %31, %33
br i1 %34, label %41, label %35
35: ; preds = %29
%36 = load i8, i8* %13
%37 = sext i8 %36 to i32
%38 = load i8, i8* %19
%39 = sext i8 %38 to i32
%40 = icmp ne i32 %37, %39
br i1 %40, label %41, label %42
41: ; preds = %35, %29, %9
store i8 0, i8* %10
br label %106
42: ; preds = %35
store i8 0, i8* %20
br label %43
43: ; preds = %101, %42
%44 = load i8, i8* %20
%45 = sext i8 %44 to i32
%46 = load i8, i8* %12
%47 = sext i8 %46 to i32
%48 = icmp slt i32 %45, %47
br i1 %48, label %49, label %104
49: ; preds = %43
store i8 0, i8* %21
br label %50
50: ; preds = %97, %49
%51 = load i8, i8* %21
%52 = sext i8 %51 to i32
%53 = load i8, i8* %13
%54 = sext i8 %53 to i32
%55 = icmp slt i32 %52, %54
br i1 %55, label %56, label %100
56: ; preds = %50
store i8 0, i8* %22
br label %57
57: ; preds = %84, %56
%58 = load i8, i8* %22
%59 = sext i8 %58 to i32
%60 = load i8, i8* %16
%61 = sext i8 %60 to i32
%62 = icmp slt i32 %59, %61
br i1 %62, label %63, label %87
63: ; preds = %57
%64 = load float, float* %23
%65 = load float**, float*** %14
%66 = load i8, i8* %20
%67 = sext i8 %66 to i64
%68 = getelementptr inbounds float*, float** %65, i64 %67
%69 = load float*, float** %68
%70 = load i8, i8* %22
%71 = sext i8 %70 to i64
%72 = getelementptr inbounds float, float* %69, i64 %71
%73 = load float, float* %72
%74 = load float**, float*** %17
%75 = load i8, i8* %22
%76 = sext i8 %75 to i64
%77 = getelementptr inbounds float*, float** %74, i64 %76
%78 = load float*, float** %77
%79 = load i8, i8* %21
%80 = sext i8 %79 to i64
%81 = getelementptr inbounds float, float* %78, i64 %80
%82 = load float, float* %81
%83 = call float @llvm.fmuladd.f32(float %73, float %82, float %64)
store float %83, float* %23
br label %84
84: ; preds = %63
%85 = load i8, i8* %22
%86 = add i8 %85, 1
store i8 %86, i8* %22
br label %57
87: ; preds = %57
%88 = load float, float* %23
%89 = load float**, float*** %11
%90 = load i8, i8* %20
%91 = sext i8 %90 to i64
%92 = getelementptr inbounds float*, float** %89, i64 %91
%93 = load float*, float** %92
%94 = load i8, i8* %21
%95 = sext i8 %94 to i64
%96 = getelementptr inbounds float, float* %93, i64 %95
store float %88, float* %96
store float 0.000000e+00, float* %23
br label %97
97: ; preds = %87
%98 = load i8, i8* %21
%99 = add i8 %98, 1
store i8 %99, i8* %21
br label %50
100: ; preds = %50
br label %101
101: ; preds = %100
%102 = load i8, i8* %20
%103 = add i8 %102, 1
store i8 %103, i8* %20
br label %43
104: ; preds = %43
br label %105
105: ; preds = %104
store i8 1, i8* %10
br label %106
106: ; preds = %105, %41
%107 = load i8, i8* %10
ret i8 %107
}

LLVM IR (-O2):
define signext i8 @matrixMultiply(float** nocapture %0, i8 signext %1, i8 signext %2, float** nocapture %3, i8 signext %4, i8 signext %5, float** nocapture %6, i8 signext %7, i8 signext %8) {
%10 = icmp eq i8 %5, %7
%11 = icmp eq i8 %1, %4
%12 = and i1 %11, %10
%13 = icmp eq i8 %2, %8
%14 = and i1 %12, %13
br i1 %14, label %15, label %102
15: ; preds = %9
%16 = icmp sgt i8 %1, 0
br i1 %16, label %17, label %102
17: ; preds = %15
%18 = icmp sgt i8 %2, 0
%19 = icmp sgt i8 %5, 0
%20 = zext i8 %1 to i64
%21 = zext i8 %2 to i64
%22 = zext i8 %5 to i64
%23 = add nsw i64 %22, -1
%24 = and i64 %22, 3
%25 = icmp ult i64 %23, 3
%26 = and i64 %22, 252
%27 = icmp eq i64 %24, 0
br label %28
28: ; preds = %17, %99
%29 = phi i64 [ 0, %17 ], [ %100, %99 ]
br i1 %18, label %30, label %99
30: ; preds = %28
%31 = getelementptr inbounds float*, float** %3, i64 %29
%32 = getelementptr inbounds float*, float** %0, i64 %29
%33 = load float*, float** %32
br label %34
34: ; preds = %30, %94
%35 = phi i64 [ 0, %30 ], [ %97, %94 ]
br i1 %19, label %36, label %94
36: ; preds = %34
%37 = load float*, float** %31
br i1 %25, label %76, label %38
38: ; preds = %36, %38
%39 = phi i64 [ %73, %38 ], [ 0, %36 ]
%40 = phi float [ %72, %38 ], [ 0.000000e+00, %36 ]
%41 = phi i64 [ %74, %38 ], [ 0, %36 ]
%42 = getelementptr inbounds float, float* %37, i64 %39
%43 = load float, float* %42
%44 = getelementptr inbounds float*, float** %6, i64 %39
%45 = load float*, float** %44
%46 = getelementptr inbounds float, float* %45, i64 %35
%47 = load float, float* %46
%48 = tail call float @llvm.fmuladd.f32(float %43, float %47, float %40)
%49 = or i64 %39, 1
%50 = getelementptr inbounds float, float* %37, i64 %49
%51 = load float, float* %50
%52 = getelementptr inbounds float*, float** %6, i64 %49
%53 = load float*, float** %52
%54 = getelementptr inbounds float, float* %53, i64 %35
%55 = load float, float* %54
%56 = tail call float @llvm.fmuladd.f32(float %51, float %55, float %48)
%57 = or i64 %39, 2
%58 = getelementptr inbounds float, float* %37, i64 %57
%59 = load float, float* %58
%60 = getelementptr inbounds float*, float** %6, i64 %57
%61 = load float*, float** %60
%62 = getelementptr inbounds float, float* %61, i64 %35
%63 = load float, float* %62
%64 = tail call float @llvm.fmuladd.f32(float %59, float %63, float %56)
%65 = or i64 %39, 3
%66 = getelementptr inbounds float, float* %37, i64 %65
%67 = load float, float* %66
%68 = getelementptr inbounds float*, float** %6, i64 %65
%69 = load float*, float** %68
%70 = getelementptr inbounds float, float* %69, i64 %35
%71 = load float, float* %70
%72 = tail call float @llvm.fmuladd.f32(float %67, float %71, float %64)
%73 = add nuw nsw i64 %39, 4
%74 = add i64 %41, 4
%75 = icmp eq i64 %74, %26
br i1 %75, label %76, label %38
76: ; preds = %38, %36
%77 = phi float [ undef, %36 ], [ %72, %38 ]
%78 = phi i64 [ 0, %36 ], [ %73, %38 ]
%79 = phi float [ 0.000000e+00, %36 ], [ %72, %38 ]
br i1 %27, label %94, label %80
80: ; preds = %76, %80
%81 = phi i64 [ %91, %80 ], [ %78, %76 ]
%82 = phi float [ %90, %80 ], [ %79, %76 ]
%83 = phi i64 [ %92, %80 ], [ 0, %76 ]
%84 = getelementptr inbounds float, float* %37, i64 %81
%85 = load float, float* %84
%86 = getelementptr inbounds float*, float** %6, i64 %81
%87 = load float*, float** %86
%88 = getelementptr inbounds float, float* %87, i64 %35
%89 = load float, float* %88
%90 = tail call float @llvm.fmuladd.f32(float %85, float %89, float %82)
%91 = add nuw nsw i64 %81, 1
%92 = add i64 %83, 1
%93 = icmp eq i64 %92, %24
br i1 %93, label %94, label %80
94: ; preds = %76, %80, %34
%95 = phi float [ 0.000000e+00, %34 ], [ %77, %76 ], [ %90, %80 ]
%96 = getelementptr inbounds float, float* %33, i64 %35
store float %95, float* %96
%97 = add nuw nsw i64 %35, 1
%98 = icmp eq i64 %97, %21
br i1 %98, label %99, label %34
99: ; preds = %94, %28
%100 = add nuw nsw i64 %29, 1
%101 = icmp eq i64 %100, %20
br i1 %101, label %102, label %28
102: ; preds = %99, %15, %9
%103 = phi i8 [ 0, %9 ], [ 1, %15 ], [ 1, %99 ]
ret i8 %103
}