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

Segmentation fault using basic SSE instructions

我刚开始使用 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] 

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

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