Commit df13b75a authored by Shiyou Yin's avatar Shiyou Yin Committed by Michael Niedermayer

avcodec/mips: [loongson] reoptimize simple idct with mmi.

Performance of mpeg4 decoding improved about 23%(from 128fps to 158fps, tested on loongson 3A3000).
Reoptimized following functions with mmi.
1. ff_simple_idct_put_8_mmi
2. ff_simple_idct_add_8_mmi
3. ff_simple_idct_8_mmi
Signed-off-by: 's avatarMichael Niedermayer <michael@niedermayer.cc>
parent 1124df03
...@@ -20,6 +20,7 @@ ...@@ -20,6 +20,7 @@
*/ */
#include "idctdsp_mips.h" #include "idctdsp_mips.h"
#include "xvididct_mips.h"
#if HAVE_MSA #if HAVE_MSA
static av_cold void idctdsp_init_msa(IDCTDSPContext *c, AVCodecContext *avctx, static av_cold void idctdsp_init_msa(IDCTDSPContext *c, AVCodecContext *avctx,
...@@ -48,8 +49,10 @@ static av_cold void idctdsp_init_mmi(IDCTDSPContext *c, AVCodecContext *avctx, ...@@ -48,8 +49,10 @@ static av_cold void idctdsp_init_mmi(IDCTDSPContext *c, AVCodecContext *avctx,
if ((avctx->lowres != 1) && (avctx->lowres != 2) && (avctx->lowres != 3) && if ((avctx->lowres != 1) && (avctx->lowres != 2) && (avctx->lowres != 3) &&
(avctx->bits_per_raw_sample != 10) && (avctx->bits_per_raw_sample != 10) &&
(avctx->bits_per_raw_sample != 12) && (avctx->bits_per_raw_sample != 12) &&
(avctx->idct_algo == FF_IDCT_AUTO)) { ((avctx->idct_algo == FF_IDCT_AUTO) || (avctx->idct_algo == FF_IDCT_SIMPLE))) {
c->idct = ff_simple_idct_mmi; c->idct_put = ff_simple_idct_put_8_mmi;
c->idct_add = ff_simple_idct_add_8_mmi;
c->idct = ff_simple_idct_8_mmi;
c->perm_type = FF_IDCT_PERM_NONE; c->perm_type = FF_IDCT_PERM_NONE;
} }
......
...@@ -46,8 +46,8 @@ void ff_put_signed_pixels_clamped_mmi(const int16_t *block, ...@@ -46,8 +46,8 @@ void ff_put_signed_pixels_clamped_mmi(const int16_t *block,
uint8_t *av_restrict pixels, ptrdiff_t line_size); uint8_t *av_restrict pixels, ptrdiff_t line_size);
void ff_add_pixels_clamped_mmi(const int16_t *block, void ff_add_pixels_clamped_mmi(const int16_t *block,
uint8_t *av_restrict pixels, ptrdiff_t line_size); uint8_t *av_restrict pixels, ptrdiff_t line_size);
void ff_simple_idct_mmi(int16_t *block); void ff_simple_idct_8_mmi(int16_t *block);
void ff_simple_idct_put_mmi(uint8_t *dest, ptrdiff_t line_size, int16_t *block); void ff_simple_idct_put_8_mmi(uint8_t *dest, ptrdiff_t line_size, int16_t *block);
void ff_simple_idct_add_mmi(uint8_t *dest, ptrdiff_t line_size, int16_t *block); void ff_simple_idct_add_8_mmi(uint8_t *dest, ptrdiff_t line_size, int16_t *block);
#endif // #ifndef AVCODEC_MIPS_IDCTDSP_MIPS_H #endif // #ifndef AVCODEC_MIPS_IDCTDSP_MIPS_H
...@@ -24,793 +24,400 @@ ...@@ -24,793 +24,400 @@
#include "idctdsp_mips.h" #include "idctdsp_mips.h"
#include "constants.h" #include "constants.h"
#include "libavutil/mips/asmdefs.h"
#include "libavutil/mips/mmiutils.h"
#define C0 23170 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5 #define W1 22725 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define C1 22725 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5 #define W2 21407 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define C2 21407 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5 #define W3 19266 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define C3 19266 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5 #define W4 16383 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define C4 16383 //cos(i*M_PI/16)*sqrt(2)*(1<<14) - 0.5 #define W5 12873 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define C5 12873 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5 #define W6 8867 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define C6 8867 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5 #define W7 4520 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define C7 4520 //cos(i*M_PI/16)*sqrt(2)*(1<<14) + 0.5
#define ROW_SHIFT 11 #define ROW_SHIFT 11
#define COL_SHIFT 20 #define COL_SHIFT 20
#define DC_SHIFT 3
DECLARE_ALIGNED(8, static const int16_t, coeffs)[]= {
1<<(ROW_SHIFT-1), 0, 1<<(ROW_SHIFT-1), 0, DECLARE_ALIGNED(8, const int16_t, W_arr)[46] = {
1<<(ROW_SHIFT-1), 1, 1<<(ROW_SHIFT-1), 0, W4, W2, W4, W6,
C4, C4, C4, C4, W1, W3, W5, W7,
C4, -C4, C4, -C4, W4, W6, -W4, -W2,
C2, C6, C2, C6, W3, -W7, -W1, -W5,
C6, -C2, C6, -C2, W4, -W6, -W4, W2,
C1, C3, C1, C3, W5, -W1, W7, W3,
C5, C7, C5, C7, W4, -W2, W4, -W6,
C3, -C7, C3, -C7, W7, -W5, W3, -W1,
-C1, -C5, -C1, -C5, 1024, 0, 1024, 0, //ff_p32_1024 = 0x0000040000000400ULL
C5, -C1, C5, -C1, 0, -1, -1, -1, //mask = 0xffffffffffff0000ULL
C7, C3, C7, C3, 32, 32, 32, 32 //ff_p16_32 = 0x0020002000200020ULL
C7, -C5, C7, -C5,
C3, -C1, C3, -C1
}; };
void ff_simple_idct_mmi(int16_t *block) void ff_simple_idct_8_mmi(int16_t *block)
{ {
DECLARE_ALIGNED(8, int64_t, align_tmp)[16]; BACKUP_REG
int16_t * const temp= (int16_t*)align_tmp;
__asm__ volatile ( __asm__ volatile (
#undef DC_COND_IDCT
#define DC_COND_IDCT(src0, src4, src1, src5, dst, rounder, rarg, shift) \ #define IDCT_ROW_COND_DC(src1, src2) \
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\ "dmfc1 $11, "#src1" \n\t" \
"ldc1 $f2, " #src4 " \n\t" /* R6 R2 r6 r2 */\ "dmfc1 $12, "#src2" \n\t" \
"ldc1 $f4, " #src1 " \n\t" /* R3 R1 r3 r1 */\ "and $11, $11, $9 \n\t" \
"ldc1 $f6, " #src5 " \n\t" /* R7 R5 r7 r5 */\ "or $10, $11, $12 \n\t" \
"ldc1 $f8, %3 \n\t" \
"and $f8, $f8, $f0 \n\t" \
"or $f8, $f8, $f2 \n\t" \
"or $f8, $f8, $f4 \n\t" \
"or $f8, $f8, $f6 \n\t" \
"packsswh $f8, $f8, $f8 \n\t" \
"li $11, " #shift " \n\t" \
"mfc1 $10, $f8 \n\t" \
"mtc1 $11, $f18 \n\t" \
"beqz $10, 1f \n\t" \ "beqz $10, 1f \n\t" \
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\ \
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ "punpcklhw $f30, "#src1", "#src2" \n\t" \
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\ "punpckhhw $f31, "#src1", "#src2" \n\t" \
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ /* s6, s4, s2, s0 */ \
"ldc1 $f10, 32(%2) \n\t" /* C6 C2 C6 C2 */\ "punpcklhw "#src1", $f30, $f31 \n\t" \
"pmaddhw $f10, $f10, $f2 \n\t" /* C6R6+C2R2 C6r6+C2r2 */\ /* s7, s5, s3, s1 */ \
"ldc1 $f12, 40(%2) \n\t" /* -C2 C6 -C2 C6 */\ "punpckhhw "#src2", $f30, $f31 \n\t" \
"pmaddhw $f2, $f2, $f12 \n\t" /* -C2R6+C6R2 -C2r6+C6r2 */\ \
"ldc1 $f14, 48(%2) \n\t" /* C3 C1 C3 C1 */\ "pmaddhw $f30, "#src1", $f18 \n\t" \
"ldc1 $f16, " #rarg " \n\t" \ "pmaddhw $f31, "#src2", $f19 \n\t" \
"pmaddhw $f14, $f14, $f4 \n\t" /* C3R3+C1R1 C3r3+C1r1 */\ "paddw $f28, $f30, $f31 \n\t" \
#rounder " $f8, $f8, $f16 \n\t" \ "psubw $f29, $f30, $f31 \n\t" \
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ "punpcklwd $f30, $f28, $f29 \n\t" \
"paddw $f8, $f8, $f10 \n\t" /* A0 a0 */\ "punpckhwd $f31, $f28, $f29 \n\t" \
"psubw $f12, $f12, $f10 \n\t" /* A3 a3 */\ "paddw $f26, $f30, $f31 \n\t" \
"ldc1 $f10, 56(%2) \n\t" /* C7 C5 C7 C5 */\ "paddw $f26, $f26, $f16 \n\t" \
"ldc1 $f16, " #rarg " \n\t" \ /* $f26: src[7], src[0] */ \
"pmaddhw $f10, $f10, $f6 \n\t" /* C7R7+C5R5 C7r7+C5r5 */\ "psraw $f26, $f26, $f17 \n\t" \
#rounder " $f0, $f0, $f16 \n\t" \ \
"paddw $f2, $f2, $f0 \n\t" /* A1 a1 */\ "pmaddhw $f30, "#src1", $f20 \n\t" \
"ldc1 $f16, 64(%2) \n\t" \ "pmaddhw $f31, "#src2", $f21 \n\t" \
"paddw $f0, $f0, $f0 \n\t" \ "paddw $f28, $f30, $f31 \n\t" \
"psubw $f0, $f0, $f2 \n\t" /* A2 a2 */\ "psubw $f29, $f30, $f31 \n\t" \
"pmaddhw $f4, $f4, $f16 \n\t" /* -C7R3+C3R1 -C7r3+C3r1 */\ "punpcklwd $f30, $f28, $f29 \n\t" \
"paddw $f14, $f14, $f10 \n\t" /* B0 b0 */\ "punpckhwd $f31, $f28, $f29 \n\t" \
"ldc1 $f10, 72(%2) \n\t" /* -C5 -C1 -C5 -C1 */\ "paddw $f27, $f30, $f31 \n\t" \
"pmaddhw $f10, $f10, $f6 \n\t" /* -C5R7-C1R5 -C5r7-C1r5 */\ "paddw $f27, $f27, $f16 \n\t" \
"paddw $f14, $f14, $f8 \n\t" /* A0+B0 a0+b0 */\ /* $f27: src[6], src[1] */ \
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\ "psraw $f27, $f27, $f17 \n\t" \
"psubw $f8, $f8, $f14 \n\t" /* A0-B0 a0-b0 */\ \
"paddw $f10, $f10, $f4 \n\t" /* B1 b1 */\ "pmaddhw $f30, "#src1", $f22 \n\t" \
"psraw $f14, $f14, $f18 \n\t" \ "pmaddhw $f31, "#src2", $f23 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \ "paddw $f28, $f30, $f31 \n\t" \
"mov.d $f4, $f2 \n\t" /* A1 a1 */\ "psubw $f29, $f30, $f31 \n\t" \
"paddw $f2, $f2, $f10 \n\t" /* A1+B1 a1+b1 */\ "punpcklwd $f30, $f28, $f29 \n\t" \
"psubw $f4, $f4, $f10 \n\t" /* A1-B1 a1-b1 */\ "punpckhwd $f31, $f28, $f29 \n\t" \
"psraw $f2, $f2, $f18 \n\t" \ "paddw $f28, $f30, $f31 \n\t" \
"psraw $f4, $f4, $f18 \n\t" \ "paddw $f28, $f28, $f16 \n\t" \
"packsswh $f14, $f14, $f2 \n\t" /* A1+B1 a1+b1 A0+B0 a0+b0 */\ /* $f28: src[5], src[2] */ \
"packsswh $f4, $f4, $f8 \n\t" /* A0-B0 a0-b0 A1-B1 a1-b1 */\ "psraw $f28, $f28, $f17 \n\t" \
"sdc1 $f14, " #dst " \n\t" \ \
"ldc1 $f2, " #src1 " \n\t" /* R3 R1 r3 r1 */\ "pmaddhw $f30, "#src1", $f24 \n\t" \
"ldc1 $f8, 80(%2) \n\t" /* -C1 C5 -C1 C5 */\ "pmaddhw $f31, "#src2", $f25 \n\t" \
"sdc1 $f4, 24+" #dst " \n\t" \ "paddw "#src1", $f30, $f31 \n\t" \
"pmaddhw $f8, $f8, $f2 \n\t" /* -C1R3+C5R1 -C1r3+C5r1 */\ "psubw "#src2", $f30, $f31 \n\t" \
"ldc1 $f16, 96(%2) \n\t" \ "punpcklwd $f30, "#src1", "#src2" \n\t" \
"ldc1 $f14, 88(%2) \n\t" /* C3 C7 C3 C7 */\ "punpckhwd $f31, "#src1", "#src2" \n\t" \
"pmaddhw $f2, $f2, $f16 \n\t" /* -C5R3+C7R1 -C5r3+C7r1 */\ "paddw $f29, $f30, $f31 \n\t" \
"pmaddhw $f14, $f14, $f6 \n\t" /* C3R7+C7R5 C3r7+C7r5 */\ "paddw $f29, $f29, $f16 \n\t" \
"ldc1 $f16, 104(%2) \n\t" \ /* $f29: src[4], src[3] */ \
"mov.d $f4, $f0 \n\t" /* A2 a2 */\ "psraw $f29, $f29, $f17 \n\t" \
"pmaddhw $f6, $f6, $f16 \n\t" /* -C1R7+C3R5 -C1r7+C3r5 */\ \
"paddw $f8, $f8, $f14 \n\t" /* B2 b2 */\ "punpcklhw "#src1", $f26, $f27 \n\t" \
"paddw $f4, $f4, $f8 \n\t" /* A2+B2 a2+b2 */\ "punpckhhw $f30, $f27, $f26 \n\t" \
"psubw $f0, $f0, $f8 \n\t" /* a2-B2 a2-b2 */\ "punpcklhw $f31, $f28, $f29 \n\t" \
"psraw $f4, $f4, $f18 \n\t" \ "punpckhhw "#src2", $f29, $f28 \n\t" \
"psraw $f0, $f0, $f18 \n\t" \ /* src[3], src[2], src[1], src[0] */ \
"mov.d $f8, $f12 \n\t" /* A3 a3 */\ "punpcklwd "#src1", "#src1", $f31 \n\t" \
"paddw $f6, $f6, $f2 \n\t" /* B3 b3 */\ /* src[7], src[6], src[5], src[4] */ \
"paddw $f12, $f12, $f6 \n\t" /* A3+B3 a3+b3 */\ "punpcklwd "#src2", "#src2", $f30 \n\t" \
"psubw $f8, $f8, $f6 \n\t" /* a3-B3 a3-b3 */\ "j 2f \n\t" \
"psraw $f12, $f12, $f18 \n\t" \ \
"packsswh $f4, $f4, $f12 \n\t" /* A3+B3 a3+b3 A2+B2 a2+b2 */\
"sdc1 $f4, 8+" #dst " \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"packsswh $f8, $f8, $f0 \n\t" /* A2-B2 a2-b2 A3-B3 a3-b3 */\
"sdc1 $f8, 16+" #dst " \n\t" \
"b 2f \n\t" \
"1: \n\t" \ "1: \n\t" \
"li $10, 16 \n\t" \ "li $10, 3 \n\t" \
"mtc1 $10, $f16 \n\t" \ "dmtc1 $10, $f30 \n\t" \
"psllw $f0, $f0, $f16 \n\t" \ "psllh $f28, "#src1", $f30 \n\t" \
"ldc1 $f16, %4 \n\t" \ "dmtc1 $9, $f31 \n\t" \
"paddw $f0, $f0, $f16 \n\t" \ "punpcklhw $f29, $f28, $f28 \n\t" \
"li $10, 13 \n\t" \ "and $f29, $f29, $f31 \n\t" \
"mtc1 $10, $f16 \n\t" \ "paddw $f28, $f28, $f29 \n\t" \
"psraw $f0, $f0, $f16 \n\t" \ "punpcklwd "#src1", $f28, $f28 \n\t" \
"packsswh $f0, $f0, $f0 \n\t" \ "punpcklwd "#src2", $f28, $f28 \n\t" \
"sdc1 $f0, " #dst " \n\t" \ "2: \n\t" \
"sdc1 $f0, 8+" #dst " \n\t" \
"sdc1 $f0, 16+" #dst " \n\t" \ /* idctRowCondDC row0~8 */
"sdc1 $f0, 24+" #dst " \n\t" \
"2: \n\t" /* load W */
"gslqc1 $f19, $f18, 0x00(%[w_arr]) \n\t"
#undef Z_COND_IDCT "gslqc1 $f21, $f20, 0x10(%[w_arr]) \n\t"
#define Z_COND_IDCT(src0, src4, src1, src5, dst, rounder, rarg, shift, bt) \ "gslqc1 $f23, $f22, 0x20(%[w_arr]) \n\t"
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\ "gslqc1 $f25, $f24, 0x30(%[w_arr]) \n\t"
"ldc1 $f2, " #src4 " \n\t" /* R6 R2 r6 r2 */\ "gslqc1 $f17, $f16, 0x40(%[w_arr]) \n\t"
"ldc1 $f4, " #src1 " \n\t" /* R3 R1 r3 r1 */\ /* load source in block */
"ldc1 $f6, " #src5 " \n\t" /* R7 R5 r7 r5 */\ "gslqc1 $f1, $f0, 0x00(%[block]) \n\t"
"mov.d $f8, $f0 \n\t" \ "gslqc1 $f3, $f2, 0x10(%[block]) \n\t"
"or $f8, $f8, $f2 \n\t" \ "gslqc1 $f5, $f4, 0x20(%[block]) \n\t"
"or $f8, $f8, $f4 \n\t" \ "gslqc1 $f7, $f6, 0x30(%[block]) \n\t"
"or $f8, $f8, $f6 \n\t" \ "gslqc1 $f9, $f8, 0x40(%[block]) \n\t"
"packsswh $f8, $f8, $f8 \n\t" \ "gslqc1 $f11, $f10, 0x50(%[block]) \n\t"
"mfc1 $10, $f8 \n\t" \ "gslqc1 $f13, $f12, 0x60(%[block]) \n\t"
"beqz $10, " #bt " \n\t" \ "gslqc1 $f15, $f14, 0x70(%[block]) \n\t"
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ /* $9: mask ; $f17: ROW_SHIFT */
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\ "dmfc1 $9, $f17 \n\t"
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ "li $10, 11 \n\t"
"ldc1 $f10, 32(%2) \n\t" /* C6 C2 C6 C2 */\ "mtc1 $10, $f17 \n\t"
"pmaddhw $f10, $f10, $f2 \n\t" /* C6R6+C2R2 C6r6+C2r2 */\ IDCT_ROW_COND_DC($f0,$f1)
"ldc1 $f12, 40(%2) \n\t" /* -C2 C6 -C2 C6 */\ IDCT_ROW_COND_DC($f2,$f3)
"pmaddhw $f2, $f2, $f12 \n\t" /* -C2R6+C6R2 -C2r6+C6r2 */\ IDCT_ROW_COND_DC($f4,$f5)
"ldc1 $f14, 48(%2) \n\t" /* C3 C1 C3 C1 */\ IDCT_ROW_COND_DC($f6,$f7)
"ldc1 $f16, " #rarg " \n\t" \ IDCT_ROW_COND_DC($f8,$f9)
"pmaddhw $f14, $f14, $f4 \n\t" /* C3R3+C1R1 C3r3+C1r1 */\ IDCT_ROW_COND_DC($f10,$f11)
#rounder " $f8, $f8, $f16 \n\t" \ IDCT_ROW_COND_DC($f12,$f13)
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ IDCT_ROW_COND_DC($f14,$f15)
"paddw $f8, $f8, $f10 \n\t" /* A0 a0 */\
"psubw $f12, $f12, $f10 \n\t" /* A3 a3 */\ #define IDCT_COL_CASE1(src, out1, out2) \
"ldc1 $f10, 56(%2) \n\t" /* C7 C5 C7 C5 */\ "pmaddhw $f26, "#src", $f18 \n\t" \
"ldc1 $f16, " #rarg " \n\t" \ "pmaddhw $f27, "#src", $f20 \n\t" \
"pmaddhw $f10, $f10, $f6 \n\t" /* C7R7+C5R5 C7r7+C5r5 */\ "pmaddhw $f28, "#src", $f22 \n\t" \
#rounder " $f0, $f0, $f16 \n\t" \ "pmaddhw $f29, "#src", $f24 \n\t" \
"paddw $f2, $f2, $f0 \n\t" /* A1 a1 */\ \
"paddw $f0, $f0, $f0 \n\t" \ "punpcklwd $f30, $f26, $f26 \n\t" \
"ldc1 $f16, 64(%2) \n\t" \ "punpckhwd $f31, $f26, $f26 \n\t" \
"psubw $f0, $f0, $f2 \n\t" /* A2 a2 */\ /* $f26: src[0], src[56] */ \
"pmaddhw $f4, $f4, $f16 \n\t" /* -C7R3+C3R1 -C7r3+C3r1 */\ "paddw $f26, $f30, $f31 \n\t" \
"paddw $f14, $f14, $f10 \n\t" /* B0 b0 */\ "punpcklwd $f30, $f27, $f27 \n\t" \
"ldc1 $f10, 72(%2) \n\t" /* -C5 -C1 -C5 -C1 */\ "punpckhwd $f31, $f27, $f27 \n\t" \
"pmaddhw $f10, $f10, $f6 \n\t" /* -C5R7-C1R5 -C5r7-C1r5 */\ /* $f27: src[8], src[48] */ \
"paddw $f14, $f14, $f8 \n\t" /* A0+B0 a0+b0 */\ "paddw $f27, $f30, $f31 \n\t" \
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\ "punpcklwd $f30, $f28, $f28 \n\t" \
"li $10, " #shift " \n\t" \ "punpckhwd $f31, $f28, $f28 \n\t" \
"psubw $f8, $f8, $f14 \n\t" /* A0-B0 a0-b0 */\ /* $f28: src[16], src[40] */ \
"mtc1 $10, $f18 \n\t" \ "paddw $f28, $f30, $f31 \n\t" \
"paddw $f10, $f10, $f4 \n\t" /* B1 b1 */\ "punpcklwd $f30, $f29, $f29 \n\t" \
"psraw $f14, $f14, $f18 \n\t" \ "punpckhwd $f31, $f29, $f29 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \ /* $f29: src[24], src[32] */ \
"mov.d $f4, $f2 \n\t" /* A1 a1 */\ "paddw $f29, $f30, $f31 \n\t" \
"paddw $f2, $f2, $f10 \n\t" /* A1+B1 a1+b1 */\ \
"psubw $f4, $f4, $f10 \n\t" /* A1-B1 a1-b1 */\ /* out1: src[24], src[16], src[8], src[0] */ \
"psraw $f2, $f2, $f18 \n\t" \ /* out2: src[56], src[48], src[40], src[32] */ \
"psraw $f4, $f4, $f18 \n\t" \ "punpcklhw $f30, $f26, $f27 \n\t" \
"packsswh $f14, $f14, $f2 \n\t" /* A1+B1 a1+b1 A0+B0 a0+b0 */\ "punpcklhw $f31, $f28, $f29 \n\t" \
"packsswh $f4, $f4, $f8 \n\t" /* A0-B0 a0-b0 A1-B1 a1-b1 */\ "punpckhwd "#out1", $f30, $f31 \n\t" \
"sdc1 $f14, " #dst " \n\t" \ "psrah "#out1", "#out1", $f16 \n\t" \
"ldc1 $f2, " #src1 " \n\t" /* R3 R1 r3 r1 */\ "punpcklhw $f30, $f27, $f26 \n\t" \
"ldc1 $f8, 80(%2) \n\t" /* -C1 C5 -C1 C5 */\ "punpcklhw $f31, $f29, $f28 \n\t" \
"sdc1 $f4, 24+" #dst " \n\t" \ "punpckhwd "#out2", $f31, $f30 \n\t" \
"pmaddhw $f8, $f8, $f2 \n\t" /* -C1R3+C5R1 -C1r3+C5r1 */\ "psrah "#out2", "#out2", $f16 \n\t"
"ldc1 $f16, 96(%2) \n\t" \
"ldc1 $f14, 88(%2) \n\t" /* C3 C7 C3 C7 */\ #define IDCT_COL_CASE2(src1, src2, out1, out2) \
"pmaddhw $f2, $f2, $f16 \n\t" /* -C5R3+C7R1 -C5r3+C7r1 */\ "pmaddhw $f28, "#src1", $f18 \n\t" \
"pmaddhw $f14, $f14, $f6 \n\t" /* C3R7+C7R5 C3r7+C7r5 */\ "pmaddhw $f29, "#src2", $f19 \n\t" \
"ldc1 $f16, 104(%2) \n\t" \ "paddw $f30, $f28, $f29 \n\t" \
"mov.d $f4, $f0 \n\t" /* A2 a2 */\ "psubw $f31, $f28, $f29 \n\t" \
"pmaddhw $f6, $f6, $f16 \n\t" /* -C1R7+C3R5 -C1r7+C3r5 */\ "punpcklwd $f28, $f30, $f31 \n\t" \
"paddw $f8, $f8, $f14 \n\t" /* B2 b2 */\ "punpckhwd $f29, $f30, $f31 \n\t" \
"paddw $f4, $f4, $f8 \n\t" /* A2+B2 a2+b2 */\ "pmaddhw $f30, "#src1", $f20 \n\t" \
"psubw $f0, $f0, $f8 \n\t" /* a2-B2 a2-b2 */\ "pmaddhw $f31, "#src2", $f21 \n\t" \
"psraw $f4, $f4, $f18 \n\t" \ /* $f26: src[0], src[56] */ \
"psraw $f0, $f0, $f18 \n\t" \ "paddw $f26, $f28, $f29 \n\t" \
"mov.d $f8, $f12 \n\t" /* A3 a3 */\ "paddw $f28, $f30, $f31 \n\t" \
"paddw $f6, $f6, $f2 \n\t" /* B3 b3 */\ "psubw $f29, $f30, $f31 \n\t" \
"paddw $f12, $f12, $f6 \n\t" /* A3+B3 a3+b3 */\ "punpcklwd $f30, $f28, $f29 \n\t" \
"psubw $f8, $f8, $f6 \n\t" /* a3-B3 a3-b3 */\ "punpckhwd $f31, $f28, $f29 \n\t" \
"psraw $f12, $f12, $f18 \n\t" \ "pmaddhw $f28, "#src1", $f22 \n\t" \
"packsswh $f4, $f4, $f12 \n\t" /* A3+B3 a3+b3 A2+B2 a2+b2 */\ "pmaddhw $f29, "#src2", $f23 \n\t" \
"sdc1 $f4, 8+" #dst " \n\t" \ /* $f27: src[8], src[48] */ \
"psraw $f8, $f8, $f18 \n\t" \ "paddw $f27, $f30, $f31 \n\t" \
"packsswh $f8, $f8, $f0 \n\t" /* A2-B2 a2-b2 A3-B3 a3-b3 */\ "paddw $f30, $f28, $f29 \n\t" \
"sdc1 $f8, 16+" #dst " \n\t" \ "psubw $f31, $f28, $f29 \n\t" \
"punpcklwd $f28, $f30, $f31 \n\t" \
//IDCT( src0, src4, src1, src5, dst, rounder, shift) "punpckhwd $f29, $f30, $f31 \n\t" \
DC_COND_IDCT(0(%0), 8(%0), 16(%0), 24(%0), 0(%1), paddw,8(%2), 11) "pmaddhw $f30, "#src1", $f24 \n\t" \
Z_COND_IDCT(32(%0), 40(%0), 48(%0), 56(%0), 32(%1), paddw,(%2), 11, 4f) "pmaddhw $f31, "#src2", $f25 \n\t" \
Z_COND_IDCT(64(%0), 72(%0), 80(%0), 88(%0), 64(%1), paddw,(%2), 11, 2f) /* $f28: src[16], src[40] */ \
Z_COND_IDCT(96(%0),104(%0),112(%0),120(%0), 96(%1), paddw,(%2), 11, 1f) "paddw $f28, $f28, $f29 \n\t" \
"paddw "#out1", $f30, $f31 \n\t" \
#undef IDCT "psubw "#out2", $f30, $f31 \n\t" \
#define IDCT(src0, src4, src1, src5, dst, shift) \ "punpcklwd $f30, "#out1", "#out2" \n\t" \
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\ "punpckhwd $f31, "#out1", "#out2" \n\t" \
"ldc1 $f2, " #src4 " \n\t" /* R6 R2 r6 r2 */\ /* $f29: src[24], src[32] */ \
"ldc1 $f4, " #src1 " \n\t" /* R3 R1 r3 r1 */\ "paddw $f29, $f30, $f31 \n\t" \
"ldc1 $f6, " #src5 " \n\t" /* R7 R5 r7 r5 */\ \
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\ /* out1: src[24], src[16], src[8], src[0] */ \
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ /* out2: src[56], src[48], src[40], src[32] */ \
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\ "punpcklhw "#out1", $f26, $f27 \n\t" \
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ "punpckhhw "#out2", $f27, $f26 \n\t" \
"ldc1 $f10, 32(%2) \n\t" /* C6 C2 C6 C2 */\ "punpcklhw $f30, $f28, $f29 \n\t" \
"pmaddhw $f10, $f10, $f2 \n\t" /* C6R6+C2R2 C6r6+C2r2 */\ "punpckhhw $f31, $f29, $f28 \n\t" \
"ldc1 $f12, 40(%2) \n\t" /* -C2 C6 -C2 C6 */\ "punpckhwd "#out1", "#out1", $f30 \n\t" \
"pmaddhw $f2, $f2, $f12 \n\t" /* -C2R6+C6R2 -C2r6+C6r2 */\ "punpckhwd "#out2", $f31, "#out2" \n\t" \
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ "psrah "#out1", "#out1", $f16 \n\t" \
"ldc1 $f14, 48(%2) \n\t" /* C3 C1 C3 C1 */\ "psrah "#out2", "#out2", $f16 \n\t"
"pmaddhw $f14, $f14, $f4 \n\t" /* C3R3+C1R1 C3r3+C1r1 */\
"paddw $f8, $f8, $f10 \n\t" /* A0 a0 */\
"psubw $f12, $f12, $f10 \n\t" /* A3 a3 */\ /* idctSparseCol col0~3 */
"mov.d $f10, $f0 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"paddw $f0, $f0, $f2 \n\t" /* A1 a1 */\ /* $f17: ff_p16_32; $f16: COL_SHIFT-16 */
"psubw $f10, $f10, $f2 \n\t" /* A2 a2 */\ "gsldlc1 $f17, 0x57(%[w_arr]) \n\t"
"ldc1 $f2, 56(%2) \n\t" /* C7 C5 C7 C5 */\ "gsldrc1 $f17, 0x50(%[w_arr]) \n\t"
"ldc1 $f16, 64(%2) \n\t" \ "li $10, 4 \n\t"
"pmaddhw $f2, $f2, $f6 \n\t" /* C7R7+C5R5 C7r7+C5r5 */\ "dmtc1 $10, $f16 \n\t"
"pmaddhw $f4, $f4, $f16 \n\t" /* -C7R3+C3R1 -C7r3+C3r1 */\ "paddh $f0, $f0, $f17 \n\t"
"li $10, " #shift " \n\t" \ /* Transpose row[0,2,4,6] */
"paddw $f14, $f14, $f2 \n\t" /* B0 b0 */\ "punpcklhw $f26, $f0, $f4 \n\t"
"ldc1 $f2, 72(%2) \n\t" /* -C5 -C1 -C5 -C1 */\ "punpckhhw $f27, $f0, $f4 \n\t"
"mtc1 $10, $f18 \n\t" \ "punpcklhw $f28, $f8, $f12 \n\t"
"pmaddhw $f2, $f2, $f6 \n\t" /* -C5R7-C1R5 -C5r7-C1r5 */\ "punpckhhw $f29, $f8, $f12 \n\t"
"paddw $f14, $f14, $f8 \n\t" /* A0+B0 a0+b0 */\ "punpcklwd $f0, $f26, $f28 \n\t"
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\ "punpckhwd $f4, $f26, $f28 \n\t"
"psubw $f8, $f8, $f14 \n\t" /* A0-B0 a0-b0 */\ "punpcklwd $f8, $f27, $f29 \n\t"
"paddw $f2, $f2, $f4 \n\t" /* B1 b1 */\ "punpckhwd $f12, $f27, $f29 \n\t"
"psraw $f14, $f14, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \ "or $f26, $f2, $f6 \n\t"
"mov.d $f4, $f0 \n\t" /* A1 a1 */\ "or $f26, $f26, $f10 \n\t"
"paddw $f0, $f0, $f2 \n\t" /* A1+B1 a1+b1 */\ "or $f26, $f26, $f14 \n\t"
"psubw $f4, $f4, $f2 \n\t" /* A1-B1 a1-b1 */\ "dmfc1 $10, $f26 \n\t"
"psraw $f0, $f0, $f18 \n\t" \ "bnez $10, 1f \n\t"
"psraw $f4, $f4, $f18 \n\t" \ /* case1: In this case, row[1,3,5,7] are all zero */
"packsswh $f14, $f14, $f14 \n\t" /* A0+B0 a0+b0 */\ /* col0: $f0: col[24,16,8,0]; $f2: col[56,48,40,32] */
"swc1 $f14, " #dst " \n\t" \ IDCT_COL_CASE1($f0, $f0, $f2)
"packsswh $f0, $f0, $f0 \n\t" /* A1+B1 a1+b1 */\ /* col1: $f4: col[25,17,9,1]; $f6: col[57,49,41,33] */
"swc1 $f0, 16+" #dst " \n\t" \ IDCT_COL_CASE1($f4, $f4, $f6)
"packsswh $f4, $f4, $f4 \n\t" /* A1-B1 a1-b1 */\ /* col2: $f8: col[26,18,10,2]; $f10: col[58,50,42,34] */
"swc1 $f4, 96+" #dst " \n\t" \ IDCT_COL_CASE1($f8, $f8, $f10)
"packsswh $f8, $f8, $f8 \n\t" /* A0-B0 a0-b0 */\ /* col3: $f12: col[27,19,11,3]; $f14: col[59,51,43,35] */
"swc1 $f8, 112+" #dst " \n\t" \ IDCT_COL_CASE1($f12, $f12, $f14)
"ldc1 $f0, " #src1 " \n\t" /* R3 R1 r3 r1 */\ "j 2f \n\t"
"ldc1 $f8, 80(%2) \n\t" /* -C1 C5 -C1 C5 */\
"pmaddhw $f8, $f8, $f0 \n\t" /* -C1R3+C5R1 -C1r3+C5r1 */\
"ldc1 $f16, 96(%2) \n\t" \
"ldc1 $f14, 88(%2) \n\t" /* C3 C7 C3 C7 */\
"pmaddhw $f0, $f0, $f16 \n\t" /* -C5R3+C7R1 -C5r3+C7r1 */\
"pmaddhw $f14, $f14, $f6 \n\t" /* C3R7+C7R5 C3r7+C7r5 */\
"ldc1 $f16, 104(%2) \n\t" \
"mov.d $f4, $f10 \n\t" /* A2 a2 */\
"pmaddhw $f6, $f6, $f16 \n\t" /* -C1R7+C3R5 -C1r7+C3r5 */\
"paddw $f8, $f8, $f14 \n\t" /* B2 b2 */\
"paddw $f4, $f4, $f8 \n\t" /* A2+B2 a2+b2 */\
"psubw $f10, $f10, $f8 \n\t" /* a2-B2 a2-b2 */\
"psraw $f4, $f4, $f18 \n\t" \
"psraw $f10, $f10, $f18 \n\t" \
"mov.d $f8, $f12 \n\t" /* A3 a3 */\
"paddw $f6, $f6, $f0 \n\t" /* B3 b3 */\
"paddw $f12, $f12, $f6 \n\t" /* A3+B3 a3+b3 */\
"psubw $f8, $f8, $f6 \n\t" /* a3-B3 a3-b3 */\
"psraw $f12, $f12, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"packsswh $f4, $f4, $f4 \n\t" /* A2+B2 a2+b2 */\
"packsswh $f12, $f12, $f12 \n\t" /* A3+B3 a3+b3 */\
"swc1 $f4, 32+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A3-B3 a3-b3 */\
"packsswh $f10, $f10, $f10 \n\t" /* A2-B2 a2-b2 */\
"swc1 $f12, 48+" #dst " \n\t" \
"swc1 $f8, 64+" #dst " \n\t" \
"swc1 $f10, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift)
IDCT( (%1), 64(%1), 32(%1), 96(%1), 0(%0), 20)
IDCT( 8(%1), 72(%1), 40(%1), 104(%1), 4(%0), 20)
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
IDCT( 24(%1), 88(%1), 56(%1), 120(%1), 12(%0), 20)
"b 9f \n\t"
"# .p2align 4 \n\t" "1: \n\t"
"4: \n\t" /* case2: row[1,3,5,7] are not all zero */
Z_COND_IDCT(64(%0), 72(%0), 80(%0), 88(%0), 64(%1),paddw,(%2), 11, 6f) /* Transpose */
Z_COND_IDCT(96(%0),104(%0),112(%0),120(%0), 96(%1),paddw,(%2), 11, 5f) "punpcklhw $f26, $f2, $f6 \n\t"
"punpckhhw $f27, $f2, $f6 \n\t"
#undef IDCT "punpcklhw $f28, $f10, $f14 \n\t"
#define IDCT(src0, src4, src1, src5, dst, shift) \ "punpckhhw $f29, $f10, $f14 \n\t"
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\ "punpcklwd $f2, $f26, $f28 \n\t"
"ldc1 $f2, " #src4 " \n\t" /* R6 R2 r6 r2 */\ "punpckhwd $f6, $f26, $f28 \n\t"
"ldc1 $f6, " #src5 " \n\t" /* R7 R5 r7 r5 */\ "punpcklwd $f10, $f27, $f29 \n\t"
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\ "punpckhwd $f14, $f27, $f29 \n\t"
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\ /* col0: $f0: col[24,16,8,0]; $f2: col[56,48,40,32] */
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ IDCT_COL_CASE2($f0, $f2, $f0, $f2)
"ldc1 $f10, 32(%2) \n\t" /* C6 C2 C6 C2 */\ /* col1: $f4: col[25,17,9,1]; $f6: col[57,49,41,33] */
"pmaddhw $f10, $f10, $f2 \n\t" /* C6R6+C2R2 C6r6+C2r2 */\ IDCT_COL_CASE2($f4, $f6, $f4, $f6)
"ldc1 $f12, 40(%2) \n\t" /* -C2 C6 -C2 C6 */\ /* col2: $f8: col[26,18,10,2]; $f10: col[58,50,42,34] */
"pmaddhw $f2, $f2, $f12 \n\t" /* -C2R6+C6R2 -C2r6+C6r2 */\ IDCT_COL_CASE2($f8, $f10, $f8, $f10)
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ /* col3: $f12: col[27,19,11,3]; $f14: col[59,51,43,35] */
"paddw $f8, $f8, $f10 \n\t" /* A0 a0 */\ IDCT_COL_CASE2($f12, $f14, $f12, $f14)
"psubw $f12, $f12, $f10 \n\t" /* A3 a3 */\
"mov.d $f10, $f0 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"paddw $f0, $f0, $f2 \n\t" /* A1 a1 */\
"psubw $f10, $f10, $f2 \n\t" /* A2 a2 */\
"ldc1 $f2, 56(%2) \n\t" /* C7 C5 C7 C5 */\
"li $10, " #shift " \n\t" \
"pmaddhw $f2, $f2, $f6 \n\t" /* C7R7+C5R5 C7r7+C5r5 */\
"ldc1 $f14, 72(%2) \n\t" /* -C5 -C1 -C5 -C1 */\
"mtc1 $10, $f18 \n\t" \
"pmaddhw $f14, $f14, $f6 \n\t" /* -C5R7-C1R5 -C5r7-C1r5 */\
"paddw $f2, $f2, $f8 \n\t" /* A0+B0 a0+b0 */\
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\
"psubw $f8, $f8, $f2 \n\t" /* A0-B0 a0-b0 */\
"psraw $f2, $f2, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"mov.d $f4, $f0 \n\t" /* A1 a1 */\
"paddw $f0, $f0, $f14 \n\t" /* A1+B1 a1+b1 */\
"psubw $f4, $f4, $f14 \n\t" /* A1-B1 a1-b1 */\
"psraw $f0, $f0, $f18 \n\t" \
"psraw $f4, $f4, $f18 \n\t" \
"packsswh $f2, $f2, $f2 \n\t" /* A0+B0 a0+b0 */\
"swc1 $f2, " #dst " \n\t" \
"packsswh $f0, $f0, $f0 \n\t" /* A1+B1 a1+b1 */\
"swc1 $f0, 16+" #dst " \n\t" \
"packsswh $f4, $f4, $f4 \n\t" /* A1-B1 a1-b1 */\
"swc1 $f4, 96+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A0-B0 a0-b0 */\
"swc1 $f8, 112+" #dst " \n\t" \
"ldc1 $f2, 88(%2) \n\t" /* C3 C7 C3 C7 */\
"ldc1 $f16, 104(%2) \n\t" \
"pmaddhw $f2, $f2, $f6 \n\t" /* C3R7+C7R5 C3r7+C7r5 */\
"mov.d $f4, $f10 \n\t" /* A2 a2 */\
"pmaddhw $f6, $f6, $f16 \n\t" /* -C1R7+C3R5 -C1r7+C3r5 */\
"paddw $f4, $f4, $f2 \n\t" /* A2+B2 a2+b2 */\
"psubw $f10, $f10, $f2 \n\t" /* a2-B2 a2-b2 */\
"psraw $f4, $f4, $f18 \n\t" \
"psraw $f10, $f10, $f18 \n\t" \
"mov.d $f2, $f12 \n\t" /* A3 a3 */\
"paddw $f12, $f12, $f6 \n\t" /* A3+B3 a3+b3 */\
"psubw $f2, $f2, $f6 \n\t" /* a3-B3 a3-b3 */\
"psraw $f12, $f12, $f18 \n\t" \
"psraw $f2, $f2, $f18 \n\t" \
"packsswh $f4, $f4, $f4 \n\t" /* A2+B2 a2+b2 */\
"packsswh $f12, $f12, $f12 \n\t" /* A3+B3 a3+b3 */\
"swc1 $f4, 32+" #dst " \n\t" \
"packsswh $f2, $f2, $f2 \n\t" /* A3-B3 a3-b3 */\
"packsswh $f10, $f10, $f10 \n\t" /* A2-B2 a2-b2 */\
"swc1 $f12, 48+" #dst " \n\t" \
"swc1 $f2, 64+" #dst " \n\t" \
"swc1 $f10, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift)
IDCT( (%1), 64(%1), 32(%1), 96(%1), 0(%0), 20)
IDCT( 8(%1), 72(%1), 40(%1), 104(%1), 4(%0), 20)
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
IDCT( 24(%1), 88(%1), 56(%1), 120(%1), 12(%0), 20)
"b 9f \n\t"
"# .p2align 4 \n\t"
"6: \n\t"
Z_COND_IDCT(96(%0),104(%0),112(%0),120(%0), 96(%1),paddw,(%2), 11, 7f)
#undef IDCT
#define IDCT(src0, src4, src1, src5, dst, shift) \
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\
"ldc1 $f6, " #src5 " \n\t" /* R7 R5 r7 r5 */\
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"mov.d $f10, $f0 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"ldc1 $f2, 56(%2) \n\t" /* C7 C5 C7 C5 */\
"pmaddhw $f2, $f2, $f6 \n\t" /* C7R7+C5R5 C7r7+C5r5 */\
"ldc1 $f14, 72(%2) \n\t" /* -C5 -C1 -C5 -C1 */\
"li $10, " #shift " \n\t" \
"pmaddhw $f14, $f14, $f6 \n\t" /* -C5R7-C1R5 -C5r7-C1r5 */\
"paddw $f2, $f2, $f8 \n\t" /* A0+B0 a0+b0 */\
"mtc1 $10, $f18 \n\t" \
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\
"psubw $f8, $f8, $f2 \n\t" /* A0-B0 a0-b0 */\
"psraw $f2, $f2, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"mov.d $f4, $f0 \n\t" /* A1 a1 */\
"paddw $f0, $f0, $f14 \n\t" /* A1+B1 a1+b1 */\
"psubw $f4, $f4, $f14 \n\t" /* A1-B1 a1-b1 */\
"psraw $f0, $f0, $f18 \n\t" \
"psraw $f4, $f4, $f18 \n\t" \
"packsswh $f2, $f2, $f2 \n\t" /* A0+B0 a0+b0 */\
"swc1 $f2, " #dst " \n\t" \
"packsswh $f0, $f0, $f0 \n\t" /* A1+B1 a1+b1 */\
"swc1 $f0, 16+" #dst " \n\t" \
"packsswh $f4, $f4, $f4 \n\t" /* A1-B1 a1-b1 */\
"swc1 $f4, 96+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A0-B0 a0-b0 */\
"swc1 $f8, 112+" #dst " \n\t" \
"ldc1 $f2, 88(%2) \n\t" /* C3 C7 C3 C7 */\
"ldc1 $f16, 104(%2) \n\t" \
"pmaddhw $f2, $f2, $f6 \n\t" /* C3R7+C7R5 C3r7+C7r5 */\
"mov.d $f4, $f10 \n\t" /* A2 a2 */\
"pmaddhw $f6, $f6, $f16 \n\t" /* -C1R7+C3R5 -C1r7+C3r5 */\
"paddw $f4, $f4, $f2 \n\t" /* A2+B2 a2+b2 */\
"psubw $f10, $f10, $f2 \n\t" /* a2-B2 a2-b2 */\
"psraw $f4, $f4, $f18 \n\t" \
"psraw $f10, $f10, $f18 \n\t" \
"mov.d $f2, $f12 \n\t" /* A3 a3 */\
"paddw $f12, $f12, $f6 \n\t" /* A3+B3 a3+b3 */\
"psubw $f2, $f2, $f6 \n\t" /* a3-B3 a3-b3 */\
"psraw $f12, $f12, $f18 \n\t" \
"psraw $f2, $f2, $f18 \n\t" \
"packsswh $f4, $f4, $f4 \n\t" /* A2+B2 a2+b2 */\
"packsswh $f12, $f12, $f12 \n\t" /* A3+B3 a3+b3 */\
"swc1 $f4, 32+" #dst " \n\t" \
"packsswh $f2, $f2, $f2 \n\t" /* A3-B3 a3-b3 */\
"packsswh $f10, $f10, $f10 \n\t" /* A2-B2 a2-b2 */\
"swc1 $f12, 48+" #dst " \n\t" \
"swc1 $f2, 64+" #dst " \n\t" \
"swc1 $f10, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift)
IDCT( (%1), 64(%1), 32(%1), 96(%1), 0(%0), 20)
IDCT( 8(%1), 72(%1), 40(%1), 104(%1), 4(%0), 20)
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
IDCT( 24(%1), 88(%1), 56(%1), 120(%1), 12(%0), 20)
"b 9f \n\t"
"# .p2align 4 \n\t"
"2: \n\t" "2: \n\t"
Z_COND_IDCT(96(%0),104(%0),112(%0),120(%0), 96(%1),paddw,(%2), 11, 3f) /* Transpose */
"punpcklhw $f26, $f0, $f4 \n\t"
#undef IDCT "punpckhhw $f27, $f0, $f4 \n\t"
#define IDCT(src0, src4, src1, src5, dst, shift) \ "punpcklhw $f28, $f8, $f12 \n\t"
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\ "punpckhhw $f29, $f8, $f12 \n\t"
"ldc1 $f4, " #src1 " \n\t" /* R3 R1 r3 r1 */\ "punpcklwd $f0, $f26, $f28 \n\t"
"ldc1 $f6, " #src5 " \n\t" /* R7 R5 r7 r5 */\ "punpckhwd $f4, $f26, $f28 \n\t"
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\ "punpcklwd $f8, $f27, $f29 \n\t"
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ "punpckhwd $f12, $f27, $f29 \n\t"
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\ /* Transpose */
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ "punpcklhw $f26, $f2, $f6 \n\t"
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ "punpckhhw $f27, $f2, $f6 \n\t"
"ldc1 $f14, 48(%2) \n\t" /* C3 C1 C3 C1 */\ "punpcklhw $f28, $f10, $f14 \n\t"
"pmaddhw $f14, $f14, $f4 \n\t" /* C3R3+C1R1 C3r3+C1r1 */\ "punpckhhw $f29, $f10, $f14 \n\t"
"mov.d $f10, $f0 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ "punpcklwd $f2, $f26, $f28 \n\t"
"ldc1 $f2, 56(%2) \n\t" /* C7 C5 C7 C5 */\ "punpckhwd $f6, $f26, $f28 \n\t"
"pmaddhw $f2, $f2, $f6 \n\t" /* C7R7+C5R5 C7r7+C5r5 */\ "punpcklwd $f10, $f27, $f29 \n\t"
"ldc1 $f16, 64(%2) \n\t" \ "punpckhwd $f14, $f27, $f29 \n\t"
"pmaddhw $f4, $f4, $f16 \n\t" /* -C7R3+C3R1 -C7r3+C3r1 */\
"paddw $f14, $f14, $f2 \n\t" /* B0 b0 */\ /* idctSparseCol col4~7 */
"ldc1 $f2, 72(%2) \n\t" /* -C5 -C1 -C5 -C1 */\
"li $10, " #shift " \n\t" \ "paddh $f1, $f1, $f17 \n\t"
"pmaddhw $f2, $f2, $f6 \n\t" /* -C5R7-C1R5 -C5r7-C1r5 */\ /* Transpose */
"paddw $f14, $f14, $f8 \n\t" /* A0+B0 a0+b0 */\ "punpcklhw $f26, $f1, $f5 \n\t"
"mtc1 $10, $f18 \n\t" \ "punpckhhw $f27, $f1, $f5 \n\t"
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\ "punpcklhw $f28, $f9, $f13 \n\t"
"psubw $f8, $f8, $f14 \n\t" /* A0-B0 a0-b0 */\ "punpckhhw $f29, $f9, $f13 \n\t"
"paddw $f2, $f2, $f4 \n\t" /* B1 b1 */\ "punpcklwd $f1, $f26, $f28 \n\t"
"psraw $f14, $f14, $f18 \n\t" \ "punpckhwd $f5, $f26, $f28 \n\t"
"psraw $f8, $f8, $f18 \n\t" \ "punpcklwd $f9, $f27, $f29 \n\t"
"mov.d $f4, $f0 \n\t" /* A1 a1 */\ "punpckhwd $f13, $f27, $f29 \n\t"
"paddw $f0, $f0, $f2 \n\t" /* A1+B1 a1+b1 */\
"psubw $f4, $f4, $f2 \n\t" /* A1-B1 a1-b1 */\ "or $f26, $f3, $f7 \n\t"
"psraw $f0, $f0, $f18 \n\t" \ "or $f26, $f26, $f11 \n\t"
"psraw $f4, $f4, $f18 \n\t" \ "or $f26, $f26, $f15 \n\t"
"packsswh $f14, $f14, $f14 \n\t" /* A0+B0 a0+b0 */\ "dmfc1 $10, $f26 \n\t"
"swc1 $f14, " #dst " \n\t" \ "bnez $10, 1f \n\t"
"packsswh $f0, $f0, $f0 \n\t" /* A1+B1 a1+b1 */\ /* case1: In this case, row[1,3,5,7] are all zero */
"swc1 $f0, 16+" #dst " \n\t" \ /* col4: $f1: col[24,16,8,0]; $f3: col[56,48,40,32] */
"packsswh $f4, $f4, $f4 \n\t" /* A1-B1 a1-b1 */\ IDCT_COL_CASE1($f1, $f1, $f3)
"swc1 $f4, 96+" #dst " \n\t" \ /* col5: $f5: col[25,17,9,1]; $f7: col[57,49,41,33] */
"packsswh $f8, $f8, $f8 \n\t" /* A0-B0 a0-b0 */\ IDCT_COL_CASE1($f5, $f5, $f7)
"swc1 $f8, 112+" #dst " \n\t" \ /* col6: $f9: col[26,18,10,2]; $f11: col[58,50,42,34] */
"ldc1 $f0, " #src1 " \n\t" /* R3 R1 r3 r1 */\ IDCT_COL_CASE1($f9, $f9, $f11)
"ldc1 $f8, 80(%2) \n\t" /* -C1 C5 -C1 C5 */\ /* col7: $f13: col[27,19,11,3]; $f15: col[59,51,43,35] */
"pmaddhw $f8, $f8, $f0 \n\t" /* -C1R3+C5R1 -C1r3+C5r1 */\ IDCT_COL_CASE1($f13, $f13, $f15)
"ldc1 $f14, 88(%2) \n\t" /* C3 C7 C3 C7 */\ "j 2f \n\t"
"ldc1 $f16, 96(%2) \n\t" \
"pmaddhw $f0, $f0, $f16 \n\t" /* -C5R3+C7R1 -C5r3+C7r1 */\
"pmaddhw $f14, $f14, $f6 \n\t" /* C3R7+C7R5 C3r7+C7r5 */\
"mov.d $f4, $f10 \n\t" /* A2 a2 */\
"ldc1 $f16, 104(%2) \n\t" \
"pmaddhw $f6, $f6, $f16 \n\t" /* -C1R7+C3R5 -C1r7+C3r5 */\
"paddw $f8, $f8, $f14 \n\t" /* B2 b2 */\
"paddw $f4, $f4, $f8 \n\t" /* A2+B2 a2+b2 */\
"psubw $f10, $f10, $f8 \n\t" /* a2-B2 a2-b2 */\
"psraw $f4, $f4, $f18 \n\t" \
"psraw $f10, $f10, $f18 \n\t" \
"mov.d $f8, $f12 \n\t" /* A3 a3 */\
"paddw $f6, $f6, $f0 \n\t" /* B3 b3 */\
"paddw $f12, $f12, $f6 \n\t" /* A3+B3 a3+b3 */\
"psubw $f8, $f8, $f6 \n\t" /* a3-B3 a3-b3 */\
"psraw $f12, $f12, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"packsswh $f4, $f4, $f4 \n\t" /* A2+B2 a2+b2 */\
"packsswh $f12, $f12, $f12 \n\t" /* A3+B3 a3+b3 */\
"swc1 $f4, 32+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A3-B3 a3-b3 */\
"packsswh $f10, $f10, $f10 \n\t" /* A2-B2 a2-b2 */\
"swc1 $f12, 48+" #dst " \n\t" \
"swc1 $f8, 64+" #dst " \n\t" \
"swc1 $f10, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift)
IDCT( (%1), 64(%1), 32(%1), 96(%1), 0(%0), 20)
IDCT( 8(%1), 72(%1), 40(%1), 104(%1), 4(%0), 20)
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
IDCT( 24(%1), 88(%1), 56(%1), 120(%1), 12(%0), 20)
"b 9f \n\t"
"# .p2align 4 \n\t"
"3: \n\t"
#undef IDCT
#define IDCT(src0, src4, src1, src5, dst, shift) \
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\
"ldc1 $f4, " #src1 " \n\t" /* R3 R1 r3 r1 */\
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"ldc1 $f14, 48(%2) \n\t" /* C3 C1 C3 C1 */\
"pmaddhw $f14, $f14, $f4 \n\t" /* C3R3+C1R1 C3r3+C1r1 */\
"mov.d $f10, $f0 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"ldc1 $f6, 64(%2) \n\t" \
"pmaddhw $f6, $f6, $f4 \n\t" /* -C7R3+C3R1 -C7r3+C3r1 */\
"li $10, " #shift " \n\t" \
"paddw $f14, $f14, $f8 \n\t" /* A0+B0 a0+b0 */\
"mtc1 $10, $f18 \n\t" \
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\
"psubw $f8, $f8, $f14 \n\t" /* A0-B0 a0-b0 */\
"psraw $f14, $f14, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"mov.d $f2, $f0 \n\t" /* A1 a1 */\
"paddw $f0, $f0, $f6 \n\t" /* A1+B1 a1+b1 */\
"psubw $f2, $f2, $f6 \n\t" /* A1-B1 a1-b1 */\
"psraw $f0, $f0, $f18 \n\t" \
"psraw $f2, $f2, $f18 \n\t" \
"packsswh $f14, $f14, $f14 \n\t" /* A0+B0 a0+b0 */\
"swc1 $f14, " #dst " \n\t" \
"packsswh $f0, $f0, $f0 \n\t" /* A1+B1 a1+b1 */\
"swc1 $f0, 16+" #dst " \n\t" \
"packsswh $f2, $f2, $f2 \n\t" /* A1-B1 a1-b1 */\
"swc1 $f2, 96+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A0-B0 a0-b0 */\
"swc1 $f8, 112+" #dst " \n\t" \
"ldc1 $f8, 80(%2) \n\t" /* -C1 C5 -C1 C5 */\
"ldc1 $f16, 96(%2) \n\t" \
"pmaddhw $f8, $f8, $f4 \n\t" /* -C1R3+C5R1 -C1r3+C5r1 */\
"pmaddhw $f4, $f4, $f16 \n\t" /* -C5R3+C7R1 -C5r3+C7r1 */\
"mov.d $f2, $f10 \n\t" /* A2 a2 */\
"paddw $f2, $f2, $f8 \n\t" /* A2+B2 a2+b2 */\
"psubw $f10, $f10, $f8 \n\t" /* a2-B2 a2-b2 */\
"psraw $f2, $f2, $f18 \n\t" \
"psraw $f10, $f10, $f18 \n\t" \
"mov.d $f8, $f12 \n\t" /* A3 a3 */\
"paddw $f12, $f12, $f4 \n\t" /* A3+B3 a3+b3 */\
"psubw $f8, $f8, $f4 \n\t" /* a3-B3 a3-b3 */\
"psraw $f12, $f12, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"packsswh $f2, $f2, $f2 \n\t" /* A2+B2 a2+b2 */\
"packsswh $f12, $f12, $f12 \n\t" /* A3+B3 a3+b3 */\
"swc1 $f2, 32+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A3-B3 a3-b3 */\
"packsswh $f10, $f10, $f10 \n\t" /* A2-B2 a2-b2 */\
"swc1 $f12, 48+" #dst " \n\t" \
"swc1 $f8, 64+" #dst " \n\t" \
"swc1 $f10, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift)
IDCT( (%1), 64(%1), 32(%1), 96(%1), 0(%0), 20)
IDCT( 8(%1), 72(%1), 40(%1), 104(%1), 4(%0), 20)
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
IDCT( 24(%1), 88(%1), 56(%1), 120(%1), 12(%0), 20)
"b 9f \n\t"
"# .p2align 4 \n\t"
"5: \n\t"
#undef IDCT
#define IDCT(src0, src4, src1, src5, dst, shift) \
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\
"ldc1 $f2, " #src4 " \n\t" /* R6 R2 r6 r2 */\
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"ldc1 $f10, 32(%2) \n\t" /* C6 C2 C6 C2 */\
"pmaddhw $f10, $f10, $f2 \n\t" /* C6R6+C2R2 C6r6+C2r2 */\
"ldc1 $f12, 40(%2) \n\t" /* -C2 C6 -C2 C6 */\
"pmaddhw $f2, $f2, $f12 \n\t" /* -C2R6+C6R2 -C2r6+C6r2 */\
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"paddw $f8, $f8, $f10 \n\t" /* A0 a0 */\
"psubw $f12, $f12, $f10 \n\t" /* A3 a3 */\
"mov.d $f10, $f0 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"paddw $f0, $f0, $f2 \n\t" /* A1 a1 */\
"psubw $f10, $f10, $f2 \n\t" /* A2 a2 */\
"ldc1 $f4, 8+" #src0 " \n\t" /* R4 R0 r4 r0 */\
"ldc1 $f6, 8+" #src4 " \n\t" /* R6 R2 r6 r2 */\
"ldc1 $f2, 16(%2) \n\t" /* C4 C4 C4 C4 */\
"pmaddhw $f2, $f2, $f4 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"ldc1 $f14, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\
"pmaddhw $f4, $f4, $f14 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"ldc1 $f14, 32(%2) \n\t" /* C6 C2 C6 C2 */\
"ldc1 $f16, 40(%2) \n\t" \
"pmaddhw $f14, $f14, $f6 \n\t" /* C6R6+C2R2 C6r6+C2r2 */\
"pmaddhw $f6, $f6, $f16 \n\t" /* -C2R6+C6R2 -C2r6+C6r2 */\
"paddw $f14, $f14, $f2 \n\t" /* A0 a0 */\
"paddw $f2, $f2, $f2 \n\t" /* 2C0 2c0 */\
"psubw $f2, $f2, $f14 \n\t" /* A3 a3 */\
"li $10, " #shift " \n\t" \
"paddw $f6, $f6, $f4 \n\t" /* A1 a1 */\
"mtc1 $10, $f18 \n\t" \
"paddw $f4, $f4, $f4 \n\t" /* 2C1 2c1 */\
"psubw $f4, $f4, $f6 \n\t" /* A2 a2 */\
"psraw $f8, $f8, $f18 \n\t" \
"psraw $f14, $f14, $f18 \n\t" \
"psraw $f6, $f6, $f18 \n\t" \
"packsswh $f8, $f8, $f14 \n\t" /* A0 a0 */\
"sdc1 $f8, " #dst " \n\t" \
"psraw $f0, $f0, $f18 \n\t" \
"packsswh $f0, $f0, $f6 \n\t" /* A1 a1 */\
"sdc1 $f0, 16+" #dst " \n\t" \
"sdc1 $f0, 96+" #dst " \n\t" \
"sdc1 $f8, 112+" #dst " \n\t" \
"psraw $f10, $f10, $f18 \n\t" \
"psraw $f12, $f12, $f18 \n\t" \
"psraw $f4, $f4, $f18 \n\t" \
"packsswh $f10, $f10, $f4 \n\t" /* A2-B2 a2-b2 */\
"sdc1 $f10, 32+" #dst " \n\t" \
"psraw $f2, $f2, $f18 \n\t" \
"packsswh $f12, $f12, $f2 \n\t" /* A3+B3 a3+b3 */\
"sdc1 $f12, 48+" #dst " \n\t" \
"sdc1 $f12, 64+" #dst " \n\t" \
"sdc1 $f10, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift)
IDCT( 0(%1), 64(%1), 32(%1), 96(%1), 0(%0), 20)
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
"b 9f \n\t"
"# .p2align 4 \n\t"
"1: \n\t" "1: \n\t"
/* case2: row[1,3,5,7] are not all zero */
/* Transpose */
"punpcklhw $f26, $f3, $f7 \n\t"
"punpckhhw $f27, $f3, $f7 \n\t"
"punpcklhw $f28, $f11, $f15 \n\t"
"punpckhhw $f29, $f11, $f15 \n\t"
"punpcklwd $f3, $f26, $f28 \n\t"
"punpckhwd $f7, $f26, $f28 \n\t"
"punpcklwd $f11, $f27, $f29 \n\t"
"punpckhwd $f15, $f27, $f29 \n\t"
/* col4: $f1: col[24,16,8,0]; $f3: col[56,48,40,32] */
IDCT_COL_CASE2($f1, $f3, $f1, $f3)
/* col5: $f5: col[25,17,9,1]; $f7: col[57,49,41,33] */
IDCT_COL_CASE2($f5, $f7, $f5, $f7)
/* col6: $f9: col[26,18,10,2]; $f11: col[58,50,42,34] */
IDCT_COL_CASE2($f9, $f11, $f9, $f11)
/* col7: $f13: col[27,19,11,3]; $f15: col[59,51,43,35] */
IDCT_COL_CASE2($f13, $f15, $f13, $f15)
#undef IDCT "2: \n\t"
#define IDCT(src0, src4, src1, src5, dst, shift) \ /* Transpose */
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\ "punpcklhw $f26, $f1, $f5 \n\t"
"ldc1 $f2, " #src4 " \n\t" /* R6 R2 r6 r2 */\ "punpckhhw $f27, $f1, $f5 \n\t"
"ldc1 $f4, " #src1 " \n\t" /* R3 R1 r3 r1 */\ "punpcklhw $f28, $f9, $f13 \n\t"
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\ "punpckhhw $f29, $f9, $f13 \n\t"
"li $10, " #shift " \n\t" \ "punpcklwd $f1, $f26, $f28 \n\t"
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ "punpckhwd $f5, $f26, $f28 \n\t"
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\ "punpcklwd $f9, $f27, $f29 \n\t"
"mtc1 $10, $f18 \n\t" \ "punpckhwd $f13, $f27, $f29 \n\t"
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ /* Transpose */
"ldc1 $f10, 32(%2) \n\t" /* C6 C2 C6 C2 */\ "punpcklhw $f26, $f3, $f7 \n\t"
"pmaddhw $f10, $f10, $f2 \n\t" /* C6R6+C2R2 C6r6+C2r2 */\ "punpckhhw $f27, $f3, $f7 \n\t"
"ldc1 $f12, 40(%2) \n\t" /* -C2 C6 -C2 C6 */\ "punpcklhw $f28, $f11, $f15 \n\t"
"pmaddhw $f2, $f2, $f12 \n\t" /* -C2R6+C6R2 -C2r6+C6r2 */\ "punpckhhw $f29, $f11, $f15 \n\t"
"mov.d $f12, $f8 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\ "punpcklwd $f3, $f26, $f28 \n\t"
"ldc1 $f14, 48(%2) \n\t" /* C3 C1 C3 C1 */\ "punpckhwd $f7, $f26, $f28 \n\t"
"pmaddhw $f14, $f14, $f4 \n\t" /* C3R3+C1R1 C3r3+C1r1 */\ "punpcklwd $f11, $f27, $f29 \n\t"
"paddw $f8, $f8, $f10 \n\t" /* A0 a0 */\ "punpckhwd $f15, $f27, $f29 \n\t"
"psubw $f12, $f12, $f10 \n\t" /* A3 a3 */\ /* Store */
"mov.d $f10, $f0 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\ "gssqc1 $f1, $f0, 0x00(%[block]) \n\t"
"paddw $f0, $f0, $f2 \n\t" /* A1 a1 */\ "gssqc1 $f5, $f4, 0x10(%[block]) \n\t"
"psubw $f10, $f10, $f2 \n\t" /* A2 a2 */\ "gssqc1 $f9, $f8, 0x20(%[block]) \n\t"
"ldc1 $f2, 64(%2) \n\t" \ "gssqc1 $f13, $f12, 0x30(%[block]) \n\t"
"pmaddhw $f2, $f2, $f4 \n\t" /* -C7R3+C3R1 -C7r3+C3r1 */\ "gssqc1 $f3, $f2, 0x40(%[block]) \n\t"
"paddw $f14, $f14, $f8 \n\t" /* A0+B0 a0+b0 */\ "gssqc1 $f7, $f6, 0x50(%[block]) \n\t"
"paddw $f8, $f8, $f8 \n\t" /* 2A0 2a0 */\ "gssqc1 $f11, $f10, 0x60(%[block]) \n\t"
"psubw $f8, $f8, $f14 \n\t" /* A0-B0 a0-b0 */\ "gssqc1 $f15, $f14, 0x70(%[block]) \n\t"
"psraw $f14, $f14, $f18 \n\t" \
"psraw $f8, $f8, $f18 \n\t" \ : [block]"+&r"(block)
"mov.d $f6, $f0 \n\t" /* A1 a1 */\ : [w_arr]"r"(W_arr)
"paddw $f0, $f0, $f2 \n\t" /* A1+B1 a1+b1 */\ : "memory"
"psubw $f6, $f6, $f2 \n\t" /* A1-B1 a1-b1 */\ );
"psraw $f0, $f0, $f18 \n\t" \
"psraw $f6, $f6, $f18 \n\t" \
"packsswh $f14, $f14, $f14 \n\t" /* A0+B0 a0+b0 */\
"swc1 $f14, " #dst " \n\t" \
"packsswh $f0, $f0, $f0 \n\t" /* A1+B1 a1+b1 */\
"swc1 $f0, 16+" #dst " \n\t" \
"packsswh $f6, $f6, $f6 \n\t" /* A1-B1 a1-b1 */\
"swc1 $f6, 96+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A0-B0 a0-b0 */\
"swc1 $f8, 112+" #dst " \n\t" \
"ldc1 $f8, 80(%2) \n\t" /* -C1 C5 -C1 C5 */\
"ldc1 $f16, 96(%2) \n\t" \
"pmaddhw $f8, $f8, $f4 \n\t" /* -C1R3+C5R1 -C1r3+C5r1 */\
"pmaddhw $f4, $f4, $f16 \n\t" /* -C5R3+C7R1 -C5r3+C7r1 */\
"mov.d $f6, $f10 \n\t" /* A2 a2 */\
"paddw $f6, $f6, $f8 \n\t" /* A2+B2 a2+b2 */\
"psubw $f10, $f10, $f8 \n\t" /* a2-B2 a2-b2 */\
"psraw $f6, $f6, $f18 \n\t" \
"psraw $f10, $f10, $f18 \n\t" \
"mov.d $f8, $f12 \n\t" /* A3 a3 */\
"paddw $f12, $f12, $f4 \n\t" /* A3+B3 a3+b3 */\
"psubw $f8, $f8, $f4 \n\t" /* a3-B3 a3-b3 */\
"psraw $f12, $f12, $f18 \n\t" \
"packsswh $f6, $f6, $f6 \n\t" /* A2+B2 a2+b2 */\
"swc1 $f6, 32+" #dst " \n\t" \
"psraw $f8, $f8, $f18 \n\t" \
"packsswh $f12, $f12, $f12 \n\t" /* A3+B3 a3+b3 */\
"swc1 $f12, 48+" #dst " \n\t" \
"packsswh $f8, $f8, $f8 \n\t" /* A3-B3 a3-b3 */\
"packsswh $f10, $f10, $f10 \n\t" /* A2-B2 a2-b2 */\
"swc1 $f8, 64+" #dst " \n\t" \
"swc1 $f10, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift)
IDCT( (%1), 64(%1), 32(%1), 96(%1), 0(%0), 20)
IDCT( 8(%1), 72(%1), 40(%1), 104(%1), 4(%0), 20)
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
IDCT( 24(%1), 88(%1), 56(%1), 120(%1), 12(%0), 20)
"b 9f \n\t"
"# .p2align 4 \n\t"
"7: \n\t"
#undef IDCT
#define IDCT(src0, src4, src1, src5, dst, shift) \
"ldc1 $f0, " #src0 " \n\t" /* R4 R0 r4 r0 */\
"ldc1 $f8, 16(%2) \n\t" /* C4 C4 C4 C4 */\
"li $10, " #shift " \n\t" \
"pmaddhw $f8, $f8, $f0 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"mtc1 $10, $f18 \n\t" \
"ldc1 $f10, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\
"pmaddhw $f0, $f0, $f10 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"psraw $f8, $f8, $f18 \n\t" \
"psraw $f0, $f0, $f18 \n\t" \
"ldc1 $f4, 8+" #src0 " \n\t" /* R4 R0 r4 r0 */\
"ldc1 $f2, 16(%2) \n\t" /* C4 C4 C4 C4 */\
"pmaddhw $f2, $f2, $f4 \n\t" /* C4R4+C4R0 C4r4+C4r0 */\
"ldc1 $f14, 24(%2) \n\t" /* -C4 C4 -C4 C4 */\
"pmaddhw $f4, $f4, $f14 \n\t" /* -C4R4+C4R0 -C4r4+C4r0 */\
"ldc1 $f14, 32(%2) \n\t" /* C6 C2 C6 C2 */\
"psraw $f2, $f2, $f18 \n\t" \
"packsswh $f8, $f8, $f2 \n\t" /* A0 a0 */\
"sdc1 $f8, " #dst " \n\t" \
"psraw $f4, $f4, $f18 \n\t" \
"packsswh $f0, $f0, $f4 \n\t" /* A1 a1 */\
"sdc1 $f0, 16+" #dst " \n\t" \
"sdc1 $f0, 96+" #dst " \n\t" \
"sdc1 $f8, 112+" #dst " \n\t" \
"sdc1 $f0, 32+" #dst " \n\t" \
"sdc1 $f8, 48+" #dst " \n\t" \
"sdc1 $f8, 64+" #dst " \n\t" \
"sdc1 $f0, 80+" #dst " \n\t"
//IDCT( src0, src4, src1, src5, dst, shift) RECOVER_REG
IDCT( 0(%1), 64(%1), 32(%1), 96(%1), 0(%0), 20) }
IDCT( 16(%1), 80(%1), 48(%1), 112(%1), 8(%0), 20)
"9: \n\t" void ff_simple_idct_put_8_mmi(uint8_t *dest, ptrdiff_t line_size, int16_t *block)
::"r"(block),"r"(temp),"r"(coeffs),"m"(ff_wm1010),"m"(ff_d40000) {
: "$10","$11" ff_simple_idct_8_mmi(block);
); ff_put_pixels_clamped_mmi(block, dest, line_size);
}
void ff_simple_idct_add_8_mmi(uint8_t *dest, ptrdiff_t line_size, int16_t *block)
{
ff_simple_idct_8_mmi(block);
ff_add_pixels_clamped_mmi(block, dest, line_size);
} }
...@@ -201,6 +201,55 @@ ...@@ -201,6 +201,55 @@
#endif /* HAVE_LOONGSON2 */ #endif /* HAVE_LOONGSON2 */
/**
* backup register
*/
#define BACKUP_REG \
double temp_backup_reg[8]; \
if (_MIPS_SIM == _ABI64) \
__asm__ volatile ( \
"gssqc1 $f25, $f24, 0x00(%[temp]) \n\t" \
"gssqc1 $f27, $f26, 0x10(%[temp]) \n\t" \
"gssqc1 $f29, $f28, 0x20(%[temp]) \n\t" \
"gssqc1 $f31, $f30, 0x30(%[temp]) \n\t" \
: \
: [temp]"r"(temp_backup_reg) \
: "memory" \
); \
else \
__asm__ volatile ( \
"gssqc1 $f22, $f20, 0x00(%[temp]) \n\t" \
"gssqc1 $f26, $f24, 0x10(%[temp]) \n\t" \
"gssqc1 $f30, $f28, 0x20(%[temp]) \n\t" \
: \
: [temp]"r"(temp_backup_reg) \
: "memory" \
);
/**
* recover register
*/
#define RECOVER_REG \
if (_MIPS_SIM == _ABI64) \
__asm__ volatile ( \
"gslqc1 $f25, $f24, 0x00(%[temp]) \n\t" \
"gslqc1 $f27, $f26, 0x10(%[temp]) \n\t" \
"gslqc1 $f29, $f28, 0x20(%[temp]) \n\t" \
"gslqc1 $f31, $f30, 0x30(%[temp]) \n\t" \
: \
: [temp]"r"(temp_backup_reg) \
: "memory" \
); \
else \
__asm__ volatile ( \
"gslqc1 $f22, $f20, 0x00(%[temp]) \n\t" \
"gslqc1 $f26, $f24, 0x10(%[temp]) \n\t" \
"gslqc1 $f30, $f28, 0x20(%[temp]) \n\t" \
: \
: [temp]"r"(temp_backup_reg) \
: "memory" \
);
#define TRANSPOSE_4H(m1, m2, m3, m4, t1, t2, t3, t4, t5, r1, zero, shift) \ #define TRANSPOSE_4H(m1, m2, m3, m4, t1, t2, t3, t4, t5, r1, zero, shift) \
"li "#r1", 0x93 \n\t" \ "li "#r1", 0x93 \n\t" \
"xor "#zero","#zero","#zero" \n\t" \ "xor "#zero","#zero","#zero" \n\t" \
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment