void MatrixTransformC44f4f(float *m, float *q, float *p) {
  register float q0 = *q, q1 = *(q + 1), q2 = *(q + 2);
  p[0] = m[0] * q0 + m[4] * q1 + m[8] * q2 + m[12];
  p[1] = m[1] * q0 + m[5] * q1 + m[9] * q2 + m[13];
  p[2] = m[2] * q0 + m[6] * q1 + m[10] * q2 + m[14];
  p[3] = m[3] * q0 + m[7] * q1 + m[11] * q2 + m[15];
}


x86 assembly (clang, -O0):
<MatrixTransformC44f4f>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
movss  -0x1c(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movss  0x10(%rax),%xmm0
mulss  -0x20(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movss  0x20(%rax),%xmm0
movss  -0x24(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
addss  0x30(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm1
movss  -0x1c(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movss  0x14(%rax),%xmm0
mulss  -0x20(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movss  0x24(%rax),%xmm0
movss  -0x24(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
addss  0x34(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm1
movss  -0x1c(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movss  0x18(%rax),%xmm0
mulss  -0x20(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movss  0x28(%rax),%xmm0
movss  -0x24(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
addss  0x38(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,0x8(%rax)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm1
movss  -0x1c(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movss  0x1c(%rax),%xmm0
mulss  -0x20(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movss  0x2c(%rax),%xmm0
movss  -0x24(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
addss  0x3c(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,0xc(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<MatrixTransformC44f4f>:
movss  (%rsi),%xmm1
movss  0x4(%rsi),%xmm2
movss  0x8(%rsi),%xmm0
movss  0x10(%rdi),%xmm3
mulss  %xmm2,%xmm3
movss  (%rdi),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm3,%xmm4
movss  0x20(%rdi),%xmm3
mulss  %xmm0,%xmm3
addss  %xmm4,%xmm3
addss  0x30(%rdi),%xmm3
movss  %xmm3,(%rdx)
movss  0x14(%rdi),%xmm3
mulss  %xmm2,%xmm3
movss  0x4(%rdi),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm3,%xmm4
movss  0x24(%rdi),%xmm3
mulss  %xmm0,%xmm3
addss  %xmm4,%xmm3
addss  0x34(%rdi),%xmm3
movss  %xmm3,0x4(%rdx)
movss  0x18(%rdi),%xmm3
mulss  %xmm2,%xmm3
movss  0x8(%rdi),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm3,%xmm4
movss  0x28(%rdi),%xmm3
mulss  %xmm0,%xmm3
addss  %xmm4,%xmm3
addss  0x38(%rdi),%xmm3
movss  %xmm3,0x8(%rdx)
mulss  0x1c(%rdi),%xmm2
mulss  0xc(%rdi),%xmm1
addss  %xmm2,%xmm1
mulss  0x2c(%rdi),%xmm0
addss  %xmm1,%xmm0
addss  0x3c(%rdi),%xmm0
movss  %xmm0,0xc(%rdx)
ret

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

LLVM IR (-O2):
define void @MatrixTransformC44f4f(float* nocapture %0, float* nocapture %1, float* nocapture writeonly %2) {
%4 = load float, float* %1
%5 = getelementptr inbounds float, float* %1, i64 1
%6 = load float, float* %5
%7 = getelementptr inbounds float, float* %1, i64 2
%8 = load float, float* %7
%9 = load float, float* %0
%10 = getelementptr inbounds float, float* %0, i64 4
%11 = load float, float* %10
%12 = fmul float %6, %11
%13 = tail call float @llvm.fmuladd.f32(float %9, float %4, float %12)
%14 = getelementptr inbounds float, float* %0, i64 8
%15 = load float, float* %14
%16 = tail call float @llvm.fmuladd.f32(float %15, float %8, float %13)
%17 = getelementptr inbounds float, float* %0, i64 12
%18 = load float, float* %17
%19 = fadd float %18, %16
store float %19, float* %2
%20 = getelementptr inbounds float, float* %0, i64 1
%21 = load float, float* %20
%22 = getelementptr inbounds float, float* %0, i64 5
%23 = load float, float* %22
%24 = fmul float %6, %23
%25 = tail call float @llvm.fmuladd.f32(float %21, float %4, float %24)
%26 = getelementptr inbounds float, float* %0, i64 9
%27 = load float, float* %26
%28 = tail call float @llvm.fmuladd.f32(float %27, float %8, float %25)
%29 = getelementptr inbounds float, float* %0, i64 13
%30 = load float, float* %29
%31 = fadd float %30, %28
%32 = getelementptr inbounds float, float* %2, i64 1
store float %31, float* %32
%33 = getelementptr inbounds float, float* %0, i64 2
%34 = load float, float* %33
%35 = getelementptr inbounds float, float* %0, i64 6
%36 = load float, float* %35
%37 = fmul float %6, %36
%38 = tail call float @llvm.fmuladd.f32(float %34, float %4, float %37)
%39 = getelementptr inbounds float, float* %0, i64 10
%40 = load float, float* %39
%41 = tail call float @llvm.fmuladd.f32(float %40, float %8, float %38)
%42 = getelementptr inbounds float, float* %0, i64 14
%43 = load float, float* %42
%44 = fadd float %43, %41
%45 = getelementptr inbounds float, float* %2, i64 2
store float %44, float* %45
%46 = getelementptr inbounds float, float* %0, i64 3
%47 = load float, float* %46
%48 = getelementptr inbounds float, float* %0, i64 7
%49 = load float, float* %48
%50 = fmul float %6, %49
%51 = tail call float @llvm.fmuladd.f32(float %47, float %4, float %50)
%52 = getelementptr inbounds float, float* %0, i64 11
%53 = load float, float* %52
%54 = tail call float @llvm.fmuladd.f32(float %53, float %8, float %51)
%55 = getelementptr inbounds float, float* %0, i64 15
%56 = load float, float* %55
%57 = fadd float %56, %54
%58 = getelementptr inbounds float, float* %2, i64 3
store float %57, float* %58
ret void
}