OpenMP的C++ GCC基本套路

Question

我运行一个非常简单的程序在C++中使用OpenMP和测量所经过的时间...代码在去读取，

#include <iostream> 
#include <math.h> 
#include "timer.h" 
#include <omp.h> 



int main() 
{ 
    double start,finish; 
    int i; 
    int n=8000; 
    double a[n]; 
    double b[n]; 
    double c[n]; 



    GET_TIME(start); 
#pragma omp parallel private(i,a) shared(b,c,n) 
    { 
#pragma omp for 
     for (i=0; i<n-1; i++) 
     b[i] += (a[i] + a[i+1])/2; 
#pragma omp for 
     for (i=0; i<n-1; i++) 
      c[i] += (a[i] + a[i+1])/2; 
    } 
    GET_TIME(finish); 
    std::cout<< "Elapsed time is" <<(finish-start)<<"seconds"; 
    return 0; 
} 

代码我用下面的编译bash脚本（观察到线程在环境中定义变量OMP_NUM_THREADS = $ N）：

#!/bin/bash 

clear 

g++ -O3 -o test test.cpp -fopenmp 

for n in $(seq 1 8); do 
    export OMP_NUM_THREADS=$n 
    ./test 
    echo threads=$n 
done 

其结果是，随着增加线程的数量降低性能的一般趋势被观察为如下：（当然数字可以改变）...

Elapsed time is0.000161886secondsthreads=1 
Elapsed time is0.00019002secondsthreads=2 
Elapsed time is0.00226498secondsthreads=3 
Elapsed time is0.000210047secondsthreads=4 
Elapsed time is0.000212908secondsthreads=5 
Elapsed time is0.00920105secondsthreads=6 
Elapsed time is0.00937104secondsthreads=7 
Elapsed time is0.000834942secondsthreads=8 

有关提高性能（而不是减少它）的任何建议？ 非常感谢！

Answer 1

你可以这样做，它会增加每个线程完成的操作。这是为了克服启动一个新线程所需的开销，通过让线程做更多的工作。而且，不需要声明b，c或n是共享的。

#pragma omp parallel private(i,a,b,c,n) 
{ 
#pragma omp for schedule(static) 
    for (i=0; i<n-1; i++){ 
     b[i] += (a[i] + a[i+1])/2; 
     c[i] += (a[i] + a[i+1])/2;} 
}