| Message ID | 20210204113259.20112-2-josh@itanimul.li |
|---|---|
| State | New |
| Headers | show |
| Series | avcodec/aarch64/hevcdsp | expand |
| Context | Check | Description |
|---|---|---|
| andriy/x86_make | success | Make finished |
| andriy/x86_make_fate | success | Make fate finished |
| andriy/PPC64_make | success | Make finished |
| andriy/PPC64_make_fate | success | Make fate finished |
On Thu, 4 Feb 2021, Josh Dekker wrote: > From: Reimar Döffinger <Reimar.Doeffinger@gmx.de> > > Makes SIMD-optimized 8x8 and 16x16 idcts for 8 and 10 bit depth > available on aarch64. > For a UHD HDR (10 bit) sample video these were consuming the most time > and this optimization reduced overall decode time from 19.4s to 16.4s, > approximately 15% speedup. > Test sample was the first 300 frames of "LG 4K HDR Demo - New York.ts", > running on Apple M1. > > Signed-off-by: Josh Dekker <josh@itanimul.li> > --- > libavcodec/aarch64/Makefile | 2 + > libavcodec/aarch64/hevcdsp_idct_neon.S | 380 ++++++++++++++++++++++ > libavcodec/aarch64/hevcdsp_init_aarch64.c | 45 +++ > libavcodec/hevcdsp.c | 2 + > libavcodec/hevcdsp.h | 1 + > 5 files changed, 430 insertions(+) > create mode 100644 libavcodec/aarch64/hevcdsp_idct_neon.S > create mode 100644 libavcodec/aarch64/hevcdsp_init_aarch64.c > > diff --git a/libavcodec/aarch64/Makefile b/libavcodec/aarch64/Makefile > index f6434e40da..2ea1d74a38 100644 > --- a/libavcodec/aarch64/Makefile > +++ b/libavcodec/aarch64/Makefile > @@ -61,3 +61,5 @@ NEON-OBJS-$(CONFIG_VP9_DECODER) += aarch64/vp9itxfm_16bpp_neon.o \ > aarch64/vp9lpf_neon.o \ > aarch64/vp9mc_16bpp_neon.o \ > aarch64/vp9mc_neon.o > +NEON-OBJS-$(CONFIG_HEVC_DECODER) += aarch64/hevcdsp_idct_neon.o \ > + aarch64/hevcdsp_init_aarch64.o > diff --git a/libavcodec/aarch64/hevcdsp_idct_neon.S b/libavcodec/aarch64/hevcdsp_idct_neon.S > new file mode 100644 > index 0000000000..c70d6a906d > --- /dev/null > +++ b/libavcodec/aarch64/hevcdsp_idct_neon.S > @@ -0,0 +1,380 @@ > +/* > + * ARM NEON optimised IDCT functions for HEVC decoding > + * Copyright (c) 2014 Seppo Tomperi <seppo.tomperi@vtt.fi> > + * Copyright (c) 2017 Alexandra Hájková > + * > + * Ported from arm/hevcdsp_idct_neon.S by > + * Copyright (c) 2020 Reimar Döffinger > + * > + * This file is part of FFmpeg. > + * > + * FFmpeg is free software; you can redistribute it and/or > + * modify it under the terms of the GNU Lesser General Public > + * License as published by the Free Software Foundation; either > + * version 2.1 of the License, or (at your option) any later version. > + * > + * FFmpeg is distributed in the hope that it will be useful, > + * but WITHOUT ANY WARRANTY; without even the implied warranty of > + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + * Lesser General Public License for more details. > + * > + * You should have received a copy of the GNU Lesser General Public > + * License along with FFmpeg; if not, write to the Free Software > + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA > + */ > + > +#include "libavutil/aarch64/asm.S" > + > +const trans, align=4 > + .short 64, 83, 64, 36 > + .short 89, 75, 50, 18 > + .short 90, 87, 80, 70 > + .short 57, 43, 25, 9 > + .short 90, 90, 88, 85 > + .short 82, 78, 73, 67 > + .short 61, 54, 46, 38 > + .short 31, 22, 13, 4 > +endconst > + > +.macro sum_sub out, in, c, op, p > + .ifc \op, + > + smlal\p \out, \in, \c > + .else > + smlsl\p \out, \in, \c > + .endif > +.endm > + > +.macro fixsqrshrn d, dt, n, m > + .ifc \dt, .8h > + sqrshrn2 \d\dt, \n\().4s, \m > + .else > + sqrshrn \n\().4h, \n\().4s, \m > + mov \d\().d[0], \n\().d[0] > + .endif > +.endm Just mentioning it again even though it was brought up in the previous review; this looks pretty clumsy, and causes stalls on in-order cores. I understand that this comes from trying to port the existing arm assembly with as little modification as possible. It would be good to look into restructuring things so we don't need contraptions like this - but as the patch itself gives good speedup as is (3.9-9x speedup over unvectorized code from GCC 7), I guess one could file that under "future work", not a requirement per se. > + > +// uses and clobbers v28-v31 as temp registers > +.macro tr_4x4_8 in0, in1, in2, in3, out0, out1, out2, out3, p1, p2 > + sshll\p1 v28.4s, \in0, #6 > + mov v29.16b, v28.16b > + smull\p1 v30.4s, \in1, v0.h[1] > + smull\p1 v31.4s, \in1, v0.h[3] Nit: Vertical alignment could be a bit more consistent, but it's hard to keep it sensible in macro-heavy code anyway. > + smlal\p2 v28.4s, \in2, v0.h[0] //e0 > + smlsl\p2 v29.4s, \in2, v0.h[0] //e1 > + smlal\p2 v30.4s, \in3, v0.h[3] //o0 > + smlsl\p2 v31.4s, \in3, v0.h[1] //o1 > + > + add \out0, v28.4s, v30.4s > + add \out1, v29.4s, v31.4s > + sub \out2, v29.4s, v31.4s > + sub \out3, v28.4s, v30.4s > +.endm > + > +.macro transpose8_4x4 r0, r1, r2, r3 > + trn1 v2.8h, \r0\().8h, \r1\().8h > + trn2 v3.8h, \r0\().8h, \r1\().8h > + trn1 v4.8h, \r2\().8h, \r3\().8h > + trn2 v5.8h, \r2\().8h, \r3\().8h > + trn1 \r0\().4s, v2.4s, v4.4s > + trn2 \r2\().4s, v2.4s, v4.4s > + trn1 \r1\().4s, v3.4s, v5.4s > + trn2 \r3\().4s, v3.4s, v5.4s > +.endm > + > +.macro transpose_8x8 r0, r1, r2, r3, r4, r5, r6, r7 > + transpose8_4x4 \r0, \r1, \r2, \r3 > + transpose8_4x4 \r4, \r5, \r6, \r7 > +.endm > + > +.macro tr_8x4 shift, in0,in0t, in1,in1t, in2,in2t, in3,in3t, in4,in4t, in5,in5t, in6,in6t, in7,in7t, p1, p2 > + tr_4x4_8 \in0\in0t, \in2\in2t, \in4\in4t, \in6\in6t, v24.4s, v25.4s, v26.4s, v27.4s, \p1, \p2 > + > + smull\p1 v30.4s, \in1\in1t, v0.h[6] > + smull\p1 v28.4s, \in1\in1t, v0.h[4] > + smull\p1 v29.4s, \in1\in1t, v0.h[5] > + sum_sub v30.4s, \in3\in3t, v0.h[4], -, \p1 > + sum_sub v28.4s, \in3\in3t, v0.h[5], +, \p1 > + sum_sub v29.4s, \in3\in3t, v0.h[7], -, \p1 > + > + sum_sub v30.4s, \in5\in5t, v0.h[7], +, \p2 > + sum_sub v28.4s, \in5\in5t, v0.h[6], +, \p2 > + sum_sub v29.4s, \in5\in5t, v0.h[4], -, \p2 > + > + sum_sub v30.4s, \in7\in7t, v0.h[5], +, \p2 > + sum_sub v28.4s, \in7\in7t, v0.h[7], +, \p2 > + sum_sub v29.4s, \in7\in7t, v0.h[6], -, \p2 > + > + add v31.4s, v26.4s, v30.4s > + sub v26.4s, v26.4s, v30.4s > + fixsqrshrn \in2,\in2t, v31, \shift > + > + > + smull\p1 v31.4s, \in1\in1t, v0.h[7] > + sum_sub v31.4s, \in3\in3t, v0.h[6], -, \p1 > + sum_sub v31.4s, \in5\in5t, v0.h[5], +, \p2 > + sum_sub v31.4s, \in7\in7t, v0.h[4], -, \p2 > + fixsqrshrn \in5,\in5t, v26, \shift > + > + > + add v26.4s, v24.4s, v28.4s > + sub v24.4s, v24.4s, v28.4s > + add v28.4s, v25.4s, v29.4s > + sub v25.4s, v25.4s, v29.4s > + add v30.4s, v27.4s, v31.4s > + sub v27.4s, v27.4s, v31.4s > + > + fixsqrshrn \in0,\in0t, v26, \shift > + fixsqrshrn \in7,\in7t, v24, \shift > + fixsqrshrn \in1,\in1t, v28, \shift > + fixsqrshrn \in6,\in6t, v25, \shift > + fixsqrshrn \in3,\in3t, v30, \shift > + fixsqrshrn \in4,\in4t, v27, \shift > +.endm > + > +.macro idct_8x8 bitdepth > +function ff_hevc_idct_8x8_\bitdepth\()_neon, export=1 > +//x0 - coeffs > + mov x1, x0 > + ld1 {v16.8h-v19.8h}, [x1], #64 > + ld1 {v20.8h-v23.8h}, [x1] > + > + movrel x1, trans > + ld1 {v0.8h}, [x1] > + > + tr_8x4 7, v16,.4h, v17,.4h, v18,.4h, v19,.4h, v20,.4h, v21,.4h, v22,.4h, v23,.4h > + tr_8x4 7, v16,.8h, v17,.8h, v18,.8h, v19,.8h, v20,.8h, v21,.8h, v22,.8h, v23,.8h, 2, 2 > + > + transpose_8x8 v16, v17, v18, v19, v20, v21, v22, v23 > + > + tr_8x4 20 - \bitdepth, v16,.4h, v17,.4h, v18,.4h, v19,.4h, v16,.8h, v17,.8h, v18,.8h, v19,.8h, , 2 > + tr_8x4 20 - \bitdepth, v20,.4h, v21,.4h, v22,.4h, v23,.4h, v20,.8h, v21,.8h, v22,.8h, v23,.8h, , 2 > + > + transpose_8x8 v16, v17, v18, v19, v20, v21, v22, v23 > + > + mov x1, x0 > + st1 {v16.8h-v19.8h}, [x1], #64 > + st1 {v20.8h-v23.8h}, [x1] > + > + ret > +endfunc > +.endm > + > +.macro butterfly e, o, tmp_p, tmp_m > + add \tmp_p, \e, \o > + sub \tmp_m, \e, \o > +.endm > + > +.macro tr16_8x4 in0, in1, in2, in3, offset > + tr_4x4_8 \in0\().4h, \in1\().4h, \in2\().4h, \in3\().4h, v24.4s, v25.4s, v26.4s, v27.4s > + > + smull2 v28.4s, \in0\().8h, v0.h[4] > + smull2 v29.4s, \in0\().8h, v0.h[5] > + smull2 v30.4s, \in0\().8h, v0.h[6] > + smull2 v31.4s, \in0\().8h, v0.h[7] > + sum_sub v28.4s, \in1\().8h, v0.h[5], +, 2 > + sum_sub v29.4s, \in1\().8h, v0.h[7], -, 2 > + sum_sub v30.4s, \in1\().8h, v0.h[4], -, 2 > + sum_sub v31.4s, \in1\().8h, v0.h[6], -, 2 > + > + sum_sub v28.4s, \in2\().8h, v0.h[6], +, 2 > + sum_sub v29.4s, \in2\().8h, v0.h[4], -, 2 > + sum_sub v30.4s, \in2\().8h, v0.h[7], +, 2 > + sum_sub v31.4s, \in2\().8h, v0.h[5], +, 2 > + > + sum_sub v28.4s, \in3\().8h, v0.h[7], +, 2 > + sum_sub v29.4s, \in3\().8h, v0.h[6], -, 2 > + sum_sub v30.4s, \in3\().8h, v0.h[5], +, 2 > + sum_sub v31.4s, \in3\().8h, v0.h[4], -, 2 > + > + butterfly v24.4s, v28.4s, v16.4s, v23.4s > + butterfly v25.4s, v29.4s, v17.4s, v22.4s > + butterfly v26.4s, v30.4s, v18.4s, v21.4s > + butterfly v27.4s, v31.4s, v19.4s, v20.4s > + add x4, sp, #\offset > + st1 {v16.4s-v19.4s}, [x4], #64 > + st1 {v20.4s-v23.4s}, [x4] > +.endm > + > +.macro load16 in0, in1, in2, in3 > + ld1 {\in0}[0], [x1], x2 > + ld1 {\in0}[1], [x3], x2 > + ld1 {\in1}[0], [x1], x2 > + ld1 {\in1}[1], [x3], x2 > + ld1 {\in2}[0], [x1], x2 > + ld1 {\in2}[1], [x3], x2 > + ld1 {\in3}[0], [x1], x2 > + ld1 {\in3}[1], [x3], x2 > +.endm > + > +.macro add_member in, t0, t1, t2, t3, t4, t5, t6, t7, op0, op1, op2, op3, op4, op5, op6, op7, p > + sum_sub v21.4s, \in, \t0, \op0, \p > + sum_sub v22.4s, \in, \t1, \op1, \p > + sum_sub v23.4s, \in, \t2, \op2, \p > + sum_sub v24.4s, \in, \t3, \op3, \p > + sum_sub v25.4s, \in, \t4, \op4, \p > + sum_sub v26.4s, \in, \t5, \op5, \p > + sum_sub v27.4s, \in, \t6, \op6, \p > + sum_sub v28.4s, \in, \t7, \op7, \p > +.endm > + > +.macro butterfly16 in0, in1, in2, in3, in4, in5, in6, in7 > + add v20.4s, \in0, \in1 > + sub \in0, \in0, \in1 > + add \in1, \in2, \in3 > + sub \in2, \in2, \in3 > + add \in3, \in4, \in5 > + sub \in4, \in4, \in5 > + add \in5, \in6, \in7 > + sub \in6, \in6, \in7 > +.endm > + > +.macro store16 in0, in1, in2, in3, rx > + st1 {\in0}[0], [x1], x2 > + st1 {\in0}[1], [x3], \rx > + st1 {\in1}[0], [x1], x2 > + st1 {\in1}[1], [x3], \rx > + st1 {\in2}[0], [x1], x2 > + st1 {\in2}[1], [x3], \rx > + st1 {\in3}[0], [x1], x2 > + st1 {\in3}[1], [x3], \rx > +.endm > + > +.macro scale out0, out1, out2, out3, in0, in1, in2, in3, in4, in5, in6, in7, shift > + sqrshrn \out0\().4h, \in0, \shift > + sqrshrn2 \out0\().8h, \in1, \shift > + sqrshrn \out1\().4h, \in2, \shift > + sqrshrn2 \out1\().8h, \in3, \shift > + sqrshrn \out2\().4h, \in4, \shift > + sqrshrn2 \out2\().8h, \in5, \shift > + sqrshrn \out3\().4h, \in6, \shift > + sqrshrn2 \out3\().8h, \in7, \shift > +.endm > + > +.macro transpose16_4x4_2 r0, r1, r2, r3 > + // lower halves > + trn1 v2.4h, \r0\().4h, \r1\().4h > + trn2 v3.4h, \r0\().4h, \r1\().4h > + trn1 v4.4h, \r2\().4h, \r3\().4h > + trn2 v5.4h, \r2\().4h, \r3\().4h > + trn1 v6.2s, v2.2s, v4.2s > + trn2 v7.2s, v2.2s, v4.2s > + trn1 v2.2s, v3.2s, v5.2s > + trn2 v4.2s, v3.2s, v5.2s > + mov \r0\().d[0], v6.d[0] > + mov \r2\().d[0], v7.d[0] > + mov \r1\().d[0], v2.d[0] > + mov \r3\().d[0], v4.d[0] > + > + // upper halves in reverse order > + trn1 v2.8h, \r3\().8h, \r2\().8h > + trn2 v3.8h, \r3\().8h, \r2\().8h > + trn1 v4.8h, \r1\().8h, \r0\().8h > + trn2 v5.8h, \r1\().8h, \r0\().8h > + trn1 v6.4s, v2.4s, v4.4s > + trn2 v7.4s, v2.4s, v4.4s > + trn1 v2.4s, v3.4s, v5.4s > + trn2 v4.4s, v3.4s, v5.4s > + mov \r3\().d[1], v6.d[1] > + mov \r1\().d[1], v7.d[1] > + mov \r2\().d[1], v2.d[1] > + mov \r0\().d[1], v4.d[1] > +.endm > + > +.macro tr_16x4 name, shift, offset, step > +function func_tr_16x4_\name > + mov x1, x5 > + add x3, x5, #(\step * 64) > + mov x2, #(\step * 128) > + load16 v16.d, v17.d, v18.d, v19.d > + movrel x1, trans > + ld1 {v0.8h}, [x1] > + > + tr16_8x4 v16, v17, v18, v19, \offset > + > + add x1, x5, #(\step * 32) > + add x3, x5, #(\step * 3 *32) Nit: The vertical alignment is a bit all over the place here as well (but that was the case in the original patch as well). For operand lists without layout specifiers (like .16b and such), I'd suggest aligning the second and third operand like all operands would be 3 chars, like x10 and so on. I.e. the "x1, x5" bit is fine here, but add another space before the #, and the mov/add above could use same alignment. > + mov x2, #(\step * 128) > + load16 v20.d, v17.d, v18.d, v19.d > + movrel x1, trans, 16 > + ld1 {v1.8h}, [x1] > + smull v21.4s, v20.4h, v1.h[0] > + smull v22.4s, v20.4h, v1.h[1] These are neatly aligned here. I normally try to keep such things aligned so there's no misaligned if one operand turns into e.g. v16.16b (i.e. all columns one char wider), but for isolated segments like this that line up nicely anyway, this can be ok too. > + smull v23.4s, v20.4h, v1.h[2] > + smull v24.4s, v20.4h, v1.h[3] > + smull v25.4s, v20.4h, v1.h[4] > + smull v26.4s, v20.4h, v1.h[5] > + smull v27.4s, v20.4h, v1.h[6] > + smull v28.4s, v20.4h, v1.h[7] > + > + add_member v20.8h, v1.h[1], v1.h[4], v1.h[7], v1.h[5], v1.h[2], v1.h[0], v1.h[3], v1.h[6], +, +, +, -, -, -, -, -, 2 > + add_member v17.4h, v1.h[2], v1.h[7], v1.h[3], v1.h[1], v1.h[6], v1.h[4], v1.h[0], v1.h[5], +, +, -, -, -, +, +, + > + add_member v17.8h, v1.h[3], v1.h[5], v1.h[1], v1.h[7], v1.h[0], v1.h[6], v1.h[2], v1.h[4], +, -, -, +, +, +, -, -, 2 > + add_member v18.4h, v1.h[4], v1.h[2], v1.h[6], v1.h[0], v1.h[7], v1.h[1], v1.h[5], v1.h[3], +, -, -, +, -, -, +, + > + add_member v18.8h, v1.h[5], v1.h[0], v1.h[4], v1.h[6], v1.h[1], v1.h[3], v1.h[7], v1.h[2], +, -, +, +, -, +, +, -, 2 > + add_member v19.4h, v1.h[6], v1.h[3], v1.h[0], v1.h[2], v1.h[5], v1.h[7], v1.h[4], v1.h[1], +, -, +, -, +, +, -, + > + add_member v19.8h, v1.h[7], v1.h[6], v1.h[5], v1.h[4], v1.h[3], v1.h[2], v1.h[1], v1.h[0], +, -, +, -, +, -, +, -, 2 > + > + add x4, sp, #\offset > + ld1 {v16.4s-v19.4s}, [x4], #64 > + > + butterfly16 v16.4s, v21.4s, v17.4s, v22.4s, v18.4s, v23.4s, v19.4s, v24.4s > + scale v29, v30, v31, v24, v20.4s, v16.4s, v21.4s, v17.4s, v22.4s, v18.4s, v23.4s, v19.4s, \shift > + transpose16_4x4_2 v29, v30, v31, v24 > + mov x1, x6 > + add x3, x6, #(24 +3*32) > + mov x2, #32 > + mov x4, #-32 > + store16 v29.d, v30.d, v31.d, v24.d, x4 > + > + add x4, sp, #(\offset + 64) > + ld1 {v16.4s-v19.4s}, [x4] > + butterfly16 v16.4s, v25.4s, v17.4s, v26.4s, v18.4s, v27.4s, v19.4s, v28.4s > + scale v29, v30, v31, v20, v20.4s, v16.4s, v25.4s, v17.4s, v26.4s, v18.4s, v27.4s, v19.4s, \shift > + transpose16_4x4_2 v29, v30, v31, v20 > + > + add x1, x6, #8 > + add x3, x6, #(16 + 3 * 32) > + mov x2, #32 > + mov x4, #-32 > + store16 v29.d, v30.d, v31.d, v20.d, x4 > + > + ret > +endfunc > +.endm > + > +.macro idct_16x16 bitdepth > +function ff_hevc_idct_16x16_\bitdepth\()_neon, export=1 > +//r0 - coeffs > + mov x15, lr Please use "x30" instead of "lr"; slightly older binutils versions (e.g. 2.28) don't support the "lr" register name alias. That's the only hard blocker I have for this patch, other than that, the rest is things that could be nicer. You had a different version of this patch of your own - did that get better or worse performance than this? I believe your code was based on the same arm code (but more loosely) - if the performance difference is significant, can we gain more by adopting changes from your version on top of this later? // Martin
diff --git a/libavcodec/aarch64/Makefile b/libavcodec/aarch64/Makefile index f6434e40da..2ea1d74a38 100644 --- a/libavcodec/aarch64/Makefile +++ b/libavcodec/aarch64/Makefile @@ -61,3 +61,5 @@ NEON-OBJS-$(CONFIG_VP9_DECODER) += aarch64/vp9itxfm_16bpp_neon.o \ aarch64/vp9lpf_neon.o \ aarch64/vp9mc_16bpp_neon.o \ aarch64/vp9mc_neon.o +NEON-OBJS-$(CONFIG_HEVC_DECODER) += aarch64/hevcdsp_idct_neon.o \ + aarch64/hevcdsp_init_aarch64.o diff --git a/libavcodec/aarch64/hevcdsp_idct_neon.S b/libavcodec/aarch64/hevcdsp_idct_neon.S new file mode 100644 index 0000000000..c70d6a906d --- /dev/null +++ b/libavcodec/aarch64/hevcdsp_idct_neon.S @@ -0,0 +1,380 @@ +/* + * ARM NEON optimised IDCT functions for HEVC decoding + * Copyright (c) 2014 Seppo Tomperi <seppo.tomperi@vtt.fi> + * Copyright (c) 2017 Alexandra Hájková + * + * Ported from arm/hevcdsp_idct_neon.S by + * Copyright (c) 2020 Reimar Döffinger + * + * This file is part of FFmpeg. + * + * FFmpeg is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; either + * version 2.1 of the License, or (at your option) any later version. + * + * FFmpeg is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + * + * You should have received a copy of the GNU Lesser General Public + * License along with FFmpeg; if not, write to the Free Software + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA + */ + +#include "libavutil/aarch64/asm.S" + +const trans, align=4 + .short 64, 83, 64, 36 + .short 89, 75, 50, 18 + .short 90, 87, 80, 70 + .short 57, 43, 25, 9 + .short 90, 90, 88, 85 + .short 82, 78, 73, 67 + .short 61, 54, 46, 38 + .short 31, 22, 13, 4 +endconst + +.macro sum_sub out, in, c, op, p + .ifc \op, + + smlal\p \out, \in, \c + .else + smlsl\p \out, \in, \c + .endif +.endm + +.macro fixsqrshrn d, dt, n, m + .ifc \dt, .8h + sqrshrn2 \d\dt, \n\().4s, \m + .else + sqrshrn \n\().4h, \n\().4s, \m + mov \d\().d[0], \n\().d[0] + .endif +.endm + +// uses and clobbers v28-v31 as temp registers +.macro tr_4x4_8 in0, in1, in2, in3, out0, out1, out2, out3, p1, p2 + sshll\p1 v28.4s, \in0, #6 + mov v29.16b, v28.16b + smull\p1 v30.4s, \in1, v0.h[1] + smull\p1 v31.4s, \in1, v0.h[3] + smlal\p2 v28.4s, \in2, v0.h[0] //e0 + smlsl\p2 v29.4s, \in2, v0.h[0] //e1 + smlal\p2 v30.4s, \in3, v0.h[3] //o0 + smlsl\p2 v31.4s, \in3, v0.h[1] //o1 + + add \out0, v28.4s, v30.4s + add \out1, v29.4s, v31.4s + sub \out2, v29.4s, v31.4s + sub \out3, v28.4s, v30.4s +.endm + +.macro transpose8_4x4 r0, r1, r2, r3 + trn1 v2.8h, \r0\().8h, \r1\().8h + trn2 v3.8h, \r0\().8h, \r1\().8h + trn1 v4.8h, \r2\().8h, \r3\().8h + trn2 v5.8h, \r2\().8h, \r3\().8h + trn1 \r0\().4s, v2.4s, v4.4s + trn2 \r2\().4s, v2.4s, v4.4s + trn1 \r1\().4s, v3.4s, v5.4s + trn2 \r3\().4s, v3.4s, v5.4s +.endm + +.macro transpose_8x8 r0, r1, r2, r3, r4, r5, r6, r7 + transpose8_4x4 \r0, \r1, \r2, \r3 + transpose8_4x4 \r4, \r5, \r6, \r7 +.endm + +.macro tr_8x4 shift, in0,in0t, in1,in1t, in2,in2t, in3,in3t, in4,in4t, in5,in5t, in6,in6t, in7,in7t, p1, p2 + tr_4x4_8 \in0\in0t, \in2\in2t, \in4\in4t, \in6\in6t, v24.4s, v25.4s, v26.4s, v27.4s, \p1, \p2 + + smull\p1 v30.4s, \in1\in1t, v0.h[6] + smull\p1 v28.4s, \in1\in1t, v0.h[4] + smull\p1 v29.4s, \in1\in1t, v0.h[5] + sum_sub v30.4s, \in3\in3t, v0.h[4], -, \p1 + sum_sub v28.4s, \in3\in3t, v0.h[5], +, \p1 + sum_sub v29.4s, \in3\in3t, v0.h[7], -, \p1 + + sum_sub v30.4s, \in5\in5t, v0.h[7], +, \p2 + sum_sub v28.4s, \in5\in5t, v0.h[6], +, \p2 + sum_sub v29.4s, \in5\in5t, v0.h[4], -, \p2 + + sum_sub v30.4s, \in7\in7t, v0.h[5], +, \p2 + sum_sub v28.4s, \in7\in7t, v0.h[7], +, \p2 + sum_sub v29.4s, \in7\in7t, v0.h[6], -, \p2 + + add v31.4s, v26.4s, v30.4s + sub v26.4s, v26.4s, v30.4s + fixsqrshrn \in2,\in2t, v31, \shift + + + smull\p1 v31.4s, \in1\in1t, v0.h[7] + sum_sub v31.4s, \in3\in3t, v0.h[6], -, \p1 + sum_sub v31.4s, \in5\in5t, v0.h[5], +, \p2 + sum_sub v31.4s, \in7\in7t, v0.h[4], -, \p2 + fixsqrshrn \in5,\in5t, v26, \shift + + + add v26.4s, v24.4s, v28.4s + sub v24.4s, v24.4s, v28.4s + add v28.4s, v25.4s, v29.4s + sub v25.4s, v25.4s, v29.4s + add v30.4s, v27.4s, v31.4s + sub v27.4s, v27.4s, v31.4s + + fixsqrshrn \in0,\in0t, v26, \shift + fixsqrshrn \in7,\in7t, v24, \shift + fixsqrshrn \in1,\in1t, v28, \shift + fixsqrshrn \in6,\in6t, v25, \shift + fixsqrshrn \in3,\in3t, v30, \shift + fixsqrshrn \in4,\in4t, v27, \shift +.endm + +.macro idct_8x8 bitdepth +function ff_hevc_idct_8x8_\bitdepth\()_neon, export=1 +//x0 - coeffs + mov x1, x0 + ld1 {v16.8h-v19.8h}, [x1], #64 + ld1 {v20.8h-v23.8h}, [x1] + + movrel x1, trans + ld1 {v0.8h}, [x1] + + tr_8x4 7, v16,.4h, v17,.4h, v18,.4h, v19,.4h, v20,.4h, v21,.4h, v22,.4h, v23,.4h + tr_8x4 7, v16,.8h, v17,.8h, v18,.8h, v19,.8h, v20,.8h, v21,.8h, v22,.8h, v23,.8h, 2, 2 + + transpose_8x8 v16, v17, v18, v19, v20, v21, v22, v23 + + tr_8x4 20 - \bitdepth, v16,.4h, v17,.4h, v18,.4h, v19,.4h, v16,.8h, v17,.8h, v18,.8h, v19,.8h, , 2 + tr_8x4 20 - \bitdepth, v20,.4h, v21,.4h, v22,.4h, v23,.4h, v20,.8h, v21,.8h, v22,.8h, v23,.8h, , 2 + + transpose_8x8 v16, v17, v18, v19, v20, v21, v22, v23 + + mov x1, x0 + st1 {v16.8h-v19.8h}, [x1], #64 + st1 {v20.8h-v23.8h}, [x1] + + ret +endfunc +.endm + +.macro butterfly e, o, tmp_p, tmp_m + add \tmp_p, \e, \o + sub \tmp_m, \e, \o +.endm + +.macro tr16_8x4 in0, in1, in2, in3, offset + tr_4x4_8 \in0\().4h, \in1\().4h, \in2\().4h, \in3\().4h, v24.4s, v25.4s, v26.4s, v27.4s + + smull2 v28.4s, \in0\().8h, v0.h[4] + smull2 v29.4s, \in0\().8h, v0.h[5] + smull2 v30.4s, \in0\().8h, v0.h[6] + smull2 v31.4s, \in0\().8h, v0.h[7] + sum_sub v28.4s, \in1\().8h, v0.h[5], +, 2 + sum_sub v29.4s, \in1\().8h, v0.h[7], -, 2 + sum_sub v30.4s, \in1\().8h, v0.h[4], -, 2 + sum_sub v31.4s, \in1\().8h, v0.h[6], -, 2 + + sum_sub v28.4s, \in2\().8h, v0.h[6], +, 2 + sum_sub v29.4s, \in2\().8h, v0.h[4], -, 2 + sum_sub v30.4s, \in2\().8h, v0.h[7], +, 2 + sum_sub v31.4s, \in2\().8h, v0.h[5], +, 2 + + sum_sub v28.4s, \in3\().8h, v0.h[7], +, 2 + sum_sub v29.4s, \in3\().8h, v0.h[6], -, 2 + sum_sub v30.4s, \in3\().8h, v0.h[5], +, 2 + sum_sub v31.4s, \in3\().8h, v0.h[4], -, 2 + + butterfly v24.4s, v28.4s, v16.4s, v23.4s + butterfly v25.4s, v29.4s, v17.4s, v22.4s + butterfly v26.4s, v30.4s, v18.4s, v21.4s + butterfly v27.4s, v31.4s, v19.4s, v20.4s + add x4, sp, #\offset + st1 {v16.4s-v19.4s}, [x4], #64 + st1 {v20.4s-v23.4s}, [x4] +.endm + +.macro load16 in0, in1, in2, in3 + ld1 {\in0}[0], [x1], x2 + ld1 {\in0}[1], [x3], x2 + ld1 {\in1}[0], [x1], x2 + ld1 {\in1}[1], [x3], x2 + ld1 {\in2}[0], [x1], x2 + ld1 {\in2}[1], [x3], x2 + ld1 {\in3}[0], [x1], x2 + ld1 {\in3}[1], [x3], x2 +.endm + +.macro add_member in, t0, t1, t2, t3, t4, t5, t6, t7, op0, op1, op2, op3, op4, op5, op6, op7, p + sum_sub v21.4s, \in, \t0, \op0, \p + sum_sub v22.4s, \in, \t1, \op1, \p + sum_sub v23.4s, \in, \t2, \op2, \p + sum_sub v24.4s, \in, \t3, \op3, \p + sum_sub v25.4s, \in, \t4, \op4, \p + sum_sub v26.4s, \in, \t5, \op5, \p + sum_sub v27.4s, \in, \t6, \op6, \p + sum_sub v28.4s, \in, \t7, \op7, \p +.endm + +.macro butterfly16 in0, in1, in2, in3, in4, in5, in6, in7 + add v20.4s, \in0, \in1 + sub \in0, \in0, \in1 + add \in1, \in2, \in3 + sub \in2, \in2, \in3 + add \in3, \in4, \in5 + sub \in4, \in4, \in5 + add \in5, \in6, \in7 + sub \in6, \in6, \in7 +.endm + +.macro store16 in0, in1, in2, in3, rx + st1 {\in0}[0], [x1], x2 + st1 {\in0}[1], [x3], \rx + st1 {\in1}[0], [x1], x2 + st1 {\in1}[1], [x3], \rx + st1 {\in2}[0], [x1], x2 + st1 {\in2}[1], [x3], \rx + st1 {\in3}[0], [x1], x2 + st1 {\in3}[1], [x3], \rx +.endm + +.macro scale out0, out1, out2, out3, in0, in1, in2, in3, in4, in5, in6, in7, shift + sqrshrn \out0\().4h, \in0, \shift + sqrshrn2 \out0\().8h, \in1, \shift + sqrshrn \out1\().4h, \in2, \shift + sqrshrn2 \out1\().8h, \in3, \shift + sqrshrn \out2\().4h, \in4, \shift + sqrshrn2 \out2\().8h, \in5, \shift + sqrshrn \out3\().4h, \in6, \shift + sqrshrn2 \out3\().8h, \in7, \shift +.endm + +.macro transpose16_4x4_2 r0, r1, r2, r3 + // lower halves + trn1 v2.4h, \r0\().4h, \r1\().4h + trn2 v3.4h, \r0\().4h, \r1\().4h + trn1 v4.4h, \r2\().4h, \r3\().4h + trn2 v5.4h, \r2\().4h, \r3\().4h + trn1 v6.2s, v2.2s, v4.2s + trn2 v7.2s, v2.2s, v4.2s + trn1 v2.2s, v3.2s, v5.2s + trn2 v4.2s, v3.2s, v5.2s + mov \r0\().d[0], v6.d[0] + mov \r2\().d[0], v7.d[0] + mov \r1\().d[0], v2.d[0] + mov \r3\().d[0], v4.d[0] + + // upper halves in reverse order + trn1 v2.8h, \r3\().8h, \r2\().8h + trn2 v3.8h, \r3\().8h, \r2\().8h + trn1 v4.8h, \r1\().8h, \r0\().8h + trn2 v5.8h, \r1\().8h, \r0\().8h + trn1 v6.4s, v2.4s, v4.4s + trn2 v7.4s, v2.4s, v4.4s + trn1 v2.4s, v3.4s, v5.4s + trn2 v4.4s, v3.4s, v5.4s + mov \r3\().d[1], v6.d[1] + mov \r1\().d[1], v7.d[1] + mov \r2\().d[1], v2.d[1] + mov \r0\().d[1], v4.d[1] +.endm + +.macro tr_16x4 name, shift, offset, step +function func_tr_16x4_\name + mov x1, x5 + add x3, x5, #(\step * 64) + mov x2, #(\step * 128) + load16 v16.d, v17.d, v18.d, v19.d + movrel x1, trans + ld1 {v0.8h}, [x1] + + tr16_8x4 v16, v17, v18, v19, \offset + + add x1, x5, #(\step * 32) + add x3, x5, #(\step * 3 *32) + mov x2, #(\step * 128) + load16 v20.d, v17.d, v18.d, v19.d + movrel x1, trans, 16 + ld1 {v1.8h}, [x1] + smull v21.4s, v20.4h, v1.h[0] + smull v22.4s, v20.4h, v1.h[1] + smull v23.4s, v20.4h, v1.h[2] + smull v24.4s, v20.4h, v1.h[3] + smull v25.4s, v20.4h, v1.h[4] + smull v26.4s, v20.4h, v1.h[5] + smull v27.4s, v20.4h, v1.h[6] + smull v28.4s, v20.4h, v1.h[7] + + add_member v20.8h, v1.h[1], v1.h[4], v1.h[7], v1.h[5], v1.h[2], v1.h[0], v1.h[3], v1.h[6], +, +, +, -, -, -, -, -, 2 + add_member v17.4h, v1.h[2], v1.h[7], v1.h[3], v1.h[1], v1.h[6], v1.h[4], v1.h[0], v1.h[5], +, +, -, -, -, +, +, + + add_member v17.8h, v1.h[3], v1.h[5], v1.h[1], v1.h[7], v1.h[0], v1.h[6], v1.h[2], v1.h[4], +, -, -, +, +, +, -, -, 2 + add_member v18.4h, v1.h[4], v1.h[2], v1.h[6], v1.h[0], v1.h[7], v1.h[1], v1.h[5], v1.h[3], +, -, -, +, -, -, +, + + add_member v18.8h, v1.h[5], v1.h[0], v1.h[4], v1.h[6], v1.h[1], v1.h[3], v1.h[7], v1.h[2], +, -, +, +, -, +, +, -, 2 + add_member v19.4h, v1.h[6], v1.h[3], v1.h[0], v1.h[2], v1.h[5], v1.h[7], v1.h[4], v1.h[1], +, -, +, -, +, +, -, + + add_member v19.8h, v1.h[7], v1.h[6], v1.h[5], v1.h[4], v1.h[3], v1.h[2], v1.h[1], v1.h[0], +, -, +, -, +, -, +, -, 2 + + add x4, sp, #\offset + ld1 {v16.4s-v19.4s}, [x4], #64 + + butterfly16 v16.4s, v21.4s, v17.4s, v22.4s, v18.4s, v23.4s, v19.4s, v24.4s + scale v29, v30, v31, v24, v20.4s, v16.4s, v21.4s, v17.4s, v22.4s, v18.4s, v23.4s, v19.4s, \shift + transpose16_4x4_2 v29, v30, v31, v24 + mov x1, x6 + add x3, x6, #(24 +3*32) + mov x2, #32 + mov x4, #-32 + store16 v29.d, v30.d, v31.d, v24.d, x4 + + add x4, sp, #(\offset + 64) + ld1 {v16.4s-v19.4s}, [x4] + butterfly16 v16.4s, v25.4s, v17.4s, v26.4s, v18.4s, v27.4s, v19.4s, v28.4s + scale v29, v30, v31, v20, v20.4s, v16.4s, v25.4s, v17.4s, v26.4s, v18.4s, v27.4s, v19.4s, \shift + transpose16_4x4_2 v29, v30, v31, v20 + + add x1, x6, #8 + add x3, x6, #(16 + 3 * 32) + mov x2, #32 + mov x4, #-32 + store16 v29.d, v30.d, v31.d, v20.d, x4 + + ret +endfunc +.endm + +.macro idct_16x16 bitdepth +function ff_hevc_idct_16x16_\bitdepth\()_neon, export=1 +//r0 - coeffs + mov x15, lr + + // allocate a temp buffer + sub sp, sp, #640 + +.irp i, 0, 1, 2, 3 + add x5, x0, #(8 * \i) + add x6, sp, #(8 * \i * 16) + bl func_tr_16x4_firstpass +.endr + +.irp i, 0, 1, 2, 3 + add x5, sp, #(8 * \i) + add x6, x0, #(8 * \i * 16) + bl func_tr_16x4_secondpass_\bitdepth +.endr + + add sp, sp, #640 + + mov lr, x15 + ret +endfunc +.endm + +idct_8x8 8 +idct_8x8 10 + +tr_16x4 firstpass, 7, 512, 1 +tr_16x4 secondpass_8, 20 - 8, 512, 1 +tr_16x4 secondpass_10, 20 - 10, 512, 1 + +idct_16x16 8 +idct_16x16 10 diff --git a/libavcodec/aarch64/hevcdsp_init_aarch64.c b/libavcodec/aarch64/hevcdsp_init_aarch64.c new file mode 100644 index 0000000000..19d9a7f9ed --- /dev/null +++ b/libavcodec/aarch64/hevcdsp_init_aarch64.c @@ -0,0 +1,45 @@ +/* + * Copyright (c) 2020 Reimar Döffinger + * + * This file is part of FFmpeg. + * + * FFmpeg is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; either + * version 2.1 of the License, or (at your option) any later version. + * + * FFmpeg is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + * + * You should have received a copy of the GNU Lesser General Public + * License along with FFmpeg; if not, write to the Free Software + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA + */ + +#include <stdint.h> + +#include "libavutil/attributes.h" +#include "libavutil/cpu.h" +#include "libavutil/aarch64/cpu.h" +#include "libavcodec/hevcdsp.h" + +void ff_hevc_idct_8x8_8_neon(int16_t *coeffs, int col_limit); +void ff_hevc_idct_8x8_10_neon(int16_t *coeffs, int col_limit); +void ff_hevc_idct_16x16_8_neon(int16_t *coeffs, int col_limit); +void ff_hevc_idct_16x16_10_neon(int16_t *coeffs, int col_limit); + +av_cold void ff_hevc_dsp_init_aarch64(HEVCDSPContext *c, const int bit_depth) +{ + if (!have_neon(av_get_cpu_flags())) return; + + if (bit_depth == 8) { + c->idct[1] = ff_hevc_idct_8x8_8_neon; + c->idct[2] = ff_hevc_idct_16x16_8_neon; + } + if (bit_depth == 10) { + c->idct[1] = ff_hevc_idct_8x8_10_neon; + c->idct[2] = ff_hevc_idct_16x16_10_neon; + } +} diff --git a/libavcodec/hevcdsp.c b/libavcodec/hevcdsp.c index 957e40d5ff..fe272ac1ce 100644 --- a/libavcodec/hevcdsp.c +++ b/libavcodec/hevcdsp.c @@ -257,6 +257,8 @@ int i = 0; break; } + if (ARCH_AARCH64) + ff_hevc_dsp_init_aarch64(hevcdsp, bit_depth); if (ARCH_ARM) ff_hevc_dsp_init_arm(hevcdsp, bit_depth); if (ARCH_PPC) diff --git a/libavcodec/hevcdsp.h b/libavcodec/hevcdsp.h index c605a343d6..0e013a8328 100644 --- a/libavcodec/hevcdsp.h +++ b/libavcodec/hevcdsp.h @@ -129,6 +129,7 @@ void ff_hevc_dsp_init(HEVCDSPContext *hpc, int bit_depth); extern const int8_t ff_hevc_epel_filters[7][4]; extern const int8_t ff_hevc_qpel_filters[3][16]; +void ff_hevc_dsp_init_aarch64(HEVCDSPContext *c, const int bit_depth); void ff_hevc_dsp_init_arm(HEVCDSPContext *c, const int bit_depth); void ff_hevc_dsp_init_ppc(HEVCDSPContext *c, const int bit_depth); void ff_hevc_dsp_init_x86(HEVCDSPContext *c, const int bit_depth);