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
|
#ifndef AV1_TXFM_SSE4_H_
#define AV1_TXFM_SSE4_H_
#include <smmintrin.h>
#ifdef __cplusplus
extern "C" {
#endif
static INLINE __m128i av1_round_shift_32_sse4_1(__m128i vec, int bit) {
__m128i tmp, round;
round = _mm_set1_epi32(1 << (bit - 1));
tmp = _mm_add_epi32(vec, round);
return _mm_srai_epi32(tmp, bit);
}
static INLINE void av1_round_shift_array_32_sse4_1(__m128i *input,
__m128i *output,
const int size,
const int bit) {
if (bit > 0) {
int i;
for (i = 0; i < size; i++) {
output[i] = av1_round_shift_32_sse4_1(input[i], bit);
}
} else {
int i;
for (i = 0; i < size; i++) {
output[i] = _mm_slli_epi32(input[i], -bit);
}
}
}
static INLINE void av1_round_shift_rect_array_32_sse4_1(__m128i *input,
__m128i *output,
const int size,
const int bit) {
const __m128i sqrt2 = _mm_set1_epi32(NewSqrt2);
if (bit > 0) {
int i;
for (i = 0; i < size; i++) {
const __m128i r0 = av1_round_shift_32_sse4_1(input[i], bit);
const __m128i r1 = _mm_mullo_epi32(sqrt2, r0);
output[i] = av1_round_shift_32_sse4_1(r1, NewSqrt2Bits);
}
} else {
int i;
for (i = 0; i < size; i++) {
const __m128i r0 = _mm_slli_epi32(input[i], -bit);
const __m128i r1 = _mm_mullo_epi32(sqrt2, r0);
output[i] = av1_round_shift_32_sse4_1(r1, NewSqrt2Bits);
}
}
}
#ifdef __cplusplus
}
#endif
#endif // AV1_TXFM_SSE4_H_
|
