shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
; jchuff-sse2.asm - Huffman entropy encoding (SSE2)
shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
; Copyright (C) 2009-2011, 2014-2017, D. R. Commander.
shun-iwasawa 82a8f5
; Copyright (C) 2015, Matthieu Darbois.
shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
; Based on the x86 SIMD extension for IJG JPEG library
shun-iwasawa 82a8f5
; Copyright (C) 1999-2006, MIYASAKA Masaru.
shun-iwasawa 82a8f5
; For conditions of distribution and use, see copyright notice in jsimdext.inc
shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
; This file should be assembled with NASM (Netwide Assembler),
shun-iwasawa 82a8f5
; can *not* be assembled with Microsoft's MASM or any compatible
shun-iwasawa 82a8f5
; assembler (including Borland's Turbo Assembler).
shun-iwasawa 82a8f5
; NASM is available from http://nasm.sourceforge.net/ or
shun-iwasawa 82a8f5
; http://sourceforge.net/project/showfiles.php?group_id=6208
shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
; This file contains an SSE2 implementation for Huffman coding of one block.
shun-iwasawa 82a8f5
; The following code is based directly on jchuff.c; see jchuff.c for more
shun-iwasawa 82a8f5
; details.
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%include "jsimdext.inc"
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
; --------------------------------------------------------------------------
shun-iwasawa 82a8f5
    SECTION     SEG_CONST
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    alignz      32
shun-iwasawa 82a8f5
    GLOBAL_DATA(jconst_huff_encode_one_block)
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
EXTN(jconst_huff_encode_one_block):
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%include "jpeg_nbits_table.inc"
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    alignz      32
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
; --------------------------------------------------------------------------
shun-iwasawa 82a8f5
    SECTION     SEG_TEXT
shun-iwasawa 82a8f5
    BITS        32
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
; These macros perform the same task as the emit_bits() function in the
shun-iwasawa 82a8f5
; original libjpeg code.  In addition to reducing overhead by explicitly
shun-iwasawa 82a8f5
; inlining the code, additional performance is achieved by taking into
shun-iwasawa 82a8f5
; account the size of the bit buffer and waiting until it is almost full
shun-iwasawa 82a8f5
; before emptying it.  This mostly benefits 64-bit platforms, since 6
shun-iwasawa 82a8f5
; bytes can be stored in a 64-bit bit buffer before it has to be emptied.
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%macro EMIT_BYTE 0
shun-iwasawa 82a8f5
    sub         put_bits, 8             ; put_bits -= 8;
shun-iwasawa 82a8f5
    mov         edx, put_buffer
shun-iwasawa 82a8f5
    mov         ecx, put_bits
shun-iwasawa 82a8f5
    shr         edx, cl                 ; c = (JOCTET)GETJOCTET(put_buffer >> put_bits);
shun-iwasawa 82a8f5
    mov         byte [eax], dl          ; *buffer++ = c;
shun-iwasawa 82a8f5
    add         eax, 1
shun-iwasawa 82a8f5
    cmp         dl, 0xFF                ; need to stuff a zero byte?
shun-iwasawa 82a8f5
    jne         %%.EMIT_BYTE_END
shun-iwasawa 82a8f5
    mov         byte [eax], 0           ; *buffer++ = 0;
shun-iwasawa 82a8f5
    add         eax, 1
shun-iwasawa 82a8f5
%%.EMIT_BYTE_END:
shun-iwasawa 82a8f5
%endmacro
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%macro PUT_BITS 1
shun-iwasawa 82a8f5
    add         put_bits, ecx           ; put_bits += size;
shun-iwasawa 82a8f5
    shl         put_buffer, cl          ; put_buffer = (put_buffer << size);
shun-iwasawa 82a8f5
    or          put_buffer, %1
shun-iwasawa 82a8f5
%endmacro
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%macro CHECKBUF15 0
shun-iwasawa 82a8f5
    cmp         put_bits, 16            ; if (put_bits > 31) {
shun-iwasawa 82a8f5
    jl          %%.CHECKBUF15_END
shun-iwasawa 82a8f5
    mov         eax, POINTER [esp+buffer]
shun-iwasawa 82a8f5
    EMIT_BYTE
shun-iwasawa 82a8f5
    EMIT_BYTE
shun-iwasawa 82a8f5
    mov         POINTER [esp+buffer], eax
shun-iwasawa 82a8f5
%%.CHECKBUF15_END:
shun-iwasawa 82a8f5
%endmacro
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%macro EMIT_BITS 1
shun-iwasawa 82a8f5
    PUT_BITS    %1
shun-iwasawa 82a8f5
    CHECKBUF15
shun-iwasawa 82a8f5
%endmacro
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%macro kloop_prepare 37                 ;(ko, jno0, ..., jno31, xmm0, xmm1, xmm2, xmm3)
shun-iwasawa 82a8f5
    pxor        xmm4, xmm4              ; __m128i neg = _mm_setzero_si128();
shun-iwasawa 82a8f5
    pxor        xmm5, xmm5              ; __m128i neg = _mm_setzero_si128();
shun-iwasawa 82a8f5
    pxor        xmm6, xmm6              ; __m128i neg = _mm_setzero_si128();
shun-iwasawa 82a8f5
    pxor        xmm7, xmm7              ; __m128i neg = _mm_setzero_si128();
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %2  * SIZEOF_WORD], 0  ; xmm_shadow[0] = block[jno0];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %10 * SIZEOF_WORD], 0  ; xmm_shadow[8] = block[jno8];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %18 * SIZEOF_WORD], 0  ; xmm_shadow[16] = block[jno16];
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %26 * SIZEOF_WORD], 0  ; xmm_shadow[24] = block[jno24];
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %3  * SIZEOF_WORD], 1  ; xmm_shadow[1] = block[jno1];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %11 * SIZEOF_WORD], 1  ; xmm_shadow[9] = block[jno9];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %19 * SIZEOF_WORD], 1  ; xmm_shadow[17] = block[jno17];
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %27 * SIZEOF_WORD], 1  ; xmm_shadow[25] = block[jno25];
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %4  * SIZEOF_WORD], 2  ; xmm_shadow[2] = block[jno2];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %12 * SIZEOF_WORD], 2  ; xmm_shadow[10] = block[jno10];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %20 * SIZEOF_WORD], 2  ; xmm_shadow[18] = block[jno18];
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %28 * SIZEOF_WORD], 2  ; xmm_shadow[26] = block[jno26];
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %5  * SIZEOF_WORD], 3  ; xmm_shadow[3] = block[jno3];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %13 * SIZEOF_WORD], 3  ; xmm_shadow[11] = block[jno11];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %21 * SIZEOF_WORD], 3  ; xmm_shadow[19] = block[jno19];
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %29 * SIZEOF_WORD], 3  ; xmm_shadow[27] = block[jno27];
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %6  * SIZEOF_WORD], 4  ; xmm_shadow[4] = block[jno4];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %14 * SIZEOF_WORD], 4  ; xmm_shadow[12] = block[jno12];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %22 * SIZEOF_WORD], 4  ; xmm_shadow[20] = block[jno20];
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %30 * SIZEOF_WORD], 4  ; xmm_shadow[28] = block[jno28];
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %7  * SIZEOF_WORD], 5  ; xmm_shadow[5] = block[jno5];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %15 * SIZEOF_WORD], 5  ; xmm_shadow[13] = block[jno13];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %23 * SIZEOF_WORD], 5  ; xmm_shadow[21] = block[jno21];
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %31 * SIZEOF_WORD], 5  ; xmm_shadow[29] = block[jno29];
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %8  * SIZEOF_WORD], 6  ; xmm_shadow[6] = block[jno6];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %16 * SIZEOF_WORD], 6  ; xmm_shadow[14] = block[jno14];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %24 * SIZEOF_WORD], 6  ; xmm_shadow[22] = block[jno22];
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %32 * SIZEOF_WORD], 6  ; xmm_shadow[30] = block[jno30];
shun-iwasawa 82a8f5
    pinsrw      %34, word [esi + %9  * SIZEOF_WORD], 7  ; xmm_shadow[7] = block[jno7];
shun-iwasawa 82a8f5
    pinsrw      %35, word [esi + %17 * SIZEOF_WORD], 7  ; xmm_shadow[15] = block[jno15];
shun-iwasawa 82a8f5
    pinsrw      %36, word [esi + %25 * SIZEOF_WORD], 7  ; xmm_shadow[23] = block[jno23];
shun-iwasawa 82a8f5
%if %1 != 32
shun-iwasawa 82a8f5
    pinsrw      %37, word [esi + %33 * SIZEOF_WORD], 7  ; xmm_shadow[31] = block[jno31];
shun-iwasawa 82a8f5
%else
shun-iwasawa 82a8f5
    pinsrw      %37, ecx, 7             ; xmm_shadow[31] = block[jno31];
shun-iwasawa 82a8f5
%endif
shun-iwasawa 82a8f5
    pcmpgtw     xmm4, %34               ; neg = _mm_cmpgt_epi16(neg, x1);
shun-iwasawa 82a8f5
    pcmpgtw     xmm5, %35               ; neg = _mm_cmpgt_epi16(neg, x1);
shun-iwasawa 82a8f5
    pcmpgtw     xmm6, %36               ; neg = _mm_cmpgt_epi16(neg, x1);
shun-iwasawa 82a8f5
    pcmpgtw     xmm7, %37               ; neg = _mm_cmpgt_epi16(neg, x1);
shun-iwasawa 82a8f5
    paddw       %34, xmm4               ; x1 = _mm_add_epi16(x1, neg);
shun-iwasawa 82a8f5
    paddw       %35, xmm5               ; x1 = _mm_add_epi16(x1, neg);
shun-iwasawa 82a8f5
    paddw       %36, xmm6               ; x1 = _mm_add_epi16(x1, neg);
shun-iwasawa 82a8f5
    paddw       %37, xmm7               ; x1 = _mm_add_epi16(x1, neg);
shun-iwasawa 82a8f5
    pxor        %34, xmm4               ; x1 = _mm_xor_si128(x1, neg);
shun-iwasawa 82a8f5
    pxor        %35, xmm5               ; x1 = _mm_xor_si128(x1, neg);
shun-iwasawa 82a8f5
    pxor        %36, xmm6               ; x1 = _mm_xor_si128(x1, neg);
shun-iwasawa 82a8f5
    pxor        %37, xmm7               ; x1 = _mm_xor_si128(x1, neg);
shun-iwasawa 82a8f5
    pxor        xmm4, %34               ; neg = _mm_xor_si128(neg, x1);
shun-iwasawa 82a8f5
    pxor        xmm5, %35               ; neg = _mm_xor_si128(neg, x1);
shun-iwasawa 82a8f5
    pxor        xmm6, %36               ; neg = _mm_xor_si128(neg, x1);
shun-iwasawa 82a8f5
    pxor        xmm7, %37               ; neg = _mm_xor_si128(neg, x1);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t1 + %1 * SIZEOF_WORD], %34          ; _mm_storeu_si128((__m128i *)(t1 + ko), x1);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t1 + (%1 + 8) * SIZEOF_WORD], %35    ; _mm_storeu_si128((__m128i *)(t1 + ko + 8), x1);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t1 + (%1 + 16) * SIZEOF_WORD], %36   ; _mm_storeu_si128((__m128i *)(t1 + ko + 16), x1);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t1 + (%1 + 24) * SIZEOF_WORD], %37   ; _mm_storeu_si128((__m128i *)(t1 + ko + 24), x1);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t2 + %1 * SIZEOF_WORD], xmm4         ; _mm_storeu_si128((__m128i *)(t2 + ko), neg);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t2 + (%1 + 8) * SIZEOF_WORD], xmm5   ; _mm_storeu_si128((__m128i *)(t2 + ko + 8), neg);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t2 + (%1 + 16) * SIZEOF_WORD], xmm6  ; _mm_storeu_si128((__m128i *)(t2 + ko + 16), neg);
shun-iwasawa 82a8f5
    movdqa      XMMWORD [esp + t2 + (%1 + 24) * SIZEOF_WORD], xmm7  ; _mm_storeu_si128((__m128i *)(t2 + ko + 24), neg);
shun-iwasawa 82a8f5
%endmacro
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
; Encode a single block's worth of coefficients.
shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
; GLOBAL(JOCTET *)
shun-iwasawa 82a8f5
; jsimd_huff_encode_one_block_sse2(working_state *state, JOCTET *buffer,
shun-iwasawa 82a8f5
;                                  JCOEFPTR block, int last_dc_val,
shun-iwasawa 82a8f5
;                                  c_derived_tbl *dctbl, c_derived_tbl *actbl)
shun-iwasawa 82a8f5
;
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
; eax + 8 = working_state *state
shun-iwasawa 82a8f5
; eax + 12 = JOCTET *buffer
shun-iwasawa 82a8f5
; eax + 16 = JCOEFPTR block
shun-iwasawa 82a8f5
; eax + 20 = int last_dc_val
shun-iwasawa 82a8f5
; eax + 24 = c_derived_tbl *dctbl
shun-iwasawa 82a8f5
; eax + 28 = c_derived_tbl *actbl
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
%define pad         6 * SIZEOF_DWORD    ; Align to 16 bytes
shun-iwasawa 82a8f5
%define t1          pad
shun-iwasawa 82a8f5
%define t2          t1 + (DCTSIZE2 * SIZEOF_WORD)
shun-iwasawa 82a8f5
%define block       t2 + (DCTSIZE2 * SIZEOF_WORD)
shun-iwasawa 82a8f5
%define actbl       block + SIZEOF_DWORD
shun-iwasawa 82a8f5
%define buffer      actbl + SIZEOF_DWORD
shun-iwasawa 82a8f5
%define temp        buffer + SIZEOF_DWORD
shun-iwasawa 82a8f5
%define temp2       temp + SIZEOF_DWORD
shun-iwasawa 82a8f5
%define temp3       temp2 + SIZEOF_DWORD
shun-iwasawa 82a8f5
%define temp4       temp3 + SIZEOF_DWORD
shun-iwasawa 82a8f5
%define temp5       temp4 + SIZEOF_DWORD
shun-iwasawa 82a8f5
%define gotptr      temp5 + SIZEOF_DWORD  ; void *gotptr
shun-iwasawa 82a8f5
%define put_buffer  ebx
shun-iwasawa 82a8f5
%define put_bits    edi
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    align       32
shun-iwasawa 82a8f5
    GLOBAL_FUNCTION(jsimd_huff_encode_one_block_sse2)
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
EXTN(jsimd_huff_encode_one_block_sse2):
shun-iwasawa 82a8f5
    push        ebp
shun-iwasawa 82a8f5
    mov         eax, esp                     ; eax = original ebp
shun-iwasawa 82a8f5
    sub         esp, byte 4
shun-iwasawa 82a8f5
    and         esp, byte (-SIZEOF_XMMWORD)  ; align to 128 bits
shun-iwasawa 82a8f5
    mov         [esp], eax
shun-iwasawa 82a8f5
    mov         ebp, esp                     ; ebp = aligned ebp
shun-iwasawa 82a8f5
    sub         esp, temp5+9*SIZEOF_DWORD-pad
shun-iwasawa 82a8f5
    push        ebx
shun-iwasawa 82a8f5
    push        ecx
shun-iwasawa 82a8f5
;   push        edx                     ; need not be preserved
shun-iwasawa 82a8f5
    push        esi
shun-iwasawa 82a8f5
    push        edi
shun-iwasawa 82a8f5
    push        ebp
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    mov         esi, POINTER [eax+8]       ; (working_state *state)
shun-iwasawa 82a8f5
    mov         put_buffer, dword [esi+8]  ; put_buffer = state->cur.put_buffer;
shun-iwasawa 82a8f5
    mov         put_bits, dword [esi+12]   ; put_bits = state->cur.put_bits;
shun-iwasawa 82a8f5
    push        esi                        ; esi is now scratch
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    get_GOT     edx                        ; get GOT address
shun-iwasawa 82a8f5
    movpic      POINTER [esp+gotptr], edx  ; save GOT address
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    mov         ecx, POINTER [eax+28]
shun-iwasawa 82a8f5
    mov         edx, POINTER [eax+16]
shun-iwasawa 82a8f5
    mov         esi, POINTER [eax+12]
shun-iwasawa 82a8f5
    mov         POINTER [esp+actbl], ecx
shun-iwasawa 82a8f5
    mov         POINTER [esp+block], edx
shun-iwasawa 82a8f5
    mov         POINTER [esp+buffer], esi
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; Encode the DC coefficient difference per section F.1.2.1
shun-iwasawa 82a8f5
    mov         esi, POINTER [esp+block]  ; block
shun-iwasawa 82a8f5
    movsx       ecx, word [esi]           ; temp = temp2 = block[0] - last_dc_val;
shun-iwasawa 82a8f5
    sub         ecx, dword [eax+20]
shun-iwasawa 82a8f5
    mov         esi, ecx
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; This is a well-known technique for obtaining the absolute value
shun-iwasawa 82a8f5
    ; with out a branch.  It is derived from an assembly language technique
shun-iwasawa 82a8f5
    ; presented in "How to Optimize for the Pentium Processors",
shun-iwasawa 82a8f5
    ; Copyright (c) 1996, 1997 by Agner Fog.
shun-iwasawa 82a8f5
    mov         edx, ecx
shun-iwasawa 82a8f5
    sar         edx, 31                 ; temp3 = temp >> (CHAR_BIT * sizeof(int) - 1);
shun-iwasawa 82a8f5
    xor         ecx, edx                ; temp ^= temp3;
shun-iwasawa 82a8f5
    sub         ecx, edx                ; temp -= temp3;
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; For a negative input, want temp2 = bitwise complement of abs(input)
shun-iwasawa 82a8f5
    ; This code assumes we are on a two's complement machine
shun-iwasawa 82a8f5
    add         esi, edx                ; temp2 += temp3;
shun-iwasawa 82a8f5
    mov         dword [esp+temp], esi   ; backup temp2 in temp
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; Find the number of bits needed for the magnitude of the coefficient
shun-iwasawa 82a8f5
    movpic      ebp, POINTER [esp+gotptr]                        ; load GOT address (ebp)
shun-iwasawa 82a8f5
    movzx       edx, byte [GOTOFF(ebp, jpeg_nbits_table + ecx)]  ; nbits = JPEG_NBITS(temp);
shun-iwasawa 82a8f5
    mov         dword [esp+temp2], edx                           ; backup nbits in temp2
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; Emit the Huffman-coded symbol for the number of bits
shun-iwasawa 82a8f5
    mov         ebp, POINTER [eax+24]         ; After this point, arguments are not accessible anymore
shun-iwasawa 82a8f5
    mov         eax,  INT [ebp + edx * 4]     ; code = dctbl->ehufco[nbits];
shun-iwasawa 82a8f5
    movzx       ecx, byte [ebp + edx + 1024]  ; size = dctbl->ehufsi[nbits];
shun-iwasawa 82a8f5
    EMIT_BITS   eax                           ; EMIT_BITS(code, size)
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    mov         ecx, dword [esp+temp2]        ; restore nbits
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; Mask off any extra bits in code
shun-iwasawa 82a8f5
    mov         eax, 1
shun-iwasawa 82a8f5
    shl         eax, cl
shun-iwasawa 82a8f5
    dec         eax
shun-iwasawa 82a8f5
    and         eax, dword [esp+temp]   ; temp2 &= (((JLONG)1)<
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; Emit that number of bits of the value, if positive,
shun-iwasawa 82a8f5
    ; or the complement of its magnitude, if negative.
shun-iwasawa 82a8f5
    EMIT_BITS   eax                     ; EMIT_BITS(temp2, nbits)
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    ; Prepare data
shun-iwasawa 82a8f5
    xor         ecx, ecx
shun-iwasawa 82a8f5
    mov         esi, POINTER [esp+block]
shun-iwasawa 82a8f5
    kloop_prepare  0,  1,  8,  16, 9,  2,  3,  10, 17, 24, 32, 25, \
shun-iwasawa 82a8f5
                   18, 11, 4,  5,  12, 19, 26, 33, 40, 48, 41, 34, \
shun-iwasawa 82a8f5
                   27, 20, 13, 6,  7,  14, 21, 28, 35, \
shun-iwasawa 82a8f5
                   xmm0, xmm1, xmm2, xmm3
shun-iwasawa 82a8f5
    kloop_prepare  32, 42, 49, 56, 57, 50, 43, 36, 29, 22, 15, 23, \
shun-iwasawa 82a8f5
                   30, 37, 44, 51, 58, 59, 52, 45, 38, 31, 39, 46, \
shun-iwasawa 82a8f5
                   53, 60, 61, 54, 47, 55, 62, 63, 63, \
shun-iwasawa 82a8f5
                   xmm0, xmm1, xmm2, xmm3
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    pxor        xmm7, xmm7
shun-iwasawa 82a8f5
    movdqa      xmm0, XMMWORD [esp + t1 + 0 * SIZEOF_WORD]   ; __m128i tmp0 = _mm_loadu_si128((__m128i *)(t1 + 0));
shun-iwasawa 82a8f5
    movdqa      xmm1, XMMWORD [esp + t1 + 8 * SIZEOF_WORD]   ; __m128i tmp1 = _mm_loadu_si128((__m128i *)(t1 + 8));
shun-iwasawa 82a8f5
    movdqa      xmm2, XMMWORD [esp + t1 + 16 * SIZEOF_WORD]  ; __m128i tmp2 = _mm_loadu_si128((__m128i *)(t1 + 16));
shun-iwasawa 82a8f5
    movdqa      xmm3, XMMWORD [esp + t1 + 24 * SIZEOF_WORD]  ; __m128i tmp3 = _mm_loadu_si128((__m128i *)(t1 + 24));
shun-iwasawa 82a8f5
    pcmpeqw     xmm0, xmm7              ; tmp0 = _mm_cmpeq_epi16(tmp0, zero);
shun-iwasawa 82a8f5
    pcmpeqw     xmm1, xmm7              ; tmp1 = _mm_cmpeq_epi16(tmp1, zero);
shun-iwasawa 82a8f5
    pcmpeqw     xmm2, xmm7              ; tmp2 = _mm_cmpeq_epi16(tmp2, zero);
shun-iwasawa 82a8f5
    pcmpeqw     xmm3, xmm7              ; tmp3 = _mm_cmpeq_epi16(tmp3, zero);
shun-iwasawa 82a8f5
    packsswb    xmm0, xmm1              ; tmp0 = _mm_packs_epi16(tmp0, tmp1);
shun-iwasawa 82a8f5
    packsswb    xmm2, xmm3              ; tmp2 = _mm_packs_epi16(tmp2, tmp3);
shun-iwasawa 82a8f5
    pmovmskb    edx, xmm0               ; index  = ((uint64_t)_mm_movemask_epi8(tmp0)) << 0;
shun-iwasawa 82a8f5
    pmovmskb    ecx, xmm2               ; index  = ((uint64_t)_mm_movemask_epi8(tmp2)) << 16;
shun-iwasawa 82a8f5
    shl         ecx, 16
shun-iwasawa 82a8f5
    or          edx, ecx
shun-iwasawa 82a8f5
    not         edx                     ; index = ~index;
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    lea         esi, [esp+t1]
shun-iwasawa 82a8f5
    mov         ebp, POINTER [esp+actbl]  ; ebp = actbl
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
.BLOOP:
shun-iwasawa 82a8f5
    bsf         ecx, edx                ; r = __builtin_ctzl(index);
shun-iwasawa 82a8f5
    jz          near .ELOOP
shun-iwasawa 82a8f5
    lea         esi, [esi+ecx*2]        ; k += r;
shun-iwasawa 82a8f5
    shr         edx, cl                 ; index >>= r;
shun-iwasawa 82a8f5
    mov         dword [esp+temp3], edx
shun-iwasawa 82a8f5
.BRLOOP:
shun-iwasawa 82a8f5
    cmp         ecx, 16                       ; while (r > 15) {
shun-iwasawa 82a8f5
    jl          near .ERLOOP
shun-iwasawa 82a8f5
    sub         ecx, 16                       ; r -= 16;
shun-iwasawa 82a8f5
    mov         dword [esp+temp], ecx
shun-iwasawa 82a8f5
    mov         eax, INT [ebp + 240 * 4]      ; code_0xf0 = actbl->ehufco[0xf0];
shun-iwasawa 82a8f5
    movzx       ecx, byte [ebp + 1024 + 240]  ; size_0xf0 = actbl->ehufsi[0xf0];
shun-iwasawa 82a8f5
    EMIT_BITS   eax                           ; EMIT_BITS(code_0xf0, size_0xf0)
shun-iwasawa 82a8f5
    mov         ecx, dword [esp+temp]
shun-iwasawa 82a8f5
    jmp         .BRLOOP
shun-iwasawa 82a8f5
.ERLOOP:
shun-iwasawa 82a8f5
    movsx       eax, word [esi]                                  ; temp = t1[k];
shun-iwasawa 82a8f5
    movpic      edx, POINTER [esp+gotptr]                        ; load GOT address (edx)
shun-iwasawa 82a8f5
    movzx       eax, byte [GOTOFF(edx, jpeg_nbits_table + eax)]  ; nbits = JPEG_NBITS(temp);
shun-iwasawa 82a8f5
    mov         dword [esp+temp2], eax
shun-iwasawa 82a8f5
    ; Emit Huffman symbol for run length / number of bits
shun-iwasawa 82a8f5
    shl         ecx, 4                        ; temp3 = (r << 4) + nbits;
shun-iwasawa 82a8f5
    add         ecx, eax
shun-iwasawa 82a8f5
    mov         eax,  INT [ebp + ecx * 4]     ; code = actbl->ehufco[temp3];
shun-iwasawa 82a8f5
    movzx       ecx, byte [ebp + ecx + 1024]  ; size = actbl->ehufsi[temp3];
shun-iwasawa 82a8f5
    EMIT_BITS   eax
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    movsx       edx, word [esi+DCTSIZE2*2]    ; temp2 = t2[k];
shun-iwasawa 82a8f5
    ; Mask off any extra bits in code
shun-iwasawa 82a8f5
    mov         ecx, dword [esp+temp2]
shun-iwasawa 82a8f5
    mov         eax, 1
shun-iwasawa 82a8f5
    shl         eax, cl
shun-iwasawa 82a8f5
    dec         eax
shun-iwasawa 82a8f5
    and         eax, edx                ; temp2 &= (((JLONG)1)<
shun-iwasawa 82a8f5
    EMIT_BITS   eax                     ; PUT_BITS(temp2, nbits)
shun-iwasawa 82a8f5
    mov         edx, dword [esp+temp3]
shun-iwasawa 82a8f5
    add         esi, 2                  ; ++k;
shun-iwasawa 82a8f5
    shr         edx, 1                  ; index >>= 1;
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    jmp         .BLOOP
shun-iwasawa 82a8f5
.ELOOP:
shun-iwasawa 82a8f5
    movdqa      xmm0, XMMWORD [esp + t1 + 32 * SIZEOF_WORD]  ; __m128i tmp0 = _mm_loadu_si128((__m128i *)(t1 + 0));
shun-iwasawa 82a8f5
    movdqa      xmm1, XMMWORD [esp + t1 + 40 * SIZEOF_WORD]  ; __m128i tmp1 = _mm_loadu_si128((__m128i *)(t1 + 8));
shun-iwasawa 82a8f5
    movdqa      xmm2, XMMWORD [esp + t1 + 48 * SIZEOF_WORD]  ; __m128i tmp2 = _mm_loadu_si128((__m128i *)(t1 + 16));
shun-iwasawa 82a8f5
    movdqa      xmm3, XMMWORD [esp + t1 + 56 * SIZEOF_WORD]  ; __m128i tmp3 = _mm_loadu_si128((__m128i *)(t1 + 24));
shun-iwasawa 82a8f5
    pcmpeqw     xmm0, xmm7              ; tmp0 = _mm_cmpeq_epi16(tmp0, zero);
shun-iwasawa 82a8f5
    pcmpeqw     xmm1, xmm7              ; tmp1 = _mm_cmpeq_epi16(tmp1, zero);
shun-iwasawa 82a8f5
    pcmpeqw     xmm2, xmm7              ; tmp2 = _mm_cmpeq_epi16(tmp2, zero);
shun-iwasawa 82a8f5
    pcmpeqw     xmm3, xmm7              ; tmp3 = _mm_cmpeq_epi16(tmp3, zero);
shun-iwasawa 82a8f5
    packsswb    xmm0, xmm1              ; tmp0 = _mm_packs_epi16(tmp0, tmp1);
shun-iwasawa 82a8f5
    packsswb    xmm2, xmm3              ; tmp2 = _mm_packs_epi16(tmp2, tmp3);
shun-iwasawa 82a8f5
    pmovmskb    edx, xmm0               ; index  = ((uint64_t)_mm_movemask_epi8(tmp0)) << 0;
shun-iwasawa 82a8f5
    pmovmskb    ecx, xmm2               ; index  = ((uint64_t)_mm_movemask_epi8(tmp2)) << 16;
shun-iwasawa 82a8f5
    shl         ecx, 16
shun-iwasawa 82a8f5
    or          edx, ecx
shun-iwasawa 82a8f5
    not         edx                     ; index = ~index;
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    lea         eax, [esp + t1 + (DCTSIZE2/2) * 2]
shun-iwasawa 82a8f5
    sub         eax, esi
shun-iwasawa 82a8f5
    shr         eax, 1
shun-iwasawa 82a8f5
    bsf         ecx, edx                ; r = __builtin_ctzl(index);
shun-iwasawa 82a8f5
    jz          near .ELOOP2
shun-iwasawa 82a8f5
    shr         edx, cl                 ; index >>= r;
shun-iwasawa 82a8f5
    add         ecx, eax
shun-iwasawa 82a8f5
    lea         esi, [esi+ecx*2]        ; k += r;
shun-iwasawa 82a8f5
    mov         dword [esp+temp3], edx
shun-iwasawa 82a8f5
    jmp         .BRLOOP2
shun-iwasawa 82a8f5
.BLOOP2:
shun-iwasawa 82a8f5
    bsf         ecx, edx                ; r = __builtin_ctzl(index);
shun-iwasawa 82a8f5
    jz          near .ELOOP2
shun-iwasawa 82a8f5
    lea         esi, [esi+ecx*2]        ; k += r;
shun-iwasawa 82a8f5
    shr         edx, cl                 ; index >>= r;
shun-iwasawa 82a8f5
    mov         dword [esp+temp3], edx
shun-iwasawa 82a8f5
.BRLOOP2:
shun-iwasawa 82a8f5
    cmp         ecx, 16                       ; while (r > 15) {
shun-iwasawa 82a8f5
    jl          near .ERLOOP2
shun-iwasawa 82a8f5
    sub         ecx, 16                       ; r -= 16;
shun-iwasawa 82a8f5
    mov         dword [esp+temp], ecx
shun-iwasawa 82a8f5
    mov         eax, INT [ebp + 240 * 4]      ; code_0xf0 = actbl->ehufco[0xf0];
shun-iwasawa 82a8f5
    movzx       ecx, byte [ebp + 1024 + 240]  ; size_0xf0 = actbl->ehufsi[0xf0];
shun-iwasawa 82a8f5
    EMIT_BITS   eax                           ; EMIT_BITS(code_0xf0, size_0xf0)
shun-iwasawa 82a8f5
    mov         ecx, dword [esp+temp]
shun-iwasawa 82a8f5
    jmp         .BRLOOP2
shun-iwasawa 82a8f5
.ERLOOP2:
shun-iwasawa 82a8f5
    movsx       eax, word [esi]         ; temp = t1[k];
shun-iwasawa 82a8f5
    bsr         eax, eax                ; nbits = 32 - __builtin_clz(temp);
shun-iwasawa 82a8f5
    inc         eax
shun-iwasawa 82a8f5
    mov         dword [esp+temp2], eax
shun-iwasawa 82a8f5
    ; Emit Huffman symbol for run length / number of bits
shun-iwasawa 82a8f5
    shl         ecx, 4                        ; temp3 = (r << 4) + nbits;
shun-iwasawa 82a8f5
    add         ecx, eax
shun-iwasawa 82a8f5
    mov         eax,  INT [ebp + ecx * 4]     ; code = actbl->ehufco[temp3];
shun-iwasawa 82a8f5
    movzx       ecx, byte [ebp + ecx + 1024]  ; size = actbl->ehufsi[temp3];
shun-iwasawa 82a8f5
    EMIT_BITS   eax
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    movsx       edx, word [esi+DCTSIZE2*2]    ; temp2 = t2[k];
shun-iwasawa 82a8f5
    ; Mask off any extra bits in code
shun-iwasawa 82a8f5
    mov         ecx, dword [esp+temp2]
shun-iwasawa 82a8f5
    mov         eax, 1
shun-iwasawa 82a8f5
    shl         eax, cl
shun-iwasawa 82a8f5
    dec         eax
shun-iwasawa 82a8f5
    and         eax, edx                ; temp2 &= (((JLONG)1)<
shun-iwasawa 82a8f5
    EMIT_BITS   eax                     ; PUT_BITS(temp2, nbits)
shun-iwasawa 82a8f5
    mov         edx, dword [esp+temp3]
shun-iwasawa 82a8f5
    add         esi, 2                  ; ++k;
shun-iwasawa 82a8f5
    shr         edx, 1                  ; index >>= 1;
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    jmp         .BLOOP2
shun-iwasawa 82a8f5
.ELOOP2:
shun-iwasawa 82a8f5
    ; If the last coef(s) were zero, emit an end-of-block code
shun-iwasawa 82a8f5
    lea         edx, [esp + t1 + (DCTSIZE2-1) * 2]  ; r = DCTSIZE2-1-k;
shun-iwasawa 82a8f5
    cmp         edx, esi                            ; if (r > 0) {
shun-iwasawa 82a8f5
    je          .EFN
shun-iwasawa 82a8f5
    mov         eax,  INT [ebp]                     ; code = actbl->ehufco[0];
shun-iwasawa 82a8f5
    movzx       ecx, byte [ebp + 1024]              ; size = actbl->ehufsi[0];
shun-iwasawa 82a8f5
    EMIT_BITS   eax
shun-iwasawa 82a8f5
.EFN:
shun-iwasawa 82a8f5
    mov         eax, [esp+buffer]
shun-iwasawa 82a8f5
    pop         esi
shun-iwasawa 82a8f5
    ; Save put_buffer & put_bits
shun-iwasawa 82a8f5
    mov         dword [esi+8], put_buffer  ; state->cur.put_buffer = put_buffer;
shun-iwasawa 82a8f5
    mov         dword [esi+12], put_bits   ; state->cur.put_bits = put_bits;
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
    pop         ebp
shun-iwasawa 82a8f5
    pop         edi
shun-iwasawa 82a8f5
    pop         esi
shun-iwasawa 82a8f5
;   pop         edx                     ; need not be preserved
shun-iwasawa 82a8f5
    pop         ecx
shun-iwasawa 82a8f5
    pop         ebx
shun-iwasawa 82a8f5
    mov         esp, ebp                ; esp <- aligned ebp
shun-iwasawa 82a8f5
    pop         esp                     ; esp <- original ebp
shun-iwasawa 82a8f5
    pop         ebp
shun-iwasawa 82a8f5
    ret
shun-iwasawa 82a8f5
shun-iwasawa 82a8f5
; For some reason, the OS X linker does not honor the request to align the
shun-iwasawa 82a8f5
; segment unless we do this.
shun-iwasawa 82a8f5
    align       32