OLD | NEW |
(Empty) | |
| 1 /* |
| 2 * Copyright (c) 2015 The WebM project authors. All Rights Reserved. |
| 3 * |
| 4 * Use of this source code is governed by a BSD-style license |
| 5 * that can be found in the LICENSE file in the root of the source |
| 6 * tree. An additional intellectual property rights grant can be found |
| 7 * in the file PATENTS. All contributing project authors may |
| 8 * be found in the AUTHORS file in the root of the source tree. |
| 9 */ |
| 10 |
| 11 #include <assert.h> |
| 12 |
| 13 #include "vpx_ports/mem.h" |
| 14 #include "vp9/common/vp9_idct.h" |
| 15 #include "vp9/common/mips/msa/vp9_macros_msa.h" |
| 16 |
| 17 #define SET_COSPI_PAIR(c0_h, c1_h) ({ \ |
| 18 v8i16 out0, r0_m, r1_m; \ |
| 19 \ |
| 20 r0_m = __msa_fill_h(c0_h); \ |
| 21 r1_m = __msa_fill_h(c1_h); \ |
| 22 out0 = __msa_ilvev_h(r1_m, r0_m); \ |
| 23 \ |
| 24 out0; \ |
| 25 }) |
| 26 |
| 27 #define DOTP_CONST_PAIR(reg0, reg1, const0, const1, out0, out1) { \ |
| 28 v8i16 k0_m = __msa_fill_h(const0); \ |
| 29 v8i16 s0_m, s1_m, s2_m, s3_m; \ |
| 30 \ |
| 31 s0_m = __msa_fill_h(const1); \ |
| 32 k0_m = __msa_ilvev_h(s0_m, k0_m); \ |
| 33 \ |
| 34 s0_m = __msa_ilvl_h(-reg1, reg0); \ |
| 35 s1_m = __msa_ilvr_h(-reg1, reg0); \ |
| 36 s2_m = __msa_ilvl_h(reg0, reg1); \ |
| 37 s3_m = __msa_ilvr_h(reg0, reg1); \ |
| 38 s1_m = (v8i16)__msa_dotp_s_w(s1_m, k0_m); \ |
| 39 s0_m = (v8i16)__msa_dotp_s_w(s0_m, k0_m); \ |
| 40 s1_m = (v8i16)__msa_srari_w((v4i32)s1_m, DCT_CONST_BITS); \ |
| 41 s0_m = (v8i16)__msa_srari_w((v4i32)s0_m, DCT_CONST_BITS); \ |
| 42 out0 = __msa_pckev_h(s0_m, s1_m); \ |
| 43 \ |
| 44 s1_m = (v8i16)__msa_dotp_s_w(s3_m, k0_m); \ |
| 45 s0_m = (v8i16)__msa_dotp_s_w(s2_m, k0_m); \ |
| 46 s1_m = (v8i16)__msa_srari_w((v4i32)s1_m, DCT_CONST_BITS); \ |
| 47 s0_m = (v8i16)__msa_srari_w((v4i32)s0_m, DCT_CONST_BITS); \ |
| 48 out1 = __msa_pckev_h(s0_m, s1_m); \ |
| 49 } |
| 50 |
| 51 #define VP9_MADD_SHORT(m0, m1, c0, c1, res0, res1) { \ |
| 52 v4i32 madd0_m, madd1_m, madd2_m, madd3_m; \ |
| 53 v8i16 madd_s0_m, madd_s1_m; \ |
| 54 \ |
| 55 ILV_H_LR_SH(m0, m1, madd_s1_m, madd_s0_m); \ |
| 56 \ |
| 57 DOTP_S_W_4VECS_SW(madd_s0_m, c0, madd_s1_m, c0, \ |
| 58 madd_s0_m, c1, madd_s1_m, c1, \ |
| 59 madd0_m, madd1_m, madd2_m, madd3_m); \ |
| 60 \ |
| 61 SRARI_W_4VECS_SW(madd0_m, madd1_m, madd2_m, madd3_m, \ |
| 62 madd0_m, madd1_m, madd2_m, madd3_m, \ |
| 63 DCT_CONST_BITS); \ |
| 64 \ |
| 65 PCKEV_H_2VECS_SH(madd1_m, madd0_m, madd3_m, madd2_m, \ |
| 66 res0, res1); \ |
| 67 } |
| 68 |
| 69 #define VP9_MADD_BF(inp0, inp1, inp2, inp3, \ |
| 70 cst0, cst1, cst2, cst3, \ |
| 71 out0, out1, out2, out3) { \ |
| 72 v8i16 madd_s0_m, madd_s1_m, madd_s2_m, madd_s3_m; \ |
| 73 v4i32 tmp0_m, tmp1_m, tmp2_m, tmp3_m; \ |
| 74 v4i32 m4_m, m5_m; \ |
| 75 \ |
| 76 ILV_H_LRLR_SH(inp0, inp1, inp2, inp3, \ |
| 77 madd_s1_m, madd_s0_m, madd_s3_m, madd_s2_m); \ |
| 78 \ |
| 79 DOTP_S_W_4VECS_SW(madd_s0_m, cst0, madd_s1_m, cst0, \ |
| 80 madd_s2_m, cst2, madd_s3_m, cst2, \ |
| 81 tmp0_m, tmp1_m, tmp2_m, tmp3_m); \ |
| 82 \ |
| 83 m4_m = tmp0_m + tmp2_m; \ |
| 84 m5_m = tmp1_m + tmp3_m; \ |
| 85 tmp3_m = tmp1_m - tmp3_m; \ |
| 86 tmp2_m = tmp0_m - tmp2_m; \ |
| 87 \ |
| 88 SRARI_W_4VECS_SW(m4_m, m5_m, tmp2_m, tmp3_m, \ |
| 89 m4_m, m5_m, tmp2_m, tmp3_m, \ |
| 90 DCT_CONST_BITS); \ |
| 91 \ |
| 92 PCKEV_H_2VECS_SH(m5_m, m4_m, tmp3_m, tmp2_m, out0, out1); \ |
| 93 \ |
| 94 DOTP_S_W_4VECS_SW(madd_s0_m, cst1, madd_s1_m, cst1, \ |
| 95 madd_s2_m, cst3, madd_s3_m, cst3, \ |
| 96 tmp0_m, tmp1_m, tmp2_m, tmp3_m); \ |
| 97 \ |
| 98 m4_m = tmp0_m + tmp2_m; \ |
| 99 m5_m = tmp1_m + tmp3_m; \ |
| 100 tmp3_m = tmp1_m - tmp3_m; \ |
| 101 tmp2_m = tmp0_m - tmp2_m; \ |
| 102 \ |
| 103 SRARI_W_4VECS_SW(m4_m, m5_m, tmp2_m, tmp3_m, \ |
| 104 m4_m, m5_m, tmp2_m, tmp3_m, \ |
| 105 DCT_CONST_BITS); \ |
| 106 \ |
| 107 PCKEV_H_2VECS_SH(m5_m, m4_m, tmp3_m, tmp2_m, out2, out3); \ |
| 108 } |
| 109 |
| 110 #define TRANSPOSE8x8_H1(in0, in1, in2, in3, \ |
| 111 in4, in5, in6, in7, \ |
| 112 out0, out1, out2, out3, \ |
| 113 out4, out5, out6, out7) { \ |
| 114 v8i16 loc0_m, loc1_m; \ |
| 115 v8i16 tmp0_m, tmp1_m, tmp2_m, tmp3_m; \ |
| 116 v8i16 tmp4_m, tmp5_m, tmp6_m, tmp7_m; \ |
| 117 \ |
| 118 loc0_m = __msa_ilvr_h((in6), (in4)); \ |
| 119 loc1_m = __msa_ilvr_h((in7), (in5)); \ |
| 120 tmp0_m = __msa_ilvr_h(loc1_m, loc0_m); \ |
| 121 tmp1_m = __msa_ilvl_h(loc1_m, loc0_m); \ |
| 122 \ |
| 123 loc0_m = __msa_ilvl_h((in6), (in4)); \ |
| 124 loc1_m = __msa_ilvl_h((in7), (in5)); \ |
| 125 tmp2_m = __msa_ilvr_h(loc1_m, loc0_m); \ |
| 126 tmp3_m = __msa_ilvl_h(loc1_m, loc0_m); \ |
| 127 \ |
| 128 loc0_m = __msa_ilvr_h((in2), (in0)); \ |
| 129 loc1_m = __msa_ilvr_h((in3), (in1)); \ |
| 130 tmp4_m = __msa_ilvr_h(loc1_m, loc0_m); \ |
| 131 tmp5_m = __msa_ilvl_h(loc1_m, loc0_m); \ |
| 132 \ |
| 133 loc0_m = __msa_ilvl_h((in2), (in0)); \ |
| 134 loc1_m = __msa_ilvl_h((in3), (in1)); \ |
| 135 tmp6_m = __msa_ilvr_h(loc1_m, loc0_m); \ |
| 136 tmp7_m = __msa_ilvl_h(loc1_m, loc0_m); \ |
| 137 \ |
| 138 out0 = (v8i16)__msa_pckev_d((v2i64)tmp0_m, (v2i64)tmp4_m); \ |
| 139 out1 = (v8i16)__msa_pckod_d((v2i64)tmp0_m, (v2i64)tmp4_m); \ |
| 140 out2 = (v8i16)__msa_pckev_d((v2i64)tmp1_m, (v2i64)tmp5_m); \ |
| 141 out3 = (v8i16)__msa_pckod_d((v2i64)tmp1_m, (v2i64)tmp5_m); \ |
| 142 out4 = (v8i16)__msa_pckev_d((v2i64)tmp2_m, (v2i64)tmp6_m); \ |
| 143 out5 = (v8i16)__msa_pckod_d((v2i64)tmp2_m, (v2i64)tmp6_m); \ |
| 144 out6 = (v8i16)__msa_pckev_d((v2i64)tmp3_m, (v2i64)tmp7_m); \ |
| 145 out7 = (v8i16)__msa_pckod_d((v2i64)tmp3_m, (v2i64)tmp7_m); \ |
| 146 } |
| 147 |
| 148 #define VP9_IADST8x16_1D(r0, r1, r2, r3, r4, r5, r6, r7, \ |
| 149 r8, r9, r10, r11, r12, r13, r14, r15, \ |
| 150 out0, out1, out2, out3, out4, out5, out6, out7, \ |
| 151 out8, out9, out10, out11, \ |
| 152 out12, out13, out14, out15) { \ |
| 153 v8i16 g0_m, g1_m, g2_m, g3_m, g4_m, g5_m, g6_m, g7_m; \ |
| 154 v8i16 g8_m, g9_m, g10_m, g11_m, g12_m, g13_m, g14_m, g15_m; \ |
| 155 v8i16 h0_m, h1_m, h2_m, h3_m, h4_m, h5_m, h6_m, h7_m; \ |
| 156 v8i16 h8_m, h9_m, h10_m, h11_m; \ |
| 157 v8i16 k0_m, k1_m, k2_m, k3_m; \ |
| 158 \ |
| 159 /* stage 1 */ \ |
| 160 k0_m = SET_COSPI_PAIR(cospi_1_64, cospi_31_64); \ |
| 161 k1_m = SET_COSPI_PAIR(cospi_31_64, -cospi_1_64); \ |
| 162 k2_m = SET_COSPI_PAIR(cospi_17_64, cospi_15_64); \ |
| 163 k3_m = SET_COSPI_PAIR(cospi_15_64, -cospi_17_64); \ |
| 164 VP9_MADD_BF(r15, r0, r7, r8, k0_m, k1_m, k2_m, k3_m, \ |
| 165 g0_m, g1_m, g2_m, g3_m); \ |
| 166 \ |
| 167 k0_m = SET_COSPI_PAIR(cospi_5_64, cospi_27_64); \ |
| 168 k1_m = SET_COSPI_PAIR(cospi_27_64, -cospi_5_64); \ |
| 169 k2_m = SET_COSPI_PAIR(cospi_21_64, cospi_11_64); \ |
| 170 k3_m = SET_COSPI_PAIR(cospi_11_64, -cospi_21_64); \ |
| 171 VP9_MADD_BF(r13, r2, r5, r10, k0_m, k1_m, k2_m, k3_m, \ |
| 172 g4_m, g5_m, g6_m, g7_m); \ |
| 173 \ |
| 174 k0_m = SET_COSPI_PAIR(cospi_9_64, cospi_23_64); \ |
| 175 k1_m = SET_COSPI_PAIR(cospi_23_64, -cospi_9_64); \ |
| 176 k2_m = SET_COSPI_PAIR(cospi_25_64, cospi_7_64); \ |
| 177 k3_m = SET_COSPI_PAIR(cospi_7_64, -cospi_25_64); \ |
| 178 VP9_MADD_BF(r11, r4, r3, r12, k0_m, k1_m, k2_m, k3_m, \ |
| 179 g8_m, g9_m, g10_m, g11_m); \ |
| 180 \ |
| 181 k0_m = SET_COSPI_PAIR(cospi_13_64, cospi_19_64); \ |
| 182 k1_m = SET_COSPI_PAIR(cospi_19_64, -cospi_13_64); \ |
| 183 k2_m = SET_COSPI_PAIR(cospi_29_64, cospi_3_64); \ |
| 184 k3_m = SET_COSPI_PAIR(cospi_3_64, -cospi_29_64); \ |
| 185 VP9_MADD_BF(r9, r6, r1, r14, k0_m, k1_m, k2_m, k3_m, \ |
| 186 g12_m, g13_m, g14_m, g15_m); \ |
| 187 \ |
| 188 /* stage 2 */ \ |
| 189 k0_m = SET_COSPI_PAIR(cospi_4_64, cospi_28_64); \ |
| 190 k1_m = SET_COSPI_PAIR(cospi_28_64, -cospi_4_64); \ |
| 191 k2_m = SET_COSPI_PAIR(-cospi_28_64, cospi_4_64); \ |
| 192 VP9_MADD_BF(g1_m, g3_m, g9_m, g11_m, k0_m, k1_m, k2_m, k0_m, \ |
| 193 h0_m, h1_m, h2_m, h3_m); \ |
| 194 \ |
| 195 k0_m = SET_COSPI_PAIR(cospi_12_64, cospi_20_64); \ |
| 196 k1_m = SET_COSPI_PAIR(-cospi_20_64, cospi_12_64); \ |
| 197 k2_m = SET_COSPI_PAIR(cospi_20_64, -cospi_12_64); \ |
| 198 VP9_MADD_BF(g7_m, g5_m, g15_m, g13_m, k0_m, k1_m, k2_m, k0_m, \ |
| 199 h4_m, h5_m, h6_m, h7_m); \ |
| 200 \ |
| 201 BUTTERFLY_4(h0_m, h2_m, h6_m, h4_m, out8, out9, out11, out10); \ |
| 202 \ |
| 203 BUTTERFLY_8(g0_m, g2_m, g4_m, g6_m, g14_m, g12_m, g10_m, g8_m, \ |
| 204 h8_m, h9_m, h10_m, h11_m, h6_m, h4_m, h2_m, h0_m); \ |
| 205 \ |
| 206 /* stage 3 */ \ |
| 207 BUTTERFLY_4(h8_m, h9_m, h11_m, h10_m, out0, out1, h11_m, h10_m); \ |
| 208 \ |
| 209 k0_m = SET_COSPI_PAIR(cospi_8_64, cospi_24_64); \ |
| 210 k1_m = SET_COSPI_PAIR(cospi_24_64, -cospi_8_64); \ |
| 211 k2_m = SET_COSPI_PAIR(-cospi_24_64, cospi_8_64); \ |
| 212 VP9_MADD_BF(h0_m, h2_m, h4_m, h6_m, k0_m, k1_m, k2_m, k0_m, \ |
| 213 out4, out6, out5, out7); \ |
| 214 VP9_MADD_BF(h1_m, h3_m, h5_m, h7_m, k0_m, k1_m, k2_m, k0_m, \ |
| 215 out12, out14, out13, out15); \ |
| 216 \ |
| 217 /* stage 4 */ \ |
| 218 k0_m = SET_COSPI_PAIR(cospi_16_64, cospi_16_64); \ |
| 219 k1_m = SET_COSPI_PAIR(-cospi_16_64, -cospi_16_64); \ |
| 220 k2_m = SET_COSPI_PAIR(cospi_16_64, -cospi_16_64); \ |
| 221 k3_m = SET_COSPI_PAIR(-cospi_16_64, cospi_16_64); \ |
| 222 VP9_MADD_SHORT(h10_m, h11_m, k1_m, k2_m, out2, out3); \ |
| 223 VP9_MADD_SHORT(out6, out7, k0_m, k3_m, out6, out7); \ |
| 224 VP9_MADD_SHORT(out10, out11, k0_m, k3_m, out10, out11); \ |
| 225 VP9_MADD_SHORT(out14, out15, k1_m, k2_m, out14, out15); \ |
| 226 } |
| 227 |
| 228 #define VP9_ADDBLK_CLIP_AND_STORE_8_BYTES_4(dest, dest_stride, \ |
| 229 in0, in1, in2, in3) { \ |
| 230 uint64_t out0_m, out1_m, out2_m, out3_m; \ |
| 231 v8i16 res0_m, res1_m, res2_m, res3_m; \ |
| 232 v16u8 dest0_m, dest1_m, dest2_m, dest3_m; \ |
| 233 v16i8 tmp0_m, tmp1_m; \ |
| 234 v16i8 zero_m = { 0 }; \ |
| 235 uint8_t *dst_m = (uint8_t *)(dest); \ |
| 236 \ |
| 237 LOAD_4VECS_UB(dst_m, (dest_stride), \ |
| 238 dest0_m, dest1_m, dest2_m, dest3_m); \ |
| 239 \ |
| 240 res0_m = (v8i16)__msa_ilvr_b(zero_m, (v16i8)dest0_m); \ |
| 241 res1_m = (v8i16)__msa_ilvr_b(zero_m, (v16i8)dest1_m); \ |
| 242 res2_m = (v8i16)__msa_ilvr_b(zero_m, (v16i8)dest2_m); \ |
| 243 res3_m = (v8i16)__msa_ilvr_b(zero_m, (v16i8)dest3_m); \ |
| 244 \ |
| 245 res0_m += (v8i16)(in0); \ |
| 246 res1_m += (v8i16)(in1); \ |
| 247 res2_m += (v8i16)(in2); \ |
| 248 res3_m += (v8i16)(in3); \ |
| 249 \ |
| 250 res0_m = CLIP_UNSIGNED_CHAR_H(res0_m); \ |
| 251 res1_m = CLIP_UNSIGNED_CHAR_H(res1_m); \ |
| 252 res2_m = CLIP_UNSIGNED_CHAR_H(res2_m); \ |
| 253 res3_m = CLIP_UNSIGNED_CHAR_H(res3_m); \ |
| 254 \ |
| 255 tmp0_m = __msa_pckev_b((v16i8)res1_m, (v16i8)res0_m); \ |
| 256 tmp1_m = __msa_pckev_b((v16i8)res3_m, (v16i8)res2_m); \ |
| 257 \ |
| 258 out0_m = __msa_copy_u_d((v2i64)tmp0_m, 0); \ |
| 259 out1_m = __msa_copy_u_d((v2i64)tmp0_m, 1); \ |
| 260 out2_m = __msa_copy_u_d((v2i64)tmp1_m, 0); \ |
| 261 out3_m = __msa_copy_u_d((v2i64)tmp1_m, 1); \ |
| 262 \ |
| 263 STORE_DWORD(dst_m, out0_m); \ |
| 264 dst_m += (dest_stride); \ |
| 265 STORE_DWORD(dst_m, out1_m); \ |
| 266 dst_m += (dest_stride); \ |
| 267 STORE_DWORD(dst_m, out2_m); \ |
| 268 dst_m += (dest_stride); \ |
| 269 STORE_DWORD(dst_m, out3_m); \ |
| 270 } |
| 271 |
| 272 void vp9_idct16_1d_rows_msa(const int16_t *input, int16_t *output) { |
| 273 v8i16 loc0, loc1, loc2, loc3; |
| 274 v8i16 reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14; |
| 275 v8i16 reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15; |
| 276 v8i16 tmp5, tmp6, tmp7; |
| 277 |
| 278 /* load left top 8x8 */ |
| 279 LOAD_8VECS_SH(input, 16, reg0, reg1, reg2, reg3, reg4, reg5, reg6, reg7); |
| 280 |
| 281 /* load right top 8x8 */ |
| 282 LOAD_8VECS_SH((input + 8), 16, |
| 283 reg8, reg9, reg10, reg11, reg12, reg13, reg14, reg15); |
| 284 |
| 285 /* transpose block */ |
| 286 TRANSPOSE8x8_H1(reg0, reg1, reg2, reg3, reg4, reg5, reg6, reg7, |
| 287 reg0, reg1, reg2, reg3, reg4, reg5, reg6, reg7); |
| 288 |
| 289 /* transpose block */ |
| 290 TRANSPOSE8x8_H1(reg8, reg9, reg10, reg11, reg12, reg13, reg14, reg15, |
| 291 reg8, reg9, reg10, reg11, reg12, reg13, reg14, reg15); |
| 292 |
| 293 DOTP_CONST_PAIR(reg2, reg14, cospi_28_64, cospi_4_64, reg2, reg14); |
| 294 DOTP_CONST_PAIR(reg10, reg6, cospi_12_64, cospi_20_64, reg10, reg6); |
| 295 |
| 296 loc0 = reg2 + reg10; |
| 297 reg2 = reg2 - reg10; |
| 298 loc1 = reg14 + reg6; |
| 299 reg14 = reg14 - reg6; |
| 300 |
| 301 DOTP_CONST_PAIR(reg14, reg2, cospi_16_64, cospi_16_64, loc2, loc3); |
| 302 DOTP_CONST_PAIR(reg0, reg8, cospi_16_64, cospi_16_64, reg0, reg8); |
| 303 DOTP_CONST_PAIR(reg4, reg12, cospi_24_64, cospi_8_64, reg4, reg12); |
| 304 |
| 305 reg14 = reg8 - reg12; |
| 306 reg2 = reg8 + reg12; |
| 307 reg10 = reg0 - reg4; |
| 308 reg6 = reg0 + reg4; |
| 309 |
| 310 reg0 = reg2 - loc1; |
| 311 reg2 = reg2 + loc1; |
| 312 reg12 = reg14 - loc0; |
| 313 reg14 = reg14 + loc0; |
| 314 reg4 = reg6 - loc3; |
| 315 reg6 = reg6 + loc3; |
| 316 reg8 = reg10 - loc2; |
| 317 reg10 = reg10 + loc2; |
| 318 |
| 319 /* stage 2 */ |
| 320 DOTP_CONST_PAIR(reg1, reg15, cospi_30_64, cospi_2_64, reg1, reg15); |
| 321 DOTP_CONST_PAIR(reg9, reg7, cospi_14_64, cospi_18_64, loc2, loc3); |
| 322 |
| 323 reg9 = reg1 - loc2; |
| 324 reg1 = reg1 + loc2; |
| 325 reg7 = reg15 - loc3; |
| 326 reg15 = reg15 + loc3; |
| 327 |
| 328 DOTP_CONST_PAIR(reg5, reg11, cospi_22_64, cospi_10_64, reg5, reg11); |
| 329 DOTP_CONST_PAIR(reg13, reg3, cospi_6_64, cospi_26_64, loc0, loc1); |
| 330 |
| 331 reg13 = loc0 + reg5; |
| 332 reg5 = loc0 - reg5; |
| 333 reg3 = loc1 + reg11; |
| 334 reg11 = loc1 - reg11; |
| 335 |
| 336 loc1 = reg15 + reg3; |
| 337 reg3 = reg15 - reg3; |
| 338 loc2 = reg2 + loc1; |
| 339 reg15 = reg2 - loc1; |
| 340 |
| 341 loc1 = reg1 + reg13; |
| 342 reg13 = reg1 - reg13; |
| 343 loc0 = reg0 + loc1; |
| 344 loc1 = reg0 - loc1; |
| 345 tmp6 = loc0; |
| 346 tmp7 = loc1; |
| 347 reg0 = loc2; |
| 348 |
| 349 DOTP_CONST_PAIR(reg7, reg9, cospi_24_64, cospi_8_64, reg7, reg9); |
| 350 DOTP_CONST_PAIR((-reg5), (-reg11), cospi_8_64, cospi_24_64, reg5, reg11); |
| 351 |
| 352 loc0 = reg9 + reg5; |
| 353 reg5 = reg9 - reg5; |
| 354 reg2 = reg6 + loc0; |
| 355 reg1 = reg6 - loc0; |
| 356 |
| 357 loc0 = reg7 + reg11; |
| 358 reg11 = reg7 - reg11; |
| 359 loc1 = reg4 + loc0; |
| 360 loc2 = reg4 - loc0; |
| 361 tmp5 = loc1; |
| 362 |
| 363 DOTP_CONST_PAIR(reg5, reg11, cospi_16_64, cospi_16_64, reg5, reg11); |
| 364 |
| 365 loc0 = reg8 + reg5; |
| 366 loc1 = reg8 - reg5; |
| 367 reg4 = reg10 + reg11; |
| 368 reg9 = reg10 - reg11; |
| 369 reg10 = loc0; |
| 370 reg11 = loc1; |
| 371 |
| 372 DOTP_CONST_PAIR(reg3, reg13, cospi_16_64, cospi_16_64, reg3, reg13); |
| 373 |
| 374 reg8 = reg12 + reg3; |
| 375 reg5 = reg12 - reg3; |
| 376 reg6 = reg14 + reg13; |
| 377 reg7 = reg14 - reg13; |
| 378 reg13 = loc2; |
| 379 |
| 380 /* Transpose and store the output */ |
| 381 reg12 = tmp5; |
| 382 reg14 = tmp6; |
| 383 reg3 = tmp7; |
| 384 |
| 385 /* transpose block */ |
| 386 TRANSPOSE8x8_H1(reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14, |
| 387 reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14); |
| 388 |
| 389 STORE_8VECS_SH(output, 16, reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14); |
| 390 |
| 391 /* transpose block */ |
| 392 TRANSPOSE8x8_H1(reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15, |
| 393 reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15); |
| 394 |
| 395 STORE_8VECS_SH((output + 8), 16, |
| 396 reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15); |
| 397 } |
| 398 |
| 399 void vp9_idct16_1d_columns_addblk_msa(int16_t *input, uint8_t *dest, |
| 400 int32_t dest_stride) { |
| 401 v8i16 loc0, loc1, loc2, loc3; |
| 402 v8i16 reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14; |
| 403 v8i16 reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15; |
| 404 v8i16 tmp5, tmp6, tmp7; |
| 405 |
| 406 /* load up 8x8 */ |
| 407 LOAD_8VECS_SH(input, 16, reg0, reg1, reg2, reg3, reg4, reg5, reg6, reg7); |
| 408 |
| 409 /* load bottom 8x8 */ |
| 410 LOAD_8VECS_SH((input + 8 * 16), 16, |
| 411 reg8, reg9, reg10, reg11, reg12, reg13, reg14, reg15); |
| 412 |
| 413 DOTP_CONST_PAIR(reg2, reg14, cospi_28_64, cospi_4_64, reg2, reg14); |
| 414 DOTP_CONST_PAIR(reg10, reg6, cospi_12_64, cospi_20_64, reg10, reg6); |
| 415 |
| 416 loc0 = reg2 + reg10; |
| 417 reg2 = reg2 - reg10; |
| 418 loc1 = reg14 + reg6; |
| 419 reg14 = reg14 - reg6; |
| 420 |
| 421 DOTP_CONST_PAIR(reg14, reg2, cospi_16_64, cospi_16_64, loc2, loc3); |
| 422 DOTP_CONST_PAIR(reg0, reg8, cospi_16_64, cospi_16_64, reg0, reg8); |
| 423 DOTP_CONST_PAIR(reg4, reg12, cospi_24_64, cospi_8_64, reg4, reg12); |
| 424 |
| 425 reg14 = reg8 - reg12; |
| 426 reg2 = reg8 + reg12; |
| 427 reg10 = reg0 - reg4; |
| 428 reg6 = reg0 + reg4; |
| 429 |
| 430 reg0 = reg2 - loc1; |
| 431 reg2 = reg2 + loc1; |
| 432 reg12 = reg14 - loc0; |
| 433 reg14 = reg14 + loc0; |
| 434 reg4 = reg6 - loc3; |
| 435 reg6 = reg6 + loc3; |
| 436 reg8 = reg10 - loc2; |
| 437 reg10 = reg10 + loc2; |
| 438 |
| 439 /* stage 2 */ |
| 440 DOTP_CONST_PAIR(reg1, reg15, cospi_30_64, cospi_2_64, reg1, reg15); |
| 441 DOTP_CONST_PAIR(reg9, reg7, cospi_14_64, cospi_18_64, loc2, loc3); |
| 442 |
| 443 reg9 = reg1 - loc2; |
| 444 reg1 = reg1 + loc2; |
| 445 reg7 = reg15 - loc3; |
| 446 reg15 = reg15 + loc3; |
| 447 |
| 448 DOTP_CONST_PAIR(reg5, reg11, cospi_22_64, cospi_10_64, reg5, reg11); |
| 449 DOTP_CONST_PAIR(reg13, reg3, cospi_6_64, cospi_26_64, loc0, loc1); |
| 450 |
| 451 reg13 = loc0 + reg5; |
| 452 reg5 = loc0 - reg5; |
| 453 reg3 = loc1 + reg11; |
| 454 reg11 = loc1 - reg11; |
| 455 |
| 456 loc1 = reg15 + reg3; |
| 457 reg3 = reg15 - reg3; |
| 458 loc2 = reg2 + loc1; |
| 459 reg15 = reg2 - loc1; |
| 460 |
| 461 loc1 = reg1 + reg13; |
| 462 reg13 = reg1 - reg13; |
| 463 loc0 = reg0 + loc1; |
| 464 loc1 = reg0 - loc1; |
| 465 tmp6 = loc0; |
| 466 tmp7 = loc1; |
| 467 reg0 = loc2; |
| 468 |
| 469 DOTP_CONST_PAIR(reg7, reg9, cospi_24_64, cospi_8_64, reg7, reg9); |
| 470 DOTP_CONST_PAIR((-reg5), (-reg11), cospi_8_64, cospi_24_64, reg5, reg11); |
| 471 |
| 472 loc0 = reg9 + reg5; |
| 473 reg5 = reg9 - reg5; |
| 474 reg2 = reg6 + loc0; |
| 475 reg1 = reg6 - loc0; |
| 476 |
| 477 loc0 = reg7 + reg11; |
| 478 reg11 = reg7 - reg11; |
| 479 loc1 = reg4 + loc0; |
| 480 loc2 = reg4 - loc0; |
| 481 tmp5 = loc1; |
| 482 |
| 483 DOTP_CONST_PAIR(reg5, reg11, cospi_16_64, cospi_16_64, reg5, reg11); |
| 484 |
| 485 loc0 = reg8 + reg5; |
| 486 loc1 = reg8 - reg5; |
| 487 reg4 = reg10 + reg11; |
| 488 reg9 = reg10 - reg11; |
| 489 reg10 = loc0; |
| 490 reg11 = loc1; |
| 491 |
| 492 DOTP_CONST_PAIR(reg3, reg13, cospi_16_64, cospi_16_64, reg3, reg13); |
| 493 |
| 494 reg8 = reg12 + reg3; |
| 495 reg5 = reg12 - reg3; |
| 496 reg6 = reg14 + reg13; |
| 497 reg7 = reg14 - reg13; |
| 498 reg13 = loc2; |
| 499 |
| 500 /* Transpose and store the output */ |
| 501 reg12 = tmp5; |
| 502 reg14 = tmp6; |
| 503 reg3 = tmp7; |
| 504 |
| 505 SRARI_H_4VECS_SH(reg0, reg2, reg4, reg6, reg0, reg2, reg4, reg6, 6); |
| 506 VP9_ADDBLK_CLIP_AND_STORE_8_BYTES_4(dest, dest_stride, |
| 507 reg0, reg2, reg4, reg6); |
| 508 SRARI_H_4VECS_SH(reg8, reg10, reg12, reg14, reg8, reg10, reg12, reg14, 6); |
| 509 VP9_ADDBLK_CLIP_AND_STORE_8_BYTES_4((dest + (4 * dest_stride)), |
| 510 dest_stride, reg8, reg10, reg12, reg14); |
| 511 SRARI_H_4VECS_SH(reg3, reg13, reg11, reg5, reg3, reg13, reg11, reg5, 6); |
| 512 VP9_ADDBLK_CLIP_AND_STORE_8_BYTES_4((dest + (8 * dest_stride)), |
| 513 dest_stride, reg3, reg13, reg11, reg5); |
| 514 SRARI_H_4VECS_SH(reg7, reg9, reg1, reg15, reg7, reg9, reg1, reg15, 6); |
| 515 VP9_ADDBLK_CLIP_AND_STORE_8_BYTES_4((dest + (12 * dest_stride)), |
| 516 dest_stride, reg7, reg9, reg1, reg15); |
| 517 } |
| 518 |
| 519 void vp9_idct16x16_256_add_msa(const int16_t *input, uint8_t *dest, |
| 520 int32_t dest_stride) { |
| 521 int32_t i; |
| 522 DECLARE_ALIGNED(32, int16_t, out_arr[16 * 16]); |
| 523 int16_t *out = out_arr; |
| 524 |
| 525 /* transform rows */ |
| 526 for (i = 0; i < 2; ++i) { |
| 527 /* process 16 * 8 block */ |
| 528 vp9_idct16_1d_rows_msa((input + (i << 7)), (out + (i << 7))); |
| 529 } |
| 530 |
| 531 /* transform columns */ |
| 532 for (i = 0; i < 2; ++i) { |
| 533 /* process 8 * 16 block */ |
| 534 vp9_idct16_1d_columns_addblk_msa((out + (i << 3)), (dest + (i << 3)), |
| 535 dest_stride); |
| 536 } |
| 537 } |
| 538 |
| 539 void vp9_idct16x16_10_add_msa(const int16_t *input, uint8_t *dest, |
| 540 int32_t dest_stride) { |
| 541 uint8_t i; |
| 542 DECLARE_ALIGNED(32, int16_t, out_arr[16 * 16]); |
| 543 int16_t *out = out_arr; |
| 544 |
| 545 /* process 16 * 8 block */ |
| 546 vp9_idct16_1d_rows_msa(input, out); |
| 547 |
| 548 /* short case just considers top 4 rows as valid output */ |
| 549 out += 4 * 16; |
| 550 for (i = 12; i--;) { |
| 551 __asm__ __volatile__ ( |
| 552 "sw $zero, 0(%[out]) \n\t" |
| 553 "sw $zero, 4(%[out]) \n\t" |
| 554 "sw $zero, 8(%[out]) \n\t" |
| 555 "sw $zero, 12(%[out]) \n\t" |
| 556 "sw $zero, 16(%[out]) \n\t" |
| 557 "sw $zero, 20(%[out]) \n\t" |
| 558 "sw $zero, 24(%[out]) \n\t" |
| 559 "sw $zero, 28(%[out]) \n\t" |
| 560 |
| 561 : |
| 562 : [out] "r" (out) |
| 563 ); |
| 564 |
| 565 out += 16; |
| 566 } |
| 567 |
| 568 out = out_arr; |
| 569 |
| 570 /* transform columns */ |
| 571 for (i = 0; i < 2; ++i) { |
| 572 /* process 8 * 16 block */ |
| 573 vp9_idct16_1d_columns_addblk_msa((out + (i << 3)), (dest + (i << 3)), |
| 574 dest_stride); |
| 575 } |
| 576 } |
| 577 |
| 578 void vp9_idct16x16_1_add_msa(const int16_t *input, uint8_t *dest, |
| 579 int32_t dest_stride) { |
| 580 uint8_t i; |
| 581 int32_t const1; |
| 582 int16_t out; |
| 583 v8i16 const2, res0, res1, res2, res3, res4, res5, res6, res7; |
| 584 v16u8 dest0, dest1, dest2, dest3; |
| 585 v16u8 tmp0, tmp1, tmp2, tmp3; |
| 586 v16i8 zero = { 0 }; |
| 587 |
| 588 out = dct_const_round_shift(input[0] * cospi_16_64); |
| 589 out = dct_const_round_shift(out * cospi_16_64); |
| 590 const1 = ROUND_POWER_OF_TWO(out, 6); |
| 591 |
| 592 const2 = __msa_fill_h(const1); |
| 593 |
| 594 for (i = 0; i < 4; ++i) { |
| 595 LOAD_4VECS_UB(dest, dest_stride, dest0, dest1, dest2, dest3); |
| 596 |
| 597 res0 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest0); |
| 598 res1 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest1); |
| 599 res2 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest2); |
| 600 res3 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest3); |
| 601 res4 = (v8i16)__msa_ilvl_b(zero, (v16i8)dest0); |
| 602 res5 = (v8i16)__msa_ilvl_b(zero, (v16i8)dest1); |
| 603 res6 = (v8i16)__msa_ilvl_b(zero, (v16i8)dest2); |
| 604 res7 = (v8i16)__msa_ilvl_b(zero, (v16i8)dest3); |
| 605 |
| 606 res0 += const2; |
| 607 res1 += const2; |
| 608 res2 += const2; |
| 609 res3 += const2; |
| 610 res4 += const2; |
| 611 res5 += const2; |
| 612 res6 += const2; |
| 613 res7 += const2; |
| 614 |
| 615 res0 = CLIP_UNSIGNED_CHAR_H(res0); |
| 616 res1 = CLIP_UNSIGNED_CHAR_H(res1); |
| 617 res2 = CLIP_UNSIGNED_CHAR_H(res2); |
| 618 res3 = CLIP_UNSIGNED_CHAR_H(res3); |
| 619 res4 = CLIP_UNSIGNED_CHAR_H(res4); |
| 620 res5 = CLIP_UNSIGNED_CHAR_H(res5); |
| 621 res6 = CLIP_UNSIGNED_CHAR_H(res6); |
| 622 res7 = CLIP_UNSIGNED_CHAR_H(res7); |
| 623 |
| 624 tmp0 = (v16u8)__msa_pckev_b((v16i8)res4, (v16i8)res0); |
| 625 tmp1 = (v16u8)__msa_pckev_b((v16i8)res5, (v16i8)res1); |
| 626 tmp2 = (v16u8)__msa_pckev_b((v16i8)res6, (v16i8)res2); |
| 627 tmp3 = (v16u8)__msa_pckev_b((v16i8)res7, (v16i8)res3); |
| 628 |
| 629 STORE_4VECS_UB(dest, dest_stride, tmp0, tmp1, tmp2, tmp3); |
| 630 dest += (4 * dest_stride); |
| 631 } |
| 632 } |
| 633 |
| 634 static void vp9_iadst16_1d_rows_msa(const int16_t *input, int16_t *output) { |
| 635 v8i16 r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11, r12, r13, r14, r15; |
| 636 v8i16 l0, l1, l2, l3, l4, l5, l6, l7, l8, l9, l10, l11, l12, l13, l14, l15; |
| 637 |
| 638 /* load input data */ |
| 639 LOAD_16VECS_SH(input, 8, |
| 640 l0, l8, l1, l9, l2, l10, l3, l11, |
| 641 l4, l12, l5, l13, l6, l14, l7, l15); |
| 642 |
| 643 TRANSPOSE8x8_H_SH(l0, l1, l2, l3, l4, l5, l6, l7, |
| 644 l0, l1, l2, l3, l4, l5, l6, l7); |
| 645 |
| 646 TRANSPOSE8x8_H_SH(l8, l9, l10, l11, l12, l13, l14, l15, |
| 647 l8, l9, l10, l11, l12, l13, l14, l15); |
| 648 |
| 649 /* ADST in horizontal */ |
| 650 VP9_IADST8x16_1D(l0, l1, l2, l3, l4, l5, l6, l7, |
| 651 l8, l9, l10, l11, l12, l13, l14, l15, |
| 652 r0, r1, r2, r3, r4, r5, r6, r7, |
| 653 r8, r9, r10, r11, r12, r13, r14, r15); |
| 654 |
| 655 l1 = -r8; |
| 656 l3 = -r4; |
| 657 l13 = -r13; |
| 658 l15 = -r1; |
| 659 |
| 660 TRANSPOSE8x8_H_SH(r0, l1, r12, l3, r6, r14, r10, r2, |
| 661 l0, l1, l2, l3, l4, l5, l6, l7); |
| 662 |
| 663 STORE_8VECS_SH(output, 16, l0, l1, l2, l3, l4, l5, l6, l7); |
| 664 |
| 665 TRANSPOSE8x8_H_SH(r3, r11, r15, r7, r5, l13, r9, l15, |
| 666 l8, l9, l10, l11, l12, l13, l14, l15); |
| 667 |
| 668 STORE_8VECS_SH((output + 8), 16, l8, l9, l10, l11, l12, l13, l14, l15); |
| 669 } |
| 670 |
| 671 static void vp9_iadst16_1d_columns_addblk_msa(int16_t *input, uint8_t *dest, |
| 672 int32_t dest_stride) { |
| 673 v8i16 v0, v2, v4, v6, k0, k1, k2, k3; |
| 674 v8i16 r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11, r12, r13, r14, r15; |
| 675 v8i16 out0, out1, out2, out3, out4, out5, out6, out7; |
| 676 v8i16 out8, out9, out10, out11, out12, out13, out14, out15; |
| 677 v8i16 g0, g1, g2, g3, g4, g5, g6, g7, g8, g9, g10, g11, g12, g13, g14, g15; |
| 678 v8i16 h0, h1, h2, h3, h4, h5, h6, h7, h8, h9, h10, h11; |
| 679 v8i16 res0, res1, res2, res3, res4, res5, res6, res7; |
| 680 v8i16 res8, res9, res10, res11, res12, res13, res14, res15; |
| 681 v16u8 dest0, dest1, dest2, dest3, dest4, dest5, dest6, dest7; |
| 682 v16u8 dest8, dest9, dest10, dest11, dest12, dest13, dest14, dest15; |
| 683 v16i8 zero = { 0 }; |
| 684 |
| 685 r0 = LOAD_SH(input + 0 * 16); |
| 686 r3 = LOAD_SH(input + 3 * 16); |
| 687 r4 = LOAD_SH(input + 4 * 16); |
| 688 r7 = LOAD_SH(input + 7 * 16); |
| 689 r8 = LOAD_SH(input + 8 * 16); |
| 690 r11 = LOAD_SH(input + 11 * 16); |
| 691 r12 = LOAD_SH(input + 12 * 16); |
| 692 r15 = LOAD_SH(input + 15 * 16); |
| 693 |
| 694 /* stage 1 */ |
| 695 k0 = SET_COSPI_PAIR(cospi_1_64, cospi_31_64); |
| 696 k1 = SET_COSPI_PAIR(cospi_31_64, -cospi_1_64); |
| 697 k2 = SET_COSPI_PAIR(cospi_17_64, cospi_15_64); |
| 698 k3 = SET_COSPI_PAIR(cospi_15_64, -cospi_17_64); |
| 699 VP9_MADD_BF(r15, r0, r7, r8, k0, k1, k2, k3, g0, g1, g2, g3); |
| 700 |
| 701 k0 = SET_COSPI_PAIR(cospi_9_64, cospi_23_64); |
| 702 k1 = SET_COSPI_PAIR(cospi_23_64, -cospi_9_64); |
| 703 k2 = SET_COSPI_PAIR(cospi_25_64, cospi_7_64); |
| 704 k3 = SET_COSPI_PAIR(cospi_7_64, -cospi_25_64); |
| 705 VP9_MADD_BF(r11, r4, r3, r12, k0, k1, k2, k3, g8, g9, g10, g11); |
| 706 |
| 707 BUTTERFLY_4(g0, g2, g10, g8, h8, h9, v2, v0); |
| 708 |
| 709 k0 = SET_COSPI_PAIR(cospi_4_64, cospi_28_64); |
| 710 k1 = SET_COSPI_PAIR(cospi_28_64, -cospi_4_64); |
| 711 k2 = SET_COSPI_PAIR(-cospi_28_64, cospi_4_64); |
| 712 VP9_MADD_BF(g1, g3, g9, g11, k0, k1, k2, k0, h0, h1, h2, h3); |
| 713 |
| 714 r1 = LOAD_SH(input + 1 * 16); |
| 715 r2 = LOAD_SH(input + 2 * 16); |
| 716 r5 = LOAD_SH(input + 5 * 16); |
| 717 r6 = LOAD_SH(input + 6 * 16); |
| 718 r9 = LOAD_SH(input + 9 * 16); |
| 719 r10 = LOAD_SH(input + 10 * 16); |
| 720 r13 = LOAD_SH(input + 13 * 16); |
| 721 r14 = LOAD_SH(input + 14 * 16); |
| 722 |
| 723 k0 = SET_COSPI_PAIR(cospi_5_64, cospi_27_64); |
| 724 k1 = SET_COSPI_PAIR(cospi_27_64, -cospi_5_64); |
| 725 k2 = SET_COSPI_PAIR(cospi_21_64, cospi_11_64); |
| 726 k3 = SET_COSPI_PAIR(cospi_11_64, -cospi_21_64); |
| 727 VP9_MADD_BF(r13, r2, r5, r10, k0, k1, k2, k3, g4, g5, g6, g7); |
| 728 |
| 729 k0 = SET_COSPI_PAIR(cospi_13_64, cospi_19_64); |
| 730 k1 = SET_COSPI_PAIR(cospi_19_64, -cospi_13_64); |
| 731 k2 = SET_COSPI_PAIR(cospi_29_64, cospi_3_64); |
| 732 k3 = SET_COSPI_PAIR(cospi_3_64, -cospi_29_64); |
| 733 VP9_MADD_BF(r9, r6, r1, r14, k0, k1, k2, k3, g12, g13, g14, g15); |
| 734 |
| 735 BUTTERFLY_4(g4, g6, g14, g12, h10, h11, v6, v4); |
| 736 |
| 737 BUTTERFLY_4(h8, h9, h11, h10, out0, out1, h11, h10); |
| 738 out1 = -out1; |
| 739 out0 = __msa_srari_h(out0, 6); |
| 740 out1 = __msa_srari_h(out1, 6); |
| 741 dest0 = LOAD_UB(dest + 0 * dest_stride); |
| 742 dest1 = LOAD_UB(dest + 15 * dest_stride); |
| 743 res0 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest0); |
| 744 res1 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest1); |
| 745 res0 += out0; |
| 746 res1 += out1; |
| 747 res0 = CLIP_UNSIGNED_CHAR_H(res0); |
| 748 res1 = CLIP_UNSIGNED_CHAR_H(res1); |
| 749 res0 = (v8i16)__msa_pckev_b((v16i8)res0, (v16i8)res0); |
| 750 res1 = (v8i16)__msa_pckev_b((v16i8)res1, (v16i8)res1); |
| 751 STORE_DWORD(dest, __msa_copy_u_d((v2i64)res0, 0)); |
| 752 STORE_DWORD(dest + 15 * dest_stride, __msa_copy_u_d((v2i64)res1, 0)); |
| 753 |
| 754 k0 = SET_COSPI_PAIR(cospi_12_64, cospi_20_64); |
| 755 k1 = SET_COSPI_PAIR(-cospi_20_64, cospi_12_64); |
| 756 k2 = SET_COSPI_PAIR(cospi_20_64, -cospi_12_64); |
| 757 VP9_MADD_BF(g7, g5, g15, g13, k0, k1, k2, k0, h4, h5, h6, h7); |
| 758 |
| 759 BUTTERFLY_4(h0, h2, h6, h4, out8, out9, out11, out10); |
| 760 out8 = -out8; |
| 761 |
| 762 out8 = __msa_srari_h(out8, 6); |
| 763 out9 = __msa_srari_h(out9, 6); |
| 764 dest8 = LOAD_UB(dest + 1 * dest_stride); |
| 765 dest9 = LOAD_UB(dest + 14 * dest_stride); |
| 766 res8 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest8); |
| 767 res9 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest9); |
| 768 res8 += out8; |
| 769 res9 += out9; |
| 770 res8 = CLIP_UNSIGNED_CHAR_H(res8); |
| 771 res9 = CLIP_UNSIGNED_CHAR_H(res9); |
| 772 res8 = (v8i16)__msa_pckev_b((v16i8)res8, (v16i8)res8); |
| 773 res9 = (v8i16)__msa_pckev_b((v16i8)res9, (v16i8)res9); |
| 774 STORE_DWORD(dest + dest_stride, __msa_copy_u_d((v2i64)res8, 0)); |
| 775 STORE_DWORD(dest + 14 * dest_stride, __msa_copy_u_d((v2i64)res9, 0)); |
| 776 |
| 777 k0 = SET_COSPI_PAIR(cospi_8_64, cospi_24_64); |
| 778 k1 = SET_COSPI_PAIR(cospi_24_64, -cospi_8_64); |
| 779 k2 = SET_COSPI_PAIR(-cospi_24_64, cospi_8_64); |
| 780 VP9_MADD_BF(v0, v2, v4, v6, k0, k1, k2, k0, out4, out6, out5, out7); |
| 781 out4 = -out4; |
| 782 out4 = __msa_srari_h(out4, 6); |
| 783 out5 = __msa_srari_h(out5, 6); |
| 784 dest4 = LOAD_UB(dest + 3 * dest_stride); |
| 785 dest5 = LOAD_UB(dest + 12 * dest_stride); |
| 786 res4 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest4); |
| 787 res5 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest5); |
| 788 res4 += out4; |
| 789 res5 += out5; |
| 790 res4 = CLIP_UNSIGNED_CHAR_H(res4); |
| 791 res5 = CLIP_UNSIGNED_CHAR_H(res5); |
| 792 res4 = (v8i16)__msa_pckev_b((v16i8)res4, (v16i8)res4); |
| 793 res5 = (v8i16)__msa_pckev_b((v16i8)res5, (v16i8)res5); |
| 794 STORE_DWORD(dest + 3 * dest_stride, __msa_copy_u_d((v2i64)res4, 0)); |
| 795 STORE_DWORD(dest + 12 * dest_stride, __msa_copy_u_d((v2i64)res5, 0)); |
| 796 |
| 797 VP9_MADD_BF(h1, h3, h5, h7, k0, k1, k2, k0, out12, out14, out13, out15); |
| 798 out13 = -out13; |
| 799 out12 = __msa_srari_h(out12, 6); |
| 800 out13 = __msa_srari_h(out13, 6); |
| 801 dest12 = LOAD_UB(dest + 2 * dest_stride); |
| 802 dest13 = LOAD_UB(dest + 13 * dest_stride); |
| 803 res12 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest12); |
| 804 res13 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest13); |
| 805 res12 += out12; |
| 806 res13 += out13; |
| 807 res12 = CLIP_UNSIGNED_CHAR_H(res12); |
| 808 res13 = CLIP_UNSIGNED_CHAR_H(res13); |
| 809 res12 = (v8i16)__msa_pckev_b((v16i8)res12, (v16i8)res12); |
| 810 res13 = (v8i16)__msa_pckev_b((v16i8)res13, (v16i8)res13); |
| 811 STORE_DWORD(dest + 2 * dest_stride, __msa_copy_u_d((v2i64)res12, 0)); |
| 812 STORE_DWORD(dest + 13 * dest_stride, __msa_copy_u_d((v2i64)res13, 0)); |
| 813 |
| 814 k0 = SET_COSPI_PAIR(cospi_16_64, cospi_16_64); |
| 815 k3 = SET_COSPI_PAIR(-cospi_16_64, cospi_16_64); |
| 816 VP9_MADD_SHORT(out6, out7, k0, k3, out6, out7); |
| 817 out6 = __msa_srari_h(out6, 6); |
| 818 out7 = __msa_srari_h(out7, 6); |
| 819 dest6 = LOAD_UB(dest + 4 * dest_stride); |
| 820 dest7 = LOAD_UB(dest + 11 * dest_stride); |
| 821 res6 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest6); |
| 822 res7 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest7); |
| 823 res6 += out6; |
| 824 res7 += out7; |
| 825 res6 = CLIP_UNSIGNED_CHAR_H(res6); |
| 826 res7 = CLIP_UNSIGNED_CHAR_H(res7); |
| 827 res6 = (v8i16)__msa_pckev_b((v16i8)res6, (v16i8)res6); |
| 828 res7 = (v8i16)__msa_pckev_b((v16i8)res7, (v16i8)res7); |
| 829 STORE_DWORD(dest + 4 * dest_stride, __msa_copy_u_d((v2i64)res6, 0)); |
| 830 STORE_DWORD(dest + 11 * dest_stride, __msa_copy_u_d((v2i64)res7, 0)); |
| 831 |
| 832 VP9_MADD_SHORT(out10, out11, k0, k3, out10, out11); |
| 833 out10 = __msa_srari_h(out10, 6); |
| 834 out11 = __msa_srari_h(out11, 6); |
| 835 dest10 = LOAD_UB(dest + 6 * dest_stride); |
| 836 dest11 = LOAD_UB(dest + 9 * dest_stride); |
| 837 res10 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest10); |
| 838 res11 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest11); |
| 839 res10 += out10; |
| 840 res11 += out11; |
| 841 res10 = CLIP_UNSIGNED_CHAR_H(res10); |
| 842 res11 = CLIP_UNSIGNED_CHAR_H(res11); |
| 843 res10 = (v8i16)__msa_pckev_b((v16i8)res10, (v16i8)res10); |
| 844 res11 = (v8i16)__msa_pckev_b((v16i8)res11, (v16i8)res11); |
| 845 STORE_DWORD(dest + 6 * dest_stride, __msa_copy_u_d((v2i64)res10, 0)); |
| 846 STORE_DWORD(dest + 9 * dest_stride, __msa_copy_u_d((v2i64)res11, 0)); |
| 847 |
| 848 k1 = SET_COSPI_PAIR(-cospi_16_64, -cospi_16_64); |
| 849 k2 = SET_COSPI_PAIR(cospi_16_64, -cospi_16_64); |
| 850 VP9_MADD_SHORT(h10, h11, k1, k2, out2, out3); |
| 851 out2 = __msa_srari_h(out2, 6); |
| 852 out3 = __msa_srari_h(out3, 6); |
| 853 dest2 = LOAD_UB(dest + 7 * dest_stride); |
| 854 dest3 = LOAD_UB(dest + 8 * dest_stride); |
| 855 res2 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest2); |
| 856 res3 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest3); |
| 857 res2 += out2; |
| 858 res3 += out3; |
| 859 res2 = CLIP_UNSIGNED_CHAR_H(res2); |
| 860 res3 = CLIP_UNSIGNED_CHAR_H(res3); |
| 861 res2 = (v8i16)__msa_pckev_b((v16i8)res2, (v16i8)res2); |
| 862 res3 = (v8i16)__msa_pckev_b((v16i8)res3, (v16i8)res3); |
| 863 STORE_DWORD(dest + 7 * dest_stride, __msa_copy_u_d((v2i64)res2, 0)); |
| 864 STORE_DWORD(dest + 8 * dest_stride, __msa_copy_u_d((v2i64)res3, 0)); |
| 865 |
| 866 VP9_MADD_SHORT(out14, out15, k1, k2, out14, out15); |
| 867 out14 = __msa_srari_h(out14, 6); |
| 868 out15 = __msa_srari_h(out15, 6); |
| 869 dest14 = LOAD_UB(dest + 5 * dest_stride); |
| 870 dest15 = LOAD_UB(dest + 10 * dest_stride); |
| 871 res14 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest14); |
| 872 res15 = (v8i16)__msa_ilvr_b(zero, (v16i8)dest15); |
| 873 res14 += out14; |
| 874 res15 += out15; |
| 875 res14 = CLIP_UNSIGNED_CHAR_H(res14); |
| 876 res15 = CLIP_UNSIGNED_CHAR_H(res15); |
| 877 res14 = (v8i16)__msa_pckev_b((v16i8)res14, (v16i8)res14); |
| 878 res15 = (v8i16)__msa_pckev_b((v16i8)res15, (v16i8)res15); |
| 879 STORE_DWORD(dest + 5 * dest_stride, __msa_copy_u_d((v2i64)res14, 0)); |
| 880 STORE_DWORD(dest + 10 * dest_stride, __msa_copy_u_d((v2i64)res15, 0)); |
| 881 } |
| 882 |
| 883 void vp9_iht16x16_256_add_msa(const int16_t *input, uint8_t *dest, |
| 884 int32_t dest_stride, int32_t tx_type) { |
| 885 int32_t i; |
| 886 DECLARE_ALIGNED(32, int16_t, out[16 * 16]); |
| 887 int16_t *out_ptr = &out[0]; |
| 888 |
| 889 switch (tx_type) { |
| 890 case DCT_DCT: |
| 891 /* transform rows */ |
| 892 for (i = 0; i < 2; ++i) { |
| 893 /* process 16 * 8 block */ |
| 894 vp9_idct16_1d_rows_msa((input + (i << 7)), (out_ptr + (i << 7))); |
| 895 } |
| 896 |
| 897 /* transform columns */ |
| 898 for (i = 0; i < 2; ++i) { |
| 899 /* process 8 * 16 block */ |
| 900 vp9_idct16_1d_columns_addblk_msa((out_ptr + (i << 3)), |
| 901 (dest + (i << 3)), dest_stride); |
| 902 } |
| 903 break; |
| 904 case ADST_DCT: |
| 905 /* transform rows */ |
| 906 for (i = 0; i < 2; ++i) { |
| 907 /* process 16 * 8 block */ |
| 908 vp9_idct16_1d_rows_msa((input + (i << 7)), (out_ptr + (i << 7))); |
| 909 } |
| 910 |
| 911 /* transform columns */ |
| 912 for (i = 0; i < 2; ++i) { |
| 913 vp9_iadst16_1d_columns_addblk_msa((out_ptr + (i << 3)), |
| 914 (dest + (i << 3)), dest_stride); |
| 915 } |
| 916 break; |
| 917 case DCT_ADST: |
| 918 /* transform rows */ |
| 919 for (i = 0; i < 2; ++i) { |
| 920 /* process 16 * 8 block */ |
| 921 vp9_iadst16_1d_rows_msa((input + (i << 7)), (out_ptr + (i << 7))); |
| 922 } |
| 923 |
| 924 /* transform columns */ |
| 925 for (i = 0; i < 2; ++i) { |
| 926 /* process 8 * 16 block */ |
| 927 vp9_idct16_1d_columns_addblk_msa((out_ptr + (i << 3)), |
| 928 (dest + (i << 3)), dest_stride); |
| 929 } |
| 930 break; |
| 931 case ADST_ADST: |
| 932 /* transform rows */ |
| 933 for (i = 0; i < 2; ++i) { |
| 934 /* process 16 * 8 block */ |
| 935 vp9_iadst16_1d_rows_msa((input + (i << 7)), (out_ptr + (i << 7))); |
| 936 } |
| 937 |
| 938 /* transform columns */ |
| 939 for (i = 0; i < 2; ++i) { |
| 940 vp9_iadst16_1d_columns_addblk_msa((out_ptr + (i << 3)), |
| 941 (dest + (i << 3)), dest_stride); |
| 942 } |
| 943 break; |
| 944 default: |
| 945 assert(0); |
| 946 break; |
| 947 } |
| 948 } |
OLD | NEW |