FFmpeg
vp8_idct_msa.c
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2015 Manojkumar Bhosale (Manojkumar.Bhosale@imgtec.com)
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20 
21 #include <string.h>
22 #include "libavcodec/vp8dsp.h"
23 #include "libavutil/mips/generic_macros_msa.h"
24 #include "vp8dsp_mips.h"
25 
26 static const int cospi8sqrt2minus1 = 20091;
27 static const int sinpi8sqrt2 = 35468;
28 
29 #define VP8_IDCT_1D_W(in0, in1, in2, in3, out0, out1, out2, out3) \
30 { \
31  v4i32 a1_m, b1_m, c1_m, d1_m; \
32  v4i32 c_tmp1_m, c_tmp2_m, d_tmp1_m, d_tmp2_m; \
33  v4i32 const_cospi8sqrt2minus1_m, sinpi8_sqrt2_m; \
34  \
35  const_cospi8sqrt2minus1_m = __msa_fill_w(cospi8sqrt2minus1); \
36  sinpi8_sqrt2_m = __msa_fill_w(sinpi8sqrt2); \
37  a1_m = in0 + in2; \
38  b1_m = in0 - in2; \
39  c_tmp1_m = ((in1) * sinpi8_sqrt2_m) >> 16; \
40  c_tmp2_m = in3 + (((in3) * const_cospi8sqrt2minus1_m) >> 16); \
41  c1_m = c_tmp1_m - c_tmp2_m; \
42  d_tmp1_m = (in1) + (((in1) * const_cospi8sqrt2minus1_m) >> 16); \
43  d_tmp2_m = ((in3) * sinpi8_sqrt2_m) >> 16; \
44  d1_m = d_tmp1_m + d_tmp2_m; \
45  BUTTERFLY_4(a1_m, b1_m, c1_m, d1_m, out0, out1, out2, out3); \
46 }
47 
48 void ff_vp8_idct_add_msa(uint8_t *dst, int16_t input[16], ptrdiff_t stride)
49 {
50  v8i16 input0, input1;
51  v4i32 in0, in1, in2, in3, hz0, hz1, hz2, hz3, vt0, vt1, vt2, vt3;
52  v4i32 res0, res1, res2, res3;
53  v16i8 zero = { 0 };
54  v16i8 pred0, pred1, pred2, pred3, dest0, dest1;
55  v16i8 mask = { 0, 4, 8, 12, 16, 20, 24, 28, 0, 0, 0, 0, 0, 0, 0, 0 };
56 
57  /* load short vector elements of 4x4 block */
58  LD_SH2(input, 8, input0, input1);
59  UNPCK_SH_SW(input0, in0, in1);
60  UNPCK_SH_SW(input1, in2, in3);
61  VP8_IDCT_1D_W(in0, in1, in2, in3, hz0, hz1, hz2, hz3);
62  /* transpose the block */
63  TRANSPOSE4x4_SW_SW(hz0, hz1, hz2, hz3, hz0, hz1, hz2, hz3);
64  VP8_IDCT_1D_W(hz0, hz1, hz2, hz3, vt0, vt1, vt2, vt3);
65  SRARI_W4_SW(vt0, vt1, vt2, vt3, 3);
66  /* transpose the block */
67  TRANSPOSE4x4_SW_SW(vt0, vt1, vt2, vt3, vt0, vt1, vt2, vt3);
68  LD_SB4(dst, stride, pred0, pred1, pred2, pred3);
69  ILVR_B4_SW(zero, pred0, zero, pred1, zero, pred2, zero, pred3,
70  res0, res1, res2, res3);
71  ILVR_H4_SW(zero, res0, zero, res1, zero, res2, zero, res3,
72  res0, res1, res2, res3);
73  ADD4(res0, vt0, res1, vt1, res2, vt2, res3, vt3, res0, res1, res2, res3);
74  res0 = CLIP_SW_0_255(res0);
75  res1 = CLIP_SW_0_255(res1);
76  res2 = CLIP_SW_0_255(res2);
77  res3 = CLIP_SW_0_255(res3);
78  VSHF_B2_SB(res0, res1, res2, res3, mask, mask, dest0, dest1);
79  ST_W2(dest0, 0, 1, dst, stride);
80  ST_W2(dest1, 0, 1, dst + 2 * stride, stride);
81 
82  memset(input, 0, 4 * 4 * sizeof(*input));
83 }
84 
85 void ff_vp8_idct_dc_add_msa(uint8_t *dst, int16_t in_dc[16], ptrdiff_t stride)
86 {
87  v8i16 vec;
88  v8i16 res0, res1, res2, res3;
89  v16i8 zero = { 0 };
90  v16i8 pred0, pred1, pred2, pred3, dest0, dest1;
91  v16i8 mask = { 0, 2, 4, 6, 16, 18, 20, 22, 0, 0, 0, 0, 0, 0, 0, 0 };
92 
93  vec = __msa_fill_h(in_dc[0]);
94  vec = __msa_srari_h(vec, 3);
95  LD_SB4(dst, stride, pred0, pred1, pred2, pred3);
96  ILVR_B4_SH(zero, pred0, zero, pred1, zero, pred2, zero, pred3,
97  res0, res1, res2, res3);
98  ADD4(res0, vec, res1, vec, res2, vec, res3, vec, res0, res1, res2, res3);
99  CLIP_SH4_0_255(res0, res1, res2, res3);
100  VSHF_B2_SB(res0, res1, res2, res3, mask, mask, dest0, dest1);
101  ST_W2(dest0, 0, 1, dst, stride);
102  ST_W2(dest1, 0, 1, dst + 2 * stride, stride);
103 
104  in_dc[0] = 0;
105 }
106 
107 void ff_vp8_luma_dc_wht_msa(int16_t block[4][4][16], int16_t input[16])
108 {
109  int16_t *mb_dq_coeff = &block[0][0][0];
110  v8i16 input0, input1;
111  v4i32 in0, in1, in2, in3, a1, b1, c1, d1;
112  v4i32 hz0, hz1, hz2, hz3, vt0, vt1, vt2, vt3;
113 
114  /* load short vector elements of 4x4 block */
115  LD_SH2(input, 8, input0, input1);
116  UNPCK_SH_SW(input0, in0, in1);
117  UNPCK_SH_SW(input1, in2, in3);
118  BUTTERFLY_4(in0, in1, in2, in3, a1, b1, c1, d1);
119  BUTTERFLY_4(a1, d1, c1, b1, hz0, hz1, hz3, hz2);
120  /* transpose the block */
121  TRANSPOSE4x4_SW_SW(hz0, hz1, hz2, hz3, hz0, hz1, hz2, hz3);
122  BUTTERFLY_4(hz0, hz1, hz2, hz3, a1, b1, c1, d1);
123  BUTTERFLY_4(a1, d1, c1, b1, vt0, vt1, vt3, vt2);
124  ADD4(vt0, 3, vt1, 3, vt2, 3, vt3, 3, vt0, vt1, vt2, vt3);
125  SRA_4V(vt0, vt1, vt2, vt3, 3);
126  mb_dq_coeff[0] = __msa_copy_s_h((v8i16) vt0, 0);
127  mb_dq_coeff[16] = __msa_copy_s_h((v8i16) vt1, 0);
128  mb_dq_coeff[32] = __msa_copy_s_h((v8i16) vt2, 0);
129  mb_dq_coeff[48] = __msa_copy_s_h((v8i16) vt3, 0);
130  mb_dq_coeff[64] = __msa_copy_s_h((v8i16) vt0, 2);
131  mb_dq_coeff[80] = __msa_copy_s_h((v8i16) vt1, 2);
132  mb_dq_coeff[96] = __msa_copy_s_h((v8i16) vt2, 2);
133  mb_dq_coeff[112] = __msa_copy_s_h((v8i16) vt3, 2);
134  mb_dq_coeff[128] = __msa_copy_s_h((v8i16) vt0, 4);
135  mb_dq_coeff[144] = __msa_copy_s_h((v8i16) vt1, 4);
136  mb_dq_coeff[160] = __msa_copy_s_h((v8i16) vt2, 4);
137  mb_dq_coeff[176] = __msa_copy_s_h((v8i16) vt3, 4);
138  mb_dq_coeff[192] = __msa_copy_s_h((v8i16) vt0, 6);
139  mb_dq_coeff[208] = __msa_copy_s_h((v8i16) vt1, 6);
140  mb_dq_coeff[224] = __msa_copy_s_h((v8i16) vt2, 6);
141  mb_dq_coeff[240] = __msa_copy_s_h((v8i16) vt3, 6);
142 
143  memset(input, 0, 4 * 4 * sizeof(int16_t));
144 }
145 
146 void ff_vp8_idct_dc_add4y_msa(uint8_t *dst, int16_t block[4][16],
147  ptrdiff_t stride)
148 {
149  ff_vp8_idct_dc_add_msa(dst, &block[0][0], stride);
150  ff_vp8_idct_dc_add_msa(dst + 4, &block[1][0], stride);
151  ff_vp8_idct_dc_add_msa(dst + 8, &block[2][0], stride);
152  ff_vp8_idct_dc_add_msa(dst + 12, &block[3][0], stride);
153 }
154 
155 void ff_vp8_idct_dc_add4uv_msa(uint8_t *dst, int16_t block[4][16],
156  ptrdiff_t stride)
157 {
158  ff_vp8_idct_dc_add_msa(dst, &block[0][0], stride);
159  ff_vp8_idct_dc_add_msa(dst + 4, &block[1][0], stride);
160  ff_vp8_idct_dc_add_msa(dst + stride * 4, &block[2][0], stride);
161  ff_vp8_idct_dc_add_msa(dst + stride * 4 + 4, &block[3][0], stride);
162 }
VSHF_B2_SB
#define VSHF_B2_SB(...)
Definition: generic_macros_msa.h:681
LD_SB4
#define LD_SB4(...)
Definition: generic_macros_msa.h:299
stride
int stride
Definition: mace.c:144
LD_SH2
#define LD_SH2(...)
Definition: generic_macros_msa.h:282
ff_vp8_idct_add_msa
void ff_vp8_idct_add_msa(uint8_t *dst, int16_t input[16], ptrdiff_t stride)
Definition: vp8_idct_msa.c:48
ILVR_H4_SW
#define ILVR_H4_SW(...)
Definition: generic_macros_msa.h:1458
ff_vp8_luma_dc_wht_msa
void ff_vp8_luma_dc_wht_msa(int16_t block[4][4][16], int16_t input[16])
Definition: vp8_idct_msa.c:107
ff_vp8_idct_dc_add4y_msa
void ff_vp8_idct_dc_add4y_msa(uint8_t *dst, int16_t block[4][16], ptrdiff_t stride)
Definition: vp8_idct_msa.c:146
c1
static const uint64_t c1
Definition: murmur3.c:49
ADD4
#define ADD4(in0, in1, in2, in3, in4, in5, in6, in7, out0, out1, out2, out3)
Definition: generic_macros_msa.h:2172
ff_vp8_idct_dc_add4uv_msa
void ff_vp8_idct_dc_add4uv_msa(uint8_t *dst, int16_t block[4][16], ptrdiff_t stride)
Definition: vp8_idct_msa.c:155
ILVR_B4_SW
#define ILVR_B4_SW(...)
Definition: generic_macros_msa.h:1413
VP8_IDCT_1D_W
#define VP8_IDCT_1D_W(in0, in1, in2, in3, out0, out1, out2, out3)
Definition: vp8_idct_msa.c:29
generic_macros_msa.h
a1
#define a1
Definition: regdef.h:47
vp8dsp.h
mask
static const uint16_t mask[17]
Definition: lzw.c:38
SRA_4V
#define SRA_4V(in0, in1, in2, in3, shift)
Definition: generic_macros_msa.h:1988
TRANSPOSE4x4_SW_SW
#define TRANSPOSE4x4_SW_SW(in0, in1, in2, in3, out0, out1, out2, out3)
Definition: generic_macros_msa.h:2586
UNPCK_SH_SW
#define UNPCK_SH_SW(in, out0, out1)
Definition: generic_macros_msa.h:2295
ST_W2
#define ST_W2(in, idx0, idx1, pdst, stride)
Definition: generic_macros_msa.h:450
input
and forward the test the status of outputs and forward it to the corresponding return FFERROR_NOT_READY If the filters stores internally one or a few frame for some input
Definition: filter_design.txt:172
CLIP_SH4_0_255
#define CLIP_SH4_0_255(in0, in1, in2, in3)
Definition: generic_macros_msa.h:971
uint8_t
uint8_t
Definition: audio_convert.c:194
vp8dsp_mips.h
sinpi8sqrt2
static const int sinpi8sqrt2
Definition: vp8_idct_msa.c:27
SRARI_W4_SW
#define SRARI_W4_SW(...)
Definition: generic_macros_msa.h:2141
BUTTERFLY_4
#define BUTTERFLY_4(in0, in1, in2, in3, out0, out1, out2, out3)
Definition: generic_macros_msa.h:2321
ILVR_B4_SH
#define ILVR_B4_SH(...)
Definition: generic_macros_msa.h:1412
zero
#define zero
Definition: regdef.h:64
ff_vp8_idct_dc_add_msa
void ff_vp8_idct_dc_add_msa(uint8_t *dst, int16_t in_dc[16], ptrdiff_t stride)
Definition: vp8_idct_msa.c:85
CLIP_SW_0_255
#define CLIP_SW_0_255(in)
Definition: generic_macros_msa.h:1002
cospi8sqrt2minus1
static const int cospi8sqrt2minus1
Definition: vp8_idct_msa.c:26
block
The exact code depends on how similar the blocks are and how related they are to the block
Definition: filter_design.txt:207
Generated on Wed Aug 24 2022 21:26:49 for FFmpeg by   doxygen 1.8.17