123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208 |
- %include "jsimdext.inc"
- %include "jdct.inc"
- SECTION SEG_TEXT
- BITS 32
- %define sample_data ebp + 8
- %define start_col ebp + 12
- %define workspace ebp + 16
- align 32
- GLOBAL_FUNCTION(jsimd_convsamp_float_sse)
- EXTN(jsimd_convsamp_float_sse):
- push ebp
- mov ebp, esp
- push ebx
- push esi
- push edi
- pcmpeqw mm7, mm7
- psllw mm7, 7
- packsswb mm7, mm7
- mov esi, JSAMPARRAY [sample_data]
- mov eax, JDIMENSION [start_col]
- mov edi, POINTER [workspace]
- mov ecx, DCTSIZE/2
- alignx 16, 7
- .convloop:
- mov ebx, JSAMPROW [esi+0*SIZEOF_JSAMPROW]
- mov edx, JSAMPROW [esi+1*SIZEOF_JSAMPROW]
- movq mm0, MMWORD [ebx+eax*SIZEOF_JSAMPLE]
- movq mm1, MMWORD [edx+eax*SIZEOF_JSAMPLE]
- psubb mm0, mm7
- psubb mm1, mm7
- punpcklbw mm2, mm0
- punpckhbw mm0, mm0
- punpcklbw mm3, mm1
- punpckhbw mm1, mm1
- punpcklwd mm4, mm2
- punpckhwd mm2, mm2
- punpcklwd mm5, mm0
- punpckhwd mm0, mm0
- psrad mm4, (DWORD_BIT-BYTE_BIT)
- psrad mm2, (DWORD_BIT-BYTE_BIT)
- cvtpi2ps xmm0, mm4
- cvtpi2ps xmm1, mm2
- psrad mm5, (DWORD_BIT-BYTE_BIT)
- psrad mm0, (DWORD_BIT-BYTE_BIT)
- cvtpi2ps xmm2, mm5
- cvtpi2ps xmm3, mm0
- punpcklwd mm6, mm3
- punpckhwd mm3, mm3
- punpcklwd mm4, mm1
- punpckhwd mm1, mm1
- psrad mm6, (DWORD_BIT-BYTE_BIT)
- psrad mm3, (DWORD_BIT-BYTE_BIT)
- cvtpi2ps xmm4, mm6
- cvtpi2ps xmm5, mm3
- psrad mm4, (DWORD_BIT-BYTE_BIT)
- psrad mm1, (DWORD_BIT-BYTE_BIT)
- cvtpi2ps xmm6, mm4
- cvtpi2ps xmm7, mm1
- movlhps xmm0, xmm1
- movlhps xmm2, xmm3
- movlhps xmm4, xmm5
- movlhps xmm6, xmm7
- movaps XMMWORD [XMMBLOCK(0,0,edi,SIZEOF_FAST_FLOAT)], xmm0
- movaps XMMWORD [XMMBLOCK(0,1,edi,SIZEOF_FAST_FLOAT)], xmm2
- movaps XMMWORD [XMMBLOCK(1,0,edi,SIZEOF_FAST_FLOAT)], xmm4
- movaps XMMWORD [XMMBLOCK(1,1,edi,SIZEOF_FAST_FLOAT)], xmm6
- add esi, byte 2*SIZEOF_JSAMPROW
- add edi, byte 2*DCTSIZE*SIZEOF_FAST_FLOAT
- dec ecx
- jnz near .convloop
- emms
- pop edi
- pop esi
- pop ebx
- pop ebp
- ret
- %define coef_block ebp + 8
- %define divisors ebp + 12
- %define workspace ebp + 16
- align 32
- GLOBAL_FUNCTION(jsimd_quantize_float_sse)
- EXTN(jsimd_quantize_float_sse):
- push ebp
- mov ebp, esp
- push esi
- push edi
- mov esi, POINTER [workspace]
- mov edx, POINTER [divisors]
- mov edi, JCOEFPTR [coef_block]
- mov eax, DCTSIZE2/16
- alignx 16, 7
- .quantloop:
- movaps xmm0, XMMWORD [XMMBLOCK(0,0,esi,SIZEOF_FAST_FLOAT)]
- movaps xmm1, XMMWORD [XMMBLOCK(0,1,esi,SIZEOF_FAST_FLOAT)]
- mulps xmm0, XMMWORD [XMMBLOCK(0,0,edx,SIZEOF_FAST_FLOAT)]
- mulps xmm1, XMMWORD [XMMBLOCK(0,1,edx,SIZEOF_FAST_FLOAT)]
- movaps xmm2, XMMWORD [XMMBLOCK(1,0,esi,SIZEOF_FAST_FLOAT)]
- movaps xmm3, XMMWORD [XMMBLOCK(1,1,esi,SIZEOF_FAST_FLOAT)]
- mulps xmm2, XMMWORD [XMMBLOCK(1,0,edx,SIZEOF_FAST_FLOAT)]
- mulps xmm3, XMMWORD [XMMBLOCK(1,1,edx,SIZEOF_FAST_FLOAT)]
- movhlps xmm4, xmm0
- movhlps xmm5, xmm1
- cvtps2pi mm0, xmm0
- cvtps2pi mm1, xmm1
- cvtps2pi mm4, xmm4
- cvtps2pi mm5, xmm5
- movhlps xmm6, xmm2
- movhlps xmm7, xmm3
- cvtps2pi mm2, xmm2
- cvtps2pi mm3, xmm3
- cvtps2pi mm6, xmm6
- cvtps2pi mm7, xmm7
- packssdw mm0, mm4
- packssdw mm1, mm5
- packssdw mm2, mm6
- packssdw mm3, mm7
- movq MMWORD [MMBLOCK(0,0,edi,SIZEOF_JCOEF)], mm0
- movq MMWORD [MMBLOCK(0,1,edi,SIZEOF_JCOEF)], mm1
- movq MMWORD [MMBLOCK(1,0,edi,SIZEOF_JCOEF)], mm2
- movq MMWORD [MMBLOCK(1,1,edi,SIZEOF_JCOEF)], mm3
- add esi, byte 16*SIZEOF_FAST_FLOAT
- add edx, byte 16*SIZEOF_FAST_FLOAT
- add edi, byte 16*SIZEOF_JCOEF
- dec eax
- jnz short .quantloop
- emms
- pop edi
- pop esi
- pop ebp
- ret
- align 32
|