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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
|
/*
* Copyright 2014 The Android Open Source Project
*
* Use of this source code is governed by a BSD-style license that can be
* found in the LICENSE file.
*/
#ifndef SkColor_opts_SSE2_DEFINED
#define SkColor_opts_SSE2_DEFINED
#include <emmintrin.h>
#define ASSERT_EQ(a,b) SkASSERT(0xffff == _mm_movemask_epi8(_mm_cmpeq_epi8((a), (b))))
// Because no _mm_mul_epi32() in SSE2, we emulate it here.
// Multiplies 4 32-bit integers from a by 4 32-bit intergers from b.
// The 4 multiplication results should be represented within 32-bit
// integers, otherwise they would be overflow.
static inline __m128i Multiply32_SSE2(const __m128i& a, const __m128i& b) {
// Calculate results of a0 * b0 and a2 * b2.
__m128i r1 = _mm_mul_epu32(a, b);
// Calculate results of a1 * b1 and a3 * b3.
__m128i r2 = _mm_mul_epu32(_mm_srli_si128(a, 4), _mm_srli_si128(b, 4));
// Shuffle results to [63..0] and interleave the results.
__m128i r = _mm_unpacklo_epi32(_mm_shuffle_epi32(r1, _MM_SHUFFLE(0,0,2,0)),
_mm_shuffle_epi32(r2, _MM_SHUFFLE(0,0,2,0)));
return r;
}
static inline __m128i SkAlpha255To256_SSE2(const __m128i& alpha) {
return _mm_add_epi32(alpha, _mm_set1_epi32(1));
}
// See #define SkAlphaMulAlpha(a, b) SkMulDiv255Round(a, b) in SkXfermode.cpp.
static inline __m128i SkAlphaMulAlpha_SSE2(const __m128i& a,
const __m128i& b) {
__m128i prod = _mm_mullo_epi16(a, b);
prod = _mm_add_epi32(prod, _mm_set1_epi32(128));
prod = _mm_add_epi32(prod, _mm_srli_epi32(prod, 8));
prod = _mm_srli_epi32(prod, 8);
return prod;
}
// Portable version SkAlphaMulQ is in SkColorPriv.h.
static inline __m128i SkAlphaMulQ_SSE2(const __m128i& c, const __m128i& scale) {
const __m128i mask = _mm_set1_epi32(0xFF00FF);
__m128i s = _mm_or_si128(_mm_slli_epi32(scale, 16), scale);
// uint32_t rb = ((c & mask) * scale) >> 8
__m128i rb = _mm_and_si128(mask, c);
rb = _mm_mullo_epi16(rb, s);
rb = _mm_srli_epi16(rb, 8);
// uint32_t ag = ((c >> 8) & mask) * scale
__m128i ag = _mm_srli_epi16(c, 8);
ASSERT_EQ(ag, _mm_and_si128(mask, ag)); // ag = _mm_srli_epi16(c, 8) did this for us.
ag = _mm_mullo_epi16(ag, s);
// (rb & mask) | (ag & ~mask)
ASSERT_EQ(rb, _mm_and_si128(mask, rb)); // rb = _mm_srli_epi16(rb, 8) did this for us.
ag = _mm_andnot_si128(mask, ag);
return _mm_or_si128(rb, ag);
}
// Fast path for SkAlphaMulQ_SSE2 with a constant scale factor.
static inline __m128i SkAlphaMulQ_SSE2(const __m128i& c, const unsigned scale) {
const __m128i mask = _mm_set1_epi32(0xFF00FF);
__m128i s = _mm_set1_epi16(scale << 8); // Move scale factor to upper byte of word.
// With mulhi, red and blue values are already in the right place and
// don't need to be divided by 256.
__m128i rb = _mm_and_si128(mask, c);
rb = _mm_mulhi_epu16(rb, s);
__m128i ag = _mm_andnot_si128(mask, c);
ag = _mm_mulhi_epu16(ag, s); // Alpha and green values are in the higher byte of each word.
ag = _mm_andnot_si128(mask, ag);
return _mm_or_si128(rb, ag);
}
static inline __m128i SkGetPackedA32_SSE2(const __m128i& src) {
#if SK_A32_SHIFT == 24 // It's very common (universal?) that alpha is the top byte.
return _mm_srli_epi32(src, 24); // You'd hope the compiler would remove the left shift then,
#else // but I've seen Clang just do a dumb left shift of zero. :(
__m128i a = _mm_slli_epi32(src, (24 - SK_A32_SHIFT));
return _mm_srli_epi32(a, 24);
#endif
}
static inline __m128i SkGetPackedR32_SSE2(const __m128i& src) {
__m128i r = _mm_slli_epi32(src, (24 - SK_R32_SHIFT));
return _mm_srli_epi32(r, 24);
}
static inline __m128i SkGetPackedG32_SSE2(const __m128i& src) {
__m128i g = _mm_slli_epi32(src, (24 - SK_G32_SHIFT));
return _mm_srli_epi32(g, 24);
}
static inline __m128i SkGetPackedB32_SSE2(const __m128i& src) {
__m128i b = _mm_slli_epi32(src, (24 - SK_B32_SHIFT));
return _mm_srli_epi32(b, 24);
}
static inline __m128i SkMul16ShiftRound_SSE2(const __m128i& a,
const __m128i& b, int shift) {
__m128i prod = _mm_mullo_epi16(a, b);
prod = _mm_add_epi16(prod, _mm_set1_epi16(1 << (shift - 1)));
prod = _mm_add_epi16(prod, _mm_srli_epi16(prod, shift));
prod = _mm_srli_epi16(prod, shift);
return prod;
}
static inline __m128i SkPackRGB16_SSE2(const __m128i& r,
const __m128i& g, const __m128i& b) {
__m128i dr = _mm_slli_epi16(r, SK_R16_SHIFT);
__m128i dg = _mm_slli_epi16(g, SK_G16_SHIFT);
__m128i db = _mm_slli_epi16(b, SK_B16_SHIFT);
__m128i c = _mm_or_si128(dr, dg);
return _mm_or_si128(c, db);
}
static inline __m128i SkPackARGB32_SSE2(const __m128i& a, const __m128i& r,
const __m128i& g, const __m128i& b) {
__m128i da = _mm_slli_epi32(a, SK_A32_SHIFT);
__m128i dr = _mm_slli_epi32(r, SK_R32_SHIFT);
__m128i dg = _mm_slli_epi32(g, SK_G32_SHIFT);
__m128i db = _mm_slli_epi32(b, SK_B32_SHIFT);
__m128i c = _mm_or_si128(da, dr);
c = _mm_or_si128(c, dg);
return _mm_or_si128(c, db);
}
static inline __m128i SkPacked16ToR32_SSE2(const __m128i& src) {
__m128i r = _mm_srli_epi32(src, SK_R16_SHIFT);
r = _mm_and_si128(r, _mm_set1_epi32(SK_R16_MASK));
r = _mm_or_si128(_mm_slli_epi32(r, (8 - SK_R16_BITS)),
_mm_srli_epi32(r, (2 * SK_R16_BITS - 8)));
return r;
}
static inline __m128i SkPacked16ToG32_SSE2(const __m128i& src) {
__m128i g = _mm_srli_epi32(src, SK_G16_SHIFT);
g = _mm_and_si128(g, _mm_set1_epi32(SK_G16_MASK));
g = _mm_or_si128(_mm_slli_epi32(g, (8 - SK_G16_BITS)),
_mm_srli_epi32(g, (2 * SK_G16_BITS - 8)));
return g;
}
static inline __m128i SkPacked16ToB32_SSE2(const __m128i& src) {
__m128i b = _mm_srli_epi32(src, SK_B16_SHIFT);
b = _mm_and_si128(b, _mm_set1_epi32(SK_B16_MASK));
b = _mm_or_si128(_mm_slli_epi32(b, (8 - SK_B16_BITS)),
_mm_srli_epi32(b, (2 * SK_B16_BITS - 8)));
return b;
}
static inline __m128i SkPixel16ToPixel32_SSE2(const __m128i& src) {
__m128i r = SkPacked16ToR32_SSE2(src);
__m128i g = SkPacked16ToG32_SSE2(src);
__m128i b = SkPacked16ToB32_SSE2(src);
return SkPackARGB32_SSE2(_mm_set1_epi32(0xFF), r, g, b);
}
static inline __m128i SkPixel32ToPixel16_ToU16_SSE2(const __m128i& src_pixel1,
const __m128i& src_pixel2) {
// Calculate result r.
__m128i r1 = _mm_srli_epi32(src_pixel1,
SK_R32_SHIFT + (8 - SK_R16_BITS));
r1 = _mm_and_si128(r1, _mm_set1_epi32(SK_R16_MASK));
__m128i r2 = _mm_srli_epi32(src_pixel2,
SK_R32_SHIFT + (8 - SK_R16_BITS));
r2 = _mm_and_si128(r2, _mm_set1_epi32(SK_R16_MASK));
__m128i r = _mm_packs_epi32(r1, r2);
// Calculate result g.
__m128i g1 = _mm_srli_epi32(src_pixel1,
SK_G32_SHIFT + (8 - SK_G16_BITS));
g1 = _mm_and_si128(g1, _mm_set1_epi32(SK_G16_MASK));
__m128i g2 = _mm_srli_epi32(src_pixel2,
SK_G32_SHIFT + (8 - SK_G16_BITS));
g2 = _mm_and_si128(g2, _mm_set1_epi32(SK_G16_MASK));
__m128i g = _mm_packs_epi32(g1, g2);
// Calculate result b.
__m128i b1 = _mm_srli_epi32(src_pixel1,
SK_B32_SHIFT + (8 - SK_B16_BITS));
b1 = _mm_and_si128(b1, _mm_set1_epi32(SK_B16_MASK));
__m128i b2 = _mm_srli_epi32(src_pixel2,
SK_B32_SHIFT + (8 - SK_B16_BITS));
b2 = _mm_and_si128(b2, _mm_set1_epi32(SK_B16_MASK));
__m128i b = _mm_packs_epi32(b1, b2);
// Store 8 16-bit colors in dst.
__m128i d_pixel = SkPackRGB16_SSE2(r, g, b);
return d_pixel;
}
// Portable version is SkPMSrcOver in SkColorPriv.h.
static inline __m128i SkPMSrcOver_SSE2(const __m128i& src, const __m128i& dst) {
return _mm_add_epi32(src,
SkAlphaMulQ_SSE2(dst, _mm_sub_epi32(_mm_set1_epi32(256),
SkGetPackedA32_SSE2(src))));
}
// Portable version is SkBlendARGB32 in SkColorPriv.h.
static inline __m128i SkBlendARGB32_SSE2(const __m128i& src, const __m128i& dst,
const __m128i& aa) {
__m128i src_scale = SkAlpha255To256_SSE2(aa);
// SkAlpha255To256(255 - SkAlphaMul(SkGetPackedA32(src), src_scale))
__m128i dst_scale = SkGetPackedA32_SSE2(src);
dst_scale = _mm_mullo_epi16(dst_scale, src_scale);
dst_scale = _mm_srli_epi16(dst_scale, 8);
dst_scale = _mm_sub_epi32(_mm_set1_epi32(256), dst_scale);
__m128i result = SkAlphaMulQ_SSE2(src, src_scale);
return _mm_add_epi8(result, SkAlphaMulQ_SSE2(dst, dst_scale));
}
// Fast path for SkBlendARGB32_SSE2 with a constant alpha factor.
static inline __m128i SkBlendARGB32_SSE2(const __m128i& src, const __m128i& dst,
const unsigned aa) {
unsigned alpha = SkAlpha255To256(aa);
__m128i src_scale = _mm_set1_epi32(alpha);
// SkAlpha255To256(255 - SkAlphaMul(SkGetPackedA32(src), src_scale))
__m128i dst_scale = SkGetPackedA32_SSE2(src);
dst_scale = _mm_mullo_epi16(dst_scale, src_scale);
dst_scale = _mm_srli_epi16(dst_scale, 8);
dst_scale = _mm_sub_epi32(_mm_set1_epi32(256), dst_scale);
__m128i result = SkAlphaMulQ_SSE2(src, alpha);
return _mm_add_epi8(result, SkAlphaMulQ_SSE2(dst, dst_scale));
}
#undef ASSERT_EQ
#endif // SkColor_opts_SSE2_DEFINED
|
