From: daichengrong <daichengr...@iscas.ac.cn>

     riscv/hevcdsp_idct_rvv: Optimize idct_32x32_8

     On Banana PI F3:

     hevc_idct_32x32_8_c:                                118945.0 ( 1.00x)
     hevc_idct_32x32_8_rvv_i64:                           28503.7 ( 4.17x)

Signed-off-by: daichengrong <daichengr...@iscas.ac.cn>
---
 libavcodec/riscv/Makefile           |   1 +
 libavcodec/riscv/hevcdsp_idct_rvv.S | 973 ++++++++++++++++++++++++++++
 libavcodec/riscv/hevcdsp_init.c     |  52 +-
 3 files changed, 1006 insertions(+), 20 deletions(-)
 create mode 100644 libavcodec/riscv/hevcdsp_idct_rvv.S

diff --git a/libavcodec/riscv/Makefile b/libavcodec/riscv/Makefile
index a80d2fa2e7..dfc33afbee 100644
--- a/libavcodec/riscv/Makefile
+++ b/libavcodec/riscv/Makefile
@@ -36,6 +36,7 @@ RVV-OBJS-$(CONFIG_H264DSP) += riscv/h264addpx_rvv.o 
riscv/h264dsp_rvv.o \
 OBJS-$(CONFIG_H264QPEL) += riscv/h264qpel_init.o
 RVV-OBJS-$(CONFIG_H264QPEL) += riscv/h264qpel_rvv.o
 OBJS-$(CONFIG_HEVC_DECODER) += riscv/hevcdsp_init.o
+OBJS-$(CONFIG_HEVC_DECODER) += riscv/hevcdsp_idct_rvv.o
 RVV-OBJS-$(CONFIG_HEVC_DECODER)  += riscv/h26x/h2656_inter_rvv.o
 OBJS-$(CONFIG_HUFFYUV_DECODER) += riscv/huffyuvdsp_init.o
 RVV-OBJS-$(CONFIG_HUFFYUV_DECODER) += riscv/huffyuvdsp_rvv.o
diff --git a/libavcodec/riscv/hevcdsp_idct_rvv.S 
b/libavcodec/riscv/hevcdsp_idct_rvv.S
new file mode 100644
index 0000000000..561b8ada47
--- /dev/null
+++ b/libavcodec/riscv/hevcdsp_idct_rvv.S
@@ -0,0 +1,973 @@
+/*
+ * Copyright (c) 2025 Institue of Software Chinese Academy of Sciences (ISCAS).
+ *
+ * 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/riscv/asm.S"
+
+const trans, align=4
+        .2byte          64, 83, 64, 36
+        .2byte          89, 75, 50, 18
+        .2byte          90, 87, 80, 70
+        .2byte          57, 43, 25, 9
+        .2byte          90, 90, 88, 85
+        .2byte          82, 78, 73, 67
+        .2byte          61, 54, 46, 38
+        .2byte          31, 22, 13, 4
+endconst
+
+.macro sum_sub out, in, c, op, p
+        mv t0, \c
+        .ifc \op, -
+                        neg    t0, t0
+        .endif
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        .ifc \p, 2  
+                vslidedown.vi    v8, \in, 4
+                vwmacc.vx    \out, t0, v8
+        .else
+                vwmacc.vx    \out, t0, \in
+        .endif
+.endm
+
+.macro add_member32 in, t0, t1, t2, t3, op0, op1, op2, op3, p
+        sum_sub         v24, \in, \t0, \op0, \p
+        sum_sub         v25, \in, \t1, \op1, \p
+        sum_sub         v26, \in, \t2, \op2, \p
+        sum_sub         v27, \in, \t3, \op3, \p
+.endm
+
+.macro butterfly e, o, tmp_p, tmp_m
+        vsetivli    zero, 4, e32, m1, tu, ma
+        vadd.vv             \tmp_p, \e, \o
+        vsub.vv             \tmp_m, \e, \o
+.endm
+
+.macro butterfly16 in0, in1, in2, in3, in4, in5, in6, in7
+        vsetivli    zero, 4, e32, m1, tu, ma
+        vadd.vv              v20, \in0, \in1
+        vsub.vv              \in0, \in0, \in1
+        vadd.vv              \in1, \in2, \in3
+        vsub.vv              \in2, \in2, \in3
+        vadd.vv              \in3, \in4, \in5
+        vsub.vv              \in4, \in4, \in5
+        vadd.vv              \in5, \in6, \in7
+        vsub.vv              \in6, \in6, \in7
+.endm
+
+.macro multiply in
+        vsetivli    zero, 4, e16, m1, tu, ma
+        vse16.v         \in, (s0)
+        ld      s2, 0*2(s0)
+        ld      s3, 1*2(s0)
+        ld      s4, 2*2(s0)
+        ld      s5, 3*2(s0)
+
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vwmul.vx        v24, v4, s2
+        vwmul.vx        v25, v4, s3
+        vwmul.vx        v26, v4, s4
+        vwmul.vx        v27, v4, s5
+.endm
+
+func tr_block1, zve64x
+        multiply        v0
+
+        addi    sp,sp,-8*16
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        sd    x\i,8*(\i - 10)(sp)
+.endr
+        vsetivli    zero, 4, e16, m1, tu, ma
+        vse16.v         v0, (s0)
+        ld      x10, 0*2(s0)
+        ld      x11, 1*2(s0)
+        ld      x12, 2*2(s0)
+        ld      x13, 3*2(s0)
+        vse16.v         v1, (s0)
+        ld      x14, 0*2(s0)
+        ld      x15, 1*2(s0)
+        ld      x16, 2*2(s0)
+        ld      x17, 3*2(s0)
+        vse16.v         v2, (s0)
+        ld      x18, 0*2(s0)
+        ld      x19, 1*2(s0)
+        ld      x20, 2*2(s0)
+        ld      x21, 3*2(s0)
+        vse16.v         v3, (s0)
+        ld      x22, 0*2(s0)
+        ld      x23, 1*2(s0)
+        ld      x24, 2*2(s0)
+        ld      x25, 3*2(s0)
+        
+        add_member32    v4,  x11, x14, x17, x20, +, +, +, +, 2
+        add_member32    v5,  x12, x17, x22, x24, +, +, +, -
+        add_member32    v5,  x13, x20, x24, x17, +, +, -, -, 2
+        add_member32    v6,  x14, x23, x19, x10, +, +, -, -
+        add_member32    v6,  x15, x25, x14, x16, +, -, -, -, 2
+        add_member32    v7,  x16, x22, x10, x23, +, -, -, -
+        add_member32    v7,  x17, x19, x15, x21, +, -, -, +, 2
+        add_member32    v16, x18, x16, x20, x14, +, -, -, +
+        add_member32    v16, x19, x13, x25, x12, +, -, -, +, 2
+        add_member32    v17, x20, x11, x21, x19, +, -, +, +
+        add_member32    v17, x21, x12, x16, x25, +, -, +, -, 2
+        add_member32    v18, x22, x15, x11, x18, +, -, +, -
+        add_member32    v18, x23, x18, x13, x11, +, -, +, -, 2
+        add_member32    v19, x24, x21, x18, x15, +, -, +, -
+        add_member32    v19, x25, x24, x23, x22, +, -, +, -, 2
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        ld    x\i, 8*(\i - 10)(sp)
+.endr
+        addi    sp, sp, 8*16
+
+        ret
+endfunc
+
+func tr_block2, zve64x
+        multiply        v1
+
+        addi    sp,sp,-8*16
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        sd    x\i,8*(\i - 10)(sp)
+.endr
+        vsetivli    zero, 4, e16, m1, tu, ma
+        vse16.v         v0, (s0)
+        ld      x10, 0*2(s0)
+        ld      x11, 1*2(s0)
+        ld      x12, 2*2(s0)
+        ld      x13, 3*2(s0)
+        vse16.v         v1, (s0)
+        ld      x14, 0*2(s0)
+        ld      x15, 1*2(s0)
+        ld      x16, 2*2(s0)
+        ld      x17, 3*2(s0)
+        vse16.v         v2, (s0)
+        ld      x18, 0*2(s0)
+        ld      x19, 1*2(s0)
+        ld      x20, 2*2(s0)
+        ld      x21, 3*2(s0)
+        vse16.v         v3, (s0)
+        ld      x22, 0*2(s0)
+        ld      x23, 1*2(s0)
+        ld      x24, 2*2(s0)
+        ld      x25, 3*2(s0)
+
+        add_member32    v4,  x23, x25, x22, x19, +, -, -, -, 2
+        add_member32    v5,  x19, x14, x10, x15, -, -, -, -
+        add_member32    v5,  x10, x16, x23, x21, -, -, -, +, 2
+        add_member32    v6,  x18, x24, x15, x13, -, +, +, +
+        add_member32    v6,  x24, x13, x17, x23, +, +, +, -, 2
+        add_member32    v7,  x15, x17, x21, x10, +, +, -, -
+        add_member32    v7,  x13, x23, x11, x25, +, -, -, +, 2
+        add_member32    v16, x22, x12, x24, x11, +, -, -, +
+        add_member32    v16, x20, x18, x14, x24, -, -, +, +, 2
+        add_member32    v17, x11, x22, x18, x12, -, +, +, -
+        add_member32    v17, x17, x11, x20, x22, -, +, -, -, 2
+        add_member32    v18, x25, x19, x12, x14, +, +, -, +
+        add_member32    v18, x16, x21, x25, x20, +, -, -, +, 2
+        add_member32    v19, x12, x11, x13, x16, +, -, +, -
+        add_member32    v19, x21, x20, x19, x18, +, -, +, -, 2
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        ld    x\i, 8*(\i - 10)(sp)
+.endr
+        addi    sp, sp, 8*16
+
+        ret
+endfunc
+
+func tr_block3, zve64x
+        multiply        v2
+        addi    sp,sp,-8*16
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        sd    x\i,8*(\i - 10)(sp)
+.endr
+        vsetivli    zero, 4, e16, m1, tu, ma
+        vse16.v         v0, (s0)
+        ld      x10, 0*2(s0)
+        ld      x11, 1*2(s0)
+        ld      x12, 2*2(s0)
+        ld      x13, 3*2(s0)
+        vse16.v         v1, (s0)
+        ld      x14, 0*2(s0)
+        ld      x15, 1*2(s0)
+        ld      x16, 2*2(s0)
+        ld      x17, 3*2(s0)
+        vse16.v         v2, (s0)
+        ld      x18, 0*2(s0)
+        ld      x19, 1*2(s0)
+        ld      x20, 2*2(s0)
+        ld      x21, 3*2(s0)
+        vse16.v         v3, (s0)
+        ld      x22, 0*2(s0)
+        ld      x23, 1*2(s0)
+        ld      x24, 2*2(s0)
+        ld      x25, 3*2(s0)
+
+        add_member32    v4,  x16, x13, x10, x12, -, -, -, -, 2
+        add_member32    v5,  x20, x25, x21, x16, -, -, +, +
+        add_member32    v5,  x14, x12, x19, x25, +, +, +, -, 2
+        add_member32    v6,  x22, x20, x11, x17, +, -, -, -
+        add_member32    v6,  x12, x18, x22, x10, -, -, +, +, 2
+        add_member32    v7,  x24, x14, x18, x20, -, +, +, -
+        add_member32    v7,  x10, x24, x12, x22, +, +, -, -, 2
+        add_member32    v16, x25, x11, x23, x13, -, -, +, +
+        add_member32    v16, x11, x21, x17, x15, -, +, +, -, 2
+        add_member32    v17, x23, x17, x13, x24, +, +, -, +
+        add_member32    v17, x13, x15, x24, x18, +, -, +, +, 2
+        add_member32    v18, x21, x23, x16, x11, -, -, +, -
+        add_member32    v18, x15, x10, x14, x19, -, +, -, +, 2
+        add_member32    v19, x19, x22, x25, x23, +, -, +, +
+        add_member32    v19, x17, x16, x15, x14, +, -, +, -, 2
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        ld    x\i, 8*(\i - 10)(sp)
+.endr
+        addi    sp, sp, 8*16
+
+        ret
+endfunc
+
+func tr_block4, zve64x
+        multiply        v3
+        addi    sp,sp,-8*16
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        sd    x\i,8*(\i - 10)(sp)
+.endr
+        vsetivli    zero, 4, e16, m1, tu, ma
+        vse16.v         v0, (s0)
+        ld      x10, 0*2(s0)
+        ld      x11, 1*2(s0)
+        ld      x12, 2*2(s0)
+        ld      x13, 3*2(s0)
+        vse16.v         v1, (s0)
+        ld      x14, 0*2(s0)
+        ld      x15, 1*2(s0)
+        ld      x16, 2*2(s0)
+        ld      x17, 3*2(s0)
+        vse16.v         v2, (s0)
+        ld      x18, 0*2(s0)
+        ld      x19, 1*2(s0)
+        ld      x20, 2*2(s0)
+        ld      x21, 3*2(s0)
+        vse16.v         v3, (s0)
+        ld      x22, 0*2(s0)
+        ld      x23, 1*2(s0)
+        ld      x24, 2*2(s0)
+        ld      x25, 3*2(s0)
+
+        add_member32    v4,  x15, x18, x21, x24, -, -, -, -, 2
+        add_member32    v5,  x10, x13, x18, x23, +, +, +, +
+        add_member32    v5,  x18, x10, x15, x22, -, -, -, -, 2
+        add_member32    v6,  x25, x16, x12, x21, +, +, +, +
+        add_member32    v6,  x19, x21, x10, x20, +, -, -, -, 2
+        add_member32    v7,  x12, x25, x13, x19, -, -, +, +
+        add_member32    v7,  x14, x20, x16, x18, +, +, -, -, 2
+        add_member32    v16, x21, x15, x19, x17, -, -, +, +
+        add_member32    v16, x23, x11, x22, x16, -, +, -, -, 2
+        add_member32    v17, x16, x14, x25, x15, +, -, +, +
+        add_member32    v17, x11, x19, x23, x14, -, +, +, -, 2
+        add_member32    v18, x17, x24, x20, x13, +, -, -, +
+        add_member32    v18, x24, x22, x17, x12, -, -, +, -, 2
+        add_member32    v19, x20, x17, x14, x11, -, +, -, +
+        add_member32    v19, x13, x12, x11, x10, +, -, +, -, 2
+
+.irp i, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
+        ld    x\i, 8*(\i - 10)(sp)
+.endr
+        addi    sp, sp, 8*16
+
+        ret
+endfunc
+
+.macro butterfly32 in0, in1, in2, in3, out
+        vsetivli    zero, 4, e32, m1, tu, ma
+        vadd.vv             \out, \in0, \in1
+        vsub.vv             \in0, \in0, \in1
+        vadd.vv             \in1, \in2, \in3
+        vsub.vv             \in2, \in2, \in3
+.endm
+
+.macro load16 in0, in1, in2, in3
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v         v8, (a3)
+        
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx    \in0, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v         \in0, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+     
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx \in1, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v \in1, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+        
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx \in2, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v \in2, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+        
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx \in3, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v \in3, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+.endm
+
+.macro store16 in0, in1, in2, in3, rx
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v \in0, (a1)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1down.vx v8, \in0, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v v8, (a3)
+        add a1, a1, a2
+        add a3, a3, \rx
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v \in1, (a1)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1down.vx v8, \in1, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v v8, (a3)
+        add a1, a1, a2
+        add a3, a3, \rx
+        
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v \in2, (a1)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1down.vx v8, \in2, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v v8, (a3)
+        add a1, a1, a2
+        add a3, a3, \rx
+        
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v \in3, (a1)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1down.vx v8, \in3, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vse64.v v8, (a3)
+        add a1, a1, a2
+        add a3, a3, \rx 
+.endm
+
+.macro load32
+        addi             a1,  a5, 64
+        addi             a3,  a1, 128
+        li               a2,  256
+        
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v4, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v4, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v5, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v5, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v6, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v6, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v7, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v7, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v16, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v16, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v17, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v17, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v18, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v18, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v8, (a3)
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx v19, v8, zero
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vle64.v v19, (a1) 
+        add a1, a1, a2
+        add a3, a3, a2
+
+.endm
+
+.macro add_member in, tt0, tt1, tt2, tt3, tt4, tt5, tt6, tt7, op0, op1, op2, 
op3, op4, op5, op6, op7, p
+        sum_sub         v21, \in, \tt0, \op0, \p
+        sum_sub         v22, \in, \tt1, \op1, \p
+        sum_sub         v23, \in, \tt2, \op2, \p
+        sum_sub         v24, \in, \tt3, \op3, \p
+        sum_sub         v25, \in, \tt4, \op4, \p
+        sum_sub         v26, \in, \tt5, \op5, \p
+        sum_sub         v27, \in, \tt6, \op6, \p
+        sum_sub         v28, \in, \tt7, \op7, \p
+.endm
+
+.macro scale_store shift
+        vsetivli    zero, 8, e16, m1, tu, ma 
+        vle16.v        v28, (a4)
+        addi            a4, a4, 2*8
+        vle16.v        v29, (a4)
+        addi            a4, a4, 2*8
+        vle16.v        v30, (a4)
+        addi            a4, a4, 2*8
+        vle16.v        v31, (a4)
+        addi            a4, a4, 2*8
+
+        butterfly32     v28, v24, v29, v25, v2
+        butterfly32     v30, v26, v31, v27, v3
+
+        scale           v20, v21, v22, v23, v2, v28, v24, v29, v3, v30, v26, 
v31, \shift
+
+        transpose16_4x4_2 20, 21, 22, 23, 24, 25, 26, 27, 28, 29
+
+        store16         v20, v21, v22, v23, t1
+
+        vsetivli    zero, 4, e16, m1, tu, ma 
+        vle16.v         v2, (t2)
+        addi            t2, t2, 8
+        vle16.v         v3, (t2)
+        addi            t2, t2, -8
+.endm
+
+.macro store_to_stack off1, off2, in0, in2, in4, in6, in7, in5, in3, in1
+        li      a7, \off1
+        add     a1, sp, a7
+        li      a7, \off2
+        add     a3, sp, a7  
+        li      a2, -16
+        li      a4, 16
+
+        vsetivli    zero, 4, e32, m1, tu, ma
+        vse32.v         \in0, (a1)
+        add             a1, a1, a4
+        vse32.v         \in1, (a3)
+        add             a3, a3, a2
+        vse32.v         \in2, (a1)
+        add             a1, a1, a4
+        vse32.v         \in3, (a3)
+        add             a3, a3, a2
+        vse32.v         \in4, (a1)
+        add             a1, a1, a4
+        vse32.v         \in5, (a3)
+        add             a3, a3, a2
+        vse32.v         \in6, (a1)
+        vse32.v         \in7, (a3)
+.endm
+
+.macro transpose16_4x4_2 r0, r1, r2, r3, tmp0, tmp1, tmp2, tmp3, tmp4, tmp5
+        vsetivli    zero, 8, e16, m1, tu, ma
+        vse16.v         v0, (s0)
+
+        vsetivli    zero, 8, e16, m1, tu, ma
+        vid.v     v0
+        vand.vi   v8, v0, 1
+        vmsne.vi  v0, v8, 0
+
+        vslideup.vi    v8, v\r1, 1
+        vsetivli    zero, 4, e16, m1, ta, ma
+        vmerge.vvm       v\tmp0\(), v\r0\(), v8, v0
+
+        vslidedown.vi    v8, v\r0, 1
+        vmerge.vvm       v\tmp1\(), v8, v\r1\(), v0
+
+        vslideup.vi    v8, v\r3, 1
+        vmerge.vvm       v\tmp2\(), v\r2\(), v8, v0
+        
+        vslidedown.vi    v8, v\r2\(), 1
+        vmerge.vvm       v\tmp3\(), v8, v\r3\(), v0
+
+        vsetivli    zero, 2, e32, m1, tu
+
+        vslideup.vi    v8, v\tmp2\(), 1
+        vmerge.vvm       v\tmp4\(), v\tmp0\(), v8, v0
+
+        vslidedown.vi    v8, v\tmp0\(), 1
+        vmerge.vvm       v\tmp5\(), v8, v\tmp2\(), v0
+
+        vslideup.vi    v8, v\tmp3\(), 1
+        vmerge.vvm       v\tmp0\(), v\tmp1\(), v8, v0
+
+        vslidedown.vi    v8, v\tmp1\(), 1
+        vmerge.vvm       v\tmp2\(), v8, v\tmp3\(), v0
+
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vmv.v.v  v\r0\(), v\tmp4\()
+        vmv.v.v  v\r2\(), v\tmp5\()
+        vmv.v.v  v\r1\(), v\tmp0\()
+        vmv.v.v  v\r3\(), v\tmp2\()
+
+        vsetivli    zero, 8, e16, m1, tu, ma
+
+        vslideup.vi    v8, v\r2\(), 1
+        vmerge.vvm       v\tmp0\(), v\r3\(), v8, v0
+
+        vslidedown.vi    v8, v\r3\(), 1
+        vmerge.vvm        v\tmp1\(), v8, v\r2\(), v0
+
+        vslideup.vi    v8, v\r0\(), 1
+        vmerge.vvm       v\tmp2\(), v\r1\(),v8, v0
+
+        vslidedown.vi    v8, v\r1\(), 1
+        vmerge.vvm        v\tmp3\(), v8, v\r0\(), v0
+
+        vsetivli    zero, 4, e32, m1, tu, ma
+
+        vslideup.vi    v8, v\tmp2\(), 1
+        vmerge.vvm       v\tmp4\(), v\tmp0\(), v8, v0
+
+        vslidedown.vi    v8, v\tmp0\(), 1
+        vmerge.vvm       v\tmp5\(), v8, v\tmp2\(), v0
+
+        vslideup.vi    v8, v\tmp3\(), 1
+        vmerge.vvm       v\tmp0\(), v\tmp1\(), v8, v0
+
+        vslidedown.vi    v8, v\tmp1\(), 1
+        vmerge.vvm       v\tmp2\(), v8, v\tmp3\(), v0
+
+        vsetivli    zero, 2, e64, m1, tu, ma
+
+        vmerge.vvm       v\r3\(), v\r3\(), v\tmp4\(), v0
+        vmerge.vvm       v\r1\(), v\r1\(), v\tmp5\(), v0
+        vmerge.vvm       v\r2\(), v\r2\(), v\tmp0\(), v0
+        vmerge.vvm       v\r0\(), v\r0\(), v\tmp2\(), v0
+
+        vsetivli    zero, 8, e16, m1, tu, ma
+        vle16.v         v0, (s0)
+.endm
+
+.macro tr16_8x4 in0, in1, in2, in3, offset
+        tr_4x4_8        \in0, \in1, \in2, \in3, v24, v25, v26, v27
+        
+        vsetivli    zero, 8, e16, m1, tu, ma
+        vse16.v     v0, (s0)
+        lh      s6, 4*2(s0)
+        lh      s7, 5*2(s0)
+        lh      s8, 6*2(s0)
+        lh      s9, 7*2(s0)
+        
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vslidedown.vi    v8, \in0, 4
+        vwmul.vx    v28, v8, s6
+        vslidedown.vi    v8, \in0, 4
+        vwmul.vx    v29, v8, s7
+        vslidedown.vi    v8, \in0, 4
+        vwmul.vx    v30, v8, s8
+        vslidedown.vi    v8, \in0, 4
+        vwmul.vx    v31, v8, s9
+
+        sum_sub         v28, \in1, s7, +, 2
+        sum_sub         v29, \in1, s9, -, 2
+        sum_sub         v30, \in1, s6, -, 2
+        sum_sub         v31, \in1, s8, -, 2
+        sum_sub         v28, \in2, s8, +, 2
+        sum_sub         v29, \in2, s6, -, 2
+        sum_sub         v30, \in2, s9, +, 2
+        sum_sub         v31, \in2, s7, +, 2
+        sum_sub         v28, \in3, s9, +, 2
+        sum_sub         v29, \in3, s8, -, 2
+        sum_sub         v30, \in3, s7, +, 2
+        sum_sub         v31, \in3, s6, -, 2
+
+        butterfly       v24, v28, v16, v23
+        butterfly       v25, v29, v17, v22
+        butterfly       v26, v30, v18, v21
+        butterfly       v27, v31, v19, v20
+        
+        li a7, \offset
+        add a4, sp, a7
+
+        vsetivli    zero, 4, e32, m1, tu, ma
+        vse32.v v16, (a4)
+        add a4, a4, 16
+        vse32.v v17, (a4)
+        add a4, a4, 16
+        vse32.v v18, (a4)
+        add a4, a4, 16
+        vse32.v v19, (a4)
+        add a4, a4, 16
+
+        vse32.v v20, (a4)
+        add a4, a4, 16
+        vse32.v v21, (a4)
+        add a4, a4, 16
+        vse32.v v22, (a4)
+        add a4, a4, 16
+        vse32.v v23, (a4)
+        add a4, a4, 16
+
+        add a4, a4, -64
+.endm
+
+.macro scale out0, out1, out2, out3, in0, in1, in2, in3, in4, in5, in6, in7, 
shift
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    \out0\(), \in0\(), \shift
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vmv.x.s    a7, \out0\()
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    v8, \in1\(), \shift
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx    \out0\(), v8, a7
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    \out1\(), \in2\(), \shift
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vmv.x.s    a7, \out1\()
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    v8, \in3\(), \shift
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx    \out1\(), v8, a7
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    \out2\(), \in4\(), \shift
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vmv.x.s    a7, \out2\()
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    v8, \in5\(), \shift
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx    \out2\(), v8, a7
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    \out3\(), \in6\(), \shift
+        vsetivli    zero, 1, e64, m1, tu, ma
+        vmv.x.s    a7, \out3\()
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vnclip.wi    v8, \in7\(), \shift
+        vsetivli    zero, 2, e64, m1, tu, ma
+        vslide1up.vx    \out3\(), v8, a7
+.endm
+
+.macro tr_4x4_8 in0, in1, in2, in3, out0, out1, out2, out3, p1, p2
+        vsetivli    zero, 4, e16, m1, tu, ma
+        vwcvt.x.x.v     v8, \in0  
+        vsetivli    zero, 4, e32, m1, tu, ma
+        vsll.vi         v28, v8, 6
+
+        vsetivli    zero, 16, e8, m1, tu, ma
+        vmv.v.v v29, v28
+
+        vsetivli    zero, 1, e16, m1, tu, ma
+        vmv.x.s    s2, v0
+        vslidedown.vi    v12, v0, 1
+        vmv.x.s    s3, v12
+        vslidedown.vi    v12, v0, 3
+        vmv.x.s    s5, v12
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vwmul.vx    v30, \in1, s3
+        vwmul.vx    v31, \in1, s5
+        vwmacc.vx    v28, s2, \in2
+        neg    s2, s2
+        vwmacc.vx    v29, s2, \in2
+        neg    s2, s2
+        vwmacc.vx    v30, s5, \in3
+        neg    s3, s3
+        vwmacc.vx    v31, s3, \in3
+        neg    s3, s3
+
+        vsetivli    zero, 4, e32, m1, tu, ma
+        vadd.vv             \out0, v28, v30
+        vadd.vv             \out1, v29, v31
+        vsub.vv             \out2, v29, v31
+        vsub.vv             \out3, v28, v30
+.endm
+
+.macro tr_16x4 name, shift, offset, step
+func func_tr_16x4_\name, zve64x
+        mv             a1,  a5
+        addi           a3,  a5, \step * 64
+        li             a2,  \step * 128
+        load16          v16, v17, v18, v19
+
+        lla              a1, trans
+
+        vsetivli    zero, 8, e16, m1, tu, ma
+        vle16.v         v0, (a1) 
+
+        tr16_8x4        v16, v17, v18, v19, \offset
+
+        addi             a1,  a5, \step * 32
+        addi             a3,  a5, \step * 3 *32
+        li             a2,  \step * 128
+        load16          v20, v17, v18, v19
+
+        lla              a1, trans
+        addi             a1, a1, 16
+
+        vsetivli    zero, 8, e16, m1, tu, ma
+        vle16.v          v1, (a1) 
+
+        lh      s2, 0*2(a1)
+        lh      s3, 1*2(a1)
+        lh      s4, 2*2(a1)
+        lh      s5, 3*2(a1)
+        lh      s6, 4*2(a1)
+        lh      s7, 5*2(a1)
+        lh      s8, 6*2(a1)
+        lh      s9, 7*2(a1)
+
+        vsetivli    zero, 4, e16, mf2, tu, ma
+        vwmul.vx            v21, v20, s2
+        vwmul.vx            v22, v20, s3
+        vwmul.vx            v23, v20, s4
+        vwmul.vx            v24, v20, s5
+        vwmul.vx            v25, v20, s6
+        vwmul.vx            v26, v20, s7
+        vwmul.vx            v27, v20, s8
+        vwmul.vx            v28, v20, s9
+
+        add_member      v20, s3, s6, s9, s7, s4, s2, s5, s8, +, +, +, -, -, -, 
-, -, 2
+        add_member      v17, s4, s9, s5, s3, s8, s6, s2, s7, +, +, -, -, -, +, 
+, +
+        add_member      v17, s5, s7, s3, s9, s2, s8, s4, s6, +, -, -, +, +, +, 
-, -, 2
+        add_member      v18, s6, s4, s8, s2, s9, s3, s7, s5, +, -, -, +, -, -, 
+, +
+        add_member      v18, s7, s2, s6, s8, s3, s5, s9, s4, +, -, +, +, -, +, 
+, -, 2
+        add_member      v19, s8, s5, s2, s4, s7, s9, s6, s3, +, -, +, -, +, +, 
-, +
+        add_member      v19, s9, s8, s7, s6, s5, s4, s3, s2, +, -, +, -, +, -, 
+, -, 2
+
+        li a7, \offset
+        add a4, sp, a7
+
+        vsetivli    zero, 4, e32, m1, tu, ma 
+        vle32.v        v16, (a4)
+        addi            a4, a4, 16
+        vle32.v        v17, (a4)
+        addi            a4, a4, 16
+        vle32.v        v18, (a4)
+        addi            a4, a4, 16
+        vle32.v        v19, (a4)
+        addi            a4, a4, 16
+
+        butterfly16       v16, v21, v17, v22, v18, v23, v19, v24
+    .if \shift > 0
+        scale             v29, v30, v31, v24, v20, v16, v21, v17, v22, v18, 
v23, v19, \shift
+
+        transpose16_4x4_2 29, 30, 31, 24, 2, 3, 4, 5, 6, 7
+
+        mv             a1,  a6
+        addi             a3,  a6, 24 +3*32
+        li             a2, 32
+        li             a4, -32
+
+        store16         v29, v30, v31, v24, a4
+    .else
+        store_to_stack  \offset, (\offset + 240), v20, v21, v22, v23, v19, 
v18, v17, v16
+    .endif
+
+        li a7, \offset+64
+        add a4, sp, a7
+
+        vsetivli    zero, 4, e32, m1, tu, ma 
+        vle32.v        v16, (a4)
+        addi            a4, a4, 16
+        vle32.v        v17, (a4)
+        addi            a4, a4, 16
+        vle32.v        v18, (a4)
+        addi            a4, a4, 16
+        vle32.v        v19, (a4)
+        addi            a4, a4, 16
+
+        butterfly16     v16, v25, v17, v26, v18, v27, v19, v28
+   .if \shift > 0
+        scale             v29, v30, v31, v20, v20, v16, v25, v17, v26, v18, 
v27, v19, \shift
+
+        transpose16_4x4_2 29, 30, 31, 20, 2, 3, 4, 5, 6, 7
+
+        add            a1,  a6, 8
+        add            a3,  a6, (16 + 3 * 32)
+        li             a2, 32
+        li             a4, -32
+        store16        v29, v30, v31, v20, a4
+   .else
+        store_to_stack  (\offset + 64), (\offset + 176), v20, v25, v26, v27, 
v19, v18, v17, v16
+   .endif
+        ret
+endfunc
+.endm
+
+tr_16x4 noscale, 0, 2048, 4
+
+.macro tr_32x4 name, shift
+func func_tr_32x4_\name, zve64x
+        mv             t3, ra        
+        jal              func_tr_16x4_noscale
+
+        load32
+
+        lla              t2, trans
+        addi             t2, t2, 32
+
+        vsetivli    zero, 4, e16, m1, tu, ma 
+        vle16.v v0, (t2)
+        addi t2, t2, 2*4
+        vle16.v v1, (t2)
+        addi t2, t2, 2*4
+        vle16.v v2, (t2)
+        addi t2, t2, 2*4
+        vle16.v v3, (t2)
+        addi t2, t2, -2*4
+
+        li a7, 2048
+        add a4, sp, a7
+
+        li               a2, 64
+        li               t1, -64
+        
+        jal               tr_block1
+        mv                a1, t4
+        addi              a3, t4, (56 + 3 * 64)
+        scale_store       \shift
+
+        jal              tr_block2
+        addi             a1, t4, 8
+        addi             a3, t4, (48 + 3 * 64)
+        scale_store      \shift
+
+        jal              tr_block3
+        addi             a1, t4, 16
+        addi             a3, t4, (40 + 3 * 64)
+        scale_store      \shift
+
+        jal              tr_block4
+        addi             a1, t4, 24
+        addi             a3, t4, (32 + 3 * 64)
+        scale_store      \shift
+
+        jr t3
+endfunc
+.endm
+
+tr_32x4 firstpass, 7
+tr_32x4 secondpass_8, 20 - 8
+
+.macro idct_32x32 bitdepth
+func ff_hevc_idct_32x32_\bitdepth\()_rvv, zve64x
+        mv t6, ra
+        addi    sp, sp, -8*13
+        sd    ra, 8*12(sp)
+        sd    s0, 8*11(sp)
+        sd    s1, 8*10(sp)
+        sd    s2, 8*9(sp)
+        sd    s3, 8*8(sp)
+        sd    s4, 8*7(sp)
+        sd    s5, 8*6(sp)
+        sd    s6, 8*5(sp)
+        sd    s7, 8*4(sp)
+        sd    s8, 8*3(sp)
+        sd    s9, 8*2(sp)
+        sd    s10, 8*1(sp)
+        sd    s11, 8*0(sp)
+
+        add sp, sp, -16
+        mv s0, sp
+
+        csrwi   vxrm, 1
+        li  a7, 2432 
+        sub sp, sp, a7
+
+.irp i, 0, 1, 2, 3, 4, 5, 6, 7
+        li               a7, 8 * \i
+        add              a5, a0, a7
+
+        li               a7, 8 * \i * 32
+        add              t4, sp, a7
+        jal              func_tr_32x4_firstpass
+.endr
+
+.irp i, 0, 1, 2, 3, 4, 5, 6, 7
+        addi             a5, sp, 8 * \i
+        addi             t4, a0, 8 * \i * 32
+        jal              func_tr_32x4_secondpass_\bitdepth
+.endr
+
+        li              a7, 2432 
+        add             sp,  sp,  a7
+        
+        add  sp, sp, 16
+
+        ld    ra, 8*12(sp)
+        ld    s0, 8*11(sp)
+        ld    s1, 8*10(sp)
+        ld    s2, 8*9(sp)
+        ld    s3, 8*8(sp)
+        ld    s4, 8*7(sp)
+        ld    s5, 8*6(sp)
+        ld    s6, 8*5(sp)
+        ld    s7, 8*4(sp)
+        ld    s8, 8*3(sp)
+        ld    s9, 8*2(sp)
+        ld    s10, 8*1(sp)
+        ld    s11, 8*0(sp)
+        addi    sp, sp, 8*13
+
+        jr t6
+endfunc
+.endm
+
+idct_32x32 8
diff --git a/libavcodec/riscv/hevcdsp_init.c b/libavcodec/riscv/hevcdsp_init.c
index 1d8326a573..6dfb889eec 100644
--- a/libavcodec/riscv/hevcdsp_init.c
+++ b/libavcodec/riscv/hevcdsp_init.c
@@ -27,6 +27,8 @@
 #include "libavcodec/hevc/dsp.h"
 #include "libavcodec/riscv/h26x/h2656dsp.h"
 
+void ff_hevc_idct_32x32_8_rvv(int16_t *coeffs, int col_limit);
+
 #define RVV_FNASSIGN(member, v, h, fn, ext) \
         member[1][v][h] = ff_h2656_put_pixels_##8_##ext;  \
         member[3][v][h] = ff_h2656_put_pixels_##8_##ext;  \
@@ -40,27 +42,37 @@ void ff_hevc_dsp_init_riscv(HEVCDSPContext *c, const int 
bit_depth)
     const int flags = av_get_cpu_flags();
     int vlenb;
 
-    if (!(flags & AV_CPU_FLAG_RVV_I32) || !(flags & AV_CPU_FLAG_RVB))
-        return;
-
     vlenb = ff_get_rv_vlenb();
-    if (vlenb >= 32) {
-        switch (bit_depth) {
-            case 8:
-                RVV_FNASSIGN(c->put_hevc_qpel, 0, 0, pel_pixels, rvv_256);
-                RVV_FNASSIGN(c->put_hevc_epel, 0, 0, pel_pixels, rvv_256);
-                break;
-            default:
-                break;
-        }
-    } else if (vlenb >= 16) {
-        switch (bit_depth) {
-            case 8:
-                RVV_FNASSIGN(c->put_hevc_qpel, 0, 0, pel_pixels, rvv_128);
-                RVV_FNASSIGN(c->put_hevc_epel, 0, 0, pel_pixels, rvv_128);
-                break;
-            default:
-                break;
+
+    if (flags & AV_CPU_FLAG_RVV_I64)
+        if (vlenb >= 16)
+            switch (bit_depth) {
+                case 8:
+                    c->idct[3]                     = ff_hevc_idct_32x32_8_rvv;
+                    break;
+                default:
+                    break;
+            }
+
+    if ((flags & AV_CPU_FLAG_RVV_I32) && (flags & AV_CPU_FLAG_RVB)){
+        if (vlenb >= 32) {
+            switch (bit_depth) {
+                case 8:
+                    RVV_FNASSIGN(c->put_hevc_qpel, 0, 0, pel_pixels, rvv_256);
+                    RVV_FNASSIGN(c->put_hevc_epel, 0, 0, pel_pixels, rvv_256);
+                    break;
+                default:
+                    break;
+            }
+        } else if (vlenb >= 16) {
+            switch (bit_depth) {
+                case 8:
+                    RVV_FNASSIGN(c->put_hevc_qpel, 0, 0, pel_pixels, rvv_128);
+                    RVV_FNASSIGN(c->put_hevc_epel, 0, 0, pel_pixels, rvv_128);
+                    break;
+                default:
+                    break;
+            }
         }
     }
 #endif
-- 
2.25.1

_______________________________________________
ffmpeg-devel mailing list
ffmpeg-devel@ffmpeg.org
https://ffmpeg.org/mailman/listinfo/ffmpeg-devel

To unsubscribe, visit link above, or email
ffmpeg-devel-requ...@ffmpeg.org with subject "unsubscribe".

Reply via email to