void _matrix_3x3_mul_3x3(float m1[3][3], float m2[3][3], float m3[3][3]) {
  int i;
  for (i = 0; i < 3; ++i) {
    m3[0][i] = m1[0][0] * m2[0][i] + m1[0][1] * m2[1][i] + m1[0][2] * m2[2][i];
    m3[1][i] = m1[1][0] * m2[0][i] + m1[1][1] * m2[1][i] + m1[1][2] * m2[2][i];
    m3[2][i] = m1[2][0] * m2[0][i] + m1[2][1] * m2[1][i] + m1[2][2] * m2[2][i];
  }
}


x86 assembly (clang, -O0):
<_matrix_3x3_mul_3x3>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x3,-0x1c(%rbp)
jge    1253 <_matrix_3x3_mul_3x3+0x153>
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mulss  0xc(%rax,%rcx,4),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  0x18(%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  0x10(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mulss  0xc(%rax,%rcx,4),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movss  0x14(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  0x18(%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,0xc(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movss  0x18(%rax),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  0x1c(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mulss  0xc(%rax,%rcx,4),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movss  0x20(%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  0x18(%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,0x18(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1117 <_matrix_3x3_mul_3x3+0x17>
pop    %rbp
ret

x86 assembly (clang, -O2):
<_matrix_3x3_mul_3x3>:
movss  (%rdi),%xmm0
movss  0x4(%rdi),%xmm1
mulss  0xc(%rsi),%xmm1
mulss  (%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rdi),%xmm1
mulss  0x18(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rdx)
movss  0xc(%rdi),%xmm0
movss  0x10(%rdi),%xmm1
mulss  0xc(%rsi),%xmm1
mulss  (%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x14(%rdi),%xmm1
mulss  0x18(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0xc(%rdx)
movss  0x18(%rdi),%xmm0
movss  0x1c(%rdi),%xmm1
mulss  0xc(%rsi),%xmm1
mulss  (%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x20(%rdi),%xmm1
mulss  0x18(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x18(%rdx)
movss  (%rdi),%xmm0
movss  0x4(%rdi),%xmm1
mulss  0x10(%rsi),%xmm1
mulss  0x4(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rdi),%xmm1
mulss  0x1c(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x4(%rdx)
movss  0xc(%rdi),%xmm0
movss  0x10(%rdi),%xmm1
mulss  0x10(%rsi),%xmm1
mulss  0x4(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x14(%rdi),%xmm1
mulss  0x1c(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x10(%rdx)
movss  0x18(%rdi),%xmm0
movss  0x1c(%rdi),%xmm1
mulss  0x10(%rsi),%xmm1
mulss  0x4(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x20(%rdi),%xmm1
mulss  0x1c(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x1c(%rdx)
movss  (%rdi),%xmm0
movss  0x4(%rdi),%xmm1
mulss  0x14(%rsi),%xmm1
mulss  0x8(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rdi),%xmm1
mulss  0x20(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x8(%rdx)
movss  0xc(%rdi),%xmm0
movss  0x10(%rdi),%xmm1
mulss  0x14(%rsi),%xmm1
mulss  0x8(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x14(%rdi),%xmm1
mulss  0x20(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x14(%rdx)
movss  0x18(%rdi),%xmm0
movss  0x1c(%rdi),%xmm1
mulss  0x14(%rsi),%xmm1
mulss  0x8(%rsi),%xmm0
addss  %xmm1,%xmm0
movss  0x20(%rdi),%xmm1
mulss  0x20(%rsi),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x20(%rdx)
ret

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

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