2017-02-21 148 views
4

我想使用SIMD重寫這段代碼:使用SIMD,我如何將8位掩碼擴展爲16位掩碼?

int16_t v; 
int32_t a[16]; 
int8_t b[32]; 
... 
((int16_t *)a[i])[0] = b[i]==1? -v:v; 
((int16_t *)a[i])[1] = b[i]==1? -v:v; 

我想用_mm256_cmpeq_epi8生成一個掩碼矢量,在那之後我可以用_mm256_and_si256_mm256_andnot_si256進行價值選擇的。
問題是b [i]是8位整數,而v是16位。
如果掩碼向量與{0xff, 0x00, 0xff, 0x00...}相似,則需要將其擴展爲{0xffff, 0x0000, 0xffff, 0x0000...}以執行16位值選擇。
我該怎麼做? (對不起,我的英文版)

編輯:
我找到了靈感來自this question的解決方案。
_mm256_shuffle_epi256只能在128位通道內執行。 因此,我打破_mm256i掩碼到2 _mm128i寄存器。然後用_mm256_broadcastsi128_si256_mm256_shuffle_epi256我得到了結果。

回答

2

有一個解決方案:

int16_t v; 
int32_t a[16]; 
int8_t b[32]; 

//((int16_t *)a[i])[0] = b[i]==1? -v:v; 
//((int16_t *)a[i])[1] = b[i]==1? -v:v; 

__m256i _1 = _mm256_set1_epi8(1); 
__m256i _b = _mm256_loadu_si256((__m256i*)b); 

__m256i mask8i = _mm256_cmpeq_epi8(_b, _1); // 8-bit compare mask 

__m256i permutedMask8i = _mm256_permute4x64_epi64(mask8i, 0xD8); 
__m256i mask16iLo = _mm256_unpacklo_epi8(permutedMask8i, permutedMask8i); // low part of 16-bit compare mask 
__m256i mask16iHi = _mm256_unpackhi_epi8(permutedMask8i, permutedMask8i); // high part of 16-bit compare mask 

__m256i positiveV = _mm256_set1_epi16(-v); //positive mask condition 
__m256i negativeV = _mm256_set1_epi16(v); //negative mask condition 

__m256i _aLo = _mm256_blendv_epi8(negativeV, positiveV, mask16iLo); 
__m256i _aHi = _mm256_blendv_epi8(negativeV, positiveV, mask16iHi); 

_mm256_storeu_si256((__m256i*)a + 0, _aLo); 
_mm256_storeu_si256((__m256i*)a + 1, _aHi);