Question

我编写了一个使用SSE指令的简单C ++ openMP程序，当线程数大于4时，我遇到了分段错误。我在Linux上使用g ++。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <sys/time.h>
#include <emmintrin.h>
#include <assert.h>
#include <stdint.h>
#include <omp.h>

unsigned **a;

void f(int input_index){
    int j;
    __m128i* t = (__m128i*) a[input_index];
    for(j=0; j<4; j++)
            t[j] = _mm_set1_epi32(input_index* lenS);
}

int main(int argc, char *argv[])
{
int i,j,nThreads,tid;
bitD = new unsigned*[4];
for(i=0; i<4; i++)
    bitD[i] = new unsigned[16];

omp_set_num_threads(8); 

#pragma omp parallel private(i,nThreads,tid)
{
    tid = omp_get_thread_num();
    nThreads = omp_get_num_threads();                
    for(i=0; i<(4/nThreads); i++){
            f(i*nThreads+tid);                        
     }              
}


for(i=0; i<4; i++)
    for(j=0; j<16; j++)
        printf("a[%d][%d]=%d\n",i,j,bitD[i][j]);
}

Answer 1

正如我在上面的评论中已经提到的，您的问题与使用SSE指令无关（至少不是您发布的代码）。原因是如果你使用超过4个线程，那么循环

for(i=0; i<(4/nThreads); i++)  /* (4/nThreads) == 0 */

永远不会输入

，永远不会调用函数f。

对此的结论是，在超过4个线程的情况下，bitD[i][j]的值未初始化。但这通常不应导致分段错误。为安全起见，您可以在分配中初始化内存：

bitD[i] = new unsigned[16]();

请注意最后的()。

openMP程序中的分段错误，带有带线程的SSE指令＆gt; 4

1 个答案: