media/libvpx/vp8/encoder/x86/quantize_sse4.asm

Thu, 22 Jan 2015 13:21:57 +0100

author
Michael Schloh von Bennewitz <michael@schloh.com>
date
Thu, 22 Jan 2015 13:21:57 +0100
branch
TOR_BUG_9701
changeset 15
b8a032363ba2
permissions
-rw-r--r--

Incorporate requested changes from Mozilla in review:
https://bugzilla.mozilla.org/show_bug.cgi?id=1123480#c6

     1 ;
     2 ;  Copyright (c) 2010 The WebM project authors. All Rights Reserved.
     3 ;
     4 ;  Use of this source code is governed by a BSD-style license and patent
     5 ;  grant that can be found in the LICENSE file in the root of the source
     6 ;  tree. All contributing project authors may be found in the AUTHORS
     7 ;  file in the root of the source tree.
     8 ;
    11 %include "vpx_ports/x86_abi_support.asm"
    12 %include "vp8_asm_enc_offsets.asm"
    15 ; void vp8_regular_quantize_b_sse4 | arg
    16 ;  (BLOCK  *b,                     |  0
    17 ;   BLOCKD *d)                     |  1
    19 global sym(vp8_regular_quantize_b_sse4) PRIVATE
    20 sym(vp8_regular_quantize_b_sse4):
    22 %if ABI_IS_32BIT
    23     push        rbp
    24     mov         rbp, rsp
    25     GET_GOT     rbx
    26     push        rdi
    27     push        rsi
    29     ALIGN_STACK 16, rax
    30     %define qcoeff      0 ; 32
    31     %define stack_size 32
    32     sub         rsp, stack_size
    33 %else
    34   %if LIBVPX_YASM_WIN64
    35     SAVE_XMM 8, u
    36     push        rdi
    37     push        rsi
    38   %endif
    39 %endif
    40     ; end prolog
    42 %if ABI_IS_32BIT
    43     mov         rdi, arg(0)                 ; BLOCK *b
    44     mov         rsi, arg(1)                 ; BLOCKD *d
    45 %else
    46   %if LIBVPX_YASM_WIN64
    47     mov         rdi, rcx                    ; BLOCK *b
    48     mov         rsi, rdx                    ; BLOCKD *d
    49   %else
    50     ;mov         rdi, rdi                    ; BLOCK *b
    51     ;mov         rsi, rsi                    ; BLOCKD *d
    52   %endif
    53 %endif
    55     mov         rax, [rdi + vp8_block_coeff]
    56     mov         rcx, [rdi + vp8_block_zbin]
    57     mov         rdx, [rdi + vp8_block_round]
    58     movd        xmm7, [rdi + vp8_block_zbin_extra]
    60     ; z
    61     movdqa      xmm0, [rax]
    62     movdqa      xmm1, [rax + 16]
    64     ; duplicate zbin_oq_value
    65     pshuflw     xmm7, xmm7, 0
    66     punpcklwd   xmm7, xmm7
    68     movdqa      xmm2, xmm0
    69     movdqa      xmm3, xmm1
    71     ; sz
    72     psraw       xmm0, 15
    73     psraw       xmm1, 15
    75     ; (z ^ sz)
    76     pxor        xmm2, xmm0
    77     pxor        xmm3, xmm1
    79     ; x = abs(z)
    80     psubw       xmm2, xmm0
    81     psubw       xmm3, xmm1
    83     ; zbin
    84     movdqa      xmm4, [rcx]
    85     movdqa      xmm5, [rcx + 16]
    87     ; *zbin_ptr + zbin_oq_value
    88     paddw       xmm4, xmm7
    89     paddw       xmm5, xmm7
    91     movdqa      xmm6, xmm2
    92     movdqa      xmm7, xmm3
    94     ; x - (*zbin_ptr + zbin_oq_value)
    95     psubw       xmm6, xmm4
    96     psubw       xmm7, xmm5
    98     ; round
    99     movdqa      xmm4, [rdx]
   100     movdqa      xmm5, [rdx + 16]
   102     mov         rax, [rdi + vp8_block_quant_shift]
   103     mov         rcx, [rdi + vp8_block_quant]
   104     mov         rdx, [rdi + vp8_block_zrun_zbin_boost]
   106     ; x + round
   107     paddw       xmm2, xmm4
   108     paddw       xmm3, xmm5
   110     ; quant
   111     movdqa      xmm4, [rcx]
   112     movdqa      xmm5, [rcx + 16]
   114     ; y = x * quant_ptr >> 16
   115     pmulhw      xmm4, xmm2
   116     pmulhw      xmm5, xmm3
   118     ; y += x
   119     paddw       xmm2, xmm4
   120     paddw       xmm3, xmm5
   122     pxor        xmm4, xmm4
   123 %if ABI_IS_32BIT
   124     movdqa      [rsp + qcoeff], xmm4
   125     movdqa      [rsp + qcoeff + 16], xmm4
   126 %else
   127     pxor        xmm8, xmm8
   128 %endif
   130     ; quant_shift
   131     movdqa      xmm5, [rax]
   133     ; zrun_zbin_boost
   134     mov         rax, rdx
   136 %macro ZIGZAG_LOOP 5
   137     ; x
   138     pextrw      ecx, %4, %2
   140     ; if (x >= zbin)
   141     sub         cx, WORD PTR[rdx]           ; x - zbin
   142     lea         rdx, [rdx + 2]              ; zbin_boost_ptr++
   143     jl          .rq_zigzag_loop_%1          ; x < zbin
   145     pextrw      edi, %3, %2                 ; y
   147     ; downshift by quant_shift[rc]
   148     pextrb      ecx, xmm5, %1               ; quant_shift[rc]
   149     sar         edi, cl                     ; also sets Z bit
   150     je          .rq_zigzag_loop_%1          ; !y
   151 %if ABI_IS_32BIT
   152     mov         WORD PTR[rsp + qcoeff + %1 *2], di
   153 %else
   154     pinsrw      %5, edi, %2                 ; qcoeff[rc]
   155 %endif
   156     mov         rdx, rax                    ; reset to b->zrun_zbin_boost
   157 .rq_zigzag_loop_%1:
   158 %endmacro
   159 ; in vp8_default_zig_zag1d order: see vp8/common/entropy.c
   160 ZIGZAG_LOOP  0, 0, xmm2, xmm6, xmm4
   161 ZIGZAG_LOOP  1, 1, xmm2, xmm6, xmm4
   162 ZIGZAG_LOOP  4, 4, xmm2, xmm6, xmm4
   163 ZIGZAG_LOOP  8, 0, xmm3, xmm7, xmm8
   164 ZIGZAG_LOOP  5, 5, xmm2, xmm6, xmm4
   165 ZIGZAG_LOOP  2, 2, xmm2, xmm6, xmm4
   166 ZIGZAG_LOOP  3, 3, xmm2, xmm6, xmm4
   167 ZIGZAG_LOOP  6, 6, xmm2, xmm6, xmm4
   168 ZIGZAG_LOOP  9, 1, xmm3, xmm7, xmm8
   169 ZIGZAG_LOOP 12, 4, xmm3, xmm7, xmm8
   170 ZIGZAG_LOOP 13, 5, xmm3, xmm7, xmm8
   171 ZIGZAG_LOOP 10, 2, xmm3, xmm7, xmm8
   172 ZIGZAG_LOOP  7, 7, xmm2, xmm6, xmm4
   173 ZIGZAG_LOOP 11, 3, xmm3, xmm7, xmm8
   174 ZIGZAG_LOOP 14, 6, xmm3, xmm7, xmm8
   175 ZIGZAG_LOOP 15, 7, xmm3, xmm7, xmm8
   177     mov         rcx, [rsi + vp8_blockd_dequant]
   178     mov         rdi, [rsi + vp8_blockd_dqcoeff]
   180 %if ABI_IS_32BIT
   181     movdqa      xmm4, [rsp + qcoeff]
   182     movdqa      xmm5, [rsp + qcoeff + 16]
   183 %else
   184     %define     xmm5 xmm8
   185 %endif
   187     ; y ^ sz
   188     pxor        xmm4, xmm0
   189     pxor        xmm5, xmm1
   190     ; x = (y ^ sz) - sz
   191     psubw       xmm4, xmm0
   192     psubw       xmm5, xmm1
   194     ; dequant
   195     movdqa      xmm0, [rcx]
   196     movdqa      xmm1, [rcx + 16]
   198     mov         rcx, [rsi + vp8_blockd_qcoeff]
   200     pmullw      xmm0, xmm4
   201     pmullw      xmm1, xmm5
   203     ; store qcoeff
   204     movdqa      [rcx], xmm4
   205     movdqa      [rcx + 16], xmm5
   207     ; store dqcoeff
   208     movdqa      [rdi], xmm0
   209     movdqa      [rdi + 16], xmm1
   211     mov         rcx, [rsi + vp8_blockd_eob]
   213     ; select the last value (in zig_zag order) for EOB
   214     pxor        xmm6, xmm6
   215     pcmpeqw     xmm4, xmm6
   216     pcmpeqw     xmm5, xmm6
   218     packsswb    xmm4, xmm5
   219     pshufb      xmm4, [GLOBAL(zig_zag1d)]
   220     pmovmskb    edx, xmm4
   221     xor         rdi, rdi
   222     mov         eax, -1
   223     xor         dx, ax
   224     bsr         eax, edx
   225     sub         edi, edx
   226     sar         edi, 31
   227     add         eax, 1
   228     and         eax, edi
   230     mov         BYTE PTR [rcx], al          ; store eob
   232     ; begin epilog
   233 %if ABI_IS_32BIT
   234     add         rsp, stack_size
   235     pop         rsp
   237     pop         rsi
   238     pop         rdi
   239     RESTORE_GOT
   240     pop         rbp
   241 %else
   242   %undef xmm5
   243   %if LIBVPX_YASM_WIN64
   244     pop         rsi
   245     pop         rdi
   246     RESTORE_XMM
   247   %endif
   248 %endif
   250     ret
   252 SECTION_RODATA
   253 align 16
   254 ; vp8/common/entropy.c: vp8_default_zig_zag1d
   255 zig_zag1d:
   256     db 0, 1, 4, 8, 5, 2, 3, 6, 9, 12, 13, 10, 7, 11, 14, 15

mercurial