media/libjpeg/simd/jdclrss2.asm

Wed, 31 Dec 2014 06:09:35 +0100

author
Michael Schloh von Bennewitz <michael@schloh.com>
date
Wed, 31 Dec 2014 06:09:35 +0100
changeset 0
6474c204b198
permissions
-rw-r--r--

Cloned upstream origin tor-browser at tor-browser-31.3.0esr-4.5-1-build1
revision ID fc1c9ff7c1b2defdbc039f12214767608f46423f for hacking purpose.

     1 ;
     2 ; jdclrss2.asm - colorspace conversion (SSE2)
     3 ;
     4 ; Copyright 2009, 2012 Pierre Ossman <ossman@cendio.se> for Cendio AB
     5 ; Copyright 2012 D. R. Commander
     6 ;
     7 ; Based on
     8 ; x86 SIMD extension for IJG JPEG library
     9 ; Copyright (C) 1999-2006, MIYASAKA Masaru.
    10 ; For conditions of distribution and use, see copyright notice in jsimdext.inc
    11 ;
    12 ; This file should be assembled with NASM (Netwide Assembler),
    13 ; can *not* be assembled with Microsoft's MASM or any compatible
    14 ; assembler (including Borland's Turbo Assembler).
    15 ; NASM is available from http://nasm.sourceforge.net/ or
    16 ; http://sourceforge.net/project/showfiles.php?group_id=6208
    17 ;
    18 ; [TAB8]
    20 %include "jcolsamp.inc"
    22 ; --------------------------------------------------------------------------
    23 ;
    24 ; Convert some rows of samples to the output colorspace.
    25 ;
    26 ; GLOBAL(void)
    27 ; jsimd_ycc_rgb_convert_sse2 (JDIMENSION out_width,
    28 ;                             JSAMPIMAGE input_buf, JDIMENSION input_row,
    29 ;                             JSAMPARRAY output_buf, int num_rows)
    30 ;
    32 %define out_width(b)	(b)+8			; JDIMENSION out_width
    33 %define input_buf(b)	(b)+12		; JSAMPIMAGE input_buf
    34 %define input_row(b)	(b)+16		; JDIMENSION input_row
    35 %define output_buf(b)	(b)+20		; JSAMPARRAY output_buf
    36 %define num_rows(b)	(b)+24		; int num_rows
    38 %define original_ebp	ebp+0
    39 %define wk(i)		ebp-(WK_NUM-(i))*SIZEOF_XMMWORD	; xmmword wk[WK_NUM]
    40 %define WK_NUM		2
    41 %define gotptr		wk(0)-SIZEOF_POINTER	; void * gotptr
    43 	align	16
    44 	global	EXTN(jsimd_ycc_rgb_convert_sse2)
    46 EXTN(jsimd_ycc_rgb_convert_sse2):
    47 	push	ebp
    48 	mov	eax,esp				; eax = original ebp
    49 	sub	esp, byte 4
    50 	and	esp, byte (-SIZEOF_XMMWORD)	; align to 128 bits
    51 	mov	[esp],eax
    52 	mov	ebp,esp				; ebp = aligned ebp
    53 	lea	esp, [wk(0)]
    54 	pushpic	eax		; make a room for GOT address
    55 	push	ebx
    56 ;	push	ecx		; need not be preserved
    57 ;	push	edx		; need not be preserved
    58 	push	esi
    59 	push	edi
    61 	get_GOT	ebx			; get GOT address
    62 	movpic	POINTER [gotptr], ebx	; save GOT address
    64 	mov	ecx, JDIMENSION [out_width(eax)]	; num_cols
    65 	test	ecx,ecx
    66 	jz	near .return
    68 	push	ecx
    70 	mov	edi, JSAMPIMAGE [input_buf(eax)]
    71 	mov	ecx, JDIMENSION [input_row(eax)]
    72 	mov	esi, JSAMPARRAY [edi+0*SIZEOF_JSAMPARRAY]
    73 	mov	ebx, JSAMPARRAY [edi+1*SIZEOF_JSAMPARRAY]
    74 	mov	edx, JSAMPARRAY [edi+2*SIZEOF_JSAMPARRAY]
    75 	lea	esi, [esi+ecx*SIZEOF_JSAMPROW]
    76 	lea	ebx, [ebx+ecx*SIZEOF_JSAMPROW]
    77 	lea	edx, [edx+ecx*SIZEOF_JSAMPROW]
    79 	pop	ecx
    81 	mov	edi, JSAMPARRAY [output_buf(eax)]
    82 	mov	eax, INT [num_rows(eax)]
    83 	test	eax,eax
    84 	jle	near .return
    85 	alignx	16,7
    86 .rowloop:
    87 	push	eax
    88 	push	edi
    89 	push	edx
    90 	push	ebx
    91 	push	esi
    92 	push	ecx			; col
    94 	mov	esi, JSAMPROW [esi]	; inptr0
    95 	mov	ebx, JSAMPROW [ebx]	; inptr1
    96 	mov	edx, JSAMPROW [edx]	; inptr2
    97 	mov	edi, JSAMPROW [edi]	; outptr
    98 	movpic	eax, POINTER [gotptr]	; load GOT address (eax)
    99 	alignx	16,7
   100 .columnloop:
   102 	movdqa	xmm5, XMMWORD [ebx]	; xmm5=Cb(0123456789ABCDEF)
   103 	movdqa	xmm1, XMMWORD [edx]	; xmm1=Cr(0123456789ABCDEF)
   105 	pcmpeqw	xmm4,xmm4
   106 	pcmpeqw	xmm7,xmm7
   107 	psrlw	xmm4,BYTE_BIT
   108 	psllw	xmm7,7			; xmm7={0xFF80 0xFF80 0xFF80 0xFF80 ..}
   109 	movdqa	xmm0,xmm4		; xmm0=xmm4={0xFF 0x00 0xFF 0x00 ..}
   111 	pand	xmm4,xmm5		; xmm4=Cb(02468ACE)=CbE
   112 	psrlw	xmm5,BYTE_BIT		; xmm5=Cb(13579BDF)=CbO
   113 	pand	xmm0,xmm1		; xmm0=Cr(02468ACE)=CrE
   114 	psrlw	xmm1,BYTE_BIT		; xmm1=Cr(13579BDF)=CrO
   116 	paddw	xmm4,xmm7
   117 	paddw	xmm5,xmm7
   118 	paddw	xmm0,xmm7
   119 	paddw	xmm1,xmm7
   121 	; (Original)
   122 	; R = Y                + 1.40200 * Cr
   123 	; G = Y - 0.34414 * Cb - 0.71414 * Cr
   124 	; B = Y + 1.77200 * Cb
   125 	;
   126 	; (This implementation)
   127 	; R = Y                + 0.40200 * Cr + Cr
   128 	; G = Y - 0.34414 * Cb + 0.28586 * Cr - Cr
   129 	; B = Y - 0.22800 * Cb + Cb + Cb
   131 	movdqa	xmm2,xmm4		; xmm2=CbE
   132 	movdqa	xmm3,xmm5		; xmm3=CbO
   133 	paddw	xmm4,xmm4		; xmm4=2*CbE
   134 	paddw	xmm5,xmm5		; xmm5=2*CbO
   135 	movdqa	xmm6,xmm0		; xmm6=CrE
   136 	movdqa	xmm7,xmm1		; xmm7=CrO
   137 	paddw	xmm0,xmm0		; xmm0=2*CrE
   138 	paddw	xmm1,xmm1		; xmm1=2*CrO
   140 	pmulhw	xmm4,[GOTOFF(eax,PW_MF0228)]	; xmm4=(2*CbE * -FIX(0.22800))
   141 	pmulhw	xmm5,[GOTOFF(eax,PW_MF0228)]	; xmm5=(2*CbO * -FIX(0.22800))
   142 	pmulhw	xmm0,[GOTOFF(eax,PW_F0402)]	; xmm0=(2*CrE * FIX(0.40200))
   143 	pmulhw	xmm1,[GOTOFF(eax,PW_F0402)]	; xmm1=(2*CrO * FIX(0.40200))
   145 	paddw	xmm4,[GOTOFF(eax,PW_ONE)]
   146 	paddw	xmm5,[GOTOFF(eax,PW_ONE)]
   147 	psraw	xmm4,1			; xmm4=(CbE * -FIX(0.22800))
   148 	psraw	xmm5,1			; xmm5=(CbO * -FIX(0.22800))
   149 	paddw	xmm0,[GOTOFF(eax,PW_ONE)]
   150 	paddw	xmm1,[GOTOFF(eax,PW_ONE)]
   151 	psraw	xmm0,1			; xmm0=(CrE * FIX(0.40200))
   152 	psraw	xmm1,1			; xmm1=(CrO * FIX(0.40200))
   154 	paddw	xmm4,xmm2
   155 	paddw	xmm5,xmm3
   156 	paddw	xmm4,xmm2		; xmm4=(CbE * FIX(1.77200))=(B-Y)E
   157 	paddw	xmm5,xmm3		; xmm5=(CbO * FIX(1.77200))=(B-Y)O
   158 	paddw	xmm0,xmm6		; xmm0=(CrE * FIX(1.40200))=(R-Y)E
   159 	paddw	xmm1,xmm7		; xmm1=(CrO * FIX(1.40200))=(R-Y)O
   161 	movdqa	XMMWORD [wk(0)], xmm4	; wk(0)=(B-Y)E
   162 	movdqa	XMMWORD [wk(1)], xmm5	; wk(1)=(B-Y)O
   164 	movdqa    xmm4,xmm2
   165 	movdqa    xmm5,xmm3
   166 	punpcklwd xmm2,xmm6
   167 	punpckhwd xmm4,xmm6
   168 	pmaddwd   xmm2,[GOTOFF(eax,PW_MF0344_F0285)]
   169 	pmaddwd   xmm4,[GOTOFF(eax,PW_MF0344_F0285)]
   170 	punpcklwd xmm3,xmm7
   171 	punpckhwd xmm5,xmm7
   172 	pmaddwd   xmm3,[GOTOFF(eax,PW_MF0344_F0285)]
   173 	pmaddwd   xmm5,[GOTOFF(eax,PW_MF0344_F0285)]
   175 	paddd     xmm2,[GOTOFF(eax,PD_ONEHALF)]
   176 	paddd     xmm4,[GOTOFF(eax,PD_ONEHALF)]
   177 	psrad     xmm2,SCALEBITS
   178 	psrad     xmm4,SCALEBITS
   179 	paddd     xmm3,[GOTOFF(eax,PD_ONEHALF)]
   180 	paddd     xmm5,[GOTOFF(eax,PD_ONEHALF)]
   181 	psrad     xmm3,SCALEBITS
   182 	psrad     xmm5,SCALEBITS
   184 	packssdw  xmm2,xmm4	; xmm2=CbE*-FIX(0.344)+CrE*FIX(0.285)
   185 	packssdw  xmm3,xmm5	; xmm3=CbO*-FIX(0.344)+CrO*FIX(0.285)
   186 	psubw     xmm2,xmm6	; xmm2=CbE*-FIX(0.344)+CrE*-FIX(0.714)=(G-Y)E
   187 	psubw     xmm3,xmm7	; xmm3=CbO*-FIX(0.344)+CrO*-FIX(0.714)=(G-Y)O
   189 	movdqa    xmm5, XMMWORD [esi]	; xmm5=Y(0123456789ABCDEF)
   191 	pcmpeqw   xmm4,xmm4
   192 	psrlw     xmm4,BYTE_BIT		; xmm4={0xFF 0x00 0xFF 0x00 ..}
   193 	pand      xmm4,xmm5		; xmm4=Y(02468ACE)=YE
   194 	psrlw     xmm5,BYTE_BIT		; xmm5=Y(13579BDF)=YO
   196 	paddw     xmm0,xmm4		; xmm0=((R-Y)E+YE)=RE=R(02468ACE)
   197 	paddw     xmm1,xmm5		; xmm1=((R-Y)O+YO)=RO=R(13579BDF)
   198 	packuswb  xmm0,xmm0		; xmm0=R(02468ACE********)
   199 	packuswb  xmm1,xmm1		; xmm1=R(13579BDF********)
   201 	paddw     xmm2,xmm4		; xmm2=((G-Y)E+YE)=GE=G(02468ACE)
   202 	paddw     xmm3,xmm5		; xmm3=((G-Y)O+YO)=GO=G(13579BDF)
   203 	packuswb  xmm2,xmm2		; xmm2=G(02468ACE********)
   204 	packuswb  xmm3,xmm3		; xmm3=G(13579BDF********)
   206 	paddw     xmm4, XMMWORD [wk(0)]	; xmm4=(YE+(B-Y)E)=BE=B(02468ACE)
   207 	paddw     xmm5, XMMWORD [wk(1)]	; xmm5=(YO+(B-Y)O)=BO=B(13579BDF)
   208 	packuswb  xmm4,xmm4		; xmm4=B(02468ACE********)
   209 	packuswb  xmm5,xmm5		; xmm5=B(13579BDF********)
   211 %if RGB_PIXELSIZE == 3 ; ---------------
   213 	; xmmA=(00 02 04 06 08 0A 0C 0E **), xmmB=(01 03 05 07 09 0B 0D 0F **)
   214 	; xmmC=(10 12 14 16 18 1A 1C 1E **), xmmD=(11 13 15 17 19 1B 1D 1F **)
   215 	; xmmE=(20 22 24 26 28 2A 2C 2E **), xmmF=(21 23 25 27 29 2B 2D 2F **)
   216 	; xmmG=(** ** ** ** ** ** ** ** **), xmmH=(** ** ** ** ** ** ** ** **)
   218 	punpcklbw xmmA,xmmC	; xmmA=(00 10 02 12 04 14 06 16 08 18 0A 1A 0C 1C 0E 1E)
   219 	punpcklbw xmmE,xmmB	; xmmE=(20 01 22 03 24 05 26 07 28 09 2A 0B 2C 0D 2E 0F)
   220 	punpcklbw xmmD,xmmF	; xmmD=(11 21 13 23 15 25 17 27 19 29 1B 2B 1D 2D 1F 2F)
   222 	movdqa    xmmG,xmmA
   223 	movdqa    xmmH,xmmA
   224 	punpcklwd xmmA,xmmE	; xmmA=(00 10 20 01 02 12 22 03 04 14 24 05 06 16 26 07)
   225 	punpckhwd xmmG,xmmE	; xmmG=(08 18 28 09 0A 1A 2A 0B 0C 1C 2C 0D 0E 1E 2E 0F)
   227 	psrldq    xmmH,2	; xmmH=(02 12 04 14 06 16 08 18 0A 1A 0C 1C 0E 1E -- --)
   228 	psrldq    xmmE,2	; xmmE=(22 03 24 05 26 07 28 09 2A 0B 2C 0D 2E 0F -- --)
   230 	movdqa    xmmC,xmmD
   231 	movdqa    xmmB,xmmD
   232 	punpcklwd xmmD,xmmH	; xmmD=(11 21 02 12 13 23 04 14 15 25 06 16 17 27 08 18)
   233 	punpckhwd xmmC,xmmH	; xmmC=(19 29 0A 1A 1B 2B 0C 1C 1D 2D 0E 1E 1F 2F -- --)
   235 	psrldq    xmmB,2	; xmmB=(13 23 15 25 17 27 19 29 1B 2B 1D 2D 1F 2F -- --)
   237 	movdqa    xmmF,xmmE
   238 	punpcklwd xmmE,xmmB	; xmmE=(22 03 13 23 24 05 15 25 26 07 17 27 28 09 19 29)
   239 	punpckhwd xmmF,xmmB	; xmmF=(2A 0B 1B 2B 2C 0D 1D 2D 2E 0F 1F 2F -- -- -- --)
   241 	pshufd    xmmH,xmmA,0x4E; xmmH=(04 14 24 05 06 16 26 07 00 10 20 01 02 12 22 03)
   242 	movdqa    xmmB,xmmE
   243 	punpckldq xmmA,xmmD	; xmmA=(00 10 20 01 11 21 02 12 02 12 22 03 13 23 04 14)
   244 	punpckldq xmmE,xmmH	; xmmE=(22 03 13 23 04 14 24 05 24 05 15 25 06 16 26 07)
   245 	punpckhdq xmmD,xmmB	; xmmD=(15 25 06 16 26 07 17 27 17 27 08 18 28 09 19 29)
   247 	pshufd    xmmH,xmmG,0x4E; xmmH=(0C 1C 2C 0D 0E 1E 2E 0F 08 18 28 09 0A 1A 2A 0B)
   248 	movdqa    xmmB,xmmF
   249 	punpckldq xmmG,xmmC	; xmmG=(08 18 28 09 19 29 0A 1A 0A 1A 2A 0B 1B 2B 0C 1C)
   250 	punpckldq xmmF,xmmH	; xmmF=(2A 0B 1B 2B 0C 1C 2C 0D 2C 0D 1D 2D 0E 1E 2E 0F)
   251 	punpckhdq xmmC,xmmB	; xmmC=(1D 2D 0E 1E 2E 0F 1F 2F 1F 2F -- -- -- -- -- --)
   253 	punpcklqdq xmmA,xmmE	; xmmA=(00 10 20 01 11 21 02 12 22 03 13 23 04 14 24 05)
   254 	punpcklqdq xmmD,xmmG	; xmmD=(15 25 06 16 26 07 17 27 08 18 28 09 19 29 0A 1A)
   255 	punpcklqdq xmmF,xmmC	; xmmF=(2A 0B 1B 2B 0C 1C 2C 0D 1D 2D 0E 1E 2E 0F 1F 2F)
   257 	cmp	ecx, byte SIZEOF_XMMWORD
   258 	jb	short .column_st32
   260 	test	edi, SIZEOF_XMMWORD-1
   261 	jnz	short .out1
   262 	; --(aligned)-------------------
   263 	movntdq	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   264 	movntdq	XMMWORD [edi+1*SIZEOF_XMMWORD], xmmD
   265 	movntdq	XMMWORD [edi+2*SIZEOF_XMMWORD], xmmF
   266 	jmp	short .out0
   267 .out1:	; --(unaligned)-----------------
   268 	movdqu	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   269 	movdqu	XMMWORD [edi+1*SIZEOF_XMMWORD], xmmD
   270 	movdqu	XMMWORD [edi+2*SIZEOF_XMMWORD], xmmF
   271 .out0:
   272 	add	edi, byte RGB_PIXELSIZE*SIZEOF_XMMWORD	; outptr
   273 	sub	ecx, byte SIZEOF_XMMWORD
   274 	jz	near .nextrow
   276 	add	esi, byte SIZEOF_XMMWORD	; inptr0
   277 	add	ebx, byte SIZEOF_XMMWORD	; inptr1
   278 	add	edx, byte SIZEOF_XMMWORD	; inptr2
   279 	jmp	near .columnloop
   280 	alignx	16,7
   282 .column_st32:
   283 	lea	ecx, [ecx+ecx*2]		; imul ecx, RGB_PIXELSIZE
   284 	cmp	ecx, byte 2*SIZEOF_XMMWORD
   285 	jb	short .column_st16
   286 	movdqu	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   287 	movdqu	XMMWORD [edi+1*SIZEOF_XMMWORD], xmmD
   288 	add	edi, byte 2*SIZEOF_XMMWORD	; outptr
   289 	movdqa	xmmA,xmmF
   290 	sub	ecx, byte 2*SIZEOF_XMMWORD
   291 	jmp	short .column_st15
   292 .column_st16:
   293 	cmp	ecx, byte SIZEOF_XMMWORD
   294 	jb	short .column_st15
   295 	movdqu	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   296 	add	edi, byte SIZEOF_XMMWORD	; outptr
   297 	movdqa	xmmA,xmmD
   298 	sub	ecx, byte SIZEOF_XMMWORD
   299 .column_st15:
   300 	; Store the lower 8 bytes of xmmA to the output when it has enough
   301 	; space.
   302 	cmp	ecx, byte SIZEOF_MMWORD
   303 	jb	short .column_st7
   304 	movq	XMM_MMWORD [edi], xmmA
   305 	add	edi, byte SIZEOF_MMWORD
   306 	sub	ecx, byte SIZEOF_MMWORD
   307 	psrldq	xmmA, SIZEOF_MMWORD
   308 .column_st7:
   309 	; Store the lower 4 bytes of xmmA to the output when it has enough
   310 	; space.
   311 	cmp	ecx, byte SIZEOF_DWORD
   312 	jb	short .column_st3
   313 	movd	XMM_DWORD [edi], xmmA
   314 	add	edi, byte SIZEOF_DWORD
   315 	sub	ecx, byte SIZEOF_DWORD
   316 	psrldq	xmmA, SIZEOF_DWORD
   317 .column_st3:
   318 	; Store the lower 2 bytes of eax to the output when it has enough
   319 	; space.
   320 	movd	eax, xmmA
   321 	cmp	ecx, byte SIZEOF_WORD
   322 	jb	short .column_st1
   323 	mov	WORD [edi], ax
   324 	add	edi, byte SIZEOF_WORD
   325 	sub	ecx, byte SIZEOF_WORD
   326 	shr	eax, 16
   327 .column_st1:
   328 	; Store the lower 1 byte of eax to the output when it has enough
   329 	; space.
   330 	test	ecx, ecx
   331 	jz	short .nextrow
   332 	mov	BYTE [edi], al
   334 %else ; RGB_PIXELSIZE == 4 ; -----------
   336 %ifdef RGBX_FILLER_0XFF
   337 	pcmpeqb   xmm6,xmm6		; xmm6=XE=X(02468ACE********)
   338 	pcmpeqb   xmm7,xmm7		; xmm7=XO=X(13579BDF********)
   339 %else
   340 	pxor      xmm6,xmm6		; xmm6=XE=X(02468ACE********)
   341 	pxor      xmm7,xmm7		; xmm7=XO=X(13579BDF********)
   342 %endif
   343 	; xmmA=(00 02 04 06 08 0A 0C 0E **), xmmB=(01 03 05 07 09 0B 0D 0F **)
   344 	; xmmC=(10 12 14 16 18 1A 1C 1E **), xmmD=(11 13 15 17 19 1B 1D 1F **)
   345 	; xmmE=(20 22 24 26 28 2A 2C 2E **), xmmF=(21 23 25 27 29 2B 2D 2F **)
   346 	; xmmG=(30 32 34 36 38 3A 3C 3E **), xmmH=(31 33 35 37 39 3B 3D 3F **)
   348 	punpcklbw xmmA,xmmC	; xmmA=(00 10 02 12 04 14 06 16 08 18 0A 1A 0C 1C 0E 1E)
   349 	punpcklbw xmmE,xmmG	; xmmE=(20 30 22 32 24 34 26 36 28 38 2A 3A 2C 3C 2E 3E)
   350 	punpcklbw xmmB,xmmD	; xmmB=(01 11 03 13 05 15 07 17 09 19 0B 1B 0D 1D 0F 1F)
   351 	punpcklbw xmmF,xmmH	; xmmF=(21 31 23 33 25 35 27 37 29 39 2B 3B 2D 3D 2F 3F)
   353 	movdqa    xmmC,xmmA
   354 	punpcklwd xmmA,xmmE	; xmmA=(00 10 20 30 02 12 22 32 04 14 24 34 06 16 26 36)
   355 	punpckhwd xmmC,xmmE	; xmmC=(08 18 28 38 0A 1A 2A 3A 0C 1C 2C 3C 0E 1E 2E 3E)
   356 	movdqa    xmmG,xmmB
   357 	punpcklwd xmmB,xmmF	; xmmB=(01 11 21 31 03 13 23 33 05 15 25 35 07 17 27 37)
   358 	punpckhwd xmmG,xmmF	; xmmG=(09 19 29 39 0B 1B 2B 3B 0D 1D 2D 3D 0F 1F 2F 3F)
   360 	movdqa    xmmD,xmmA
   361 	punpckldq xmmA,xmmB	; xmmA=(00 10 20 30 01 11 21 31 02 12 22 32 03 13 23 33)
   362 	punpckhdq xmmD,xmmB	; xmmD=(04 14 24 34 05 15 25 35 06 16 26 36 07 17 27 37)
   363 	movdqa    xmmH,xmmC
   364 	punpckldq xmmC,xmmG	; xmmC=(08 18 28 38 09 19 29 39 0A 1A 2A 3A 0B 1B 2B 3B)
   365 	punpckhdq xmmH,xmmG	; xmmH=(0C 1C 2C 3C 0D 1D 2D 3D 0E 1E 2E 3E 0F 1F 2F 3F)
   367 	cmp	ecx, byte SIZEOF_XMMWORD
   368 	jb	short .column_st32
   370 	test	edi, SIZEOF_XMMWORD-1
   371 	jnz	short .out1
   372 	; --(aligned)-------------------
   373 	movntdq	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   374 	movntdq	XMMWORD [edi+1*SIZEOF_XMMWORD], xmmD
   375 	movntdq	XMMWORD [edi+2*SIZEOF_XMMWORD], xmmC
   376 	movntdq	XMMWORD [edi+3*SIZEOF_XMMWORD], xmmH
   377 	jmp	short .out0
   378 .out1:	; --(unaligned)-----------------
   379 	movdqu	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   380 	movdqu	XMMWORD [edi+1*SIZEOF_XMMWORD], xmmD
   381 	movdqu	XMMWORD [edi+2*SIZEOF_XMMWORD], xmmC
   382 	movdqu	XMMWORD [edi+3*SIZEOF_XMMWORD], xmmH
   383 .out0:
   384 	add	edi, byte RGB_PIXELSIZE*SIZEOF_XMMWORD	; outptr
   385 	sub	ecx, byte SIZEOF_XMMWORD
   386 	jz	near .nextrow
   388 	add	esi, byte SIZEOF_XMMWORD	; inptr0
   389 	add	ebx, byte SIZEOF_XMMWORD	; inptr1
   390 	add	edx, byte SIZEOF_XMMWORD	; inptr2
   391 	jmp	near .columnloop
   392 	alignx	16,7
   394 .column_st32:
   395 	cmp	ecx, byte SIZEOF_XMMWORD/2
   396 	jb	short .column_st16
   397 	movdqu	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   398 	movdqu	XMMWORD [edi+1*SIZEOF_XMMWORD], xmmD
   399 	add	edi, byte 2*SIZEOF_XMMWORD	; outptr
   400 	movdqa	xmmA,xmmC
   401 	movdqa	xmmD,xmmH
   402 	sub	ecx, byte SIZEOF_XMMWORD/2
   403 .column_st16:
   404 	cmp	ecx, byte SIZEOF_XMMWORD/4
   405 	jb	short .column_st15
   406 	movdqu	XMMWORD [edi+0*SIZEOF_XMMWORD], xmmA
   407 	add	edi, byte SIZEOF_XMMWORD	; outptr
   408 	movdqa	xmmA,xmmD
   409 	sub	ecx, byte SIZEOF_XMMWORD/4
   410 .column_st15:
   411 	; Store two pixels (8 bytes) of xmmA to the output when it has enough
   412 	; space.
   413 	cmp	ecx, byte SIZEOF_XMMWORD/8
   414 	jb	short .column_st7
   415 	movq	XMM_MMWORD [edi], xmmA
   416 	add	edi, byte SIZEOF_XMMWORD/8*4
   417 	sub	ecx, byte SIZEOF_XMMWORD/8
   418 	psrldq	xmmA, SIZEOF_XMMWORD/8*4
   419 .column_st7:
   420 	; Store one pixel (4 bytes) of xmmA to the output when it has enough
   421 	; space.
   422 	test	ecx, ecx
   423 	jz	short .nextrow
   424 	movd	XMM_DWORD [edi], xmmA
   426 %endif ; RGB_PIXELSIZE ; ---------------
   428 	alignx	16,7
   430 .nextrow:
   431 	pop	ecx
   432 	pop	esi
   433 	pop	ebx
   434 	pop	edx
   435 	pop	edi
   436 	pop	eax
   438 	add	esi, byte SIZEOF_JSAMPROW
   439 	add	ebx, byte SIZEOF_JSAMPROW
   440 	add	edx, byte SIZEOF_JSAMPROW
   441 	add	edi, byte SIZEOF_JSAMPROW	; output_buf
   442 	dec	eax				; num_rows
   443 	jg	near .rowloop
   445 	sfence		; flush the write buffer
   447 .return:
   448 	pop	edi
   449 	pop	esi
   450 ;	pop	edx		; need not be preserved
   451 ;	pop	ecx		; need not be preserved
   452 	pop	ebx
   453 	mov	esp,ebp		; esp <- aligned ebp
   454 	pop	esp		; esp <- original ebp
   455 	pop	ebp
   456 	ret
   458 ; For some reason, the OS X linker does not honor the request to align the
   459 ; segment unless we do this.
   460 	align	16

mercurial