shithub: openh264

ref: e65525d7b45f0b7ad07dd45b8259baedf3314160
dir: /processing/src/asm/denoisefilter.asm/

View raw version
;*!
;* \copy
;*     Copyright (c)  2010-2013, Cisco Systems
;*     All rights reserved.
;*
;*     Redistribution and use in source and binary forms, with or without
;*     modification, are permitted provided that the following conditions
;*     are met:
;*
;*        * Redistributions of source code must retain the above copyright
;*          notice, this list of conditions and the following disclaimer.
;*
;*        * Redistributions in binary form must reproduce the above copyright
;*          notice, this list of conditions and the following disclaimer in
;*          the documentation and/or other materials provided with the
;*          distribution.
;*
;*     THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
;*     "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
;*     LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
;*     FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
;*     COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
;*     INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
;*     BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
;*     LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
;*     CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
;*     LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
;*     ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
;*     POSSIBILITY OF SUCH DAMAGE.
;*
;*
;*  predenoise.asm
;*
;*  Abstract
;*      denoise for SVC2.1
;*  History
;*      4/13/2010 Created
;*      7/30/2010 Modified
;*
;*
;*************************************************************************/
%include "asm_inc.asm"

;***********************************************************************
; Constant
;***********************************************************************
SECTION .rodata align=16

sse2_32 times 8 dw 32
sse2_20 times 8 dw 20


BITS 32
;***********************************************************************
; Code
;***********************************************************************
SECTION .text

%macro	WEIGHT_LINE	9
		movq		%2,	%9
		punpcklbw	%2,	%7
		movdqa		%8,	%2

		movdqa		%1,	%6
		psubusb		%1,	%8
		psubusb		%8,	%6
		por			%8,	%1		; ABS(curPixel - centerPixel);

		movdqa		%1,	%3
		psubusb		%1,	%8

		pmullw		%1,	%1
		psrlw		%1,	5
		pmullw		%2,	%1
		paddusw		%4,	%1
		paddusw		%5,	%2
%endmacro

%macro	WEIGHT_LINE1_UV	4
		movdqa		%2,	%1
		punpcklbw	%2,	%4
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	1
		punpcklbw	%2,	%4
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	2
		punpcklbw	%2,	%4
		psllw		%2,	1
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	3
		punpcklbw	%2,	%4
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	4
		punpcklbw	%2,	%4
		paddw		%3,	%2
%endmacro

%macro	WEIGHT_LINE2_UV	4
		movdqa		%2,	%1
		punpcklbw	%2,	%4
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	1
		punpcklbw	%2,	%4
		psllw		%2,	1
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	2
		punpcklbw	%2,	%4
		psllw		%2,	2
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	3
		punpcklbw	%2,	%4
		psllw		%2,	1
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	4
		punpcklbw	%2,	%4
		paddw		%3,	%2
%endmacro

%macro	WEIGHT_LINE3_UV	4
		movdqa		%2,	%1
		punpcklbw	%2,	%4
		psllw		%2,	1
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	1
		punpcklbw	%2,	%4
		psllw		%2,	2
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	2
		punpcklbw	%2,	%4
		pmullw		%2,	[sse2_20]
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	3
		punpcklbw	%2,	%4
		psllw		%2,	2
		paddw		%3,	%2

		movdqa		%2,	%1
		psrldq		%2,	4
		punpcklbw	%2,	%4
		psllw		%2,	1
		paddw		%3,	%2
%endmacro

ALIGN 16
WELS_EXTERN BilateralLumaFilter8_sse2
;***********************************************************************
;  BilateralLumaFilter8_sse2(uint8_t *pixels, int stride);
;***********************************************************************
;	1	2	3
;	4	0	5
;	6	7	8
;	0:	the center point
%define		pushsize	4
%define		pixel		esp + pushsize + 4
%define		stride		esp + pushsize + 8
BilateralLumaFilter8_sse2:
		push		ebx

		pxor		xmm7,	xmm7
		mov			eax,	[pixel]
		mov			ebx,	eax
		movq		xmm6,	[eax]
		punpcklbw	xmm6,	xmm7
		movdqa		xmm3,	[sse2_32]
		pxor		xmm4,	xmm4		; nTotWeight
		pxor		xmm5,	xmm5		; nSum

		dec			eax
		mov			ecx,	[stride]

		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax]			; pixel 4
		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax + 2]		; pixel 5

		sub			eax,	ecx
		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax]			; pixel 1
		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax + 1]		; pixel 2
		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax + 2]		; pixel 3

		lea			eax,	[eax + ecx * 2]
		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax]			; pixel 6
		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax + 1]		; pixel 7
		WEIGHT_LINE	xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7, xmm0,  [eax + 2]		; pixel 8

		pcmpeqw		xmm0,	xmm0
		psrlw		xmm0,	15
		psllw		xmm0,	8
		psubusw		xmm0,	xmm4
		pmullw		xmm0,	xmm6
		paddusw		xmm5,	xmm0
		psrlw		xmm5,	8
		packuswb	xmm5,	xmm5
		movq		[ebx],	xmm5

		pop ebx
		ret

WELS_EXTERN WaverageChromaFilter8_sse2
;***********************************************************************
; void		WaverageChromaFilter8_sse2(uint8_t *pixels, int stride);
;***********************************************************************
;5x5 filter:
;1	1	2	1	1
;1	2	4	2	1
;2	4	20	4	2
;1	2	4	2	1
;1	1	2	1	1

ALIGN 16
WaverageChromaFilter8_sse2:
		mov		edx,	[esp + 4]	; pixels
		mov		ecx,	[esp + 8]	; stride

		mov		eax,	ecx
		add		eax,	eax
		sub		edx,	eax			; pixels - 2 * stride
		sub		edx,	2

		pxor	xmm0,	xmm0
		pxor	xmm3,	xmm3

		movdqu		xmm1,	[edx]
		WEIGHT_LINE1_UV	xmm1,	xmm2,	xmm3,	xmm0

		movdqu		xmm1,	[edx + ecx]
		WEIGHT_LINE2_UV	xmm1,	xmm2,	xmm3,	xmm0

		add		edx,	eax
		movdqu		xmm1,	[edx]
		WEIGHT_LINE3_UV	xmm1,	xmm2,	xmm3,	xmm0

		movdqu		xmm1,	[edx + ecx]
		WEIGHT_LINE2_UV	xmm1,	xmm2,	xmm3,	xmm0

		movdqu		xmm1,	[edx + ecx * 2]
		WEIGHT_LINE1_UV	xmm1,	xmm2,	xmm3,	xmm0

		psrlw		xmm3,		6
		packuswb	xmm3,		xmm3
		movq		[edx + 2],		xmm3

		ret