使用基本 SSE 指令的分段错误

Question

我刚开始使用 SSE 指令： 程序将两个矩阵读入一个连续的缓冲区，然后尝试从那里对它们进行操作，代码如下：

void multMatrix(void * buffer,int n){
    int i,j;
    float *p = (float*) buffer;
    float *q = (float*) buffer + (n*n);
    float mr[n][n];

    __m128 va = _mm_load_ps(&p[0]);
    __m128 vb = _mm_load_ps(&q[0]);
    __m128 mm_r = _mm_add_ps( va, vb );
    _mm_store_ps (mr[0],mm_r);

    return;
}

缓冲区定义为：

void * buffer = malloc(sizeof(float)*n*n*2 + sizeof(float));

我得到段。这条线的故障：__m128 mm_r = _mm_add_ps( va, vb );

但是如果我评论_mm_store_ps (mr[0],mm_r);它会正常退出。

我只是在测试和学习 SSE 指令，但我对这个程序的目标是：

A[n][n] * B[n][n] = C[n][n] 

Answer 1

首先尝试使用未对齐的加载和存储：loadu、storeu。然后使缓冲区/指针对齐到 16 字节边界。

此外，指针 q 访问的矩阵超出范围——当最后一条存储指令被注释掉时，编译器可能足够聪明，不会产生任何这些指令。