Source code
Revision control
Copy as Markdown
Other Tools
; Copyright © 2021, VideoLAN and dav1d authors
; Copyright © 2021, Two Orioles, LLC
; All rights reserved.
;
; Redistribution and use in source and binary forms, with or without
; modification, are permitted provided that the following conditions are met:
;
; 1. Redistributions of source code must retain the above copyright notice, this
; list of conditions and the following disclaimer.
;
; 2. Redistributions in binary form must reproduce the above copyright notice,
; this list of conditions and the following disclaimer in the documentation
; and/or other materials provided with the distribution.
;
; THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
; ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
; WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
; DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR
; ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
; (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
; ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
; (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
; SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
%include "config.asm"
%include "ext/x86/x86inc.asm"
%if ARCH_X86_64
SECTION_RODATA 16
wiener_shufA: db 2, 3, 4, 5, 4, 5, 6, 7, 6, 7, 8, 9, 8, 9, 10, 11
wiener_shufB: db 6, 7, 4, 5, 8, 9, 6, 7, 10, 11, 8, 9, 12, 13, 10, 11
wiener_shufC: db 6, 7, 8, 9, 8, 9, 10, 11, 10, 11, 12, 13, 12, 13, 14, 15
wiener_shufD: db 2, 3, -1, -1, 4, 5, -1, -1, 6, 7, -1, -1, 8, 9, -1, -1
wiener_shufE: db 0, 1, 8, 9, 2, 3, 10, 11, 4, 5, 12, 13, 6, 7, 14, 15
r_ext_mask: times 72 db -1
times 8 db 0
wiener_hshift: dw 4, 4, 1, 1
wiener_vshift: dw 1024, 1024, 4096, 4096
wiener_round: dd 1049600, 1048832
pw_164_455: dw 164, 455
pw_1023: times 2 dw 1023
pw_61448: times 2 dw 61448
pd_m262128: dd -262128
pd_m34816: dd -34816
pd_m25: dd -25
pd_m9: dd -9
pd_8: dd 8
pd_2147483648: dd 2147483648
cextern sgr_x_by_x
SECTION .text
DECLARE_REG_TMP 8, 7, 9, 11, 12, 13, 14 ; ring buffer pointers
INIT_ZMM avx512icl
cglobal wiener_filter7_16bpc, 4, 15, 17, -384*12-16, dst, stride, left, lpf, \
w, h, edge, flt
%define base t4-wiener_hshift
mov fltq, r6mp
movifnidn wd, wm
movifnidn hd, hm
mov edged, r7m
mov t3d, r8m ; pixel_max
vbroadcasti128 m6, [wiener_shufA]
vpbroadcastd m12, [fltq+ 0] ; x0 x1
lea t4, [wiener_hshift]
vbroadcasti128 m7, [wiener_shufB]
add wd, wd
vpbroadcastd m13, [fltq+ 4] ; x2 x3
shr t3d, 11
vpbroadcastd m14, [fltq+16] ; y0 y1
add lpfq, wq
vpbroadcastd m15, [fltq+20] ; y2 y3
add dstq, wq
vbroadcasti128 m8, [wiener_shufC]
lea t1, [rsp+wq+16]
vbroadcasti128 m9, [wiener_shufD]
neg wq
vpbroadcastd m0, [base+wiener_hshift+t3*4]
mov r10d, 0xfe
vpbroadcastd m10, [base+wiener_round+t3*4]
kmovb k1, r10d
vpbroadcastd m11, [base+wiener_vshift+t3*4]
pmullw m12, m0 ; upshift filter coefs to make the
vpbroadcastd m16, [pd_m262128]
pmullw m13, m0 ; horizontal downshift constant
test edgeb, 4 ; LR_HAVE_TOP
jz .no_top
call .h_top
add lpfq, strideq
mov t6, t1
mov t5, t1
add t1, 384*2
call .h_top
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
mov t4, t1
add t1, 384*2
add r10, strideq
mov [rsp], r10 ; below
call .h
mov t3, t1
mov t2, t1
dec hd
jz .v1
add lpfq, strideq
add t1, 384*2
call .h
mov t2, t1
dec hd
jz .v2
add lpfq, strideq
add t1, 384*2
call .h
dec hd
jz .v3
.main:
lea t0, [t1+384*2]
.main_loop:
call .hv
dec hd
jnz .main_loop
test edgeb, 8 ; LR_HAVE_BOTTOM
jz .v3
mov lpfq, [rsp]
call .hv_bottom
add lpfq, strideq
call .hv_bottom
.v1:
call .v
RET
.no_top:
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
lea r10, [r10+strideq*2]
mov [rsp], r10
call .h
mov t6, t1
mov t5, t1
mov t4, t1
mov t3, t1
mov t2, t1
dec hd
jz .v1
add lpfq, strideq
add t1, 384*2
call .h
mov t2, t1
dec hd
jz .v2
add lpfq, strideq
add t1, 384*2
call .h
dec hd
jz .v3
lea t0, [t1+384*2]
call .hv
dec hd
jz .v3
add t0, 384*8
call .hv
dec hd
jnz .main
.v3:
call .v
.v2:
call .v
jmp .v1
.h:
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
movq xm3, [leftq]
vmovdqu64 m3{k1}, [lpfq+r10-8]
add leftq, 8
jmp .h_main
.h_extend_left:
mova m4, [lpfq+r10+0]
vpbroadcastw xm3, xm4
vmovdqu64 m3{k1}, [lpfq+r10-8]
jmp .h_main2
.h_top:
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
.h_loop:
movu m3, [lpfq+r10-8]
.h_main:
mova m4, [lpfq+r10+0]
.h_main2:
movu m5, [lpfq+r10+8]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .h_have_right
cmp r10d, -68
jl .h_have_right
push r0
lea r0, [r_ext_mask+66]
vpbroadcastw m0, [lpfq-2]
vpternlogd m3, m0, [r0+r10+ 0], 0xe4 ; c ? a : b
vpternlogd m4, m0, [r0+r10+ 8], 0xe4
vpternlogd m5, m0, [r0+r10+16], 0xe4
pop r0
.h_have_right:
pshufb m2, m3, m6
pshufb m1, m4, m7
paddw m2, m1
pshufb m3, m8
mova m0, m16
vpdpwssd m0, m2, m12
pshufb m1, m4, m9
paddw m3, m1
pshufb m1, m4, m6
vpdpwssd m0, m3, m13
pshufb m2, m5, m7
paddw m2, m1
mova m1, m16
pshufb m4, m8
vpdpwssd m1, m2, m12
pshufb m5, m9
paddw m4, m5
vpdpwssd m1, m4, m13
psrad m0, 4
psrad m1, 4
packssdw m0, m1
psraw m0, 1
mova [t1+r10], m0
add r10, 64
jl .h_loop
ret
ALIGN function_align
.hv:
add lpfq, strideq
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv_extend_left
movq xm3, [leftq]
vmovdqu64 m3{k1}, [lpfq+r10-8]
add leftq, 8
jmp .hv_main
.hv_extend_left:
mova m4, [lpfq+r10+0]
vpbroadcastw xm3, xm4
vmovdqu64 m3{k1}, [lpfq+r10-8]
jmp .hv_main2
.hv_bottom:
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv_extend_left
.hv_loop:
movu m3, [lpfq+r10-8]
.hv_main:
mova m4, [lpfq+r10+0]
.hv_main2:
movu m5, [lpfq+r10+8]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .hv_have_right
cmp r10d, -68
jl .hv_have_right
push r0
lea r0, [r_ext_mask+66]
vpbroadcastw m0, [lpfq-2]
vpternlogd m3, m0, [r0+r10+ 0], 0xe4
vpternlogd m4, m0, [r0+r10+ 8], 0xe4
vpternlogd m5, m0, [r0+r10+16], 0xe4
pop r0
.hv_have_right:
pshufb m2, m3, m6
pshufb m1, m4, m7
paddw m2, m1
pshufb m3, m8
mova m0, m16
vpdpwssd m0, m2, m12
pshufb m1, m4, m9
paddw m3, m1
pshufb m1, m4, m6
vpdpwssd m0, m3, m13
pshufb m2, m5, m7
paddw m2, m1
pshufb m4, m8
mova m1, m16
vpdpwssd m1, m2, m12
pshufb m5, m9
paddw m4, m5
vpdpwssd m1, m4, m13
mova m2, [t4+r10]
paddw m2, [t2+r10]
mova m5, [t3+r10]
psrad m0, 4
psrad m1, 4
packssdw m0, m1
mova m4, [t5+r10]
paddw m4, [t1+r10]
psraw m0, 1
paddw m3, m0, [t6+r10]
mova [t0+r10], m0
punpcklwd m1, m2, m5
mova m0, m10
vpdpwssd m0, m1, m15
punpckhwd m2, m5
mova m1, m10
vpdpwssd m1, m2, m15
punpcklwd m2, m3, m4
vpdpwssd m0, m2, m14
punpckhwd m3, m4
vpdpwssd m1, m3, m14
psrad m0, 5
psrad m1, 5
packusdw m0, m1
pmulhuw m0, m11
mova [dstq+r10], m0
add r10, 64
jl .hv_loop
mov t6, t5
mov t5, t4
mov t4, t3
mov t3, t2
mov t2, t1
mov t1, t0
mov t0, t6
add dstq, strideq
ret
.v:
mov r10, wq
.v_loop:
mova m2, [t4+r10]
paddw m2, [t2+r10]
mova m3, [t3+r10]
punpcklwd m1, m2, m3
mova m0, m10
vpdpwssd m0, m1, m15
punpckhwd m2, m3
mova m1, m10
vpdpwssd m1, m2, m15
mova m4, [t1+r10]
paddw m3, m4, [t6+r10]
paddw m4, [t5+r10]
punpcklwd m2, m3, m4
vpdpwssd m0, m2, m14
punpckhwd m3, m4
vpdpwssd m1, m3, m14
psrad m0, 5
psrad m1, 5
packusdw m0, m1
pmulhuw m0, m11
mova [dstq+r10], m0
add r10, 64
jl .v_loop
mov t6, t5
mov t5, t4
mov t4, t3
mov t3, t2
mov t2, t1
add dstq, strideq
ret
cglobal wiener_filter5_16bpc, 4, 14, 15, 384*8+16, dst, stride, left, lpf, \
w, h, edge, flt
%define base r13-r_ext_mask-70
mov fltq, r6mp
movifnidn wd, wm
movifnidn hd, hm
mov edged, r7m
mov t3d, r8m ; pixel_max
vbroadcasti128 m5, [wiener_shufE]
vpbroadcastw m11, [fltq+ 2] ; x1
vbroadcasti128 m6, [wiener_shufB]
lea r13, [r_ext_mask+70]
vbroadcasti128 m7, [wiener_shufD]
add wd, wd
vpbroadcastd m12, [fltq+ 4] ; x2 x3
shr t3d, 11
vpbroadcastd m8, [pd_m262128] ; (1 << 4) - (1 << 18)
add lpfq, wq
vpbroadcastw m13, [fltq+18] ; y1
add dstq, wq
vpbroadcastd m14, [fltq+20] ; y2 y3
lea t1, [rsp+wq+16]
vpbroadcastd m0, [base+wiener_hshift+t3*4]
neg wq
vpbroadcastd m9, [base+wiener_round+t3*4]
mov r10d, 0xfffe
vpbroadcastd m10, [base+wiener_vshift+t3*4]
kmovw k1, r10d
pmullw m11, m0
pmullw m12, m0
test edgeb, 4 ; LR_HAVE_TOP
jz .no_top
call .h_top
add lpfq, strideq
mov t4, t1
add t1, 384*2
call .h_top
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
mov t3, t1
add t1, 384*2
add r10, strideq
mov [rsp], r10 ; below
call .h
mov t2, t1
dec hd
jz .v1
add lpfq, strideq
add t1, 384*2
call .h
dec hd
jz .v2
.main:
mov t0, t4
.main_loop:
call .hv
dec hd
jnz .main_loop
test edgeb, 8 ; LR_HAVE_BOTTOM
jz .v2
mov lpfq, [rsp]
call .hv_bottom
add lpfq, strideq
call .hv_bottom
.end:
RET
.no_top:
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
lea r10, [r10+strideq*2]
mov [rsp], r10
call .h
mov t4, t1
mov t3, t1
mov t2, t1
dec hd
jz .v1
add lpfq, strideq
add t1, 384*2
call .h
dec hd
jz .v2
lea t0, [t1+384*2]
call .hv
dec hd
jz .v2
add t0, 384*6
call .hv
dec hd
jnz .main
.v2:
call .v
mov t4, t3
mov t3, t2
mov t2, t1
add dstq, strideq
.v1:
call .v
jmp .end
.h:
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
movd xm3, [leftq+4]
vmovdqu32 m3{k1}, [lpfq+r10-4]
add leftq, 8
jmp .h_main
.h_extend_left:
vpbroadcastw xm3, [lpfq+r10]
vmovdqu32 m3{k1}, [lpfq+r10-4]
jmp .h_main
.h_top:
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
.h_loop:
movu m3, [lpfq+r10-4]
.h_main:
movu m4, [lpfq+r10+4]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .h_have_right
cmp r10d, -66
jl .h_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m3, m0, [r13+r10+0], 0xe4 ; c ? a : b
vpternlogd m4, m0, [r13+r10+8], 0xe4
.h_have_right:
pshufb m1, m3, m5
mova m0, m8
vpdpwssd m0, m1, m11
pshufb m2, m4, m5
mova m1, m8
vpdpwssd m1, m2, m11
pshufb m2, m3, m6
pshufb m3, m7
paddw m2, m3
pshufb m3, m4, m6
vpdpwssd m0, m2, m12
pshufb m4, m7
paddw m3, m4
vpdpwssd m1, m3, m12
psrad m0, 4
psrad m1, 4
packssdw m0, m1
psraw m0, 1
mova [t1+r10], m0
add r10, 64
jl .h_loop
ret
ALIGN function_align
.hv:
add lpfq, strideq
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv_extend_left
movd xm3, [leftq+4]
vmovdqu32 m3{k1}, [lpfq+r10-4]
add leftq, 8
jmp .hv_main
.hv_extend_left:
vpbroadcastw xm3, [lpfq+r10]
vmovdqu32 m3{k1}, [lpfq+r10-4]
jmp .hv_main
.hv_bottom:
mov r10, wq
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv_extend_left
.hv_loop:
movu m3, [lpfq+r10-4]
.hv_main:
movu m4, [lpfq+r10+4]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .hv_have_right
cmp r10d, -66
jl .hv_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m3, m0, [r13+r10+0], 0xe4
vpternlogd m4, m0, [r13+r10+8], 0xe4
.hv_have_right:
pshufb m1, m3, m5
mova m0, m8
vpdpwssd m0, m1, m11
pshufb m2, m4, m5
mova m1, m8
vpdpwssd m1, m2, m11
pshufb m2, m3, m6
pshufb m3, m7
paddw m2, m3
pshufb m3, m4, m6
vpdpwssd m0, m2, m12
pshufb m4, m7
paddw m4, m3
vpdpwssd m1, m4, m12
mova m2, [t3+r10]
paddw m2, [t1+r10]
mova m3, [t2+r10]
punpcklwd m4, m2, m3
punpckhwd m2, m3
mova m3, m9
vpdpwssd m3, m2, m14
mova m2, m9
vpdpwssd m2, m4, m14
mova m4, [t4+r10]
psrad m0, 4
psrad m1, 4
packssdw m0, m1
psraw m0, 1
mova [t0+r10], m0
punpcklwd m1, m0, m4
vpdpwssd m2, m1, m13
punpckhwd m0, m4
vpdpwssd m3, m0, m13
psrad m2, 5
psrad m3, 5
packusdw m2, m3
pmulhuw m2, m10
mova [dstq+r10], m2
add r10, 64
jl .hv_loop
mov t4, t3
mov t3, t2
mov t2, t1
mov t1, t0
mov t0, t4
add dstq, strideq
ret
.v:
mov r10, wq
.v_loop:
mova m0, [t1+r10]
paddw m2, m0, [t3+r10]
mova m1, [t2+r10]
mova m4, [t4+r10]
punpckhwd m3, m2, m1
pmaddwd m3, m14
punpcklwd m2, m1
pmaddwd m2, m14
punpckhwd m1, m0, m4
pmaddwd m1, m13
punpcklwd m0, m4
pmaddwd m0, m13
paddd m3, m9
paddd m2, m9
paddd m1, m3
paddd m0, m2
psrad m1, 5
psrad m0, 5
packusdw m0, m1
pmulhuw m0, m10
mova [dstq+r10], m0
add r10, 64
jl .v_loop
ret
cglobal sgr_filter_5x5_16bpc, 4, 14, 22, 416*24+8, dst, stride, left, lpf, \
w, h, edge, params
%define base r13-r_ext_mask-72
movifnidn wd, wm
mov paramsq, r6mp
lea r13, [r_ext_mask+72]
mov edged, r7m
movifnidn hd, hm
pxor m6, m6
vpbroadcastw m7, [paramsq+8] ; w0
add wd, wd
vpbroadcastd m8, [base+pd_8]
add lpfq, wq
vpbroadcastd m9, [base+pd_m25]
add dstq, wq
vpsubd m10, m6, [paramsq+0] {1to16} ; -s0
lea t3, [rsp+wq*2+416*12+8]
vpbroadcastd m11, [base+pw_164_455]
lea t4, [rsp+wq+416*20+8]
vpbroadcastd m12, [base+pw_61448] ; (15 << 12) + (1 << 3)
lea t1, [rsp+wq+12]
vpbroadcastd m13, [base+pd_m34816] ; -((1 << 11) + (1 << 15))
neg wq
vpbroadcastd m14, [base+pw_1023]
psllw m7, 4
mova m18, [sgr_x_by_x+64*0]
mov r10d, 0xfffffff8
mova m19, [sgr_x_by_x+64*1]
kmovd k1, r10d
mova m20, [sgr_x_by_x+64*2]
mov r10, 0x3333333333333333
mova m21, [sgr_x_by_x+64*3]
kmovq k2, r10
test edgeb, 4 ; LR_HAVE_TOP
jz .no_top
call .h_top
add lpfq, strideq
mov t2, t1
call .top_fixup
add t1, 416*6
call .h_top
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
add r10, strideq
mov [rsp], r10 ; below
mov t0, t2
dec hd
jz .height1
or edged, 16
call .h
.main:
add lpfq, strideq
call .hv
call .prep_n
sub hd, 2
jl .extend_bottom
.main_loop:
add lpfq, strideq
test hd, hd
jz .odd_height
call .h
add lpfq, strideq
call .hv
call .n0
call .n1
sub hd, 2
jge .main_loop
test edgeb, 8 ; LR_HAVE_BOTTOM
jz .extend_bottom
mov lpfq, [rsp]
call .h_top
add lpfq, strideq
call .hv_bottom
.end:
call .n0
call .n1
.end2:
RET
.height1:
call .hv
call .prep_n
jmp .odd_height_end
.odd_height:
call .hv
call .n0
call .n1
.odd_height_end:
call .v
call .n0
jmp .end2
.extend_bottom:
call .v
jmp .end
.no_top:
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
lea r10, [r10+strideq*2]
mov [rsp], r10
call .h
lea t2, [t1+416*6]
call .top_fixup
dec hd
jz .no_top_height1
or edged, 16
mov t0, t1
mov t1, t2
jmp .main
.no_top_height1:
call .v
call .prep_n
jmp .odd_height_end
.h: ; horizontal boxsum
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
movq xm16, [leftq+2]
vmovdqu16 m16{k1}, [lpfq+wq-6]
add leftq, 8
jmp .h_main
.h_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-6]
jmp .h_main
.h_top:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
.h_loop:
movu m16, [lpfq+r10- 2]
.h_main:
movu m17, [lpfq+r10+14]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .h_have_right
cmp r10d, -68
jl .h_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4 ; c ? a : b
vpternlogd m17, m0, [r13+r10+16], 0xe4
.h_have_right:
palignr m2, m17, m16, 2
paddw m0, m16, m2
palignr m3, m17, m16, 6
paddw m0, m3
punpcklwd m1, m2, m3
pmaddwd m1, m1
punpckhwd m2, m3
pmaddwd m2, m2
shufpd m17, m16, m17, 0x55
paddw m0, m17
punpcklwd m3, m16, m17
vpdpwssd m1, m3, m3
punpckhwd m3, m16, m17
vpdpwssd m2, m3, m3
shufps m16, m17, q2121
paddw m0, m16 ; sum
test edgeb, 16 ; y > 0
jz .h_loop_end
paddw m0, [t1+r10+416*0]
paddd m1, [t1+r10+416*2]
paddd m2, [t1+r10+416*4]
.h_loop_end:
punpcklwd m17, m16, m6
vpdpwssd m1, m17, m17 ; sumsq
punpckhwd m16, m6
vpdpwssd m2, m16, m16
mova [t1+r10+416*0], m0
mova [t1+r10+416*2], m1
mova [t1+r10+416*4], m2
add r10, 64
jl .h_loop
ret
.top_fixup:
lea r10, [wq-4]
.top_fixup_loop: ; the sums of the first row needs to be doubled
mova m0, [t1+r10+416*0]
mova m1, [t1+r10+416*2]
mova m2, [t1+r10+416*4]
paddw m0, m0
paddd m1, m1
paddd m2, m2
mova [t2+r10+416*0], m0
mova [t2+r10+416*2], m1
mova [t2+r10+416*4], m2
add r10, 64
jl .top_fixup_loop
ret
ALIGN function_align
.hv: ; horizontal boxsum + vertical boxsum + ab
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv_extend_left
movq xm16, [leftq+2]
vmovdqu16 m16{k1}, [lpfq+wq-6]
add leftq, 8
jmp .hv_main
.hv_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-6]
jmp .hv_main
.hv_bottom:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv_extend_left
.hv_loop:
movu m16, [lpfq+r10- 2]
.hv_main:
movu m17, [lpfq+r10+14]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .hv_have_right
cmp r10d, -68
jl .hv_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4
vpternlogd m17, m0, [r13+r10+16], 0xe4
.hv_have_right:
palignr m3, m17, m16, 2
paddw m0, m16, m3
palignr m1, m17, m16, 6
paddw m0, m1
punpcklwd m2, m3, m1
pmaddwd m2, m2
punpckhwd m3, m1
pmaddwd m3, m3
shufpd m17, m16, m17, 0x55
paddw m0, m17
punpcklwd m1, m16, m17
vpdpwssd m2, m1, m1
punpckhwd m1, m16, m17
vpdpwssd m3, m1, m1
shufps m16, m17, q2121
paddw m0, m16 ; h sum
punpcklwd m17, m16, m6
vpdpwssd m2, m17, m17 ; h sumsq
punpckhwd m16, m6
vpdpwssd m3, m16, m16
paddw m1, m0, [t1+r10+416*0]
paddd m16, m2, [t1+r10+416*2]
paddd m17, m3, [t1+r10+416*4]
test hd, hd
jz .hv_last_row
.hv_main2:
paddw m1, [t2+r10+416*0] ; hv sum
paddd m16, [t2+r10+416*2] ; hv sumsq
paddd m17, [t2+r10+416*4]
mova [t0+r10+416*0], m0
mova [t0+r10+416*2], m2
mova [t0+r10+416*4], m3
psrlw m3, m1, 1
paddd m16, m8
pavgw m3, m6 ; (b + 2) >> 2
paddd m17, m8
psrld m16, 4 ; (a + 8) >> 4
psrld m17, 4
pmulld m16, m9 ; -a * 25
pmulld m17, m9
punpcklwd m2, m3, m6
vpdpwssd m16, m2, m2 ; -p
punpckhwd m3, m6
vpdpwssd m17, m3, m3
punpcklwd m0, m1, m6 ; b
punpckhwd m1, m6
pmulld m16, m10 ; p * s
pmulld m17, m10
pmaddwd m0, m11 ; b * 164
pmaddwd m1, m11
vpalignr m17{k2}, m16, m16, 2
mova m16, m20
pmaxsw m17, m6
paddusw m17, m12
psraw m17, 4 ; min(z, 255) - 256
vpermt2b m16, m17, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m17
vpermi2b m17, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m17{k3}, m16 ; x
pandn m16, m13, m17
psrld m17, 16
pmulld m0, m16
pmulld m1, m17
packssdw m16, m17
psubd m0, m13 ; x * b * 164 + (1 << 11) + (1 << 15)
psubd m1, m13
mova [t4+r10+4], m16
psrld m16, m0, 12 ; b
psrld m17, m1, 12
mova [t3+r10*2+ 8], xm16
mova [t3+r10*2+ 24], xm17
vextracti128 [t3+r10*2+ 40], ym16, 1
vextracti128 [t3+r10*2+ 56], ym17, 1
vextracti32x4 [t3+r10*2+ 72], m16, 2
vextracti32x4 [t3+r10*2+ 88], m17, 2
vextracti32x4 [t3+r10*2+104], m16, 3
vextracti32x4 [t3+r10*2+120], m17, 3
add r10, 64
jl .hv_loop
mov t2, t1
mov t1, t0
mov t0, t2
ret
.hv_last_row: ; esoteric edge case for odd heights
mova [t1+r10+416*0], m1
paddw m1, m0
mova [t1+r10+416*2], m16
paddd m16, m2
mova [t1+r10+416*4], m17
paddd m17, m3
jmp .hv_main2
.v: ; vertical boxsum + ab
lea r10, [wq-4]
.v_loop:
mova m2, [t1+r10+416*2]
mova m3, [t1+r10+416*4]
mova m0, [t1+r10+416*0]
paddd m16, m2, [t2+r10+416*2]
paddd m17, m3, [t2+r10+416*4]
paddw m1, m0, [t2+r10+416*0]
paddd m2, m2
paddd m3, m3
paddd m16, m2 ; hv sumsq
paddd m17, m3
paddd m16, m8
paddd m17, m8
psrld m16, 4 ; (a + 8) >> 4
psrld m17, 4
pmulld m16, m9 ; -a * 25
pmulld m17, m9
paddw m0, m0
paddw m1, m0 ; hv sum
psrlw m3, m1, 1
pavgw m3, m6 ; (b + 2) >> 2
punpcklwd m2, m3, m6
vpdpwssd m16, m2, m2 ; -p
punpckhwd m3, m6
vpdpwssd m17, m3, m3
punpcklwd m0, m1, m6 ; b
punpckhwd m1, m6
pmulld m16, m10 ; p * s
pmulld m17, m10
pmaddwd m0, m11 ; b * 164
pmaddwd m1, m11
vpalignr m17{k2}, m16, m16, 2
mova m16, m20
pmaxsw m17, m6
paddusw m17, m12
psraw m17, 4 ; min(z, 255) - 256
vpermt2b m16, m17, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m17
vpermi2b m17, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m17{k3}, m16 ; x
pandn m16, m13, m17
psrld m17, 16
pmulld m0, m16
pmulld m1, m17
packssdw m16, m17
psubd m0, m13 ; x * b * 164 + (1 << 11) + (1 << 15)
psubd m1, m13
mova [t4+r10+4], m16
psrld m16, m0, 12 ; b
psrld m17, m1, 12
mova [t3+r10*2+ 8], xm16
mova [t3+r10*2+ 24], xm17
vextracti128 [t3+r10*2+ 40], ym16, 1
vextracti128 [t3+r10*2+ 56], ym17, 1
vextracti32x4 [t3+r10*2+ 72], m16, 2
vextracti32x4 [t3+r10*2+ 88], m17, 2
vextracti32x4 [t3+r10*2+104], m16, 3
vextracti32x4 [t3+r10*2+120], m17, 3
add r10, 64
jl .v_loop
ret
.prep_n: ; initial neighbor setup
mov r10, wq
.prep_n_loop:
movu m0, [t4+r10*1+ 2]
movu m1, [t3+r10*2+ 4]
movu m2, [t3+r10*2+68]
paddw m3, m0, [t4+r10*1+ 0]
paddd m16, m1, [t3+r10*2+ 0]
paddd m17, m2, [t3+r10*2+64]
paddw m3, [t4+r10*1+ 4]
paddd m16, [t3+r10*2+ 8]
paddd m17, [t3+r10*2+72]
paddw m0, m3
psllw m3, 2
paddd m1, m16
pslld m16, 2
paddd m2, m17
pslld m17, 2
paddw m0, m3 ; a 565
paddd m1, m16 ; b 565
paddd m2, m17
mova [t4+r10*1+416*2+ 0], m0
mova [t3+r10*2+416*4+ 0], m1
mova [t3+r10*2+416*4+64], m2
add r10, 64
jl .prep_n_loop
ret
ALIGN function_align
.n0: ; neighbor + output (even rows)
mov r10, wq
.n0_loop:
movu m0, [t4+r10*1+ 2]
movu m1, [t3+r10*2+ 4]
movu m2, [t3+r10*2+68]
paddw m3, m0, [t4+r10*1+ 0]
paddd m16, m1, [t3+r10*2+ 0]
paddd m17, m2, [t3+r10*2+64]
paddw m3, [t4+r10*1+ 4]
paddd m16, [t3+r10*2+ 8]
paddd m17, [t3+r10*2+72]
paddw m0, m3
psllw m3, 2
paddd m1, m16
pslld m16, 2
paddd m2, m17
pslld m17, 2
paddw m0, m3 ; a 565
paddd m1, m16 ; b 565
paddd m2, m17
paddw m3, m0, [t4+r10*1+416*2+ 0]
paddd m16, m1, [t3+r10*2+416*4+ 0]
paddd m17, m2, [t3+r10*2+416*4+64]
mova [t4+r10*1+416*2+ 0], m0
mova [t3+r10*2+416*4+ 0], m1
mova [t3+r10*2+416*4+64], m2
mova m0, [dstq+r10]
punpcklwd m1, m0, m6 ; src
punpcklwd m2, m3, m6 ; a
pmaddwd m2, m1 ; a * src
punpckhwd m1, m0, m6
punpckhwd m3, m6
pmaddwd m3, m1
vshufi32x4 m1, m16, m17, q2020
vshufi32x4 m16, m17, q3131
psubd m1, m2 ; b - a * src + (1 << 8)
psubd m16, m3
psrad m1, 9
psrad m16, 9
packssdw m1, m16
pmulhrsw m1, m7
paddw m0, m1
pmaxsw m0, m6
pminsw m0, m14
mova [dstq+r10], m0
add r10, 64
jl .n0_loop
add dstq, strideq
ret
ALIGN function_align
.n1: ; neighbor + output (odd rows)
mov r10, wq
.n1_loop:
mova m0, [dstq+r10]
mova m3, [t4+r10*1+416*2+ 0]
mova m16, [t3+r10*2+416*4+ 0]
mova m17, [t3+r10*2+416*4+64]
punpcklwd m1, m0, m6 ; src
punpcklwd m2, m3, m6 ; a
pmaddwd m2, m1
punpckhwd m1, m0, m6
punpckhwd m3, m6
pmaddwd m3, m1
vshufi32x4 m1, m16, m17, q2020
vshufi32x4 m16, m17, q3131
psubd m1, m2 ; b - a * src + (1 << 7)
psubd m16, m3
psrad m1, 8
psrad m16, 8
packssdw m1, m16
pmulhrsw m1, m7
paddw m0, m1
pmaxsw m0, m6
pminsw m0, m14
mova [dstq+r10], m0
add r10, 64
jl .n1_loop
add dstq, strideq
ret
cglobal sgr_filter_3x3_16bpc, 4, 14, 22, 416*42+8, dst, stride, left, lpf, \
w, h, edge, params
movifnidn wd, wm
mov paramsq, r6mp
lea r13, [r_ext_mask+72]
mov edged, r7m
movifnidn hd, hm
pxor m6, m6
vpbroadcastw m7, [paramsq+10] ; w1
add wd, wd
vpbroadcastd m8, [base+pd_8]
add lpfq, wq
vpbroadcastd m9, [base+pd_m9]
add dstq, wq
vpsubd m10, m6, [paramsq+4] {1to16} ; -s1
lea t3, [rsp+wq*2+416*12+8]
vpbroadcastd m11, [base+pw_164_455]
lea t4, [rsp+wq+416*32+8]
vpbroadcastd m12, [base+pw_61448]
lea t1, [rsp+wq+12]
vpbroadcastd m13, [base+pd_m34816]
neg wq
vpbroadcastd m14, [base+pw_1023]
psllw m7, 4
mova m18, [sgr_x_by_x+64*0]
mov r10d, 0xfffffffc
mova m19, [sgr_x_by_x+64*1]
kmovd k1, r10d
mova m20, [sgr_x_by_x+64*2]
mov r10, 0x3333333333333333
mova m21, [sgr_x_by_x+64*3]
kmovq k2, r10
test edgeb, 4 ; LR_HAVE_TOP
jz .no_top
call .h_top
add lpfq, strideq
mov t2, t1
add t1, 416*6
call .h_top
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
add r10, strideq
mov [rsp], r10 ; below
call .hv0
.main:
dec hd
jz .height1
add lpfq, strideq
call .hv1
call .prep_n
sub hd, 2
jl .extend_bottom
.main_loop:
add lpfq, strideq
call .hv0
test hd, hd
jz .odd_height
add lpfq, strideq
call .hv1
call .n0
call .n1
sub hd, 2
jge .main_loop
test edgeb, 8 ; LR_HAVE_BOTTOM
jz .extend_bottom
mov lpfq, [rsp]
call .hv0_bottom
add lpfq, strideq
call .hv1_bottom
.end:
call .n0
call .n1
.end2:
RET
.height1:
call .v1
call .prep_n
jmp .odd_height_end
.odd_height:
call .v1
call .n0
call .n1
.odd_height_end:
call .v0
call .v1
call .n0
jmp .end2
.extend_bottom:
call .v0
call .v1
jmp .end
.no_top:
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
lea r10, [r10+strideq*2]
mov [rsp], r10
call .h
lea r10, [wq-4]
lea t2, [t1+416*6]
.top_fixup_loop:
mova m0, [t1+r10+416*0]
mova m1, [t1+r10+416*2]
mova m2, [t1+r10+416*4]
mova [t2+r10+416*0], m0
mova [t2+r10+416*2], m1
mova [t2+r10+416*4], m2
add r10, 64
jl .top_fixup_loop
call .v0
jmp .main
.h: ; horizontal boxsum
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
movd xm16, [leftq+4]
vmovdqu16 m16{k1}, [lpfq+wq-4]
add leftq, 8
jmp .h_main
.h_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-4]
jmp .h_main
.h_top:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
.h_loop:
movu m16, [lpfq+r10+ 0]
.h_main:
movu m17, [lpfq+r10+16]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .h_have_right
cmp r10d, -66
jl .h_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4
vpternlogd m17, m0, [r13+r10+16], 0xe4
.h_have_right:
palignr m0, m17, m16, 2
paddw m1, m16, m0
punpcklwd m2, m16, m0
pmaddwd m2, m2
punpckhwd m3, m16, m0
pmaddwd m3, m3
palignr m17, m16, 4
paddw m1, m17 ; sum
punpcklwd m16, m17, m6
vpdpwssd m2, m16, m16 ; sumsq
punpckhwd m17, m6
vpdpwssd m3, m17, m17
mova [t1+r10+416*0], m1
mova [t1+r10+416*2], m2
mova [t1+r10+416*4], m3
add r10, 64
jl .h_loop
ret
ALIGN function_align
.hv0: ; horizontal boxsum + vertical boxsum + ab (even rows)
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv0_extend_left
movd xm16, [leftq+4]
vmovdqu16 m16{k1}, [lpfq+wq-4]
add leftq, 8
jmp .hv0_main
.hv0_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-4]
jmp .hv0_main
.hv0_bottom:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv0_extend_left
.hv0_loop:
movu m16, [lpfq+r10+ 0]
.hv0_main:
movu m17, [lpfq+r10+16]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .hv0_have_right
cmp r10d, -66
jl .hv0_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4
vpternlogd m17, m0, [r13+r10+16], 0xe4
.hv0_have_right:
palignr m0, m17, m16, 2
paddw m1, m16, m0
punpcklwd m2, m16, m0
pmaddwd m2, m2
punpckhwd m3, m16, m0
pmaddwd m3, m3
palignr m17, m16, 4
paddw m1, m17 ; sum
punpcklwd m16, m17, m6
vpdpwssd m2, m16, m16 ; sumsq
punpckhwd m17, m6
vpdpwssd m3, m17, m17
paddw m0, m1, [t1+r10+416*0]
paddd m16, m2, [t1+r10+416*2]
paddd m17, m3, [t1+r10+416*4]
mova [t1+r10+416*0], m1
mova [t1+r10+416*2], m2
mova [t1+r10+416*4], m3
paddw m1, m0, [t2+r10+416*0]
paddd m2, m16, [t2+r10+416*2]
paddd m3, m17, [t2+r10+416*4]
mova [t2+r10+416*0], m0
mova [t2+r10+416*2], m16
mova [t2+r10+416*4], m17
paddd m2, m8
paddd m3, m8
psrld m2, 4 ; (a + 8) >> 4
psrld m3, 4
pmulld m2, m9 ; -((a + 8) >> 4) * 9
pmulld m3, m9
psrlw m17, m1, 1
pavgw m17, m6 ; (b + 2) >> 2
punpcklwd m16, m17, m6
vpdpwssd m2, m16, m16 ; -p
punpckhwd m17, m6
vpdpwssd m3, m17, m17
punpcklwd m16, m6, m1 ; b
punpckhwd m17, m6, m1
pminsd m2, m6
pminsd m3, m6
pmulld m2, m10 ; p * s
pmulld m3, m10
pmaddwd m16, m11 ; b * 455
pmaddwd m17, m11
vpalignr m3{k2}, m2, m2, 2
mova m2, m20
paddusw m3, m12
psraw m3, 4 ; min(z, 255) - 256
vpermt2b m2, m3, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m3
vpermi2b m3, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m3{k3}, m2 ; x
pandn m2, m13, m3
psrld m3, 16
pmulld m16, m2
pmulld m17, m3
packssdw m2, m3
psubd m16, m13 ; x * b * 455 + (1 << 11) + (1 << 15)
psubd m17, m13
mova [t4+r10*1+416*0+4], m2
psrld m16, 12
psrld m17, 12
mova [t3+r10*2+416*0+ 8], xm16
mova [t3+r10*2+416*0+ 24], xm17
vextracti128 [t3+r10*2+416*0+ 40], ym16, 1
vextracti128 [t3+r10*2+416*0+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*0+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*0+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*0+104], m16, 3
vextracti32x4 [t3+r10*2+416*0+120], m17, 3
add r10, 64
jl .hv0_loop
ret
ALIGN function_align
.hv1: ; horizontal boxsums + vertical boxsums + ab (odd rows)
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv1_extend_left
movd xm16, [leftq+4]
vmovdqu16 m16{k1}, [lpfq+wq-4]
add leftq, 8
jmp .hv1_main
.hv1_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-4]
jmp .hv1_main
.hv1_bottom:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv1_extend_left
.hv1_loop:
movu m16, [lpfq+r10+ 0]
.hv1_main:
movu m17, [lpfq+r10+16]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .hv1_have_right
cmp r10d, -66
jl .hv1_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4
vpternlogd m17, m0, [r13+r10+16], 0xe4
.hv1_have_right:
palignr m1, m17, m16, 2
paddw m0, m16, m1
punpcklwd m2, m16, m1
pmaddwd m2, m2
punpckhwd m3, m16, m1
pmaddwd m3, m3
palignr m17, m16, 4
paddw m0, m17 ; h sum
punpcklwd m1, m17, m6
vpdpwssd m2, m1, m1 ; h sumsq
punpckhwd m17, m6
vpdpwssd m3, m17, m17
paddw m1, m0, [t2+r10+416*0]
paddd m16, m2, [t2+r10+416*2]
paddd m17, m3, [t2+r10+416*4]
mova [t2+r10+416*0], m0
mova [t2+r10+416*2], m2
mova [t2+r10+416*4], m3
paddd m16, m8
paddd m17, m8
psrld m16, 4 ; (a + 8) >> 4
psrld m17, 4
pmulld m16, m9 ; -((a + 8) >> 4) * 9
pmulld m17, m9
psrlw m3, m1, 1
pavgw m3, m6 ; (b + 2) >> 2
punpcklwd m2, m3, m6
vpdpwssd m16, m2, m2 ; -p
punpckhwd m3, m6
vpdpwssd m17, m3, m3
punpcklwd m0, m6, m1 ; b
punpckhwd m1, m6, m1
pminsd m16, m6
pminsd m17, m6
pmulld m16, m10 ; p * s
pmulld m17, m10
pmaddwd m0, m11 ; b * 455
pmaddwd m1, m11
vpalignr m17{k2}, m16, m16, 2
mova m16, m20
paddusw m17, m12
psraw m17, 4 ; min(z, 255) - 256
vpermt2b m16, m17, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m17
vpermi2b m17, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m17{k3}, m16 ; x
pandn m16, m13, m17
psrld m17, 16
pmulld m0, m16
pmulld m1, m17
packssdw m16, m17
psubd m0, m13 ; x * b * 455 + (1 << 11) + (1 << 15)
psubd m1, m13
mova [t4+r10*1+416*2+4], m16
psrld m16, m0, 12
psrld m17, m1, 12
mova [t3+r10*2+416*4+ 8], xm16
mova [t3+r10*2+416*4+ 24], xm17
vextracti128 [t3+r10*2+416*4+ 40], ym16, 1
vextracti128 [t3+r10*2+416*4+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*4+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*4+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*4+104], m16, 3
vextracti32x4 [t3+r10*2+416*4+120], m17, 3
add r10, 64
jl .hv1_loop
mov r10, t2
mov t2, t1
mov t1, r10
ret
.v0: ; vertical boxsums + ab (even rows)
lea r10, [wq-4]
.v0_loop:
mova m0, [t1+r10+416*0]
mova m16, [t1+r10+416*2]
mova m17, [t1+r10+416*4]
paddw m0, m0
paddd m16, m16
paddd m17, m17
paddw m1, m0, [t2+r10+416*0]
paddd m2, m16, [t2+r10+416*2]
paddd m3, m17, [t2+r10+416*4]
mova [t2+r10+416*0], m0
mova [t2+r10+416*2], m16
mova [t2+r10+416*4], m17
paddd m2, m8
paddd m3, m8
psrld m2, 4 ; (a + 8) >> 4
psrld m3, 4
pmulld m2, m9 ; -((a + 8) >> 4) * 9
pmulld m3, m9
psrlw m17, m1, 1
pavgw m17, m6 ; (b + 2) >> 2
punpcklwd m16, m17, m6
vpdpwssd m2, m16, m16 ; -p
punpckhwd m17, m6
vpdpwssd m3, m17, m17
punpcklwd m16, m6, m1 ; b
punpckhwd m17, m6, m1
pminsd m2, m6
pminsd m3, m6
pmulld m2, m10 ; p * s
pmulld m3, m10
pmaddwd m16, m11 ; b * 455
pmaddwd m17, m11
vpalignr m3{k2}, m2, m2, 2
mova m2, m20
paddusw m3, m12
psraw m3, 4 ; min(z, 255) - 256
vpermt2b m2, m3, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m3
vpermi2b m3, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m3{k3}, m2 ; x
pandn m2, m13, m3
psrld m3, 16
pmulld m16, m2
pmulld m17, m3
packssdw m2, m3
psubd m16, m13 ; x * b * 455 + (1 << 11) + (1 << 15)
psubd m17, m13
mova [t4+r10*1+416*0+4], m2
psrld m16, 12
psrld m17, 12
mova [t3+r10*2+416*0+ 8], xm16
mova [t3+r10*2+416*0+ 24], xm17
vextracti128 [t3+r10*2+416*0+ 40], ym16, 1
vextracti128 [t3+r10*2+416*0+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*0+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*0+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*0+104], m16, 3
vextracti32x4 [t3+r10*2+416*0+120], m17, 3
add r10, 64
jl .v0_loop
ret
.v1: ; vertical boxsums + ab (odd rows)
lea r10, [wq-4]
.v1_loop:
mova m0, [t1+r10+416*0]
mova m16, [t1+r10+416*2]
mova m17, [t1+r10+416*4]
paddw m1, m0, [t2+r10+416*0]
paddd m2, m16, [t2+r10+416*2]
paddd m3, m17, [t2+r10+416*4]
mova [t2+r10+416*0], m0
mova [t2+r10+416*2], m16
mova [t2+r10+416*4], m17
paddd m2, m8
paddd m3, m8
psrld m2, 4 ; (a + 8) >> 4
psrld m3, 4
pmulld m2, m9 ; -((a + 8) >> 4) * 9
pmulld m3, m9
psrlw m17, m1, 1
pavgw m17, m6 ; (b + 2) >> 2
punpcklwd m16, m17, m6
vpdpwssd m2, m16, m16 ; -p
punpckhwd m17, m6
vpdpwssd m3, m17, m17
punpcklwd m16, m6, m1 ; b
punpckhwd m17, m6, m1
pminsd m2, m6
pminsd m3, m6
pmulld m2, m10 ; p * s
pmulld m3, m10
pmaddwd m16, m11 ; b * 455
pmaddwd m17, m11
vpalignr m3{k2}, m2, m2, 2
mova m2, m20
paddusw m3, m12
psraw m3, 4 ; min(z, 255) - 256
vpermt2b m2, m3, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m3
vpermi2b m3, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m3{k3}, m2 ; x
pandn m2, m13, m3
psrld m3, 16
pmulld m16, m2
pmulld m17, m3
packssdw m2, m3
psubd m16, m13 ; x * b * 455 + (1 << 11) + (1 << 15)
psubd m17, m13
mova [t4+r10*1+416*2+4], m2
psrld m16, 12
psrld m17, 12
mova [t3+r10*2+416*4+ 8], xm16
mova [t3+r10*2+416*4+ 24], xm17
vextracti128 [t3+r10*2+416*4+ 40], ym16, 1
vextracti128 [t3+r10*2+416*4+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*4+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*4+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*4+104], m16, 3
vextracti32x4 [t3+r10*2+416*4+120], m17, 3
add r10, 64
jl .v1_loop
mov r10, t2
mov t2, t1
mov t1, r10
ret
.prep_n: ; initial neighbor setup
mov r10, wq
.prep_n_loop:
mova ym16, [t4+r10*1+416*0+0]
paddw ym16, [t4+r10*1+416*0+4]
paddw ym17, ym16, [t4+r10*1+416*0+2]
mova m0, [t3+r10*2+416*0+0]
paddd m0, [t3+r10*2+416*0+8]
paddd m1, m0, [t3+r10*2+416*0+4]
psllw ym17, 2 ; a[-1] 444
pslld m1, 2 ; b[-1] 444
psubw ym17, ym16 ; a[-1] 343
psubd m1, m0 ; b[-1] 343
vmovdqa32 [t4+r10*1+416* 4], ym17
vmovdqa32 [t3+r10*2+416* 8], m1
mova ym16, [t4+r10*1+416*2+0]
paddw ym16, [t4+r10*1+416*2+4]
paddw ym17, ym16, [t4+r10*1+416*2+2]
mova m0, [t3+r10*2+416*4+0]
paddd m0, [t3+r10*2+416*4+8]
paddd m1, m0, [t3+r10*2+416*4+4]
psllw ym17, 2 ; a[ 0] 444
pslld m1, 2 ; b[ 0] 444
vmovdqa32 [t4+r10*1+416* 6], ym17
vmovdqa32 [t3+r10*2+416*12], m1
psubw ym17, ym16 ; a[ 0] 343
psubd m1, m0 ; b[ 0] 343
vmovdqa32 [t4+r10*1+416* 8], ym17
vmovdqa32 [t3+r10*2+416*16], m1
add r10, 32
jl .prep_n_loop
ret
ALIGN function_align
.n0: ; neighbor + output (even rows)
mov r10, wq
.n0_loop:
mova m3, [t4+r10*1+416*0+0]
paddw m3, [t4+r10*1+416*0+4]
paddw m1, m3, [t4+r10*1+416*0+2]
psllw m1, 2 ; a[ 1] 444
psubw m2, m1, m3 ; a[ 1] 343
paddw m3, m2, [t4+r10*1+416*4]
paddw m3, [t4+r10*1+416*6]
mova [t4+r10*1+416*4], m2
mova [t4+r10*1+416*6], m1
mova m16, [t3+r10*2+416*0+0]
paddd m16, [t3+r10*2+416*0+8]
paddd m1, m16, [t3+r10*2+416*0+4]
pslld m1, 2 ; b[ 1] 444
psubd m2, m1, m16 ; b[ 1] 343
paddd m16, m2, [t3+r10*2+416* 8+ 0]
paddd m16, [t3+r10*2+416*12+ 0]
mova [t3+r10*2+416* 8+ 0], m2
mova [t3+r10*2+416*12+ 0], m1
mova m17, [t3+r10*2+416*0+64]
paddd m17, [t3+r10*2+416*0+72]
paddd m1, m17, [t3+r10*2+416*0+68]
pslld m1, 2
psubd m2, m1, m17
paddd m17, m2, [t3+r10*2+416* 8+64]
paddd m17, [t3+r10*2+416*12+64]
mova [t3+r10*2+416* 8+64], m2
mova [t3+r10*2+416*12+64], m1
mova m0, [dstq+r10]
punpcklwd m1, m0, m6
punpcklwd m2, m3, m6
pmaddwd m2, m1 ; a * src
punpckhwd m1, m0, m6
punpckhwd m3, m6
pmaddwd m3, m1
vshufi32x4 m1, m16, m17, q2020
vshufi32x4 m16, m17, q3131
psubd m1, m2 ; b - a * src + (1 << 8)
psubd m16, m3
psrad m1, 9
psrad m16, 9
packssdw m1, m16
pmulhrsw m1, m7
paddw m0, m1
pmaxsw m0, m6
pminsw m0, m14
mova [dstq+r10], m0
add r10, 64
jl .n0_loop
add dstq, strideq
ret
ALIGN function_align
.n1: ; neighbor + output (odd rows)
mov r10, wq
.n1_loop:
mova m3, [t4+r10*1+416*2+0]
paddw m3, [t4+r10*1+416*2+4]
paddw m1, m3, [t4+r10*1+416*2+2]
psllw m1, 2 ; a[ 1] 444
psubw m2, m1, m3 ; a[ 1] 343
paddw m3, m2, [t4+r10*1+416*6]
paddw m3, [t4+r10*1+416*8]
mova [t4+r10*1+416*6], m1
mova [t4+r10*1+416*8], m2
mova m16, [t3+r10*2+416*4+0]
paddd m16, [t3+r10*2+416*4+8]
paddd m1, m16, [t3+r10*2+416*4+4]
pslld m1, 2 ; b[ 1] 444
psubd m2, m1, m16 ; b[ 1] 343
paddd m16, m2, [t3+r10*2+416*12+ 0]
paddd m16, [t3+r10*2+416*16+ 0]
mova [t3+r10*2+416*12+ 0], m1
mova [t3+r10*2+416*16+ 0], m2
mova m17, [t3+r10*2+416*4+64]
paddd m17, [t3+r10*2+416*4+72]
paddd m1, m17, [t3+r10*2+416*4+68]
pslld m1, 2
psubd m2, m1, m17
paddd m17, m2, [t3+r10*2+416*12+64]
paddd m17, [t3+r10*2+416*16+64]
mova [t3+r10*2+416*12+64], m1
mova [t3+r10*2+416*16+64], m2
mova m0, [dstq+r10]
punpcklwd m1, m0, m6
punpcklwd m2, m3, m6
pmaddwd m2, m1 ; a * src
punpckhwd m1, m0, m6
punpckhwd m3, m6
pmaddwd m3, m1
vshufi32x4 m1, m16, m17, q2020
vshufi32x4 m16, m17, q3131
psubd m1, m2 ; b - a * src + (1 << 8)
psubd m16, m3
psrad m1, 9
psrad m16, 9
packssdw m1, m16
pmulhrsw m1, m7
paddw m0, m1
pmaxsw m0, m6
pminsw m0, m14
mova [dstq+r10], m0
add r10, 64
jl .n1_loop
add dstq, strideq
ret
cglobal sgr_filter_mix_16bpc, 4, 14, 23, 416*66+8, dst, stride, left, lpf, \
w, h, edge, params
movifnidn wd, wm
mov paramsq, r6mp
lea r13, [r_ext_mask+72]
mov edged, r7m
movifnidn hd, hm
vpbroadcastd m7, [paramsq+8] ; w0 w1
pxor m6, m6
vpbroadcastd m8, [base+pd_8]
add wd, wd
vpbroadcastd m9, [base+pd_m9]
add lpfq, wq
vpbroadcastd m10, [base+pd_m25]
add dstq, wq
vpsubd m11, m6, [paramsq+0] {1to16} ; -s0
lea t3, [rsp+wq*2+416*24+8]
vpsubd m12, m6, [paramsq+4] {1to16} ; -s1
lea t4, [rsp+wq+416*52+8]
vpbroadcastd m13, [base+pw_164_455]
lea t1, [rsp+wq+12]
vpbroadcastd m14, [base+pw_61448]
neg wq
vpbroadcastd m15, [base+pd_m34816]
psllw m7, 2
vpbroadcastd m22, [base+pd_2147483648]
mov r10d, 0xfffffff8
mova m18, [sgr_x_by_x+64*0]
kmovd k1, r10d
mova m19, [sgr_x_by_x+64*1]
mov r10, 0x3333333333333333
mova m20, [sgr_x_by_x+64*2]
kmovq k2, r10
mova m21, [sgr_x_by_x+64*3]
test edgeb, 4 ; LR_HAVE_TOP
jz .no_top
call .h_top
add lpfq, strideq
mov t2, t1
call mangle(private_prefix %+ _sgr_filter_5x5_16bpc_avx512icl).top_fixup
add t1, 416*12
call .h_top
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
add r10, strideq
mov [rsp], r10 ; below
call .hv0
.main:
dec hd
jz .height1
add lpfq, strideq
call .hv1
call .prep_n
sub hd, 2
jl .extend_bottom
.main_loop:
add lpfq, strideq
call .hv0
test hd, hd
jz .odd_height
add lpfq, strideq
call .hv1
call .n0
call .n1
sub hd, 2
jge .main_loop
test edgeb, 8 ; LR_HAVE_BOTTOM
jz .extend_bottom
mov lpfq, [rsp]
call .hv0_bottom
add lpfq, strideq
call .hv1_bottom
.end:
call .n0
call .n1
.end2:
RET
.height1:
call .v1
call .prep_n
jmp .odd_height_end
.odd_height:
call .v1
call .n0
call .n1
.odd_height_end:
call .v0
call .v1
call .n0
jmp .end2
.extend_bottom:
call .v0
call .v1
jmp .end
.no_top:
lea r10, [lpfq+strideq*4]
mov lpfq, dstq
lea r10, [r10+strideq*2]
mov [rsp], r10
call .h
lea r10, [wq-4]
lea t2, [t1+416*12]
.top_fixup_loop:
mova m0, [t1+r10+416* 0]
mova m1, [t1+r10+416* 2]
mova m2, [t1+r10+416* 4]
paddw m0, m0
mova m3, [t1+r10+416* 6]
paddd m1, m1
mova m4, [t1+r10+416* 8]
paddd m2, m2
mova m5, [t1+r10+416*10]
mova [t2+r10+416* 0], m0
mova [t2+r10+416* 2], m1
mova [t2+r10+416* 4], m2
mova [t2+r10+416* 6], m3
mova [t2+r10+416* 8], m4
mova [t2+r10+416*10], m5
add r10, 64
jl .top_fixup_loop
call .v0
jmp .main
.h: ; horizontal boxsum
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
movq xm16, [leftq+2]
vmovdqu16 m16{k1}, [lpfq+wq-6]
add leftq, 8
jmp .h_main
.h_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-6]
jmp .h_main
.h_top:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .h_extend_left
.h_loop:
movu m16, [lpfq+r10- 2]
.h_main:
movu m17, [lpfq+r10+14]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .h_have_right
cmp r10d, -68
jl .h_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4
vpternlogd m17, m0, [r13+r10+16], 0xe4
.h_have_right:
palignr m3, m17, m16, 2
palignr m0, m17, m16, 4
paddw m1, m3, m0
punpcklwd m2, m3, m0
pmaddwd m2, m2
punpckhwd m3, m0
pmaddwd m3, m3
palignr m0, m17, m16, 6
paddw m1, m0 ; sum3
punpcklwd m4, m0, m6
vpdpwssd m2, m4, m4 ; sumsq3
punpckhwd m0, m6
vpdpwssd m3, m0, m0
shufpd m4, m16, m17, 0x55
punpcklwd m17, m4, m16
paddw m0, m16, m4
punpckhwd m4, m16
mova [t1+r10+416* 6], m1
mova [t1+r10+416* 8], m2
mova [t1+r10+416*10], m3
paddw m1, m0 ; sum5
vpdpwssd m2, m17, m17 ; sumsq5
vpdpwssd m3, m4, m4
mova [t1+r10+416* 0], m1
mova [t1+r10+416* 2], m2
mova [t1+r10+416* 4], m3
add r10, 64
jl .h_loop
ret
ALIGN function_align
.hv0: ; horizontal boxsum + vertical boxsum + ab3 (even rows)
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv0_extend_left
movq xm16, [leftq+2]
vmovdqu16 m16{k1}, [lpfq+wq-6]
add leftq, 8
jmp .hv0_main
.hv0_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-6]
jmp .hv0_main
.hv0_bottom:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv0_extend_left
.hv0_loop:
movu m16, [lpfq+r10- 2]
.hv0_main:
movu m17, [lpfq+r10+14]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .hv0_have_right
cmp r10d, -68
jl .hv0_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4
vpternlogd m17, m0, [r13+r10+16], 0xe4
.hv0_have_right:
palignr m3, m17, m16, 2
palignr m0, m17, m16, 4
paddw m1, m3, m0
punpcklwd m2, m3, m0
pmaddwd m2, m2
punpckhwd m3, m0
pmaddwd m3, m3
palignr m0, m17, m16, 6
paddw m1, m0 ; h sum3
punpcklwd m4, m0, m6
vpdpwssd m2, m4, m4 ; h sumsq3
punpckhwd m0, m6
vpdpwssd m3, m0, m0
shufpd m17, m16, m17, 0x55
paddw m4, m1, [t1+r10+416* 6]
paddd m5, m2, [t1+r10+416* 8]
mova [t1+r10+416* 6], m1
mova [t1+r10+416* 8], m2
paddw m1, m16
paddw m1, m17 ; h sum5
punpcklwd m0, m17, m16
vpdpwssd m2, m0, m0 ; h sumsq5
paddd m0, m3, [t1+r10+416*10]
mova [t1+r10+416*10], m3
punpckhwd m17, m16
vpdpwssd m3, m17, m17
mova [t3+r10*2+416*8+ 8], m1 ; we need a clean copy of the last row
mova [t3+r10*2+416*0+ 8], m2 ; in case height is odd
mova [t3+r10*2+416*0+72], m3
paddw m1, [t1+r10+416* 0]
paddd m2, [t1+r10+416* 2]
paddd m3, [t1+r10+416* 4]
mova [t1+r10+416* 0], m1
mova [t1+r10+416* 2], m2
mova [t1+r10+416* 4], m3
paddw m17, m4, [t2+r10+416* 6]
paddd m2, m5, [t2+r10+416* 8]
paddd m3, m0, [t2+r10+416*10]
mova [t2+r10+416* 6], m4
mova [t2+r10+416* 8], m5
mova [t2+r10+416*10], m0
paddd m2, m8
paddd m3, m8
psrld m2, 4 ; (a3 + 8) >> 4
psrld m3, 4
pmulld m2, m9 ; -((a3 + 8) >> 4) * 9
pmulld m3, m9
psrlw m5, m17, 1
pavgw m5, m6 ; (b3 + 2) >> 2
punpcklwd m4, m5, m6
vpdpwssd m2, m4, m4 ; -p3
punpckhwd m5, m6
vpdpwssd m3, m5, m5
punpcklwd m16, m6, m17 ; b3
punpckhwd m17, m6, m17
pminsd m2, m6
pminsd m3, m6
pmulld m2, m12 ; p3 * s1
pmulld m3, m12
pmaddwd m16, m13 ; b3 * 455
pmaddwd m17, m13
vpalignr m3{k2}, m2, m2, 2
mova m2, m20
paddusw m3, m14
psraw m3, 4 ; min(z3, 255) - 256
vpermt2b m2, m3, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m3
vpermi2b m3, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m3{k3}, m2 ; x3
pandn m2, m15, m3
psrld m3, 16
pmulld m16, m2
pmulld m17, m3
packssdw m2, m3
mova [t4+r10*1+416*2+4], m2
psubd m16, m15 ; x3 * b3 * 455 + (1 << 11) + (1 << 15)
psubd m17, m15
psrld m16, 12
psrld m17, 12
mova [t3+r10*2+416*4+ 8], xm16
mova [t3+r10*2+416*4+ 24], xm17
vextracti128 [t3+r10*2+416*4+ 40], ym16, 1
vextracti128 [t3+r10*2+416*4+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*4+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*4+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*4+104], m16, 3
vextracti32x4 [t3+r10*2+416*4+120], m17, 3
add r10, 64
jl .hv0_loop
ret
ALIGN function_align
.hv1: ; horizontal boxsums + vertical boxsums + ab (odd rows)
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv1_extend_left
movq xm16, [leftq+2]
vmovdqu16 m16{k1}, [lpfq+wq-6]
add leftq, 8
jmp .hv1_main
.hv1_extend_left:
vpbroadcastw xm16, [lpfq+wq]
vmovdqu16 m16{k1}, [lpfq+wq-6]
jmp .hv1_main
.hv1_bottom:
lea r10, [wq-4]
test edgeb, 1 ; LR_HAVE_LEFT
jz .hv1_extend_left
.hv1_loop:
movu m16, [lpfq+r10- 2]
.hv1_main:
movu m17, [lpfq+r10+14]
test edgeb, 2 ; LR_HAVE_RIGHT
jnz .hv1_have_right
cmp r10d, -68
jl .hv1_have_right
vpbroadcastw m0, [lpfq-2]
vpternlogd m16, m0, [r13+r10+ 0], 0xe4
vpternlogd m17, m0, [r13+r10+16], 0xe4
.hv1_have_right:
palignr m1, m17, m16, 2
palignr m3, m17, m16, 4
paddw m2, m1, m3
punpcklwd m0, m1, m3
pmaddwd m0, m0
punpckhwd m1, m3
pmaddwd m1, m1
palignr m3, m17, m16, 6
paddw m2, m3 ; h sum3
punpcklwd m5, m3, m6
vpdpwssd m0, m5, m5 ; h sumsq3
punpckhwd m3, m6
vpdpwssd m1, m3, m3
shufpd m3, m16, m17, 0x55
punpcklwd m5, m16, m3
paddw m4, m16, m3
punpckhwd m16, m3
paddw m17, m2, [t2+r10+416* 6]
mova [t2+r10+416* 6], m2
paddw m4, m2 ; h sum5
paddd m2, m0, [t2+r10+416* 8]
paddd m3, m1, [t2+r10+416*10]
mova [t2+r10+416* 8], m0
mova [t2+r10+416*10], m1
vpdpwssd m0, m5, m5 ; h sumsq5
vpdpwssd m1, m16, m16
paddd m2, m8
paddd m3, m8
psrld m2, 4 ; (a3 + 8) >> 4
psrld m3, 4
pmulld m2, m9 ; -((a3 + 8) >> 4) * 9
pmulld m3, m9
psrlw m16, m17, 1
pavgw m16, m6 ; (b3 + 2) >> 2
punpcklwd m5, m16, m6
vpdpwssd m2, m5, m5 ; -p3
punpckhwd m16, m6
vpdpwssd m3, m16, m16
punpcklwd m16, m6, m17 ; b3
punpckhwd m17, m6, m17
pminsd m2, m6
pminsd m3, m6
pmulld m2, m12 ; p3 * s1
pmulld m3, m12
pmaddwd m16, m13 ; b3 * 455
pmaddwd m17, m13
vpalignr m3{k2}, m2, m2, 2
mova m2, m20
paddusw m3, m14
psraw m3, 4 ; min(z3, 255) - 256
vpermt2b m2, m3, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m3
vpermi2b m3, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m3{k3}, m2 ; x3
pandn m2, m15, m3
psrld m3, 16
pmulld m16, m2
pmulld m17, m3
packssdw m2, m3
mova [t4+r10*1+416*4+4], m2
psubd m16, m15 ; x3 * b3 * 455 + (1 << 11) + (1 << 15)
psubd m17, m15
psrld m16, 12
psrld m17, 12
paddw m5, m4, [t2+r10+416*0]
paddd m2, m0, [t2+r10+416*2]
paddd m3, m1, [t2+r10+416*4]
paddw m5, [t1+r10+416*0]
paddd m2, [t1+r10+416*2]
paddd m3, [t1+r10+416*4]
mova [t2+r10+416*0], m4
mova [t2+r10+416*2], m0
mova [t2+r10+416*4], m1
mova [t3+r10*2+416*8+ 8], xm16
mova [t3+r10*2+416*8+ 24], xm17
vextracti128 [t3+r10*2+416*8+ 40], ym16, 1
vextracti128 [t3+r10*2+416*8+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*8+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*8+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*8+104], m16, 3
vextracti32x4 [t3+r10*2+416*8+120], m17, 3
paddd m2, m8
paddd m3, m8
psrld m2, 4 ; (a5 + 8) >> 4
psrld m3, 4
pmulld m2, m10 ; -((a5 + 8) >> 4) * 25
pmulld m3, m10
psrlw m17, m5, 1
pavgw m17, m6 ; (b5 + 2) >> 2
punpcklwd m16, m17, m6
vpdpwssd m2, m16, m16 ; -p5
punpckhwd m17, m6
vpdpwssd m3, m17, m17
punpcklwd m16, m5, m6 ; b5
punpckhwd m17, m5, m6
pmulld m2, m11 ; p5 * s0
pmulld m3, m11
pmaddwd m16, m13 ; b5 * 164
pmaddwd m17, m13
vpalignr m3{k2}, m2, m2, 2
mova m2, m20
pmaxsw m3, m6
paddusw m3, m14
psraw m3, 4 ; min(z5, 255) - 256
vpermt2b m2, m3, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m3
vpermi2b m3, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m3{k3}, m2 ; x5
pandn m2, m15, m3
psrld m3, 16
pmulld m16, m2
pmulld m17, m3
packssdw m2, m3
mova [t4+r10*1+416*0+4], m2
psubd m16, m15 ; x5 * b5 * 164 + (1 << 11) + (1 << 15)
psubd m17, m15
psrld m16, 12
psrld m17, 12
mova [t3+r10*2+416*0+ 8], xm16
mova [t3+r10*2+416*0+ 24], xm17
vextracti128 [t3+r10*2+416*0+ 40], ym16, 1
vextracti128 [t3+r10*2+416*0+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*0+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*0+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*0+104], m16, 3
vextracti32x4 [t3+r10*2+416*0+120], m17, 3
add r10, 64
jl .hv1_loop
mov r10, t2
mov t2, t1
mov t1, r10
ret
.v0: ; vertical boxsums + ab3 (even rows)
lea r10, [wq-4]
.v0_loop:
mova m16, [t1+r10+416* 6]
mova m2, [t1+r10+416* 8]
mova m3, [t1+r10+416*10]
paddw m16, m16
paddd m2, m2
paddd m3, m3
paddw m17, m16, [t2+r10+416* 6]
paddd m4, m2, [t2+r10+416* 8]
paddd m5, m3, [t2+r10+416*10]
mova [t2+r10+416* 6], m16
mova [t2+r10+416* 8], m2
mova [t2+r10+416*10], m3
paddd m4, m8
paddd m5, m8
psrld m4, 4 ; (a3 + 8) >> 4
psrld m5, 4
pmulld m4, m9 ; -((a3 + 8) >> 4) * 9
pmulld m5, m9
psrlw m3, m17, 1
pavgw m3, m6 ; (b3 + 2) >> 2
punpcklwd m2, m3, m6
vpdpwssd m4, m2, m2 ; -p3
punpckhwd m3, m6
vpdpwssd m5, m3, m3
punpcklwd m16, m6, m17 ; b3
punpckhwd m17, m6, m17
pminsd m4, m6
pminsd m5, m6
pmulld m4, m12 ; p3 * s1
pmulld m5, m12
pmaddwd m16, m13 ; b3 * 455
pmaddwd m17, m13
vpalignr m5{k2}, m4, m4, 2
mova m4, m20
paddusw m5, m14
psraw m5, 4 ; min(z3, 255) - 256
vpermt2b m4, m5, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m5
vpermi2b m5, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m5{k3}, m4 ; x3
pandn m4, m15, m5
psrld m5, 16
pmulld m16, m4
pmulld m17, m5
packssdw m4, m5
mova [t4+r10*1+416*2+4], m4
psubd m16, m15 ; x3 * b3 * 455 + (1 << 11) + (1 << 15)
psubd m17, m15
psrld m16, 12
psrld m17, 12
mova m3, [t1+r10+416*0]
mova m4, [t1+r10+416*2]
mova m5, [t1+r10+416*4]
mova [t3+r10*2+416*8+ 8], m3
mova [t3+r10*2+416*0+ 8], m4
mova [t3+r10*2+416*0+72], m5
paddw m3, m3 ; cc5
paddd m4, m4
paddd m5, m5
mova [t1+r10+416*0], m3
mova [t1+r10+416*2], m4
mova [t1+r10+416*4], m5
mova [t3+r10*2+416*4+ 8], xm16
mova [t3+r10*2+416*4+ 24], xm17
vextracti128 [t3+r10*2+416*4+ 40], ym16, 1
vextracti128 [t3+r10*2+416*4+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*4+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*4+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*4+104], m16, 3
vextracti32x4 [t3+r10*2+416*4+120], m17, 3
add r10, 64
jl .v0_loop
ret
.v1: ; vertical boxsums + ab (odd rows)
lea r10, [wq-4]
.v1_loop:
mova m16, [t1+r10+416* 6]
mova m2, [t1+r10+416* 8]
mova m3, [t1+r10+416*10]
paddw m17, m16, [t2+r10+416* 6]
paddd m4, m2, [t2+r10+416* 8]
paddd m5, m3, [t2+r10+416*10]
mova [t2+r10+416* 6], m16
mova [t2+r10+416* 8], m2
mova [t2+r10+416*10], m3
paddd m4, m8
paddd m5, m8
psrld m4, 4 ; (a3 + 8) >> 4
psrld m5, 4
pmulld m4, m9 ; -((a3 + 8) >> 4) * 9
pmulld m5, m9
psrlw m3, m17, 1
pavgw m3, m6 ; (b3 + 2) >> 2
punpcklwd m2, m3, m6
vpdpwssd m4, m2, m2 ; -p3
punpckhwd m3, m6
vpdpwssd m5, m3, m3
punpcklwd m16, m6, m17 ; b3
punpckhwd m17, m6, m17
pminsd m4, m6
pminsd m5, m6
pmulld m4, m12 ; p3 * s1
pmulld m5, m12
pmaddwd m16, m13 ; b3 * 455
pmaddwd m17, m13
vpalignr m5{k2}, m4, m4, 2
mova m4, m20
paddusw m5, m14
psraw m5, 4 ; min(z3, 255) - 256
vpermt2b m4, m5, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m5
vpermi2b m5, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m5{k3}, m4 ; x3
pandn m4, m15, m5
psrld m5, 16
pmulld m16, m4
pmulld m17, m5
packssdw m4, m5
mova [t4+r10*1+416*4+4], m4
psubd m16, m15 ; x3 * b3 * 455 + (1 << 11) + (1 << 15)
psubd m17, m15
psrld m16, 12
psrld m17, 12
mova m0, [t3+r10*2+416*8+ 8]
mova m4, [t3+r10*2+416*0+ 8]
mova m5, [t3+r10*2+416*0+72]
paddw m1, m0, [t2+r10+416*0]
paddd m2, m4, [t2+r10+416*2]
paddd m3, m5, [t2+r10+416*4]
paddw m1, [t1+r10+416*0]
paddd m2, [t1+r10+416*2]
paddd m3, [t1+r10+416*4]
mova [t2+r10+416*0], m0
mova [t2+r10+416*2], m4
mova [t2+r10+416*4], m5
mova [t3+r10*2+416*8+ 8], xm16
mova [t3+r10*2+416*8+ 24], xm17
vextracti128 [t3+r10*2+416*8+ 40], ym16, 1
vextracti128 [t3+r10*2+416*8+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*8+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*8+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*8+104], m16, 3
vextracti32x4 [t3+r10*2+416*8+120], m17, 3
paddd m2, m8
paddd m3, m8
psrld m2, 4 ; (a5 + 8) >> 4
psrld m3, 4
pmulld m2, m10 ; -((a5 + 8) >> 4) * 25
pmulld m3, m10
psrlw m5, m1, 1
pavgw m5, m6 ; (b5 + 2) >> 2
punpcklwd m4, m5, m6
vpdpwssd m2, m4, m4 ; -p5
punpckhwd m5, m6
vpdpwssd m3, m5, m5
punpcklwd m16, m1, m6 ; b5
punpckhwd m17, m1, m6
pmulld m2, m11 ; p5 * s0
pmulld m3, m11
pmaddwd m16, m13 ; b5 * 164
pmaddwd m17, m13
vpalignr m3{k2}, m2, m2, 2
mova m2, m20
pmaxsw m3, m6
paddusw m3, m14
psraw m3, 4 ; min(z5, 255) - 256
vpermt2b m2, m3, m21 ; sgr_x_by_x[128..255]
vpmovb2m k3, m3
vpermi2b m3, m18, m19 ; sgr_x_by_x[ 0..127]
vmovdqu8 m3{k3}, m2 ; x5
pandn m2, m15, m3
psrld m3, 16
pmulld m16, m2
pmulld m17, m3
packssdw m2, m3
mova [t4+r10*1+416*0+4], m2
psubd m16, m15 ; x5 * b5 * 164 + (1 << 11) + (1 << 15)
psubd m17, m15
psrld m16, 12
psrld m17, 12
mova [t3+r10*2+416*0+ 8], xm16
mova [t3+r10*2+416*0+ 24], xm17
vextracti128 [t3+r10*2+416*0+ 40], ym16, 1
vextracti128 [t3+r10*2+416*0+ 56], ym17, 1
vextracti32x4 [t3+r10*2+416*0+ 72], m16, 2
vextracti32x4 [t3+r10*2+416*0+ 88], m17, 2
vextracti32x4 [t3+r10*2+416*0+104], m16, 3
vextracti32x4 [t3+r10*2+416*0+120], m17, 3
add r10, 64
jl .v1_loop
mov r10, t2
mov t2, t1
mov t1, r10
ret
.prep_n: ; initial neighbor setup
mov r10, wq
.prep_n_loop:
movu ym0, [t4+r10*1+416*0+2]
paddw ym2, ym0, [t4+r10*1+416*0+0]
paddw ym2, [t4+r10*1+416*0+4]
movu m1, [t3+r10*2+416*0+4]
paddd m3, m1, [t3+r10*2+416*0+0]
paddd m3, [t3+r10*2+416*0+8]
paddw ym0, ym2
paddd m1, m3
psllw ym2, 2
pslld m3, 2
paddw ym0, ym2 ; a5 565
paddd m1, m3 ; b5 565
mova [t4+r10*1+416* 6], ym0
mova [t3+r10*2+416*12], m1
mova ym0, [t4+r10*1+416*2+0]
paddw ym0, [t4+r10*1+416*2+4]
paddw ym2, ym0, [t4+r10*1+416*2+2]
mova m1, [t3+r10*2+416*4+0]
paddd m1, [t3+r10*2+416*4+8]
paddd m3, m1, [t3+r10*2+416*4+4]
psllw ym2, 2 ; a3[-1] 444
pslld m3, 2 ; b3[-1] 444
psubw ym2, ym0 ; a3[-1] 343
psubd m3, m1 ; b3[-1] 343
mova [t4+r10*1+416* 8], ym2
mova [t3+r10*2+416*16], m3
mova ym0, [t4+r10*1+416*4+0]
paddw ym0, [t4+r10*1+416*4+4]
paddw ym2, ym0, [t4+r10*1+416*4+2]
mova m1, [t3+r10*2+416*8+0]
paddd m1, [t3+r10*2+416*8+8]
paddd m3, m1, [t3+r10*2+416*8+4]
psllw ym2, 2 ; a3[ 0] 444
pslld m3, 2 ; b3[ 0] 444
mova [t4+r10*1+416*10], ym2
mova [t3+r10*2+416*20], m3
psubw ym2, ym0 ; a3[ 0] 343
psubd m3, m1 ; b3[ 0] 343
mova [t4+r10*1+416*12], ym2
mova [t3+r10*2+416*24], m3
add r10, 32
jl .prep_n_loop
ret
ALIGN function_align
.n0: ; neighbor + output (even rows)
mov r10, wq
.n0_loop:
movu ym2, [t4+r10*1+2]
paddw ym0, ym2, [t4+r10*1+0]
paddw ym0, [t4+r10*1+4]
paddw ym2, ym0
psllw ym0, 2
paddw ym0, ym2 ; a5
movu m1, [t3+r10*2+4]
paddd m4, m1, [t3+r10*2+0]
paddd m4, [t3+r10*2+8]
paddd m1, m4
pslld m4, 2
paddd m4, m1 ; b5
paddw ym2, ym0, [t4+r10*1+416* 6]
mova [t4+r10*1+416* 6], ym0
paddd m0, m4, [t3+r10*2+416*12]
mova [t3+r10*2+416*12], m4
mova ym3, [t4+r10*1+416*2+0]
paddw ym3, [t4+r10*1+416*2+4]
paddw ym5, ym3, [t4+r10*1+416*2+2]
psllw ym5, 2 ; a3[ 1] 444
psubw ym4, ym5, ym3 ; a3[ 1] 343
paddw ym3, ym4, [t4+r10*1+416* 8]
paddw ym3, [t4+r10*1+416*10]
mova [t4+r10*1+416* 8], ym4
mova [t4+r10*1+416*10], ym5
mova m1, [t3+r10*2+416*4+0]
paddd m1, [t3+r10*2+416*4+8]
paddd m5, m1, [t3+r10*2+416*4+4]
pslld m5, 2 ; b3[ 1] 444
psubd m4, m5, m1 ; b3[ 1] 343
paddd m1, m4, [t3+r10*2+416*16]
paddd m1, [t3+r10*2+416*20]
mova [t3+r10*2+416*16], m4
mova [t3+r10*2+416*20], m5
pmovzxwd m4, [dstq+r10]
pmovzxwd m2, ym2 ; a5
pmovzxwd m3, ym3 ; a3
pmaddwd m2, m4 ; a5 * src
pmaddwd m3, m4 ; a3 * src
vpshldd m4, m22, 13
psubd m0, m2 ; b5 - a5 * src + (1 << 8)
psubd m1, m3 ; b3 - a3 * src + (1 << 8)
psrld m0, 9
pslld m1, 7
vpblendmb m0{k2}, m1, m0
vpdpwssd m4, m0, m7
psrad m4, 7
pmaxsd m4, m6
vpmovusdw ym16, m4 ; clip
psrlw ym16, 6
mova [dstq+r10], ym16
add r10, 32
jl .n0_loop
add dstq, strideq
ret
ALIGN function_align
.n1: ; neighbor + output (odd rows)
mov r10, wq
.n1_loop:
mova ym3, [t4+r10*1+416*4+0]
paddw ym3, [t4+r10*1+416*4+4]
paddw ym5, ym3, [t4+r10*1+416*4+2]
psllw ym5, 2 ; a3[ 1] 444
psubw ym4, ym5, ym3 ; a3[ 1] 343
paddw ym3, ym4, [t4+r10*1+416*12]
paddw ym3, [t4+r10*1+416*10]
mova [t4+r10*1+416*10], ym5
mova [t4+r10*1+416*12], ym4
mova m0, [t3+r10*2+416*8+0]
paddd m0, [t3+r10*2+416*8+8]
paddd m5, m0, [t3+r10*2+416*8+4]
pslld m5, 2 ; b3[ 1] 444
psubd m4, m5, m0 ; b3[ 1] 343
paddd m0, m4, [t3+r10*2+416*24]
paddd m0, [t3+r10*2+416*20]
mova [t3+r10*2+416*20], m5
mova [t3+r10*2+416*24], m4
pmovzxwd m4, [dstq+r10]
pmovzxwd m2, [t4+r10*1+416* 6]
pmovzxwd m3, ym3
mova m1, [t3+r10*2+416*12]
pmaddwd m2, m4 ; a5 * src
pmaddwd m3, m4 ; a3 * src
vpshldd m4, m22, 13
psubd m1, m2 ; b5 - a5 * src + (1 << 8)
psubd m0, m3 ; b3 - a3 * src + (1 << 8)
pslld m0, 7
vpalignr m0{k2}, m1, m1, 1
vpdpwssd m4, m0, m7
psrad m4, 7
pmaxsd m4, m6
vpmovusdw ym16, m4 ; clip
psrlw ym16, 6
mova [dstq+r10], ym16
add r10, 32
jl .n1_loop
add dstq, strideq
ret
%endif ; ARCH_X86_64