我正在尝试用 openmp 并行化这个递归函数:
#include <stdio.h>
#include <omp.h>
void rec(int from, int to){
int len=to-from;
printf("%X %x %X %d\n", from, to, len, omp_get_thread_num());
if (len > 1){
int mid = (from+to)/2;
#pragma omp task
rec(from, mid);
#pragma omp task
rec(mid, to);
}
}
int main(int argc, char *argv[]){
long len=1024;
#pragma omp parallel
#pragma omp single
rec(0, len);
return 0;
}
但是当我运行它时,我得到了段错误:
$g++ -fopenmp -Wall -pedantic -lefence -g -O0 test.cpp && ./a.out
0 400 400 0
0 200 200 1
200 400 200 0
Segmentation fault
当我在 valgrind 中运行它时,它没有显示任何错误。没有 -lefence 它也可以工作。
我尝试了所有可能的#pragma omp 子句组合,它要么是单线程的,要么是段错误的。
怎么了?
非常感谢。