Index: third_party/libjpeg_turbo/simd/jcsamss2.asm |
=================================================================== |
--- third_party/libjpeg_turbo/simd/jcsamss2.asm (revision 0) |
+++ third_party/libjpeg_turbo/simd/jcsamss2.asm (revision 0) |
@@ -0,0 +1,351 @@ |
+; |
+; jcsamss2.asm - downsampling (SSE2) |
+; |
+; Copyright 2009 Pierre Ossman <ossman@cendio.se> for Cendio AB |
+; |
+; Based on |
+; x86 SIMD extension for IJG JPEG library |
+; Copyright (C) 1999-2006, MIYASAKA Masaru. |
+; For conditions of distribution and use, see copyright notice in jsimdext.inc |
+; |
+; This file should be assembled with NASM (Netwide Assembler), |
+; can *not* be assembled with Microsoft's MASM or any compatible |
+; assembler (including Borland's Turbo Assembler). |
+; NASM is available from http://nasm.sourceforge.net/ or |
+; http://sourceforge.net/project/showfiles.php?group_id=6208 |
+; |
+; [TAB8] |
+ |
+%include "jsimdext.inc" |
+ |
+; -------------------------------------------------------------------------- |
+ SECTION SEG_TEXT |
+ BITS 32 |
+; |
+; Downsample pixel values of a single component. |
+; This version handles the common case of 2:1 horizontal and 1:1 vertical, |
+; without smoothing. |
+; |
+; GLOBAL(void) |
+; jsimd_h2v1_downsample_sse2 (JDIMENSION image_width, int max_v_samp_factor, |
+; JDIMENSION v_samp_factor, JDIMENSION width_blocks, |
+; JSAMPARRAY input_data, JSAMPARRAY output_data); |
+; |
+ |
+%define img_width(b) (b)+8 ; JDIMENSION image_width |
+%define max_v_samp(b) (b)+12 ; int max_v_samp_factor |
+%define v_samp(b) (b)+16 ; JDIMENSION v_samp_factor |
+%define width_blks(b) (b)+20 ; JDIMENSION width_blocks |
+%define input_data(b) (b)+24 ; JSAMPARRAY input_data |
+%define output_data(b) (b)+28 ; JSAMPARRAY output_data |
+ |
+ align 16 |
+ global EXTN(jsimd_h2v1_downsample_sse2) |
+ |
+EXTN(jsimd_h2v1_downsample_sse2): |
+ push ebp |
+ mov ebp,esp |
+; push ebx ; unused |
+; push ecx ; need not be preserved |
+; push edx ; need not be preserved |
+ push esi |
+ push edi |
+ |
+ mov ecx, JDIMENSION [width_blks(ebp)] |
+ shl ecx,3 ; imul ecx,DCTSIZE (ecx = output_cols) |
+ jz near .return |
+ |
+ mov edx, JDIMENSION [img_width(ebp)] |
+ |
+ ; -- expand_right_edge |
+ |
+ push ecx |
+ shl ecx,1 ; output_cols * 2 |
+ sub ecx,edx |
+ jle short .expand_end |
+ |
+ mov eax, INT [max_v_samp(ebp)] |
+ test eax,eax |
+ jle short .expand_end |
+ |
+ cld |
+ mov esi, JSAMPARRAY [input_data(ebp)] ; input_data |
+ alignx 16,7 |
+.expandloop: |
+ push eax |
+ push ecx |
+ |
+ mov edi, JSAMPROW [esi] |
+ add edi,edx |
+ mov al, JSAMPLE [edi-1] |
+ |
+ rep stosb |
+ |
+ pop ecx |
+ pop eax |
+ |
+ add esi, byte SIZEOF_JSAMPROW |
+ dec eax |
+ jg short .expandloop |
+ |
+.expand_end: |
+ pop ecx ; output_cols |
+ |
+ ; -- h2v1_downsample |
+ |
+ mov eax, JDIMENSION [v_samp(ebp)] ; rowctr |
+ test eax,eax |
+ jle near .return |
+ |
+ mov edx, 0x00010000 ; bias pattern |
+ movd xmm7,edx |
+ pcmpeqw xmm6,xmm6 |
+ pshufd xmm7,xmm7,0x00 ; xmm7={0, 1, 0, 1, 0, 1, 0, 1} |
+ psrlw xmm6,BYTE_BIT ; xmm6={0xFF 0x00 0xFF 0x00 ..} |
+ |
+ mov esi, JSAMPARRAY [input_data(ebp)] ; input_data |
+ mov edi, JSAMPARRAY [output_data(ebp)] ; output_data |
+ alignx 16,7 |
+.rowloop: |
+ push ecx |
+ push edi |
+ push esi |
+ |
+ mov esi, JSAMPROW [esi] ; inptr |
+ mov edi, JSAMPROW [edi] ; outptr |
+ |
+ cmp ecx, byte SIZEOF_XMMWORD |
+ jae short .columnloop |
+ alignx 16,7 |
+ |
+.columnloop_r8: |
+ movdqa xmm0, XMMWORD [esi+0*SIZEOF_XMMWORD] |
+ pxor xmm1,xmm1 |
+ mov ecx, SIZEOF_XMMWORD |
+ jmp short .downsample |
+ alignx 16,7 |
+ |
+.columnloop: |
+ movdqa xmm0, XMMWORD [esi+0*SIZEOF_XMMWORD] |
+ movdqa xmm1, XMMWORD [esi+1*SIZEOF_XMMWORD] |
+ |
+.downsample: |
+ movdqa xmm2,xmm0 |
+ movdqa xmm3,xmm1 |
+ |
+ pand xmm0,xmm6 |
+ psrlw xmm2,BYTE_BIT |
+ pand xmm1,xmm6 |
+ psrlw xmm3,BYTE_BIT |
+ |
+ paddw xmm0,xmm2 |
+ paddw xmm1,xmm3 |
+ paddw xmm0,xmm7 |
+ paddw xmm1,xmm7 |
+ psrlw xmm0,1 |
+ psrlw xmm1,1 |
+ |
+ packuswb xmm0,xmm1 |
+ |
+ movdqa XMMWORD [edi+0*SIZEOF_XMMWORD], xmm0 |
+ |
+ sub ecx, byte SIZEOF_XMMWORD ; outcol |
+ add esi, byte 2*SIZEOF_XMMWORD ; inptr |
+ add edi, byte 1*SIZEOF_XMMWORD ; outptr |
+ cmp ecx, byte SIZEOF_XMMWORD |
+ jae short .columnloop |
+ test ecx,ecx |
+ jnz short .columnloop_r8 |
+ |
+ pop esi |
+ pop edi |
+ pop ecx |
+ |
+ add esi, byte SIZEOF_JSAMPROW ; input_data |
+ add edi, byte SIZEOF_JSAMPROW ; output_data |
+ dec eax ; rowctr |
+ jg near .rowloop |
+ |
+.return: |
+ pop edi |
+ pop esi |
+; pop edx ; need not be preserved |
+; pop ecx ; need not be preserved |
+; pop ebx ; unused |
+ pop ebp |
+ ret |
+ |
+; -------------------------------------------------------------------------- |
+; |
+; Downsample pixel values of a single component. |
+; This version handles the standard case of 2:1 horizontal and 2:1 vertical, |
+; without smoothing. |
+; |
+; GLOBAL(void) |
+; jsimd_h2v2_downsample_sse2 (JDIMENSION image_width, int max_v_samp_factor, |
+; JDIMENSION v_samp_factor, JDIMENSION width_blocks, |
+; JSAMPARRAY input_data, JSAMPARRAY output_data); |
+; |
+ |
+%define img_width(b) (b)+8 ; JDIMENSION image_width |
+%define max_v_samp(b) (b)+12 ; int max_v_samp_factor |
+%define v_samp(b) (b)+16 ; JDIMENSION v_samp_factor |
+%define width_blks(b) (b)+20 ; JDIMENSION width_blocks |
+%define input_data(b) (b)+24 ; JSAMPARRAY input_data |
+%define output_data(b) (b)+28 ; JSAMPARRAY output_data |
+ |
+ align 16 |
+ global EXTN(jsimd_h2v2_downsample_sse2) |
+ |
+EXTN(jsimd_h2v2_downsample_sse2): |
+ push ebp |
+ mov ebp,esp |
+; push ebx ; unused |
+; push ecx ; need not be preserved |
+; push edx ; need not be preserved |
+ push esi |
+ push edi |
+ |
+ mov ecx, JDIMENSION [width_blks(ebp)] |
+ shl ecx,3 ; imul ecx,DCTSIZE (ecx = output_cols) |
+ jz near .return |
+ |
+ mov edx, JDIMENSION [img_width(ebp)] |
+ |
+ ; -- expand_right_edge |
+ |
+ push ecx |
+ shl ecx,1 ; output_cols * 2 |
+ sub ecx,edx |
+ jle short .expand_end |
+ |
+ mov eax, INT [max_v_samp(ebp)] |
+ test eax,eax |
+ jle short .expand_end |
+ |
+ cld |
+ mov esi, JSAMPARRAY [input_data(ebp)] ; input_data |
+ alignx 16,7 |
+.expandloop: |
+ push eax |
+ push ecx |
+ |
+ mov edi, JSAMPROW [esi] |
+ add edi,edx |
+ mov al, JSAMPLE [edi-1] |
+ |
+ rep stosb |
+ |
+ pop ecx |
+ pop eax |
+ |
+ add esi, byte SIZEOF_JSAMPROW |
+ dec eax |
+ jg short .expandloop |
+ |
+.expand_end: |
+ pop ecx ; output_cols |
+ |
+ ; -- h2v2_downsample |
+ |
+ mov eax, JDIMENSION [v_samp(ebp)] ; rowctr |
+ test eax,eax |
+ jle near .return |
+ |
+ mov edx, 0x00020001 ; bias pattern |
+ movd xmm7,edx |
+ pcmpeqw xmm6,xmm6 |
+ pshufd xmm7,xmm7,0x00 ; xmm7={1, 2, 1, 2, 1, 2, 1, 2} |
+ psrlw xmm6,BYTE_BIT ; xmm6={0xFF 0x00 0xFF 0x00 ..} |
+ |
+ mov esi, JSAMPARRAY [input_data(ebp)] ; input_data |
+ mov edi, JSAMPARRAY [output_data(ebp)] ; output_data |
+ alignx 16,7 |
+.rowloop: |
+ push ecx |
+ push edi |
+ push esi |
+ |
+ mov edx, JSAMPROW [esi+0*SIZEOF_JSAMPROW] ; inptr0 |
+ mov esi, JSAMPROW [esi+1*SIZEOF_JSAMPROW] ; inptr1 |
+ mov edi, JSAMPROW [edi] ; outptr |
+ |
+ cmp ecx, byte SIZEOF_XMMWORD |
+ jae short .columnloop |
+ alignx 16,7 |
+ |
+.columnloop_r8: |
+ movdqa xmm0, XMMWORD [edx+0*SIZEOF_XMMWORD] |
+ movdqa xmm1, XMMWORD [esi+0*SIZEOF_XMMWORD] |
+ pxor xmm2,xmm2 |
+ pxor xmm3,xmm3 |
+ mov ecx, SIZEOF_XMMWORD |
+ jmp short .downsample |
+ alignx 16,7 |
+ |
+.columnloop: |
+ movdqa xmm0, XMMWORD [edx+0*SIZEOF_XMMWORD] |
+ movdqa xmm1, XMMWORD [esi+0*SIZEOF_XMMWORD] |
+ movdqa xmm2, XMMWORD [edx+1*SIZEOF_XMMWORD] |
+ movdqa xmm3, XMMWORD [esi+1*SIZEOF_XMMWORD] |
+ |
+.downsample: |
+ movdqa xmm4,xmm0 |
+ movdqa xmm5,xmm1 |
+ pand xmm0,xmm6 |
+ psrlw xmm4,BYTE_BIT |
+ pand xmm1,xmm6 |
+ psrlw xmm5,BYTE_BIT |
+ paddw xmm0,xmm4 |
+ paddw xmm1,xmm5 |
+ |
+ movdqa xmm4,xmm2 |
+ movdqa xmm5,xmm3 |
+ pand xmm2,xmm6 |
+ psrlw xmm4,BYTE_BIT |
+ pand xmm3,xmm6 |
+ psrlw xmm5,BYTE_BIT |
+ paddw xmm2,xmm4 |
+ paddw xmm3,xmm5 |
+ |
+ paddw xmm0,xmm1 |
+ paddw xmm2,xmm3 |
+ paddw xmm0,xmm7 |
+ paddw xmm2,xmm7 |
+ psrlw xmm0,2 |
+ psrlw xmm2,2 |
+ |
+ packuswb xmm0,xmm2 |
+ |
+ movdqa XMMWORD [edi+0*SIZEOF_XMMWORD], xmm0 |
+ |
+ sub ecx, byte SIZEOF_XMMWORD ; outcol |
+ add edx, byte 2*SIZEOF_XMMWORD ; inptr0 |
+ add esi, byte 2*SIZEOF_XMMWORD ; inptr1 |
+ add edi, byte 1*SIZEOF_XMMWORD ; outptr |
+ cmp ecx, byte SIZEOF_XMMWORD |
+ jae near .columnloop |
+ test ecx,ecx |
+ jnz near .columnloop_r8 |
+ |
+ pop esi |
+ pop edi |
+ pop ecx |
+ |
+ add esi, byte 2*SIZEOF_JSAMPROW ; input_data |
+ add edi, byte 1*SIZEOF_JSAMPROW ; output_data |
+ dec eax ; rowctr |
+ jg near .rowloop |
+ |
+.return: |
+ pop edi |
+ pop esi |
+; pop edx ; need not be preserved |
+; pop ecx ; need not be preserved |
+; pop ebx ; unused |
+ pop ebp |
+ ret |
+ |
+; For some reason, the OS X linker does not honor the request to align the |
+; segment unless we do this. |
+ align 16 |