void affine_xform(const float m[2][3], float x_ts, float y_ts, float *x_disp,
                  float *y_disp) {
  *x_disp = m[0][0] * x_ts + m[0][1] * y_ts + m[0][2];
  *y_disp = m[1][0] * x_ts + m[1][1] * y_ts + m[1][2];
}


x86 assembly (clang, -O0):
<affine_xform>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
movss  %xmm1,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  -0xc(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm1
mulss  -0x10(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
addss  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm0
movss  -0xc(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movss  0x10(%rax),%xmm1
mulss  -0x10(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
addss  0x14(%rax),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<affine_xform>:
movss  0x4(%rdi),%xmm2
mulss  %xmm1,%xmm2
movss  (%rdi),%xmm3
mulss  %xmm0,%xmm3
addss  %xmm2,%xmm3
addss  0x8(%rdi),%xmm3
movss  %xmm3,(%rsi)
mulss  0x10(%rdi),%xmm1
mulss  0xc(%rdi),%xmm0
addss  %xmm1,%xmm0
addss  0x14(%rdi),%xmm0
movss  %xmm0,(%rdx)
ret

LLVM IR (-O0):
define void @affine_xform([3 x float]* %0, float %1, float %2, float* %3, float* %4) {
%6 = alloca [3 x float]*
%7 = alloca float
%8 = alloca float
%9 = alloca float*
%10 = alloca float*
store [3 x float]* %0, [3 x float]** %6
store float %1, float* %7
store float %2, float* %8
store float* %3, float** %9
store float* %4, float** %10
%11 = load [3 x float]*, [3 x float]** %6
%12 = getelementptr inbounds [3 x float], [3 x float]* %11, i64 0
%13 = getelementptr inbounds [3 x float], [3 x float]* %12, i64 0, i64 0
%14 = load float, float* %13
%15 = load float, float* %7
%16 = load [3 x float]*, [3 x float]** %6
%17 = getelementptr inbounds [3 x float], [3 x float]* %16, i64 0
%18 = getelementptr inbounds [3 x float], [3 x float]* %17, i64 0, i64 1
%19 = load float, float* %18
%20 = load float, float* %8
%21 = fmul float %19, %20
%22 = call float @llvm.fmuladd.f32(float %14, float %15, float %21)
%23 = load [3 x float]*, [3 x float]** %6
%24 = getelementptr inbounds [3 x float], [3 x float]* %23, i64 0
%25 = getelementptr inbounds [3 x float], [3 x float]* %24, i64 0, i64 2
%26 = load float, float* %25
%27 = fadd float %22, %26
%28 = load float*, float** %9
store float %27, float* %28
%29 = load [3 x float]*, [3 x float]** %6
%30 = getelementptr inbounds [3 x float], [3 x float]* %29, i64 1
%31 = getelementptr inbounds [3 x float], [3 x float]* %30, i64 0, i64 0
%32 = load float, float* %31
%33 = load float, float* %7
%34 = load [3 x float]*, [3 x float]** %6
%35 = getelementptr inbounds [3 x float], [3 x float]* %34, i64 1
%36 = getelementptr inbounds [3 x float], [3 x float]* %35, i64 0, i64 1
%37 = load float, float* %36
%38 = load float, float* %8
%39 = fmul float %37, %38
%40 = call float @llvm.fmuladd.f32(float %32, float %33, float %39)
%41 = load [3 x float]*, [3 x float]** %6
%42 = getelementptr inbounds [3 x float], [3 x float]* %41, i64 1
%43 = getelementptr inbounds [3 x float], [3 x float]* %42, i64 0, i64 2
%44 = load float, float* %43
%45 = fadd float %40, %44
%46 = load float*, float** %10
store float %45, float* %46
ret void
}

LLVM IR (-O2):
define void @affine_xform([3 x float]* nocapture %0, float %1, float %2, float* nocapture writeonly %3, float* nocapture writeonly %4) {
%6 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 0, i64 0
%7 = load float, float* %6
%8 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 0, i64 1
%9 = load float, float* %8
%10 = fmul float %9, %2
%11 = tail call float @llvm.fmuladd.f32(float %7, float %1, float %10)
%12 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 0, i64 2
%13 = load float, float* %12
%14 = fadd float %13, %11
store float %14, float* %3
%15 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 1, i64 0
%16 = load float, float* %15
%17 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 1, i64 1
%18 = load float, float* %17
%19 = fmul float %18, %2
%20 = tail call float @llvm.fmuladd.f32(float %16, float %1, float %19)
%21 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 1, i64 2
%22 = load float, float* %21
%23 = fadd float %22, %20
store float %23, float* %4
ret void
}