void MatrixMultiplyC44f(const float *b, float *m) {
  int i;
  for (i = 0; i < 4; i++) {
    float ai0 = m[(0 << 2) + i], ai1 = m[(1 << 2) + i], ai2 = m[(2 << 2) + i],
          ai3 = m[(3 << 2) + i];
    m[(0 << 2) + i] = ai0 * b[(0 << 2) + 0] + ai1 * b[(0 << 2) + 1] +
                      ai2 * b[(0 << 2) + 2] + ai3 * b[(0 << 2) + 3];
    m[(1 << 2) + i] = ai0 * b[(1 << 2) + 0] + ai1 * b[(1 << 2) + 1] +
                      ai2 * b[(1 << 2) + 2] + ai3 * b[(1 << 2) + 3];
    m[(2 << 2) + i] = ai0 * b[(2 << 2) + 0] + ai1 * b[(2 << 2) + 1] +
                      ai2 * b[(2 << 2) + 2] + ai3 * b[(2 << 2) + 3];
    m[(3 << 2) + i] = ai0 * b[(3 << 2) + 0] + ai1 * b[(3 << 2) + 1] +
                      ai2 * b[(3 << 2) + 2] + ai3 * b[(3 << 2) + 3];
  }
}


x86 assembly (clang, -O0):
<MatrixMultiplyC44f>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
cmpl   $0x4,-0x14(%rbp)
jge    130e <MatrixMultiplyC44f+0x20e>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x8,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0xc,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm2
movss  -0x1c(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulss  0x4(%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x10(%rax),%xmm2
movss  -0x1c(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulss  0x14(%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  0x18(%rax),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x1c(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x20(%rax),%xmm2
movss  -0x1c(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulss  0x24(%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  0x28(%rax),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x2c(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x8,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x30(%rax),%xmm2
movss  -0x1c(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulss  0x34(%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  0x38(%rax),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x3c(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0xc,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1113 <MatrixMultiplyC44f+0x13>
pop    %rbp
ret

x86 assembly (clang, -O2):
<MatrixMultiplyC44f>:
mov    $0xfffffffffffffffc,%rax
nopw   0x0(%rax,%rax,1)
movss  0x10(%rsi,%rax,4),%xmm0
movss  0x20(%rsi,%rax,4),%xmm3
movss  0x30(%rsi,%rax,4),%xmm1
movss  0x40(%rsi,%rax,4),%xmm2
movss  0x4(%rdi),%xmm4
mulss  %xmm3,%xmm4
movss  (%rdi),%xmm5
mulss  %xmm0,%xmm5
addss  %xmm4,%xmm5
movss  0x8(%rdi),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm5,%xmm4
movss  0xc(%rdi),%xmm5
mulss  %xmm2,%xmm5
addss  %xmm4,%xmm5
movss  %xmm5,0x10(%rsi,%rax,4)
movss  0x14(%rdi),%xmm4
mulss  %xmm3,%xmm4
movss  0x10(%rdi),%xmm5
mulss  %xmm0,%xmm5
addss  %xmm4,%xmm5
movss  0x18(%rdi),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm5,%xmm4
movss  0x1c(%rdi),%xmm5
mulss  %xmm2,%xmm5
addss  %xmm4,%xmm5
movss  %xmm5,0x20(%rsi,%rax,4)
movss  0x24(%rdi),%xmm4
mulss  %xmm3,%xmm4
movss  0x20(%rdi),%xmm5
mulss  %xmm0,%xmm5
addss  %xmm4,%xmm5
movss  0x28(%rdi),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm5,%xmm4
movss  0x2c(%rdi),%xmm5
mulss  %xmm2,%xmm5
addss  %xmm4,%xmm5
movss  %xmm5,0x30(%rsi,%rax,4)
mulss  0x34(%rdi),%xmm3
mulss  0x30(%rdi),%xmm0
addss  %xmm3,%xmm0
mulss  0x38(%rdi),%xmm1
mulss  0x3c(%rdi),%xmm2
addss  %xmm0,%xmm1
addss  %xmm1,%xmm2
movss  %xmm2,0x40(%rsi,%rax,4)
inc    %rax
jne    1110 <MatrixMultiplyC44f+0x10>
ret

LLVM IR (-O0):
define void @MatrixMultiplyC44f(float* %0, float* %1) {
%3 = alloca float*
%4 = alloca float*
%5 = alloca i32
%6 = alloca float
%7 = alloca float
%8 = alloca float
%9 = alloca float
store float* %0, float** %3
store float* %1, float** %4
store i32 0, i32* %5
br label %10
10: ; preds = %138, %2
%11 = load i32, i32* %5
%12 = icmp slt i32 %11, 4
br i1 %12, label %13, label %141
13: ; preds = %10
%14 = load float*, float** %4
%15 = load i32, i32* %5
%16 = add nsw i32 0, %15
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds float, float* %14, i64 %17
%19 = load float, float* %18
store float %19, float* %6
%20 = load float*, float** %4
%21 = load i32, i32* %5
%22 = add nsw i32 4, %21
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds float, float* %20, i64 %23
%25 = load float, float* %24
store float %25, float* %7
%26 = load float*, float** %4
%27 = load i32, i32* %5
%28 = add nsw i32 8, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %26, i64 %29
%31 = load float, float* %30
store float %31, float* %8
%32 = load float*, float** %4
%33 = load i32, i32* %5
%34 = add nsw i32 12, %33
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds float, float* %32, i64 %35
%37 = load float, float* %36
store float %37, float* %9
%38 = load float, float* %6
%39 = load float*, float** %3
%40 = getelementptr inbounds float, float* %39, i64 0
%41 = load float, float* %40
%42 = load float, float* %7
%43 = load float*, float** %3
%44 = getelementptr inbounds float, float* %43, i64 1
%45 = load float, float* %44
%46 = fmul float %42, %45
%47 = call float @llvm.fmuladd.f32(float %38, float %41, float %46)
%48 = load float, float* %8
%49 = load float*, float** %3
%50 = getelementptr inbounds float, float* %49, i64 2
%51 = load float, float* %50
%52 = call float @llvm.fmuladd.f32(float %48, float %51, float %47)
%53 = load float, float* %9
%54 = load float*, float** %3
%55 = getelementptr inbounds float, float* %54, i64 3
%56 = load float, float* %55
%57 = call float @llvm.fmuladd.f32(float %53, float %56, float %52)
%58 = load float*, float** %4
%59 = load i32, i32* %5
%60 = add nsw i32 0, %59
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds float, float* %58, i64 %61
store float %57, float* %62
%63 = load float, float* %6
%64 = load float*, float** %3
%65 = getelementptr inbounds float, float* %64, i64 4
%66 = load float, float* %65
%67 = load float, float* %7
%68 = load float*, float** %3
%69 = getelementptr inbounds float, float* %68, i64 5
%70 = load float, float* %69
%71 = fmul float %67, %70
%72 = call float @llvm.fmuladd.f32(float %63, float %66, float %71)
%73 = load float, float* %8
%74 = load float*, float** %3
%75 = getelementptr inbounds float, float* %74, i64 6
%76 = load float, float* %75
%77 = call float @llvm.fmuladd.f32(float %73, float %76, float %72)
%78 = load float, float* %9
%79 = load float*, float** %3
%80 = getelementptr inbounds float, float* %79, i64 7
%81 = load float, float* %80
%82 = call float @llvm.fmuladd.f32(float %78, float %81, float %77)
%83 = load float*, float** %4
%84 = load i32, i32* %5
%85 = add nsw i32 4, %84
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds float, float* %83, i64 %86
store float %82, float* %87
%88 = load float, float* %6
%89 = load float*, float** %3
%90 = getelementptr inbounds float, float* %89, i64 8
%91 = load float, float* %90
%92 = load float, float* %7
%93 = load float*, float** %3
%94 = getelementptr inbounds float, float* %93, i64 9
%95 = load float, float* %94
%96 = fmul float %92, %95
%97 = call float @llvm.fmuladd.f32(float %88, float %91, float %96)
%98 = load float, float* %8
%99 = load float*, float** %3
%100 = getelementptr inbounds float, float* %99, i64 10
%101 = load float, float* %100
%102 = call float @llvm.fmuladd.f32(float %98, float %101, float %97)
%103 = load float, float* %9
%104 = load float*, float** %3
%105 = getelementptr inbounds float, float* %104, i64 11
%106 = load float, float* %105
%107 = call float @llvm.fmuladd.f32(float %103, float %106, float %102)
%108 = load float*, float** %4
%109 = load i32, i32* %5
%110 = add nsw i32 8, %109
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds float, float* %108, i64 %111
store float %107, float* %112
%113 = load float, float* %6
%114 = load float*, float** %3
%115 = getelementptr inbounds float, float* %114, i64 12
%116 = load float, float* %115
%117 = load float, float* %7
%118 = load float*, float** %3
%119 = getelementptr inbounds float, float* %118, i64 13
%120 = load float, float* %119
%121 = fmul float %117, %120
%122 = call float @llvm.fmuladd.f32(float %113, float %116, float %121)
%123 = load float, float* %8
%124 = load float*, float** %3
%125 = getelementptr inbounds float, float* %124, i64 14
%126 = load float, float* %125
%127 = call float @llvm.fmuladd.f32(float %123, float %126, float %122)
%128 = load float, float* %9
%129 = load float*, float** %3
%130 = getelementptr inbounds float, float* %129, i64 15
%131 = load float, float* %130
%132 = call float @llvm.fmuladd.f32(float %128, float %131, float %127)
%133 = load float*, float** %4
%134 = load i32, i32* %5
%135 = add nsw i32 12, %134
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds float, float* %133, i64 %136
store float %132, float* %137
br label %138
138: ; preds = %13
%139 = load i32, i32* %5
%140 = add nsw i32 %139, 1
store i32 %140, i32* %5
br label %10
141: ; preds = %10
ret void
}

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