applymatrix(lptr, mat, n) unsigned long *lptr;
float mat[4][4];
int n;
{
  int ir, ig, ib, r, g, b;
  unsigned char *cptr;
  cptr = (unsigned char *)lptr;
  while (n--) {
    ir = cptr[0];
    ig = cptr[1];
    ib = cptr[2];
    r = ir * mat[0][0] + ig * mat[1][0] + ib * mat[2][0] + mat[3][0];
    g = ir * mat[0][1] + ig * mat[1][1] + ib * mat[2][1] + mat[3][1];
    b = ir * mat[0][2] + ig * mat[1][2] + ib * mat[2][2] + mat[3][2];
    cptr[0] = CLAMP0255(r);
    cptr[1] = CLAMP0255(g);
    cptr[2] = CLAMP0255(b);
    cptr += 4;
  }
}


x86 assembly (clang, -O0):
<applymatrix>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x1c(%rbp)
cmp    $0x0,%eax
je     1285 <applymatrix+0x175>
mov    -0x40(%rbp),%rax
movzbl (%rax),%eax
mov    %eax,-0x20(%rbp)
mov    -0x40(%rbp),%rax
movzbl 0x1(%rax),%eax
mov    %eax,-0x24(%rbp)
mov    -0x40(%rbp),%rax
movzbl 0x2(%rax),%eax
mov    %eax,-0x28(%rbp)
cvtsi2ssl -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm2
cvtsi2ssl -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x10(%rax),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
cvtsi2ssl -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movss  0x20(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addss  0x30(%rax),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x2c(%rbp)
cvtsi2ssl -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm2
cvtsi2ssl -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x14(%rax),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
cvtsi2ssl -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movss  0x24(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addss  0x34(%rax),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
cvtsi2ssl -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm2
cvtsi2ssl -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x18(%rax),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
cvtsi2ssl -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movss  0x28(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addss  0x38(%rax),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%edi
mov    $0x0,%al
call   1030 <CLAMP0255@plt>
mov    %al,%cl
mov    -0x40(%rbp),%rax
mov    %cl,(%rax)
mov    -0x30(%rbp),%edi
mov    $0x0,%al
call   1030 <CLAMP0255@plt>
mov    %al,%cl
mov    -0x40(%rbp),%rax
mov    %cl,0x1(%rax)
mov    -0x34(%rbp),%edi
mov    $0x0,%al
call   1030 <CLAMP0255@plt>
mov    %al,%cl
mov    -0x40(%rbp),%rax
mov    %cl,0x2(%rax)
mov    -0x40(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x40(%rbp)
jmp    112b <applymatrix+0x1b>
mov    -0x4(%rbp),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<applymatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
test   %edx,%edx
je     1205 <applymatrix+0xf5>
mov    %rsi,%rbx
mov    %rdi,%r14
mov    %edx,%r12d
xor    %ebp,%ebp
xchg   %ax,%ax
movzbl (%r14,%rbp,4),%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movzbl 0x1(%r14,%rbp,4),%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
movzbl 0x2(%r14,%rbp,4),%eax
movss  0x10(%rbx),%xmm3
mulss  %xmm1,%xmm3
movss  (%rbx),%xmm4
mulss  %xmm0,%xmm4
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
addss  %xmm3,%xmm4
movss  0x20(%rbx),%xmm3
mulss  %xmm2,%xmm3
addss  %xmm4,%xmm3
addss  0x30(%rbx),%xmm3
cvttss2si %xmm3,%edi
movss  0x14(%rbx),%xmm3
mulss  %xmm1,%xmm3
movss  0x4(%rbx),%xmm4
mulss  %xmm0,%xmm4
addss  %xmm3,%xmm4
movss  0x24(%rbx),%xmm3
mulss  %xmm2,%xmm3
addss  %xmm4,%xmm3
addss  0x34(%rbx),%xmm3
cvttss2si %xmm3,%r15d
mulss  0x18(%rbx),%xmm1
mulss  0x8(%rbx),%xmm0
addss  %xmm1,%xmm0
mulss  0x28(%rbx),%xmm2
addss  %xmm0,%xmm2
addss  0x38(%rbx),%xmm2
cvttss2si %xmm2,%r13d
xor    %eax,%eax
call   1030 <CLAMP0255@plt>
mov    %al,(%r14,%rbp,4)
mov    %r15d,%edi
xor    %eax,%eax
call   1030 <CLAMP0255@plt>
mov    %al,0x1(%r14,%rbp,4)
mov    %r13d,%edi
xor    %eax,%eax
call   1030 <CLAMP0255@plt>
mov    %al,0x2(%r14,%rbp,4)
add    $0x1,%rbp
cmp    %ebp,%r12d
jne    1130 <applymatrix+0x20>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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