c - 使用基本 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 超出范围访问矩阵 - 当最后一条存储指令被注释掉时，编译器可能足够聪明，不会产生任何这些指令。