void invert_mtx44(float *m) {
  float t[16];
  int h, i, v;
  float tx = -m[3], ty = -m[7], tz = -m[11];
  for (h = 0; h < 3; h++)
    for (v = 0; v < 3; v++)
      t[h + (v << 2)] = m[v + (h << 2)];
  for (i = 0; i < 11; i++)
    m[i] = t[i];
  m[3] = tx * m[0] + ty * m[1] + tz * m[2];
  m[7] = tx * m[4] + ty * m[5] + tz * m[6];
  m[11] = tx * m[8] + ty * m[9] + tz * m[10];
}


x86 assembly (clang, -O0):
<invert_mtx44>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
movss  0x1c(%rax),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x64(%rbp)
mov    -0x8(%rbp),%rax
movss  0x2c(%rax),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x68(%rbp)
movl   $0x0,-0x54(%rbp)
cmpl   $0x3,-0x54(%rbp)
jge    11c6 <invert_mtx44+0xc6>
movl   $0x0,-0x5c(%rbp)
cmpl   $0x3,-0x5c(%rbp)
jge    11b3 <invert_mtx44+0xb3>
mov    -0x8(%rbp),%rax
mov    -0x5c(%rbp),%ecx
mov    -0x54(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x54(%rbp),%eax
mov    -0x5c(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
cltq
movss  %xmm0,-0x50(%rbp,%rax,4)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1171 <invert_mtx44+0x71>
jmp    11b8 <invert_mtx44+0xb8>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1160 <invert_mtx44+0x60>
movl   $0x0,-0x58(%rbp)
cmpl   $0xb,-0x58(%rbp)
jge    11fc <invert_mtx44+0xfc>
movslq -0x58(%rbp),%rax
movss  -0x50(%rbp,%rax,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x58(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11cd <invert_mtx44+0xcd>
movss  -0x60(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm2
movss  -0x64(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x4(%rax),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0xc(%rax)
movss  -0x60(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  0x10(%rax),%xmm2
movss  -0x64(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x14(%rax),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x18(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x1c(%rax)
movss  -0x60(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movss  0x20(%rax),%xmm2
movss  -0x64(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x24(%rax),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x28(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x2c(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<invert_mtx44>:
movss  0x1c(%rdi),%xmm4
movss  0x2c(%rdi),%xmm0
movss  0x4(%rdi),%xmm3
movss  0x8(%rdi),%xmm8
movss  0x10(%rdi),%xmm7
movss  0x18(%rdi),%xmm9
movss  0x20(%rdi),%xmm1
movss  0x24(%rdi),%xmm6
movss  %xmm7,0x4(%rdi)
movss  %xmm1,0x8(%rdi)
movss  %xmm3,0x10(%rdi)
movss  %xmm6,0x18(%rdi)
movss  %xmm8,0x20(%rdi)
movss  %xmm9,0x24(%rdi)
movss  0xc(%rdi),%xmm5
xorps  0xeaa(%rip),%xmm5
mulss  %xmm4,%xmm7
movss  (%rdi),%xmm2
mulss  %xmm5,%xmm2
subss  %xmm7,%xmm2
mulss  %xmm0,%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0xc(%rdi)
mulss  %xmm4,%xmm9
mulss  0x14(%rdi),%xmm4
mulss  %xmm5,%xmm3
subss  %xmm4,%xmm3
mulss  %xmm0,%xmm6
subss  %xmm6,%xmm3
movss  %xmm3,0x1c(%rdi)
mulss  %xmm8,%xmm5
subss  %xmm9,%xmm5
mulss  0x28(%rdi),%xmm0
subss  %xmm0,%xmm5
movss  %xmm5,0x2c(%rdi)
ret

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

LLVM IR (-O2):
define void @invert_mtx44(float* nocapture %0) {
%2 = getelementptr inbounds float, float* %0, i64 3
%3 = load float, float* %2
%4 = getelementptr inbounds float, float* %0, i64 7
%5 = load float, float* %4
%6 = getelementptr inbounds float, float* %0, i64 11
%7 = load float, float* %6
%8 = load float, float* %0
%9 = getelementptr inbounds float, float* %0, i64 1
%10 = load float, float* %9
%11 = getelementptr inbounds float, float* %0, i64 2
%12 = load float, float* %11
%13 = getelementptr inbounds float, float* %0, i64 4
%14 = load float, float* %13
%15 = getelementptr inbounds float, float* %0, i64 5
%16 = load float, float* %15
%17 = getelementptr inbounds float, float* %0, i64 6
%18 = load float, float* %17
%19 = getelementptr inbounds float, float* %0, i64 8
%20 = load float, float* %19
%21 = getelementptr inbounds float, float* %0, i64 9
%22 = load float, float* %21
%23 = getelementptr inbounds float, float* %0, i64 10
%24 = load float, float* %23
%25 = fneg float %5
store float %14, float* %9
store float %20, float* %11
store float %10, float* %13
store float %22, float* %17
store float %12, float* %19
store float %18, float* %21
%26 = fneg float %3
%27 = fneg float %7
%28 = fmul float %14, %25
%29 = tail call float @llvm.fmuladd.f32(float %26, float %8, float %28)
%30 = tail call float @llvm.fmuladd.f32(float %27, float %20, float %29)
store float %30, float* %2
%31 = fmul float %16, %25
%32 = tail call float @llvm.fmuladd.f32(float %26, float %10, float %31)
%33 = tail call float @llvm.fmuladd.f32(float %27, float %22, float %32)
store float %33, float* %4
%34 = fmul float %18, %25
%35 = tail call float @llvm.fmuladd.f32(float %26, float %12, float %34)
%36 = tail call float @llvm.fmuladd.f32(float %27, float %24, float %35)
store float %36, float* %6
ret void
}