cnmgrMatrixMultiply4x4(float matrixOne[4][4], float matrixTwo[4][4],
                       float outputMatrix[4][4]) {
  int notDone = 1;
  int x = 0;
  int y = 0;
  while (notDone) {
    outputMatrix[y][x] =
        matrixOne[y][0] * matrixTwo[0][x] + matrixOne[y][1] * matrixTwo[1][x] +
        matrixOne[y][2] * matrixTwo[2][x] + matrixOne[y][3] * matrixTwo[3][x];
    x++;
    if (x == 4) {
      x = 0;
      y++;
      if (y == 4) {
        notDone = 0;
      }
    }
  }
}


x86 assembly (clang, -O0):
<cnmgrMatrixMultiply4x4>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movl   $0x1,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     1223 <cnmgrMatrixMultiply4x4+0x123>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulss  0x10(%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movss  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  0x20(%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movss  0xc(%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  0x30(%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x4,-0x28(%rbp)
jne    121e <cnmgrMatrixMultiply4x4+0x11e>
movl   $0x0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x4,-0x2c(%rbp)
jne    1219 <cnmgrMatrixMultiply4x4+0x119>
movl   $0x0,-0x24(%rbp)
jmp    121e <cnmgrMatrixMultiply4x4+0x11e>
jmp    1125 <cnmgrMatrixMultiply4x4+0x25>
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<cnmgrMatrixMultiply4x4>:
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rdi,%rax,1),%xmm0
movss  0x4(%rdi,%rax,1),%xmm1
mulss  0x10(%rsi),%xmm1
mulss  (%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rdi,%rax,1),%xmm1
mulss  0x20(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  0xc(%rdi,%rax,1),%xmm0
mulss  0x30(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rdx,%rax,1)
movss  (%rdi,%rax,1),%xmm0
movss  0x4(%rdi,%rax,1),%xmm1
mulss  0x14(%rsi),%xmm1
mulss  0x4(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rdi,%rax,1),%xmm1
mulss  0x24(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  0xc(%rdi,%rax,1),%xmm0
mulss  0x34(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x4(%rdx,%rax,1)
movss  (%rdi,%rax,1),%xmm0
movss  0x4(%rdi,%rax,1),%xmm1
mulss  0x18(%rsi),%xmm1
mulss  0x8(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rdi,%rax,1),%xmm1
mulss  0x28(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  0xc(%rdi,%rax,1),%xmm0
mulss  0x38(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x8(%rdx,%rax,1)
movss  (%rdi,%rax,1),%xmm0
movss  0x4(%rdi,%rax,1),%xmm1
mulss  0x1c(%rsi),%xmm1
mulss  0xc(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rdi,%rax,1),%xmm1
mulss  0x2c(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  0xc(%rdi,%rax,1),%xmm0
mulss  0x3c(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0xc(%rdx,%rax,1)
add    $0x10,%rax
cmp    $0x40,%rax
jne    1110 <cnmgrMatrixMultiply4x4+0x10>
ret

LLVM IR (-O0):
define i32 @cnmgrMatrixMultiply4x4([4 x float]* %0, [4 x float]* %1, [4 x float]* %2) {
%4 = alloca i32
%5 = alloca [4 x float]*
%6 = alloca [4 x float]*
%7 = alloca [4 x float]*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store [4 x float]* %0, [4 x float]** %5
store [4 x float]* %1, [4 x float]** %6
store [4 x float]* %2, [4 x float]** %7
store i32 1, i32* %8
store i32 0, i32* %9
store i32 0, i32* %10
br label %11
11: ; preds = %85, %3
%12 = load i32, i32* %8
%13 = icmp ne i32 %12, 0
br i1 %13, label %14, label %86
14: ; preds = %11
%15 = load [4 x float]*, [4 x float]** %5
%16 = load i32, i32* %10
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds [4 x float], [4 x float]* %15, i64 %17
%19 = getelementptr inbounds [4 x float], [4 x float]* %18, i64 0, i64 0
%20 = load float, float* %19
%21 = load [4 x float]*, [4 x float]** %6
%22 = getelementptr inbounds [4 x float], [4 x float]* %21, i64 0
%23 = load i32, i32* %9
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds [4 x float], [4 x float]* %22, i64 0, i64 %24
%26 = load float, float* %25
%27 = load [4 x float]*, [4 x float]** %5
%28 = load i32, i32* %10
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds [4 x float], [4 x float]* %27, i64 %29
%31 = getelementptr inbounds [4 x float], [4 x float]* %30, i64 0, i64 1
%32 = load float, float* %31
%33 = load [4 x float]*, [4 x float]** %6
%34 = getelementptr inbounds [4 x float], [4 x float]* %33, i64 1
%35 = load i32, i32* %9
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [4 x float], [4 x float]* %34, i64 0, i64 %36
%38 = load float, float* %37
%39 = fmul float %32, %38
%40 = call float @llvm.fmuladd.f32(float %20, float %26, float %39)
%41 = load [4 x float]*, [4 x float]** %5
%42 = load i32, i32* %10
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds [4 x float], [4 x float]* %41, i64 %43
%45 = getelementptr inbounds [4 x float], [4 x float]* %44, i64 0, i64 2
%46 = load float, float* %45
%47 = load [4 x float]*, [4 x float]** %6
%48 = getelementptr inbounds [4 x float], [4 x float]* %47, i64 2
%49 = load i32, i32* %9
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds [4 x float], [4 x float]* %48, i64 0, i64 %50
%52 = load float, float* %51
%53 = call float @llvm.fmuladd.f32(float %46, float %52, float %40)
%54 = load [4 x float]*, [4 x float]** %5
%55 = load i32, i32* %10
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds [4 x float], [4 x float]* %54, i64 %56
%58 = getelementptr inbounds [4 x float], [4 x float]* %57, i64 0, i64 3
%59 = load float, float* %58
%60 = load [4 x float]*, [4 x float]** %6
%61 = getelementptr inbounds [4 x float], [4 x float]* %60, i64 3
%62 = load i32, i32* %9
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds [4 x float], [4 x float]* %61, i64 0, i64 %63
%65 = load float, float* %64
%66 = call float @llvm.fmuladd.f32(float %59, float %65, float %53)
%67 = load [4 x float]*, [4 x float]** %7
%68 = load i32, i32* %10
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds [4 x float], [4 x float]* %67, i64 %69
%71 = load i32, i32* %9
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds [4 x float], [4 x float]* %70, i64 0, i64 %72
store float %66, float* %73
%74 = load i32, i32* %9
%75 = add nsw i32 %74, 1
store i32 %75, i32* %9
%76 = load i32, i32* %9
%77 = icmp eq i32 %76, 4
br i1 %77, label %78, label %85
78: ; preds = %14
store i32 0, i32* %9
%79 = load i32, i32* %10
%80 = add nsw i32 %79, 1
store i32 %80, i32* %10
%81 = load i32, i32* %10
%82 = icmp eq i32 %81, 4
br i1 %82, label %83, label %84
83: ; preds = %78
store i32 0, i32* %8
br label %84
84: ; preds = %83, %78
br label %85
85: ; preds = %84, %14
br label %11
86: ; preds = %11
%87 = load i32, i32* %4
ret i32 %87
}

LLVM IR (-O2):
define i32 @cnmgrMatrixMultiply4x4([4 x float]* nocapture %0, [4 x float]* nocapture %1, [4 x float]* nocapture writeonly %2) {
%4 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 0
%5 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 0
%6 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 0
%7 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 0
%8 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 1
%9 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 1
%10 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 1
%11 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 1
%12 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 2
%13 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 2
%14 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 2
%15 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 2
%16 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 3
%17 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 3
%18 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 3
%19 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 3
%20 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 0
%21 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 1
%22 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 2
%23 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 0, i64 3
br label %24
24: ; preds = %3, %24
%25 = phi float* [ %23, %3 ], [ %86, %24 ]
%26 = phi float* [ %22, %3 ], [ %85, %24 ]
%27 = phi float* [ %21, %3 ], [ %84, %24 ]
%28 = phi float* [ %20, %3 ], [ %83, %24 ]
%29 = phi i64 [ 0, %3 ], [ %82, %24 ]
%30 = load float, float* %28
%31 = load float, float* %4
%32 = load float, float* %27
%33 = load float, float* %5
%34 = fmul float %32, %33
%35 = tail call float @llvm.fmuladd.f32(float %30, float %31, float %34)
%36 = load float, float* %26
%37 = load float, float* %6
%38 = tail call float @llvm.fmuladd.f32(float %36, float %37, float %35)
%39 = load float, float* %25
%40 = load float, float* %7
%41 = tail call float @llvm.fmuladd.f32(float %39, float %40, float %38)
%42 = getelementptr inbounds [4 x float], [4 x float]* %2, i64 %29, i64 0
store float %41, float* %42
%43 = load float, float* %28
%44 = load float, float* %8
%45 = load float, float* %27
%46 = load float, float* %9
%47 = fmul float %45, %46
%48 = tail call float @llvm.fmuladd.f32(float %43, float %44, float %47)
%49 = load float, float* %26
%50 = load float, float* %10
%51 = tail call float @llvm.fmuladd.f32(float %49, float %50, float %48)
%52 = load float, float* %25
%53 = load float, float* %11
%54 = tail call float @llvm.fmuladd.f32(float %52, float %53, float %51)
%55 = getelementptr inbounds [4 x float], [4 x float]* %2, i64 %29, i64 1
store float %54, float* %55
%56 = load float, float* %28
%57 = load float, float* %12
%58 = load float, float* %27
%59 = load float, float* %13
%60 = fmul float %58, %59
%61 = tail call float @llvm.fmuladd.f32(float %56, float %57, float %60)
%62 = load float, float* %26
%63 = load float, float* %14
%64 = tail call float @llvm.fmuladd.f32(float %62, float %63, float %61)
%65 = load float, float* %25
%66 = load float, float* %15
%67 = tail call float @llvm.fmuladd.f32(float %65, float %66, float %64)
%68 = getelementptr inbounds [4 x float], [4 x float]* %2, i64 %29, i64 2
store float %67, float* %68
%69 = load float, float* %28
%70 = load float, float* %16
%71 = load float, float* %27
%72 = load float, float* %17
%73 = fmul float %71, %72
%74 = tail call float @llvm.fmuladd.f32(float %69, float %70, float %73)
%75 = load float, float* %26
%76 = load float, float* %18
%77 = tail call float @llvm.fmuladd.f32(float %75, float %76, float %74)
%78 = load float, float* %25
%79 = load float, float* %19
%80 = tail call float @llvm.fmuladd.f32(float %78, float %79, float %77)
%81 = getelementptr inbounds [4 x float], [4 x float]* %2, i64 %29, i64 3
store float %80, float* %81
%82 = add nuw nsw i64 %29, 1
%83 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 %82, i64 0
%84 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 %82, i64 1
%85 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 %82, i64 2
%86 = getelementptr inbounds [4 x float], [4 x float]* %0, i64 %82, i64 3
%87 = icmp eq i64 %82, 4
br i1 %87, label %88, label %24
88: ; preds = %24
ret i32 undef
}