// File: /workspace/github_cpp_makefile_star2_O3/kylegranger[P]ImageScaling/filter16.cpp
// Function: DownfilterRight(color_t*, color_t*, color_t*, color_t*, color_t*)
void DownfilterRight(color_t *src1, color_t *src2, color_t *src3, color_t *src4, color_t *dst)
{
	// we mix 16 source pixels into one destination pixel
	int r = (src1[0].r)     + 3 * (src1[1].r + src1[2].r);
	r    += (src2[0].r) * 3 + 9 * (src2[1].r + src2[2].r);
	r    += (src3[0].r) * 3 + 9 * (src3[1].r + src3[2].r);
	r    += (src4[0].r)     + 3 * (src4[1].r + src4[2].r);
	dst->r = ((r+28)/56);
	int g = (src1[0].g)     + 3 * (src1[1].g + src1[2].g);
	g    += (src2[0].g) * 3 + 9 * (src2[1].g + src2[2].g);
	g    += (src3[0].g) * 3 + 9 * (src3[1].g + src3[2].g);
	g    += (src4[0].g)     + 3 * (src4[1].g + src4[2].g);
	dst->g = ((g+28)/56);
	int b = (src1[0].b)     + 3 * (src1[1].b + src1[2].b);
	b    += (src2[0].b) * 3 + 9 * (src2[1].b + src2[2].b);
	b    += (src3[0].b) * 3 + 9 * (src3[1].b + src3[2].b);
	b    += (src4[0].b)     + 3 * (src4[1].b + src4[2].b);
	dst->b = ((b+28)/56);
	int a = (src1[0].a)     + 3 * (src1[1].a + src1[2].a);
	a    += (src2[0].a) * 3 + 9 * (src2[1].a + src2[2].a);
	a    += (src3[0].a) * 3 + 9 * (src3[1].a + src3[2].a);
	a    += (src4[0].a)     + 3 * (src4[1].a + src4[2].a);
	dst->a = ((a+28)/56);
}

// Assembly:
movq (%rdi), %xmm1
pxor %xmm0, %xmm0
punpcklwd %xmm0, %xmm1    # xmm1 = xmm1[0],xmm0[0],xmm1[1],xmm0[1],xmm1[2],xmm0[2],xmm1[3],xmm0[3]
movq 0x8(%rdi), %xmm2
punpcklwd %xmm0, %xmm2    # xmm2 = xmm2[0],xmm0[0],xmm2[1],xmm0[1],xmm2[2],xmm0[2],xmm2[3],xmm0[3]
movq 0x10(%rdi), %xmm3
punpcklwd %xmm0, %xmm3    # xmm3 = xmm3[0],xmm0[0],xmm3[1],xmm0[1],xmm3[2],xmm0[2],xmm3[3],xmm0[3]
paddd %xmm2, %xmm3
movq (%rsi), %xmm2
punpcklwd %xmm0, %xmm2    # xmm2 = xmm2[0],xmm0[0],xmm2[1],xmm0[1],xmm2[2],xmm0[2],xmm2[3],xmm0[3]
movq 0x8(%rsi), %xmm4
punpcklwd %xmm0, %xmm4    # xmm4 = xmm4[0],xmm0[0],xmm4[1],xmm0[1],xmm4[2],xmm0[2],xmm4[3],xmm0[3]
movq 0x10(%rsi), %xmm5
punpcklwd %xmm0, %xmm5    # xmm5 = xmm5[0],xmm0[0],xmm5[1],xmm0[1],xmm5[2],xmm0[2],xmm5[3],xmm0[3]
paddd %xmm4, %xmm5
movq (%rdx), %xmm4
punpcklwd %xmm0, %xmm4    # xmm4 = xmm4[0],xmm0[0],xmm4[1],xmm0[1],xmm4[2],xmm0[2],xmm4[3],xmm0[3]
paddd %xmm2, %xmm4
paddd %xmm3, %xmm4
movq 0x8(%rdx), %xmm2
punpcklwd %xmm0, %xmm2    # xmm2 = xmm2[0],xmm0[0],xmm2[1],xmm0[1],xmm2[2],xmm0[2],xmm2[3],xmm0[3]
movq 0x10(%rdx), %xmm3
punpcklwd %xmm0, %xmm3    # xmm3 = xmm3[0],xmm0[0],xmm3[1],xmm0[1],xmm3[2],xmm0[2],xmm3[3],xmm0[3]
paddd %xmm2, %xmm3
paddd %xmm5, %xmm3
movq (%rcx), %xmm2
punpcklwd %xmm0, %xmm2    # xmm2 = xmm2[0],xmm0[0],xmm2[1],xmm0[1],xmm2[2],xmm0[2],xmm2[3],xmm0[3]
paddd %xmm1, %xmm2
movq 0x8(%rcx), %xmm1
punpcklwd %xmm0, %xmm1    # xmm1 = xmm1[0],xmm0[0],xmm1[1],xmm0[1],xmm1[2],xmm0[2],xmm1[3],xmm0[3]
movq 0x10(%rcx), %xmm5
punpcklwd %xmm0, %xmm5    # xmm5 = xmm5[0],xmm0[0],xmm5[1],xmm0[1],xmm5[2],xmm0[2],xmm5[3],xmm0[3]
paddd %xmm1, %xmm5
paddd %xmm4, %xmm5
movdqa %xmm5, %xmm0
paddd %xmm5, %xmm0
paddd %xmm5, %xmm0
movdqa %xmm3, %xmm1
pslld $0x3, %xmm1
paddd %xmm3, %xmm1
paddd %xmm2, %xmm1
paddd %xmm0, %xmm1
paddd 0x18504(%rip), %xmm1    # 0x1e1f0
psrld $0x3, %xmm1
movdqa 0x184d7(%rip), %xmm0    # 0x1e1d0
pshufd $0xf5, %xmm1, %xmm2     # xmm2 = xmm1[1,1,3,3]
pmuludq %xmm0, %xmm1
pshufd $0xed, %xmm1, %xmm1     # xmm1 = xmm1[1,3,2,3]
pmuludq %xmm0, %xmm2
pshufd $0xed, %xmm2, %xmm0     # xmm0 = xmm2[1,3,2,3]
punpcklwd %xmm0, %xmm1    # xmm1 = xmm1[0],xmm0[0],xmm1[1],xmm0[1],xmm1[2],xmm0[2],xmm1[3],xmm0[3]
pshufd $0xe8, %xmm1, %xmm0     # xmm0 = xmm1[0,2,2,3]
movq %xmm0, (%r8)
retq