28

我正在尝试学习使用内在函数进行编码，下面是添加的代码

compiler used: icc

#include<stdio.h>
#include<emmintrin.h>
int main()
{
        __m128i a = _mm_set_epi32(1,2,3,4);
        __m128i b = _mm_set_epi32(1,2,3,4);
        __m128i c;
        c = _mm_add_epi32(a,b);
        printf("%d\n",c[2]);
        return 0;
}

我收到以下错误：

test.c(9): error: expression must have pointer-to-object type
        printf("%d\n",c[2]);

如何打印c类型变量中的值__m128i

4

4 回答 4

26

使用此功能打印它们：

#include <stdint.h>
#include <string.h>

void print128_num(__m128i var)
{
    uint16_t val[8];
    memcpy(val, &var, sizeof(val));
    printf("Numerical: %i %i %i %i %i %i %i %i \n", 
           val[0], val[1], val[2], val[3], val[4], val[5], 
           val[6], val[7]);
}

您在打印之前将 128 位拆分为 16 位（或 32 位）。

如果您有 64 位支持，这是 64 位拆分和打印的一种方式：

#include <inttypes.h>

void print128_num(__m128i var) 
{
    int64_t v64val[2];
    memcpy(v64val, &var, sizeof(v64val));
    printf("%.16llx %.16llx\n", v64val[1], v64val[0]);
}

注意：将&var直接转换为int*oruint16_t*也可以使用 MSVC，但这违反了严格的别名并且是未定义的行为。使用memcpy是执行相同操作的标准兼容方式，并且通过最少的优化，编译器将生成完全相同的二进制代码。

于 2012-11-06T18:42:48.537 回答
18
  • 可跨 gcc/clang/ICC/MSVC、C 和 C++ 移植。
  • 在所有优化级别下完全安全：无严格混叠违规 UB
  • 以十六进制打印为 u8、u16、u32 或 u64 元素（基于@AG1 的回答）
  • 按内存顺序打印（最不重要的元素在前，如_mm_setr_epiX）。如果您希望以英特尔手册使用的相同顺序打印，则反转数组索引，其中最重要的元素位于左侧（如_mm_set_epiX）。相关：显示向量寄存器的约定

使用 a__m128i*从数组加载int是安全的，因为__m128类型被定义为允许别名，就像 ISO C 一样unsigned char*。（例如在 gcc 的标题中，定义包括__attribute__((may_alias)).）

反过来是不安全的（指向对象的int*一部分__m128i）。MSVC 保证这是安全的，但 GCC/clang 不保证。（-fstrict-aliasing默认开启）。它有时适用于 GCC/clang，但为什么要冒险呢？它有时甚至会干扰优化；看到这个问答。另请参阅硬件 SIMD 向量指针和相应类型之间的“重新解释转换”是否是未定义的行为？

(uint32_t*) &my_vector违反了 C 和 C++ 别名规则，并且不能保证按您期望的方式工作。存储到本地数组然后访问它可以保证是安全的。它甚至可以优化大多数编译器，例如，您可以直接从 xmm 获取movq/pextrq到整数寄存器，而不是实际的存储/重新加载。

Godbolt 编译器资源管理器上的 Source + asm 输出：证明它可以使用 MSVC 等进行编译。

#include <immintrin.h>
#include <stdint.h>
#include <stdio.h>

#ifndef __cplusplus
#include <stdalign.h>   // C11 defines _Alignas().  This header defines alignas()
#endif

void p128_hex_u8(__m128i in) {
    alignas(16) uint8_t v[16];
    _mm_store_si128((__m128i*)v, in);
    printf("v16_u8: %x %x %x %x | %x %x %x %x | %x %x %x %x | %x %x %x %x\n",
           v[0], v[1],  v[2],  v[3],  v[4],  v[5],  v[6],  v[7],
           v[8], v[9], v[10], v[11], v[12], v[13], v[14], v[15]);
}

void p128_hex_u16(__m128i in) {
    alignas(16) uint16_t v[8];
    _mm_store_si128((__m128i*)v, in);
    printf("v8_u16: %x %x %x %x,  %x %x %x %x\n", v[0], v[1], v[2], v[3], v[4], v[5], v[6], v[7]);
}

void p128_hex_u32(__m128i in) {
    alignas(16) uint32_t v[4];
    _mm_store_si128((__m128i*)v, in);
    printf("v4_u32: %x %x %x %x\n", v[0], v[1], v[2], v[3]);
}

void p128_hex_u64(__m128i in) {
    alignas(16) unsigned long long v[2];  // uint64_t might give format-string warnings with %llx; it's just long in some ABIs
    _mm_store_si128((__m128i*)v, in);
    printf("v2_u64: %llx %llx\n", v[0], v[1]);
}

如果您需要移植到 C99 或 C++03 或更早版本（即没有 C11/C++11），请删除并alignas()使用storeu. store或使用__attribute__((aligned(16)))or__declspec( align(16) )代替。

（如果您使用内在函数编写代码，则应该使用最新的编译器版本。较旧的编译器通常比较旧的编译器生成更好的 asm，包括 SSE/AVX 内在函数。但也许您想将 gcc-6.3 与-std=gnu++03C++03一起使用未准备好 C++11 或其他东西的代码库的模式。）

调用所有 4 个函数的示例输出

// source used:
__m128i vec = _mm_setr_epi8(1, 2, 3, 4, 5, 6, 7,
                            8, 9, 10, 11, 12, 13, 14, 15, 16);

// output:

v2_u64: 0x807060504030201 0x100f0e0d0c0b0a09
v4_u32: 0x4030201 0x8070605 0xc0b0a09 0x100f0e0d
v8_u16: 0x201 0x403 0x605 0x807  | 0xa09 0xc0b 0xe0d 0x100f
v16_u8: 0x1 0x2 0x3 0x4 | 0x5 0x6 0x7 0x8 | 0x9 0xa 0xb 0xc | 0xd 0xe 0xf 0x10

如果要使用前导零填充以获得一致的输出宽度，请调整格式字符串。见printf(3)。

于 2017-10-15T07:01:02.657 回答
4

我知道这个问题被标记为 C，但在寻找相同问题的 C++ 解决方案时，它也是最好的搜索结果。

所以，这可能是一个 C++ 实现：

#include <string>
#include <cstring>
#include <sstream>

#if defined(__SSE2__)
template <typename T>
std::string __m128i_toString(const __m128i var) {
    std::stringstream sstr;
    T values[16/sizeof(T)];
    std::memcpy(values,&var,sizeof(values)); //See discussion below
    if (sizeof(T) == 1) {
        for (unsigned int i = 0; i < sizeof(__m128i); i++) { //C++11: Range for also possible
            sstr << (int) values[i] << " ";
        }
    } else {
        for (unsigned int i = 0; i < sizeof(__m128i) / sizeof(T); i++) { //C++11: Range for also possible
            sstr << values[i] << " ";
        }
    }
    return sstr.str();
}
#endif

用法：

#include <iostream>
[..]
__m128i x
[..]
std::cout << __m128i_toString<uint8_t>(x) << std::endl;
std::cout << __m128i_toString<uint16_t>(x) << std::endl;
std::cout << __m128i_toString<uint32_t>(x) << std::endl;
std::cout << __m128i_toString<uint64_t>(x) << std::endl;

结果：

141 114 0 0 0 0 0 0 151 104 0 0 0 0 0 0
29325 0 0 0 26775 0 0 0
29325 0 26775 0
29325 26775

注意：有一种简单的方法可以避免if (size(T)==1)，请参阅https://stackoverflow.com/a/28414758/2436175

于 2014-09-24T08:32:37.077 回答
2 
#include<stdio.h>
#include<emmintrin.h>
int main()
{
    __m128i a = _mm_set_epi32(1,2,3,4);
    __m128i b = _mm_set_epi32(1,2,3,4);
    __m128i c;

    const int32_t* q; 
    //add a pointer 
    c = _mm_add_epi32(a,b);

    q = (const int32_t*) &c;
    printf("%d\n",q[2]);
    //printf("%d\n",c[2]);
    return 0;
}

试试这个代码。

于 2013-12-25T01:57:20.430 回答