void colmult(float *tgt, float *m1, const float *m2, int n) {
  while (n--)
    *(tgt++) = (*(m1++)) * (*(m2++));
}


x86 assembly (clang, -O0):
<colmult>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x1c(%rbp)
cmp    $0x0,%eax
je     1165 <colmult+0x65>
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x10(%rbp)
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
mulss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  %xmm0,(%rax)
jmp    1113 <colmult+0x13>
pop    %rbp
ret

x86 assembly (clang, -O2):
<colmult>:
push   %rbp
push   %rbx
test   %ecx,%ecx
je     12ab <colmult+0x1ab>
lea    -0x1(%rcx),%r9d
cmp    $0x7,%r9d
jae    11b2 <colmult+0xb2>
lea    -0x1(%rcx),%r8d
mov    %ecx,%r9d
and    $0x3,%r9d
je     1152 <colmult+0x52>
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movss  (%rsi),%xmm0
add    $0x4,%rsi
mulss  (%rdx),%xmm0
add    $0x4,%rdx
movss  %xmm0,(%rdi)
add    $0x4,%rdi
add    $0x1,%eax
cmp    %eax,%r9d
jne    1130 <colmult+0x30>
sub    %eax,%ecx
cmp    $0x3,%r8d
jb     12ab <colmult+0x1ab>
mov    %ecx,%eax
xor    %ecx,%ecx
movss  (%rsi,%rcx,4),%xmm0
mulss  (%rdx,%rcx,4),%xmm0
movss  %xmm0,(%rdi,%rcx,4)
movss  0x4(%rsi,%rcx,4),%xmm0
mulss  0x4(%rdx,%rcx,4),%xmm0
movss  %xmm0,0x4(%rdi,%rcx,4)
movss  0x8(%rsi,%rcx,4),%xmm0
mulss  0x8(%rdx,%rcx,4),%xmm0
movss  %xmm0,0x8(%rdi,%rcx,4)
movss  0xc(%rsi,%rcx,4),%xmm0
mulss  0xc(%rdx,%rcx,4),%xmm0
movss  %xmm0,0xc(%rdi,%rcx,4)
add    $0x4,%rcx
cmp    %ecx,%eax
jne    1160 <colmult+0x60>
jmp    12ab <colmult+0x1ab>
lea    (%rsi,%r9,4),%rax
add    $0x4,%rax
cmp    %rdi,%rax
seta   %r10b
lea    (%rdi,%r9,4),%rax
add    $0x4,%rax
cmp    %rsi,%rax
seta   %bpl
lea    (%rdx,%r9,4),%rbx
add    $0x4,%rbx
cmp    %rdi,%rbx
seta   %r11b
cmp    %rdx,%rax
seta   %r8b
test   %bpl,%r10b
jne    1118 <colmult+0x18>
and    %r8b,%r11b
jne    1118 <colmult+0x18>
add    $0x1,%r9
mov    %r9,%r8
and    $0xfffffffffffffff8,%r8
lea    -0x8(%r8),%rax
mov    %rax,%r10
shr    $0x3,%r10
add    $0x1,%r10
test   %rax,%rax
je     12c2 <colmult+0x1c2>
mov    %r10,%r11
and    $0xfffffffffffffffe,%r11
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movups (%rsi,%rax,4),%xmm0
movups 0x10(%rsi,%rax,4),%xmm1
movups (%rdx,%rax,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x10(%rdx,%rax,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,(%rdi,%rax,4)
movups %xmm0,0x10(%rdi,%rax,4)
movups 0x20(%rsi,%rax,4),%xmm0
movups 0x30(%rsi,%rax,4),%xmm1
movups 0x20(%rdx,%rax,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x30(%rdx,%rax,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,0x20(%rdi,%rax,4)
movups %xmm0,0x30(%rdi,%rax,4)
add    $0x10,%rax
add    $0xfffffffffffffffe,%r11
jne    1230 <colmult+0x130>
test   $0x1,%r10b
je     12a6 <colmult+0x1a6>
movups (%rsi,%rax,4),%xmm0
movups 0x10(%rsi,%rax,4),%xmm1
movups (%rdx,%rax,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x10(%rdx,%rax,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,(%rdi,%rax,4)
movups %xmm0,0x10(%rdi,%rax,4)
cmp    %r8,%r9
jne    12ae <colmult+0x1ae>
pop    %rbx
pop    %rbp
ret
sub    %r8d,%ecx
lea    (%rdx,%r8,4),%rdx
lea    (%rsi,%r8,4),%rsi
lea    (%rdi,%r8,4),%rdi
jmp    1118 <colmult+0x18>
xor    %eax,%eax
test   $0x1,%r10b
jne    1285 <colmult+0x185>
jmp    12a6 <colmult+0x1a6>

LLVM IR (-O0):
define void @colmult(float* %0, float* %1, float* %2, i32 %3) {
%5 = alloca float*
%6 = alloca float*
%7 = alloca float*
%8 = alloca i32
store float* %0, float** %5
store float* %1, float** %6
store float* %2, float** %7
store i32 %3, i32* %8
br label %9
9: ; preds = %13, %4
%10 = load i32, i32* %8
%11 = add nsw i32 %10, -1
store i32 %11, i32* %8
%12 = icmp ne i32 %10, 0
br i1 %12, label %13, label %23
13: ; preds = %9
%14 = load float*, float** %6
%15 = getelementptr inbounds float, float* %14, i32 1
store float* %15, float** %6
%16 = load float, float* %14
%17 = load float*, float** %7
%18 = getelementptr inbounds float, float* %17, i32 1
store float* %18, float** %7
%19 = load float, float* %17
%20 = fmul float %16, %19
%21 = load float*, float** %5
%22 = getelementptr inbounds float, float* %21, i32 1
store float* %22, float** %5
store float %20, float* %21
br label %9
23: ; preds = %9
ret void
}

LLVM IR (-O2):
define void @colmult(float* nocapture writeonly %0, float* nocapture %1, float* nocapture %2, i32 %3) {
%5 = icmp eq i32 %3, 0
br i1 %5, label %166, label %6
6: ; preds = %4
%7 = add i32 %3, -1
%8 = zext i32 %7 to i64
%9 = add nuw nsw i64 %8, 1
%10 = icmp ult i32 %7, 7
br i1 %10, label %106, label %11
11: ; preds = %6
%12 = add i32 %3, -1
%13 = zext i32 %12 to i64
%14 = add nuw nsw i64 %13, 1
%15 = getelementptr float, float* %0, i64 %14
%16 = getelementptr float, float* %1, i64 %14
%17 = getelementptr float, float* %2, i64 %14
%18 = icmp ugt float* %16, %0
%19 = icmp ugt float* %15, %1
%20 = and i1 %18, %19
%21 = icmp ugt float* %17, %0
%22 = icmp ugt float* %15, %2
%23 = and i1 %21, %22
%24 = or i1 %20, %23
br i1 %24, label %106, label %25
25: ; preds = %11
%26 = and i64 %9, 8589934584
%27 = trunc i64 %26 to i32
%28 = sub i32 %3, %27
%29 = getelementptr float, float* %2, i64 %26
%30 = getelementptr float, float* %1, i64 %26
%31 = getelementptr float, float* %0, i64 %26
%32 = add nsw i64 %26, -8
%33 = lshr exact i64 %32, 3
%34 = add nuw nsw i64 %33, 1
%35 = and i64 %34, 1
%36 = icmp eq i64 %32, 0
br i1 %36, label %82, label %37
37: ; preds = %25
%38 = and i64 %34, 4611686018427387902
br label %39
39: ; preds = %39, %37
%40 = phi i64 [ 0, %37 ], [ %79, %39 ]
%41 = phi i64 [ 0, %37 ], [ %80, %39 ]
%42 = getelementptr float, float* %2, i64 %40
%43 = getelementptr float, float* %1, i64 %40
%44 = getelementptr float, float* %0, i64 %40
%45 = bitcast float* %43 to <4 x float>*
%46 = load <4 x float>, <4 x float>* %45
%47 = getelementptr float, float* %43, i64 4
%48 = bitcast float* %47 to <4 x float>*
%49 = load <4 x float>, <4 x float>* %48
%50 = bitcast float* %42 to <4 x float>*
%51 = load <4 x float>, <4 x float>* %50
%52 = getelementptr float, float* %42, i64 4
%53 = bitcast float* %52 to <4 x float>*
%54 = load <4 x float>, <4 x float>* %53
%55 = fmul <4 x float> %46, %51
%56 = fmul <4 x float> %49, %54
%57 = bitcast float* %44 to <4 x float>*
store <4 x float> %55, <4 x float>* %57
%58 = getelementptr float, float* %44, i64 4
%59 = bitcast float* %58 to <4 x float>*
store <4 x float> %56, <4 x float>* %59
%60 = or i64 %40, 8
%61 = getelementptr float, float* %2, i64 %60
%62 = getelementptr float, float* %1, i64 %60
%63 = getelementptr float, float* %0, i64 %60
%64 = bitcast float* %62 to <4 x float>*
%65 = load <4 x float>, <4 x float>* %64
%66 = getelementptr float, float* %62, i64 4
%67 = bitcast float* %66 to <4 x float>*
%68 = load <4 x float>, <4 x float>* %67
%69 = bitcast float* %61 to <4 x float>*
%70 = load <4 x float>, <4 x float>* %69
%71 = getelementptr float, float* %61, i64 4
%72 = bitcast float* %71 to <4 x float>*
%73 = load <4 x float>, <4 x float>* %72
%74 = fmul <4 x float> %65, %70
%75 = fmul <4 x float> %68, %73
%76 = bitcast float* %63 to <4 x float>*
store <4 x float> %74, <4 x float>* %76
%77 = getelementptr float, float* %63, i64 4
%78 = bitcast float* %77 to <4 x float>*
store <4 x float> %75, <4 x float>* %78
%79 = add nuw i64 %40, 16
%80 = add i64 %41, 2
%81 = icmp eq i64 %80, %38
br i1 %81, label %82, label %39
82: ; preds = %39, %25
%83 = phi i64 [ 0, %25 ], [ %79, %39 ]
%84 = icmp eq i64 %35, 0
br i1 %84, label %104, label %85
85: ; preds = %82
%86 = getelementptr float, float* %2, i64 %83
%87 = getelementptr float, float* %1, i64 %83
%88 = getelementptr float, float* %0, i64 %83
%89 = bitcast float* %87 to <4 x float>*
%90 = load <4 x float>, <4 x float>* %89
%91 = getelementptr float, float* %87, i64 4
%92 = bitcast float* %91 to <4 x float>*
%93 = load <4 x float>, <4 x float>* %92
%94 = bitcast float* %86 to <4 x float>*
%95 = load <4 x float>, <4 x float>* %94
%96 = getelementptr float, float* %86, i64 4
%97 = bitcast float* %96 to <4 x float>*
%98 = load <4 x float>, <4 x float>* %97
%99 = fmul <4 x float> %90, %95
%100 = fmul <4 x float> %93, %98
%101 = bitcast float* %88 to <4 x float>*
store <4 x float> %99, <4 x float>* %101
%102 = getelementptr float, float* %88, i64 4
%103 = bitcast float* %102 to <4 x float>*
store <4 x float> %100, <4 x float>* %103
br label %104
104: ; preds = %82, %85
%105 = icmp eq i64 %9, %26
br i1 %105, label %166, label %106
106: ; preds = %11, %6, %104
%107 = phi i32 [ %3, %11 ], [ %3, %6 ], [ %28, %104 ]
%108 = phi float* [ %2, %11 ], [ %2, %6 ], [ %29, %104 ]
%109 = phi float* [ %1, %11 ], [ %1, %6 ], [ %30, %104 ]
%110 = phi float* [ %0, %11 ], [ %0, %6 ], [ %31, %104 ]
%111 = add i32 %107, -1
%112 = and i32 %107, 3
%113 = icmp eq i32 %112, 0
br i1 %113, label %129, label %114
114: ; preds = %106, %114
%115 = phi i32 [ %120, %114 ], [ %107, %106 ]
%116 = phi float* [ %123, %114 ], [ %108, %106 ]
%117 = phi float* [ %121, %114 ], [ %109, %106 ]
%118 = phi float* [ %126, %114 ], [ %110, %106 ]
%119 = phi i32 [ %127, %114 ], [ 0, %106 ]
%120 = add nsw i32 %115, -1
%121 = getelementptr inbounds float, float* %117, i64 1
%122 = load float, float* %117
%123 = getelementptr inbounds float, float* %116, i64 1
%124 = load float, float* %116
%125 = fmul float %122, %124
%126 = getelementptr inbounds float, float* %118, i64 1
store float %125, float* %118
%127 = add i32 %119, 1
%128 = icmp eq i32 %127, %112
br i1 %128, label %129, label %114
129: ; preds = %114, %106
%130 = phi i32 [ %107, %106 ], [ %120, %114 ]
%131 = phi float* [ %108, %106 ], [ %123, %114 ]
%132 = phi float* [ %109, %106 ], [ %121, %114 ]
%133 = phi float* [ %110, %106 ], [ %126, %114 ]
%134 = icmp ult i32 %111, 3
br i1 %134, label %166, label %135
135: ; preds = %129, %135
%136 = phi i32 [ %158, %135 ], [ %130, %129 ]
%137 = phi float* [ %161, %135 ], [ %131, %129 ]
%138 = phi float* [ %159, %135 ], [ %132, %129 ]
%139 = phi float* [ %164, %135 ], [ %133, %129 ]
%140 = getelementptr inbounds float, float* %138, i64 1
%141 = load float, float* %138
%142 = getelementptr inbounds float, float* %137, i64 1
%143 = load float, float* %137
%144 = fmul float %141, %143
%145 = getelementptr inbounds float, float* %139, i64 1
store float %144, float* %139
%146 = getelementptr inbounds float, float* %138, i64 2
%147 = load float, float* %140
%148 = getelementptr inbounds float, float* %137, i64 2
%149 = load float, float* %142
%150 = fmul float %147, %149
%151 = getelementptr inbounds float, float* %139, i64 2
store float %150, float* %145
%152 = getelementptr inbounds float, float* %138, i64 3
%153 = load float, float* %146
%154 = getelementptr inbounds float, float* %137, i64 3
%155 = load float, float* %148
%156 = fmul float %153, %155
%157 = getelementptr inbounds float, float* %139, i64 3
store float %156, float* %151
%158 = add nsw i32 %136, -4
%159 = getelementptr inbounds float, float* %138, i64 4
%160 = load float, float* %152
%161 = getelementptr inbounds float, float* %137, i64 4
%162 = load float, float* %154
%163 = fmul float %160, %162
%164 = getelementptr inbounds float, float* %139, i64 4
store float %163, float* %157
%165 = icmp eq i32 %158, 0
br i1 %165, label %166, label %135
166: ; preds = %129, %135, %104, %4
ret void
}