| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296 |
- ;
- ; Copyright (c) 2010 The WebM project authors. All Rights Reserved.
- ;
- ; Use of this source code is governed by a BSD-style license
- ; that can be found in the LICENSE file in the root of the source
- ; tree. An additional intellectual property rights grant can be found
- ; in the file PATENTS. All contributing project authors may
- ; be found in the AUTHORS file in the root of the source tree.
- ;
- %include "vpx_ports/x86_abi_support.asm"
- ; /****************************************************************************
- ; * Notes:
- ; *
- ; * This implementation makes use of 16 bit fixed point version of two multiply
- ; * constants:
- ; * 1. sqrt(2) * cos (pi/8)
- ; * 2. sqrt(2) * sin (pi/8)
- ; * Because the first constant is bigger than 1, to maintain the same 16 bit
- ; * fixed point precision as the second one, we use a trick of
- ; * x * a = x + x*(a-1)
- ; * so
- ; * x * sqrt(2) * cos (pi/8) = x + x * (sqrt(2) *cos(pi/8)-1).
- ; *
- ; * For the second constant, because of the 16bit version is 35468, which
- ; * is bigger than 32768, in signed 16 bit multiply, it becomes a negative
- ; * number.
- ; * (x * (unsigned)35468 >> 16) = x * (signed)35468 >> 16 + x
- ; *
- ; **************************************************************************/
- SECTION .text
- ;void vp8_short_idct4x4llm_mmx(short *input, unsigned char *pred,
- ;int pitch, unsigned char *dest,int stride)
- global sym(vp8_short_idct4x4llm_mmx) PRIVATE
- sym(vp8_short_idct4x4llm_mmx):
- push rbp
- mov rbp, rsp
- SHADOW_ARGS_TO_STACK 5
- GET_GOT rbx
- push rsi
- push rdi
- ; end prolog
- mov rax, arg(0) ;input
- mov rsi, arg(1) ;pred
- movq mm0, [rax ]
- movq mm1, [rax+ 8]
- movq mm2, [rax+16]
- movq mm3, [rax+24]
- %if 0
- pxor mm7, mm7
- movq [rax], mm7
- movq [rax+8], mm7
- movq [rax+16],mm7
- movq [rax+24],mm7
- %endif
- movsxd rax, dword ptr arg(2) ;pitch
- mov rdx, arg(3) ;dest
- movsxd rdi, dword ptr arg(4) ;stride
- psubw mm0, mm2 ; b1= 0-2
- paddw mm2, mm2 ;
- movq mm5, mm1
- paddw mm2, mm0 ; a1 =0+2
- pmulhw mm5, [GLOBAL(x_s1sqr2)];
- paddw mm5, mm1 ; ip1 * sin(pi/8) * sqrt(2)
- movq mm7, mm3 ;
- pmulhw mm7, [GLOBAL(x_c1sqr2less1)];
- paddw mm7, mm3 ; ip3 * cos(pi/8) * sqrt(2)
- psubw mm7, mm5 ; c1
- movq mm5, mm1
- movq mm4, mm3
- pmulhw mm5, [GLOBAL(x_c1sqr2less1)]
- paddw mm5, mm1
- pmulhw mm3, [GLOBAL(x_s1sqr2)]
- paddw mm3, mm4
- paddw mm3, mm5 ; d1
- movq mm6, mm2 ; a1
- movq mm4, mm0 ; b1
- paddw mm2, mm3 ;0
- paddw mm4, mm7 ;1
- psubw mm0, mm7 ;2
- psubw mm6, mm3 ;3
- movq mm1, mm2 ; 03 02 01 00
- movq mm3, mm4 ; 23 22 21 20
- punpcklwd mm1, mm0 ; 11 01 10 00
- punpckhwd mm2, mm0 ; 13 03 12 02
- punpcklwd mm3, mm6 ; 31 21 30 20
- punpckhwd mm4, mm6 ; 33 23 32 22
- movq mm0, mm1 ; 11 01 10 00
- movq mm5, mm2 ; 13 03 12 02
- punpckldq mm0, mm3 ; 30 20 10 00
- punpckhdq mm1, mm3 ; 31 21 11 01
- punpckldq mm2, mm4 ; 32 22 12 02
- punpckhdq mm5, mm4 ; 33 23 13 03
- movq mm3, mm5 ; 33 23 13 03
- psubw mm0, mm2 ; b1= 0-2
- paddw mm2, mm2 ;
- movq mm5, mm1
- paddw mm2, mm0 ; a1 =0+2
- pmulhw mm5, [GLOBAL(x_s1sqr2)];
- paddw mm5, mm1 ; ip1 * sin(pi/8) * sqrt(2)
- movq mm7, mm3 ;
- pmulhw mm7, [GLOBAL(x_c1sqr2less1)];
- paddw mm7, mm3 ; ip3 * cos(pi/8) * sqrt(2)
- psubw mm7, mm5 ; c1
- movq mm5, mm1
- movq mm4, mm3
- pmulhw mm5, [GLOBAL(x_c1sqr2less1)]
- paddw mm5, mm1
- pmulhw mm3, [GLOBAL(x_s1sqr2)]
- paddw mm3, mm4
- paddw mm3, mm5 ; d1
- paddw mm0, [GLOBAL(fours)]
- paddw mm2, [GLOBAL(fours)]
- movq mm6, mm2 ; a1
- movq mm4, mm0 ; b1
- paddw mm2, mm3 ;0
- paddw mm4, mm7 ;1
- psubw mm0, mm7 ;2
- psubw mm6, mm3 ;3
- psraw mm2, 3
- psraw mm0, 3
- psraw mm4, 3
- psraw mm6, 3
- movq mm1, mm2 ; 03 02 01 00
- movq mm3, mm4 ; 23 22 21 20
- punpcklwd mm1, mm0 ; 11 01 10 00
- punpckhwd mm2, mm0 ; 13 03 12 02
- punpcklwd mm3, mm6 ; 31 21 30 20
- punpckhwd mm4, mm6 ; 33 23 32 22
- movq mm0, mm1 ; 11 01 10 00
- movq mm5, mm2 ; 13 03 12 02
- punpckldq mm0, mm3 ; 30 20 10 00
- punpckhdq mm1, mm3 ; 31 21 11 01
- punpckldq mm2, mm4 ; 32 22 12 02
- punpckhdq mm5, mm4 ; 33 23 13 03
- pxor mm7, mm7
- movd mm4, [rsi]
- punpcklbw mm4, mm7
- paddsw mm0, mm4
- packuswb mm0, mm7
- movd [rdx], mm0
- movd mm4, [rsi+rax]
- punpcklbw mm4, mm7
- paddsw mm1, mm4
- packuswb mm1, mm7
- movd [rdx+rdi], mm1
- movd mm4, [rsi+2*rax]
- punpcklbw mm4, mm7
- paddsw mm2, mm4
- packuswb mm2, mm7
- movd [rdx+rdi*2], mm2
- add rdx, rdi
- add rsi, rax
- movd mm4, [rsi+2*rax]
- punpcklbw mm4, mm7
- paddsw mm5, mm4
- packuswb mm5, mm7
- movd [rdx+rdi*2], mm5
- ; begin epilog
- pop rdi
- pop rsi
- RESTORE_GOT
- UNSHADOW_ARGS
- pop rbp
- ret
- ;void vp8_dc_only_idct_add_mmx(
- ;short input_dc,
- ;unsigned char *pred_ptr,
- ;int pred_stride,
- ;unsigned char *dst_ptr,
- ;int stride)
- global sym(vp8_dc_only_idct_add_mmx) PRIVATE
- sym(vp8_dc_only_idct_add_mmx):
- push rbp
- mov rbp, rsp
- SHADOW_ARGS_TO_STACK 5
- GET_GOT rbx
- ; end prolog
- movd mm5, arg(0) ;input_dc
- mov rax, arg(1) ;pred_ptr
- movsxd rdx, dword ptr arg(2) ;pred_stride
- pxor mm0, mm0
- paddw mm5, [GLOBAL(fours)]
- lea rcx, [rdx + rdx*2]
- psraw mm5, 3
- punpcklwd mm5, mm5
- punpckldq mm5, mm5
- movd mm1, [rax]
- movd mm2, [rax+rdx]
- movd mm3, [rax+2*rdx]
- movd mm4, [rax+rcx]
- mov rax, arg(3) ;d -- destination
- movsxd rdx, dword ptr arg(4) ;dst_stride
- punpcklbw mm1, mm0
- paddsw mm1, mm5
- packuswb mm1, mm0 ; pack and unpack to saturate
- lea rcx, [rdx + rdx*2]
- punpcklbw mm2, mm0
- paddsw mm2, mm5
- packuswb mm2, mm0 ; pack and unpack to saturate
- punpcklbw mm3, mm0
- paddsw mm3, mm5
- packuswb mm3, mm0 ; pack and unpack to saturate
- punpcklbw mm4, mm0
- paddsw mm4, mm5
- packuswb mm4, mm0 ; pack and unpack to saturate
- movd [rax], mm1
- movd [rax+rdx], mm2
- movd [rax+2*rdx], mm3
- movd [rax+rcx], mm4
- ; begin epilog
- RESTORE_GOT
- UNSHADOW_ARGS
- pop rbp
- ret
- SECTION_RODATA
- align 16
- x_s1sqr2:
- times 4 dw 0x8A8C
- align 16
- x_c1sqr2less1:
- times 4 dw 0x4E7B
- align 16
- fours:
- times 4 dw 0x0004
|
