Question

我正在玩GCC中的矢量内含物，特别是对于AVX，我很想写这样的东西在两个数组之间进行向量乘法运算：

#include <unistd.h>

void __attribute__((target("avx"))) vmul(float* __restrict__ cc, const float* __restrict__ aa, const float* __restrict__ bb, ssize_t size) {
    const ssize_t VECSIZE=8;
    typedef float vfloat __attribute__((vector_size(sizeof(float)*VECSIZE)));

    // duff's device, process any remainder up front
    ssize_t rem = size % VECSIZE;
    switch (rem) {
        case 7: cc[6] = aa[6]*bb[6]; /* FALLTHRU */
        case 6: cc[5] = aa[5]*bb[5]; /* FALLTHRU */
        case 5: cc[4] = aa[4]*bb[4]; /* FALLTHRU */
        case 4: cc[3] = aa[3]*bb[3]; /* FALLTHRU */
        case 3: cc[2] = aa[2]*bb[2]; /* FALLTHRU */
        case 2: cc[1] = aa[1]*bb[1]; /* FALLTHRU */
        case 1: cc[0] = aa[0]*bb[0]; /* FALLTHRU */
        case 0: break;
    }
    size -= rem;

    // process rest of array
    const vfloat *va = (const vfloat*)(aa+rem);
    const vfloat *vb = (const vfloat*)(bb+rem);
    vfloat *vc = (vfloat*)(cc+rem);

    for (ssize_t ii=0; ii < size; ii++) {
        vc[ii] = va[ii]*vb[ii];
    }    
}

int main() {
}

问题是将数据导入向量类型所需的指针别名。 GCC很乐意让你这样做（没有-Wall -Wextra -ansi -pedantic的警告），但随后假设底层内存对齐是合适的。因此它在内部循环中生成vmovaps指令：

   0x0000000000400660 <+176>:   vmovaps (%rsi,%rax,1),%ymm0
   0x0000000000400665 <+181>:   vmulps (%rdx,%rax,1),%ymm0,%ymm0
   0x000000000040066a <+186>:   vmovaps %ymm0,(%rdi,%rax,1)
   0x000000000040066f <+191>:   add    $0x20,%rax
   0x0000000000400673 <+195>:   cmp    %r8,%rax
   0x0000000000400676 <+198>:   jne    0x400660 <_Z4vmulPfPKfS1_l+176>

哪个好，直到你传入一些非对齐的内存（或者在我的情况下大小不是8的倍数），然后很高兴看到程序试图用对齐的指令加载未对齐的内存。

使用向量扩展是否有正确的方法？

Answer 1

你可以reduce the alignment这样：

vmovups

通过这项更改，我收到<winsock2.h>条指示。

使用gcc中的向量内部函数对常规数组进行别名

1 个答案: