123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501 |
- %include "jsimdext.inc"
- %include "jdct.inc"
- %define CONST_BITS 8
- %define PASS1_BITS 2
- %if IFAST_SCALE_BITS != PASS1_BITS
- %error "'IFAST_SCALE_BITS' must be equal to 'PASS1_BITS'."
- %endif
- %if CONST_BITS == 8
- F_1_082 equ 277
- F_1_414 equ 362
- F_1_847 equ 473
- F_2_613 equ 669
- F_1_613 equ (F_2_613 - 256)
- %else
- %define DESCALE(x, n) (((x) + (1 << ((n) - 1))) >> (n))
- F_1_082 equ DESCALE(1162209775, 30 - CONST_BITS)
- F_1_414 equ DESCALE(1518500249, 30 - CONST_BITS)
- F_1_847 equ DESCALE(1984016188, 30 - CONST_BITS)
- F_2_613 equ DESCALE(2805822602, 30 - CONST_BITS)
- F_1_613 equ (F_2_613 - (1 << CONST_BITS))
- %endif
- SECTION SEG_CONST
- %define PRE_MULTIPLY_SCALE_BITS 2
- %define CONST_SHIFT (16 - PRE_MULTIPLY_SCALE_BITS - CONST_BITS)
- alignz 32
- GLOBAL_DATA(jconst_idct_ifast_sse2)
- EXTN(jconst_idct_ifast_sse2):
- PW_F1414 times 8 dw F_1_414 << CONST_SHIFT
- PW_F1847 times 8 dw F_1_847 << CONST_SHIFT
- PW_MF1613 times 8 dw -F_1_613 << CONST_SHIFT
- PW_F1082 times 8 dw F_1_082 << CONST_SHIFT
- PB_CENTERJSAMP times 16 db CENTERJSAMPLE
- alignz 32
- SECTION SEG_TEXT
- BITS 32
- %define dct_table(b) (b) + 8
- %define coef_block(b) (b) + 12
- %define output_buf(b) (b) + 16
- %define output_col(b) (b) + 20
- %define original_ebp ebp + 0
- %define wk(i) ebp - (WK_NUM - (i)) * SIZEOF_XMMWORD
-
- %define WK_NUM 2
- align 32
- GLOBAL_FUNCTION(jsimd_idct_ifast_sse2)
- EXTN(jsimd_idct_ifast_sse2):
- push ebp
- mov eax, esp
- sub esp, byte 4
- and esp, byte (-SIZEOF_XMMWORD)
- mov [esp], eax
- mov ebp, esp
- lea esp, [wk(0)]
- pushpic ebx
- push esi
- push edi
- get_GOT ebx
-
- mov edx, POINTER [dct_table(eax)]
- mov esi, JCOEFPTR [coef_block(eax)]
- %ifndef NO_ZERO_COLUMN_TEST_IFAST_SSE2
- mov eax, dword [DWBLOCK(1,0,esi,SIZEOF_JCOEF)]
- or eax, dword [DWBLOCK(2,0,esi,SIZEOF_JCOEF)]
- jnz near .columnDCT
- movdqa xmm0, XMMWORD [XMMBLOCK(1,0,esi,SIZEOF_JCOEF)]
- movdqa xmm1, XMMWORD [XMMBLOCK(2,0,esi,SIZEOF_JCOEF)]
- por xmm0, XMMWORD [XMMBLOCK(3,0,esi,SIZEOF_JCOEF)]
- por xmm1, XMMWORD [XMMBLOCK(4,0,esi,SIZEOF_JCOEF)]
- por xmm0, XMMWORD [XMMBLOCK(5,0,esi,SIZEOF_JCOEF)]
- por xmm1, XMMWORD [XMMBLOCK(6,0,esi,SIZEOF_JCOEF)]
- por xmm0, XMMWORD [XMMBLOCK(7,0,esi,SIZEOF_JCOEF)]
- por xmm1, xmm0
- packsswb xmm1, xmm1
- packsswb xmm1, xmm1
- movd eax, xmm1
- test eax, eax
- jnz short .columnDCT
-
- movdqa xmm0, XMMWORD [XMMBLOCK(0,0,esi,SIZEOF_JCOEF)]
- pmullw xmm0, XMMWORD [XMMBLOCK(0,0,edx,SIZEOF_ISLOW_MULT_TYPE)]
- movdqa xmm7, xmm0
- punpcklwd xmm0, xmm0
- punpckhwd xmm7, xmm7
- pshufd xmm6, xmm0, 0x00
- pshufd xmm2, xmm0, 0x55
- pshufd xmm5, xmm0, 0xAA
- pshufd xmm0, xmm0, 0xFF
- pshufd xmm1, xmm7, 0x00
- pshufd xmm4, xmm7, 0x55
- pshufd xmm3, xmm7, 0xAA
- pshufd xmm7, xmm7, 0xFF
- movdqa XMMWORD [wk(0)], xmm2
- movdqa XMMWORD [wk(1)], xmm0
- jmp near .column_end
- alignx 16, 7
- %endif
- .columnDCT:
-
- movdqa xmm0, XMMWORD [XMMBLOCK(0,0,esi,SIZEOF_JCOEF)]
- movdqa xmm1, XMMWORD [XMMBLOCK(2,0,esi,SIZEOF_JCOEF)]
- pmullw xmm0, XMMWORD [XMMBLOCK(0,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- pmullw xmm1, XMMWORD [XMMBLOCK(2,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- movdqa xmm2, XMMWORD [XMMBLOCK(4,0,esi,SIZEOF_JCOEF)]
- movdqa xmm3, XMMWORD [XMMBLOCK(6,0,esi,SIZEOF_JCOEF)]
- pmullw xmm2, XMMWORD [XMMBLOCK(4,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- pmullw xmm3, XMMWORD [XMMBLOCK(6,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- movdqa xmm4, xmm0
- movdqa xmm5, xmm1
- psubw xmm0, xmm2
- psubw xmm1, xmm3
- paddw xmm4, xmm2
- paddw xmm5, xmm3
- psllw xmm1, PRE_MULTIPLY_SCALE_BITS
- pmulhw xmm1, [GOTOFF(ebx,PW_F1414)]
- psubw xmm1, xmm5
- movdqa xmm6, xmm4
- movdqa xmm7, xmm0
- psubw xmm4, xmm5
- psubw xmm0, xmm1
- paddw xmm6, xmm5
- paddw xmm7, xmm1
- movdqa XMMWORD [wk(1)], xmm4
- movdqa XMMWORD [wk(0)], xmm0
-
- movdqa xmm2, XMMWORD [XMMBLOCK(1,0,esi,SIZEOF_JCOEF)]
- movdqa xmm3, XMMWORD [XMMBLOCK(3,0,esi,SIZEOF_JCOEF)]
- pmullw xmm2, XMMWORD [XMMBLOCK(1,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- pmullw xmm3, XMMWORD [XMMBLOCK(3,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- movdqa xmm5, XMMWORD [XMMBLOCK(5,0,esi,SIZEOF_JCOEF)]
- movdqa xmm1, XMMWORD [XMMBLOCK(7,0,esi,SIZEOF_JCOEF)]
- pmullw xmm5, XMMWORD [XMMBLOCK(5,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- pmullw xmm1, XMMWORD [XMMBLOCK(7,0,edx,SIZEOF_IFAST_MULT_TYPE)]
- movdqa xmm4, xmm2
- movdqa xmm0, xmm5
- psubw xmm2, xmm1
- psubw xmm5, xmm3
- paddw xmm4, xmm1
- paddw xmm0, xmm3
- movdqa xmm1, xmm5
- psllw xmm2, PRE_MULTIPLY_SCALE_BITS
- psllw xmm5, PRE_MULTIPLY_SCALE_BITS
- movdqa xmm3, xmm4
- psubw xmm4, xmm0
- paddw xmm3, xmm0
- psllw xmm4, PRE_MULTIPLY_SCALE_BITS
- pmulhw xmm4, [GOTOFF(ebx,PW_F1414)]
-
-
-
-
-
-
-
-
- movdqa xmm0, xmm5
- paddw xmm5, xmm2
- pmulhw xmm5, [GOTOFF(ebx,PW_F1847)]
- pmulhw xmm0, [GOTOFF(ebx,PW_MF1613)]
- pmulhw xmm2, [GOTOFF(ebx,PW_F1082)]
- psubw xmm0, xmm1
- psubw xmm2, xmm5
- paddw xmm0, xmm5
-
- psubw xmm0, xmm3
- movdqa xmm1, xmm6
- movdqa xmm5, xmm7
- paddw xmm6, xmm3
- paddw xmm7, xmm0
- psubw xmm1, xmm3
- psubw xmm5, xmm0
- psubw xmm4, xmm0
- movdqa xmm3, xmm6
- punpcklwd xmm6, xmm7
- punpckhwd xmm3, xmm7
- movdqa xmm0, xmm5
- punpcklwd xmm5, xmm1
- punpckhwd xmm0, xmm1
- movdqa xmm7, XMMWORD [wk(0)]
- movdqa xmm1, XMMWORD [wk(1)]
- movdqa XMMWORD [wk(0)], xmm5
- movdqa XMMWORD [wk(1)], xmm0
- paddw xmm2, xmm4
- movdqa xmm5, xmm7
- movdqa xmm0, xmm1
- paddw xmm7, xmm4
- paddw xmm1, xmm2
- psubw xmm5, xmm4
- psubw xmm0, xmm2
- movdqa xmm4, xmm7
- punpcklwd xmm7, xmm0
- punpckhwd xmm4, xmm0
- movdqa xmm2, xmm1
- punpcklwd xmm1, xmm5
- punpckhwd xmm2, xmm5
- movdqa xmm0, xmm3
- punpckldq xmm3, xmm4
- punpckhdq xmm0, xmm4
- movdqa xmm5, xmm6
- punpckldq xmm6, xmm7
- punpckhdq xmm5, xmm7
- movdqa xmm4, XMMWORD [wk(0)]
- movdqa xmm7, XMMWORD [wk(1)]
- movdqa XMMWORD [wk(0)], xmm3
- movdqa XMMWORD [wk(1)], xmm0
- movdqa xmm3, xmm1
- punpckldq xmm1, xmm4
- punpckhdq xmm3, xmm4
- movdqa xmm0, xmm2
- punpckldq xmm2, xmm7
- punpckhdq xmm0, xmm7
- movdqa xmm4, xmm6
- punpcklqdq xmm6, xmm1
- punpckhqdq xmm4, xmm1
- movdqa xmm7, xmm5
- punpcklqdq xmm5, xmm3
- punpckhqdq xmm7, xmm3
- movdqa xmm1, XMMWORD [wk(0)]
- movdqa xmm3, XMMWORD [wk(1)]
- movdqa XMMWORD [wk(0)], xmm4
- movdqa XMMWORD [wk(1)], xmm7
- movdqa xmm4, xmm1
- punpcklqdq xmm1, xmm2
- punpckhqdq xmm4, xmm2
- movdqa xmm7, xmm3
- punpcklqdq xmm3, xmm0
- punpckhqdq xmm7, xmm0
- .column_end:
-
- prefetchnta [esi + DCTSIZE2*SIZEOF_JCOEF + 0*32]
- prefetchnta [esi + DCTSIZE2*SIZEOF_JCOEF + 1*32]
- prefetchnta [esi + DCTSIZE2*SIZEOF_JCOEF + 2*32]
- prefetchnta [esi + DCTSIZE2*SIZEOF_JCOEF + 3*32]
-
- mov eax, [original_ebp]
- mov edi, JSAMPARRAY [output_buf(eax)]
- mov eax, JDIMENSION [output_col(eax)]
-
-
- movdqa xmm2, xmm6
- movdqa xmm0, xmm5
- psubw xmm6, xmm1
- psubw xmm5, xmm3
- paddw xmm2, xmm1
- paddw xmm0, xmm3
- psllw xmm5, PRE_MULTIPLY_SCALE_BITS
- pmulhw xmm5, [GOTOFF(ebx,PW_F1414)]
- psubw xmm5, xmm0
- movdqa xmm1, xmm2
- movdqa xmm3, xmm6
- psubw xmm2, xmm0
- psubw xmm6, xmm5
- paddw xmm1, xmm0
- paddw xmm3, xmm5
- movdqa xmm0, XMMWORD [wk(0)]
- movdqa xmm5, XMMWORD [wk(1)]
- movdqa XMMWORD [wk(0)], xmm2
- movdqa XMMWORD [wk(1)], xmm6
-
-
- movdqa xmm2, xmm0
- movdqa xmm6, xmm4
- psubw xmm0, xmm7
- psubw xmm4, xmm5
- paddw xmm2, xmm7
- paddw xmm6, xmm5
- movdqa xmm7, xmm4
- psllw xmm0, PRE_MULTIPLY_SCALE_BITS
- psllw xmm4, PRE_MULTIPLY_SCALE_BITS
- movdqa xmm5, xmm2
- psubw xmm2, xmm6
- paddw xmm5, xmm6
- psllw xmm2, PRE_MULTIPLY_SCALE_BITS
- pmulhw xmm2, [GOTOFF(ebx,PW_F1414)]
-
-
-
-
-
-
-
-
- movdqa xmm6, xmm4
- paddw xmm4, xmm0
- pmulhw xmm4, [GOTOFF(ebx,PW_F1847)]
- pmulhw xmm6, [GOTOFF(ebx,PW_MF1613)]
- pmulhw xmm0, [GOTOFF(ebx,PW_F1082)]
- psubw xmm6, xmm7
- psubw xmm0, xmm4
- paddw xmm6, xmm4
-
- psubw xmm6, xmm5
- movdqa xmm7, xmm1
- movdqa xmm4, xmm3
- paddw xmm1, xmm5
- paddw xmm3, xmm6
- psraw xmm1, (PASS1_BITS+3)
- psraw xmm3, (PASS1_BITS+3)
- psubw xmm7, xmm5
- psubw xmm4, xmm6
- psraw xmm7, (PASS1_BITS+3)
- psraw xmm4, (PASS1_BITS+3)
- psubw xmm2, xmm6
- packsswb xmm1, xmm4
- packsswb xmm3, xmm7
- movdqa xmm5, XMMWORD [wk(1)]
- movdqa xmm6, XMMWORD [wk(0)]
- paddw xmm0, xmm2
- movdqa xmm4, xmm5
- movdqa xmm7, xmm6
- paddw xmm5, xmm2
- paddw xmm6, xmm0
- psraw xmm5, (PASS1_BITS+3)
- psraw xmm6, (PASS1_BITS+3)
- psubw xmm4, xmm2
- psubw xmm7, xmm0
- psraw xmm4, (PASS1_BITS+3)
- psraw xmm7, (PASS1_BITS+3)
- movdqa xmm2, [GOTOFF(ebx,PB_CENTERJSAMP)]
- packsswb xmm5, xmm6
- packsswb xmm7, xmm4
- paddb xmm1, xmm2
- paddb xmm3, xmm2
- paddb xmm5, xmm2
- paddb xmm7, xmm2
- movdqa xmm0, xmm1
- punpcklbw xmm1, xmm3
- punpckhbw xmm0, xmm3
- movdqa xmm6, xmm5
- punpcklbw xmm5, xmm7
- punpckhbw xmm6, xmm7
- movdqa xmm4, xmm1
- punpcklwd xmm1, xmm5
- punpckhwd xmm4, xmm5
- movdqa xmm2, xmm6
- punpcklwd xmm6, xmm0
- punpckhwd xmm2, xmm0
- movdqa xmm3, xmm1
- punpckldq xmm1, xmm6
- punpckhdq xmm3, xmm6
- movdqa xmm7, xmm4
- punpckldq xmm4, xmm2
- punpckhdq xmm7, xmm2
- pshufd xmm5, xmm1, 0x4E
- pshufd xmm0, xmm3, 0x4E
- pshufd xmm6, xmm4, 0x4E
- pshufd xmm2, xmm7, 0x4E
- mov edx, JSAMPROW [edi+0*SIZEOF_JSAMPROW]
- mov esi, JSAMPROW [edi+2*SIZEOF_JSAMPROW]
- movq XMM_MMWORD [edx+eax*SIZEOF_JSAMPLE], xmm1
- movq XMM_MMWORD [esi+eax*SIZEOF_JSAMPLE], xmm3
- mov edx, JSAMPROW [edi+4*SIZEOF_JSAMPROW]
- mov esi, JSAMPROW [edi+6*SIZEOF_JSAMPROW]
- movq XMM_MMWORD [edx+eax*SIZEOF_JSAMPLE], xmm4
- movq XMM_MMWORD [esi+eax*SIZEOF_JSAMPLE], xmm7
- mov edx, JSAMPROW [edi+1*SIZEOF_JSAMPROW]
- mov esi, JSAMPROW [edi+3*SIZEOF_JSAMPROW]
- movq XMM_MMWORD [edx+eax*SIZEOF_JSAMPLE], xmm5
- movq XMM_MMWORD [esi+eax*SIZEOF_JSAMPLE], xmm0
- mov edx, JSAMPROW [edi+5*SIZEOF_JSAMPROW]
- mov esi, JSAMPROW [edi+7*SIZEOF_JSAMPROW]
- movq XMM_MMWORD [edx+eax*SIZEOF_JSAMPLE], xmm6
- movq XMM_MMWORD [esi+eax*SIZEOF_JSAMPLE], xmm2
- pop edi
- pop esi
- poppic ebx
- mov esp, ebp
- pop esp
- pop ebp
- ret
- align 32
|