mirror of
https://github.com/FFmpeg/FFmpeg.git
synced 2024-12-18 03:19:31 +02:00
c0683dce89
* commit '0b9a237b2386ff84a6f99716bd58fa27a1b767e7': hevc: Add NEON 4x4 and 8x8 IDCT [15:12:59] <@ubitux> hevc_idct_4x4_8_c: 389.1 [15:13:00] <@ubitux> hevc_idct_4x4_8_neon: 126.6 [15:13:02] <@ubitux> our ^ [15:13:06] <@ubitux> hevc_idct_4x4_8_c: 389.3 [15:13:08] <@ubitux> hevc_idct_4x4_8_neon: 107.8 [15:13:10] <@ubitux> hevc_idct_4x4_10_c: 418.6 [15:13:12] <@ubitux> hevc_idct_4x4_10_neon: 108.1 [15:13:14] <@ubitux> libav ^ [15:13:30] <@ubitux> so yeah, we can probably trash our versions here Merged-by: James Almer <jamrial@gmail.com>
436 lines
13 KiB
ArmAsm
436 lines
13 KiB
ArmAsm
/*
|
|
* ARM NEON optimised IDCT functions for HEVC decoding
|
|
* Copyright (c) 2014 Seppo Tomperi <seppo.tomperi@vtt.fi>
|
|
* Copyright (c) 2017 Alexandra Hájková
|
|
*
|
|
* 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/arm/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
|
|
endconst
|
|
|
|
function ff_hevc_idct_4x4_dc_neon_8, export=1
|
|
ldrsh r1, [r0]
|
|
ldr r2, =0x20
|
|
add r1, #1
|
|
asr r1, #1
|
|
add r1, r2
|
|
asr r1, #6
|
|
vdup.16 q0, r1
|
|
vdup.16 q1, r1
|
|
vst1.16 {q0, q1}, [r0]
|
|
bx lr
|
|
endfunc
|
|
|
|
function ff_hevc_idct_8x8_dc_neon_8, export=1
|
|
ldrsh r1, [r0]
|
|
ldr r2, =0x20
|
|
add r1, #1
|
|
asr r1, #1
|
|
add r1, r2
|
|
asr r1, #6
|
|
vdup.16 q8, r1
|
|
vdup.16 q9, r1
|
|
vmov.16 q10, q8
|
|
vmov.16 q11, q8
|
|
vmov.16 q12, q8
|
|
vmov.16 q13, q8
|
|
vmov.16 q14, q8
|
|
vmov.16 q15, q8
|
|
vstm r0, {q8-q15}
|
|
bx lr
|
|
endfunc
|
|
|
|
function ff_hevc_idct_16x16_dc_neon_8, export=1
|
|
ldrsh r1, [r0]
|
|
ldr r2, =0x20
|
|
add r1, #1
|
|
asr r1, #1
|
|
add r1, r2
|
|
asr r1, #6
|
|
vdup.16 q8, r1
|
|
vdup.16 q9, r1
|
|
vmov.16 q10, q8
|
|
vmov.16 q11, q8
|
|
vmov.16 q12, q8
|
|
vmov.16 q13, q8
|
|
vmov.16 q14, q8
|
|
vmov.16 q15, q8
|
|
vstm r0!, {q8-q15}
|
|
vstm r0!, {q8-q15}
|
|
vstm r0!, {q8-q15}
|
|
vstm r0, {q8-q15}
|
|
bx lr
|
|
endfunc
|
|
|
|
function ff_hevc_idct_32x32_dc_neon_8, export=1
|
|
ldrsh r1, [r0]
|
|
ldr r2, =0x20
|
|
add r1, #1
|
|
asr r1, #1
|
|
add r1, r2
|
|
asr r1, #6
|
|
mov r3, #16
|
|
vdup.16 q8, r1
|
|
vdup.16 q9, r1
|
|
vmov.16 q10, q8
|
|
vmov.16 q11, q8
|
|
vmov.16 q12, q8
|
|
vmov.16 q13, q8
|
|
vmov.16 q14, q8
|
|
vmov.16 q15, q8
|
|
1: subs r3, #1
|
|
vstm r0!, {q8-q15}
|
|
bne 1b
|
|
bx lr
|
|
endfunc
|
|
|
|
function ff_hevc_add_residual_4x4_neon_8, export=1
|
|
vldm r1, {q0-q1}
|
|
vld1.32 d4[0], [r0], r2
|
|
vld1.32 d4[1], [r0], r2
|
|
vld1.32 d5[0], [r0], r2
|
|
vld1.32 d5[1], [r0], r2
|
|
sub r0, r0, r2, lsl #2
|
|
vmovl.u8 q8, d4
|
|
vmovl.u8 q9, d5
|
|
vqadd.s16 q0, q0, q8
|
|
vqadd.s16 q1, q1, q9
|
|
vqmovun.s16 d0, q0
|
|
vqmovun.s16 d1, q1
|
|
vst1.32 d0[0], [r0], r2
|
|
vst1.32 d0[1], [r0], r2
|
|
vst1.32 d1[0], [r0], r2
|
|
vst1.32 d1[1], [r0], r2
|
|
bx lr
|
|
endfunc
|
|
|
|
function ff_hevc_add_residual_8x8_neon_8, export=1
|
|
mov r3, #8
|
|
1: subs r3, #1
|
|
vld1.16 {q0}, [r1]!
|
|
vld1.8 d16, [r0]
|
|
vmovl.u8 q8, d16
|
|
vqadd.s16 q0, q8
|
|
vqmovun.s16 d0, q0
|
|
vst1.32 d0, [r0], r2
|
|
bne 1b
|
|
bx lr
|
|
endfunc
|
|
|
|
function ff_hevc_add_residual_16x16_neon_8, export=1
|
|
mov r3, #16
|
|
1: subs r3, #1
|
|
vld1.16 {q0, q1}, [r1]!
|
|
vld1.8 {q8}, [r0]
|
|
vmovl.u8 q9, d16
|
|
vmovl.u8 q10, d17
|
|
vqadd.s16 q0, q9
|
|
vqadd.s16 q1, q10
|
|
vqmovun.s16 d0, q0
|
|
vqmovun.s16 d1, q1
|
|
vst1.8 {q0}, [r0], r2
|
|
bne 1b
|
|
bx lr
|
|
endfunc
|
|
|
|
function ff_hevc_add_residual_32x32_neon_8, export=1
|
|
mov r3, #32
|
|
1: subs r3, #1
|
|
vldm r1!, {q0-q3}
|
|
vld1.8 {q8, q9}, [r0]
|
|
vmovl.u8 q10, d16
|
|
vmovl.u8 q11, d17
|
|
vmovl.u8 q12, d18
|
|
vmovl.u8 q13, d19
|
|
vqadd.s16 q0, q10
|
|
vqadd.s16 q1, q11
|
|
vqadd.s16 q2, q12
|
|
vqadd.s16 q3, q13
|
|
vqmovun.s16 d0, q0
|
|
vqmovun.s16 d1, q1
|
|
vqmovun.s16 d2, q2
|
|
vqmovun.s16 d3, q3
|
|
vst1.8 {q0, q1}, [r0], r2
|
|
bne 1b
|
|
bx lr
|
|
endfunc
|
|
|
|
/* uses registers q2 - q9 for temp values */
|
|
/* TODO: reorder */
|
|
.macro tr4_luma_shift r0, r1, r2, r3, shift
|
|
vaddl.s16 q5, \r0, \r2 // c0 = src0 + src2
|
|
vaddl.s16 q2, \r2, \r3 // c1 = src2 + src3
|
|
vsubl.s16 q4, \r0, \r3 // c2 = src0 - src3
|
|
vmull.s16 q6, \r1, d0[0] // c3 = 74 * src1
|
|
|
|
vaddl.s16 q7, \r0, \r3 // src0 + src3
|
|
vsubw.s16 q7, q7, \r2 // src0 - src2 + src3
|
|
vmul.s32 q7, q7, d0[0] // dst2 = 74 * (src0 - src2 + src3)
|
|
|
|
vmul.s32 q8, q5, d0[1] // 29 * c0
|
|
vmul.s32 q9, q2, d1[0] // 55 * c1
|
|
vadd.s32 q8, q9 // 29 * c0 + 55 * c1
|
|
vadd.s32 q8, q6 // dst0 = 29 * c0 + 55 * c1 + c3
|
|
|
|
vmul.s32 q2, q2, d0[1] // 29 * c1
|
|
vmul.s32 q9, q4, d1[0] // 55 * c2
|
|
vsub.s32 q9, q2 // 55 * c2 - 29 * c1
|
|
vadd.s32 q9, q6 // dst1 = 55 * c2 - 29 * c1 + c3
|
|
|
|
vmul.s32 q5, q5, d1[0] // 55 * c0
|
|
vmul.s32 q4, q4, d0[1] // 29 * c2
|
|
vadd.s32 q5, q4 // 55 * c0 + 29 * c2
|
|
vsub.s32 q5, q6 // dst3 = 55 * c0 + 29 * c2 - c3
|
|
|
|
vqrshrn.s32 \r0, q8, \shift
|
|
vqrshrn.s32 \r1, q9, \shift
|
|
vqrshrn.s32 \r2, q7, \shift
|
|
vqrshrn.s32 \r3, q5, \shift
|
|
.endm
|
|
|
|
function ff_hevc_transform_luma_4x4_neon_8, export=1
|
|
vpush {d8-d15}
|
|
vld1.16 {q14, q15}, [r0] // coeffs
|
|
ldr r3, =0x4a // 74
|
|
vmov.32 d0[0], r3
|
|
ldr r3, =0x1d // 29
|
|
vmov.32 d0[1], r3
|
|
ldr r3, =0x37 // 55
|
|
vmov.32 d1[0], r3
|
|
|
|
tr4_luma_shift d28, d29, d30, d31, #7
|
|
|
|
vtrn.16 d28, d29
|
|
vtrn.16 d30, d31
|
|
vtrn.32 q14, q15
|
|
|
|
tr4_luma_shift d28, d29, d30, d31, #12
|
|
|
|
vtrn.16 d28, d29
|
|
vtrn.16 d30, d31
|
|
vtrn.32 q14, q15
|
|
vst1.16 {q14, q15}, [r0]
|
|
vpop {d8-d15}
|
|
bx lr
|
|
endfunc
|
|
|
|
.macro sum_sub out, in, c, op
|
|
.ifc \op, +
|
|
vmlal.s16 \out, \in, \c
|
|
.else
|
|
vmlsl.s16 \out, \in, \c
|
|
.endif
|
|
.endm
|
|
|
|
.macro tr_4x4 in0, in1, in2, in3, out0, out1, out2, out3, shift, tmp0, tmp1, tmp2, tmp3, tmp4
|
|
vshll.s16 \tmp0, \in0, #6
|
|
vmull.s16 \tmp2, \in1, d4[1]
|
|
vmov \tmp1, \tmp0
|
|
vmull.s16 \tmp3, \in1, d4[3]
|
|
vmlal.s16 \tmp0, \in2, d4[0] @e0
|
|
vmlsl.s16 \tmp1, \in2, d4[0] @e1
|
|
vmlal.s16 \tmp2, \in3, d4[3] @o0
|
|
vmlsl.s16 \tmp3, \in3, d4[1] @o1
|
|
|
|
vadd.s32 \tmp4, \tmp0, \tmp2
|
|
vsub.s32 \tmp0, \tmp0, \tmp2
|
|
vadd.s32 \tmp2, \tmp1, \tmp3
|
|
vsub.s32 \tmp1, \tmp1, \tmp3
|
|
vqrshrn.s32 \out0, \tmp4, #\shift
|
|
vqrshrn.s32 \out3, \tmp0, #\shift
|
|
vqrshrn.s32 \out1, \tmp2, #\shift
|
|
vqrshrn.s32 \out2, \tmp1, #\shift
|
|
.endm
|
|
|
|
.macro tr_4x4_8 in0, in1, in2, in3, out0, out1, out2, out3, tmp0, tmp1, tmp2, tmp3
|
|
vshll.s16 \tmp0, \in0, #6
|
|
vld1.s16 {\in0}, [r1, :64]!
|
|
vmov \tmp1, \tmp0
|
|
vmull.s16 \tmp2, \in1, \in0[1]
|
|
vmull.s16 \tmp3, \in1, \in0[3]
|
|
vmlal.s16 \tmp0, \in2, \in0[0] @e0
|
|
vmlsl.s16 \tmp1, \in2, \in0[0] @e1
|
|
vmlal.s16 \tmp2, \in3, \in0[3] @o0
|
|
vmlsl.s16 \tmp3, \in3, \in0[1] @o1
|
|
|
|
vld1.s16 {\in0}, [r1, :64]
|
|
|
|
vadd.s32 \out0, \tmp0, \tmp2
|
|
vadd.s32 \out1, \tmp1, \tmp3
|
|
vsub.s32 \out2, \tmp1, \tmp3
|
|
vsub.s32 \out3, \tmp0, \tmp2
|
|
|
|
sub r1, r1, #8
|
|
.endm
|
|
|
|
@ Do a 4x4 transpose, using q registers for the subtransposes that don't
|
|
@ need to address the indiviudal d registers.
|
|
@ r0,r1 == rq0, r2,r3 == rq1
|
|
.macro transpose_4x4 rq0, rq1, r0, r1, r2, r3
|
|
vtrn.32 \rq0, \rq1
|
|
vtrn.16 \r0, \r1
|
|
vtrn.16 \r2, \r3
|
|
.endm
|
|
|
|
.macro idct_4x4 bitdepth
|
|
function ff_hevc_idct_4x4_\bitdepth\()_neon, export=1
|
|
@r0 - coeffs
|
|
vld1.s16 {q0-q1}, [r0, :128]
|
|
|
|
movrel r1, trans
|
|
vld1.s16 {d4}, [r1, :64]
|
|
|
|
tr_4x4 d0, d1, d2, d3, d16, d17, d18, d19, 7, q10, q11, q12, q13, q0
|
|
transpose_4x4 q8, q9, d16, d17, d18, d19
|
|
|
|
tr_4x4 d16, d17, d18, d19, d0, d1, d2, d3, 20 - \bitdepth, q10, q11, q12, q13, q0
|
|
transpose_4x4 q0, q1, d0, d1, d2, d3
|
|
vst1.s16 {d0-d3}, [r0, :128]
|
|
bx lr
|
|
endfunc
|
|
.endm
|
|
|
|
.macro transpose8_4x4 r0, r1, r2, r3
|
|
vtrn.16 \r0, \r1
|
|
vtrn.16 \r2, \r3
|
|
vtrn.32 \r0, \r2
|
|
vtrn.32 \r1, \r3
|
|
.endm
|
|
|
|
.macro transpose_8x8 r0, r1, r2, r3, r4, r5, r6, r7, l0, l1, l2, l3, l4, l5, l6, l7
|
|
transpose8_4x4 \r0, \r1, \r2, \r3
|
|
transpose8_4x4 \r4, \r5, \r6, \r7
|
|
|
|
transpose8_4x4 \l0, \l1, \l2, \l3
|
|
transpose8_4x4 \l4, \l5, \l6, \l7
|
|
.endm
|
|
|
|
.macro tr_8x4 shift, in0, in1, in2, in3, in4, in5, in6, in7
|
|
tr_4x4_8 \in0, \in2, \in4, \in6, q8, q9, q10, q11, q12, q13, q14, q15
|
|
|
|
vmull.s16 q14, \in1, \in0[2]
|
|
vmull.s16 q12, \in1, \in0[0]
|
|
vmull.s16 q13, \in1, \in0[1]
|
|
sum_sub q14, \in3, \in0[0], -
|
|
sum_sub q12, \in3, \in0[1], +
|
|
sum_sub q13, \in3, \in0[3], -
|
|
|
|
sum_sub q14, \in5, \in0[3], +
|
|
sum_sub q12, \in5, \in0[2], +
|
|
sum_sub q13, \in5, \in0[0], -
|
|
|
|
sum_sub q14, \in7, \in0[1], +
|
|
sum_sub q12, \in7, \in0[3], +
|
|
sum_sub q13, \in7, \in0[2], -
|
|
|
|
vadd.s32 q15, q10, q14
|
|
vsub.s32 q10, q10, q14
|
|
vqrshrn.s32 \in2, q15, \shift
|
|
|
|
vmull.s16 q15, \in1, \in0[3]
|
|
sum_sub q15, \in3, \in0[2], -
|
|
sum_sub q15, \in5, \in0[1], +
|
|
sum_sub q15, \in7, \in0[0], -
|
|
|
|
vqrshrn.s32 \in5, q10, \shift
|
|
|
|
vadd.s32 q10, q8, q12
|
|
vsub.s32 q8, q8, q12
|
|
vadd.s32 q12, q9, q13
|
|
vsub.s32 q9, q9, q13
|
|
vadd.s32 q14, q11, q15
|
|
vsub.s32 q11, q11, q15
|
|
|
|
vqrshrn.s32 \in0, q10, \shift
|
|
vqrshrn.s32 \in7, q8, \shift
|
|
vqrshrn.s32 \in1, q12, \shift
|
|
vqrshrn.s32 \in6, q9, \shift
|
|
vqrshrn.s32 \in3, q14, \shift
|
|
vqrshrn.s32 \in4, q11, \shift
|
|
.endm
|
|
|
|
.macro idct_8x8 bitdepth
|
|
function ff_hevc_idct_8x8_\bitdepth\()_neon, export=1
|
|
@r0 - coeffs
|
|
vpush {q4-q7}
|
|
|
|
mov r1, r0
|
|
mov r2, #64
|
|
add r3, r0, #32
|
|
vld1.s16 {q0-q1}, [r1,:128], r2
|
|
vld1.s16 {q2-q3}, [r3,:128], r2
|
|
vld1.s16 {q4-q5}, [r1,:128], r2
|
|
vld1.s16 {q6-q7}, [r3,:128], r2
|
|
|
|
movrel r1, trans
|
|
|
|
tr_8x4 7, d0, d2, d4, d6, d8, d10, d12, d14
|
|
tr_8x4 7, d1, d3, d5, d7, d9, d11, d13, d15
|
|
|
|
@ Transpose each 4x4 block, and swap how d4-d7 and d8-d11 are used.
|
|
@ Layout before:
|
|
@ d0 d1
|
|
@ d2 d3
|
|
@ d4 d5
|
|
@ d6 d7
|
|
@ d8 d9
|
|
@ d10 d11
|
|
@ d12 d13
|
|
@ d14 d15
|
|
transpose_8x8 d0, d2, d4, d6, d8, d10, d12, d14, d1, d3, d5, d7, d9, d11, d13, d15
|
|
@ Now the layout is:
|
|
@ d0 d8
|
|
@ d2 d10
|
|
@ d4 d12
|
|
@ d6 d14
|
|
@ d1 d9
|
|
@ d3 d11
|
|
@ d5 d13
|
|
@ d7 d15
|
|
|
|
tr_8x4 20 - \bitdepth, d0, d2, d4, d6, d1, d3, d5, d7
|
|
vswp d0, d8
|
|
tr_8x4 20 - \bitdepth, d0, d10, d12, d14, d9, d11, d13, d15
|
|
vswp d0, d8
|
|
|
|
transpose_8x8 d0, d2, d4, d6, d8, d10, d12, d14, d1, d3, d5, d7, d9, d11, d13, d15
|
|
|
|
mov r1, r0
|
|
mov r2, #64
|
|
add r3, r0, #32
|
|
vst1.s16 {q0-q1}, [r1,:128], r2
|
|
vst1.s16 {q2-q3}, [r3,:128], r2
|
|
vst1.s16 {q4-q5}, [r1,:128], r2
|
|
vst1.s16 {q6-q7}, [r3,:128], r2
|
|
|
|
vpop {q4-q7}
|
|
bx lr
|
|
endfunc
|
|
.endm
|
|
|
|
idct_4x4 8
|
|
idct_4x4 10
|
|
idct_8x8 8
|
|
idct_8x8 10
|