Используя встроенные функции SSE, я получил вектор из четырех 32-битных чисел с плавающей запятой, ограниченный диапазоном 0–255 и округленный до ближайшего целого числа. Теперь я хотел бы записать эти четыре байта.
Существует внутренний _mm_cvtps_pi8, который преобразует 32-битное в 8-битное подписанное int, но проблема в том, что любое значение, превышающее 127, ограничивается до 127. Я не могу найти никаких инструкций, которые будут зажимать к беззнаковым 8-битным значениям.
У меня есть интуиция, что я могу захотеть сделать какую-то комбинацию _mm_cvtps_pi16 и _mm_shuffle_pi8, за которой следует инструкция перемещения, чтобы получить четыре байта, которые мне нужны, в память. Это лучший способ сделать это? Я собираюсь посмотреть, смогу ли я выяснить, как кодировать маску управления перемешиванием.
ОБНОВЛЕНИЕ: кажется, что следующее делает именно то, что я хочу. Есть ли способ лучше?
#include <tmmintrin.h>
#include <stdio.h>
unsigned char out[8];
unsigned char shuf[8] = { 0, 2, 4, 6, 128, 128, 128, 128 };
float ins[4] = {500, 0, 120, 240};
int main()
{
__m128 x = _mm_load_ps(ins); // Load the floats
__m64 y = _mm_cvtps_pi16(x); // Convert them to 16-bit ints
__m64 sh = *(__m64*)shuf; // Get the shuffle mask into a register
y = _mm_shuffle_pi8(y, sh); // Shuffle the lower byte of each into the first four bytes
*(int*)out = _mm_cvtsi64_si32(y); // Store the lower 32 bits
printf("%d\n", out[0]);
printf("%d\n", out[1]);
printf("%d\n", out[2]);
printf("%d\n", out[3]);
return 0;
}
ОБНОВЛЕНИЕ 2: Вот еще лучшее решение, основанное на ответе Гарольда:
#include <smmintrin.h>
#include <stdio.h>
unsigned char out[8];
float ins[4] = {10.4, 10.6, 120, 100000};
int main()
{
__m128 x = _mm_load_ps(ins); // Load the floats
__m128i y = _mm_cvtps_epi32(x); // Convert them to 32-bit ints
y = _mm_packus_epi32(y, y); // Pack down to 16 bits
y = _mm_packus_epi16(y, y); // Pack down to 8 bits
*(int*)out = _mm_cvtsi128_si32(y); // Store the lower 32 bits
printf("%d\n", out[0]);
printf("%d\n", out[1]);
printf("%d\n", out[2]);
printf("%d\n", out[3]);
return 0;
}
_mm_shuffle_pi8это версия с регистром мм, верно? Не забывай свой_mm_empty- person harold schedule 24.04.2015-mfpmath=sseв командной строке компилятора. - person Timothy Miller schedule 24.04.2015_mm_packus_epi32на_mm_packs_epi32? Как сказал Питер, он отлично работает и требует только SSE2. Ваш (на основе Гарольда) требует SSE4.1 - person user1593842 schedule 08.09.2017