1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
/*
* Copyright (c) 2022 The WebM project authors. All Rights Reserved.
*
* Use of this source code is governed by a BSD-style license
* that can be found in the LICENSE file in the root of the source
* tree. An additional intellectual property rights grant can be found
* in the file PATENTS. All contributing project authors may
* be found in the AUTHORS file in the root of the source tree.
*/
#include <stdlib.h>
#include "./vpx_dsp_rtcd.h"
#include "vpx_dsp/loongarch/bitdepth_conversion_lsx.h"
void vpx_hadamard_8x8_lsx(const int16_t *src, ptrdiff_t src_stride,
tran_low_t *dst) {
__m128i src0, src1, src2, src3, src4, src5, src6, src7;
__m128i tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
ptrdiff_t src_stride2 = src_stride << 1;
ptrdiff_t src_stride3 = src_stride2 + src_stride;
ptrdiff_t src_stride4 = src_stride2 << 1;
ptrdiff_t src_stride6 = src_stride3 << 1;
int16_t *src_tmp = (int16_t *)src;
src0 = __lsx_vld(src_tmp, 0);
DUP2_ARG2(__lsx_vldx, src_tmp, src_stride2, src_tmp, src_stride4, src1, src2);
src3 = __lsx_vldx(src_tmp, src_stride6);
src_tmp += src_stride4;
src4 = __lsx_vld(src_tmp, 0);
DUP2_ARG2(__lsx_vldx, src_tmp, src_stride2, src_tmp, src_stride4, src5, src6);
src7 = __lsx_vldx(src_tmp, src_stride6);
LSX_BUTTERFLY_8_H(src0, src2, src4, src6, src7, src5, src3, src1, tmp0, tmp2,
tmp4, tmp6, tmp7, tmp5, tmp3, tmp1);
LSX_BUTTERFLY_8_H(tmp0, tmp1, tmp4, tmp5, tmp7, tmp6, tmp3, tmp2, src0, src1,
src4, src5, src7, src6, src3, src2);
LSX_BUTTERFLY_8_H(src0, src1, src2, src3, src7, src6, src5, src4, tmp0, tmp7,
tmp3, tmp4, tmp5, tmp1, tmp6, tmp2);
LSX_TRANSPOSE8x8_H(tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7, src0, src1,
src2, src3, src4, src5, src6, src7);
LSX_BUTTERFLY_8_H(src0, src2, src4, src6, src7, src5, src3, src1, tmp0, tmp2,
tmp4, tmp6, tmp7, tmp5, tmp3, tmp1);
LSX_BUTTERFLY_8_H(tmp0, tmp1, tmp4, tmp5, tmp7, tmp6, tmp3, tmp2, src0, src1,
src4, src5, src7, src6, src3, src2);
LSX_BUTTERFLY_8_H(src0, src1, src2, src3, src7, src6, src5, src4, tmp0, tmp7,
tmp3, tmp4, tmp5, tmp1, tmp6, tmp2);
store_tran_low(tmp0, dst, 0);
store_tran_low(tmp1, dst, 8);
store_tran_low(tmp2, dst, 16);
store_tran_low(tmp3, dst, 24);
store_tran_low(tmp4, dst, 32);
store_tran_low(tmp5, dst, 40);
store_tran_low(tmp6, dst, 48);
store_tran_low(tmp7, dst, 56);
}
void vpx_hadamard_16x16_lsx(const int16_t *src, ptrdiff_t src_stride,
tran_low_t *dst) {
int i;
__m128i a0, a1, a2, a3, b0, b1, b2, b3;
/* Rearrange 16x16 to 8x32 and remove stride.
* Top left first. */
vpx_hadamard_8x8_lsx(src + 0 + 0 * src_stride, src_stride, dst + 0);
/* Top right. */
vpx_hadamard_8x8_lsx(src + 8 + 0 * src_stride, src_stride, dst + 64);
/* Bottom left. */
vpx_hadamard_8x8_lsx(src + 0 + 8 * src_stride, src_stride, dst + 128);
/* Bottom right. */
vpx_hadamard_8x8_lsx(src + 8 + 8 * src_stride, src_stride, dst + 192);
for (i = 0; i < 64; i += 8) {
a0 = load_tran_low(dst);
a1 = load_tran_low(dst + 64);
a2 = load_tran_low(dst + 128);
a3 = load_tran_low(dst + 192);
LSX_BUTTERFLY_4_H(a0, a2, a3, a1, b0, b2, b3, b1);
DUP4_ARG2(__lsx_vsrai_h, b0, 1, b1, 1, b2, 1, b3, 1, b0, b1, b2, b3);
LSX_BUTTERFLY_4_H(b0, b1, b3, b2, a0, a1, a3, a2);
store_tran_low(a0, dst, 0);
store_tran_low(a1, dst, 64);
store_tran_low(a2, dst, 128);
store_tran_low(a3, dst, 192);
dst += 8;
}
}
|
