我正在关注Tim Mattson关于 OpenMP 的讲座,以了解一些并行编程概念的实现方式。
我试图观察使用3x10^8步计算PI值的并行程序的运行时间行为。
这是代码,
#include <omp.h>
#include <stadio.h>
static long num_steps = 300000000;
double step;
#define PAD 8 // tried 50 too
#define NUM_THREADS 4
int main()
{
int i, nthreads;
double pi, sum[NUM_THREADS][PAD];
double ts, te;
ts = omp_get_wtime();
step = 1.0/(double) num_steps;
omp_set_num_threads(NUM_THREADS);
#pragma omp parallel
{
int i, id,nthrds;
double x;
id = omp_get_thread_num();
nthrds = omp_get_num_threads();
if (id == 0) nthreads = nthrds;
for (i=id, sum[id]=0.0;i< num_steps; i=i+nthrds) {
x = (i+0.5)*step;
sum[id][0] += 4.0/(1.0+x*x);
}
}
for(i=0, pi=0.0;i<nthreads;i++)
pi += sum[i][0] * step;
te = omp_get_wtime();
printf("%.10f\n", pi);
printf("%.f\n", te-ts);
}
现在我在双核机器上运行 Ubuntu 14.04 LTS。调用omp_get_num_procs()
返回2。运行时间完全随机,从 1.31 秒到 4.46 秒不等。而串行程序几乎总是花费 2.31 秒。
我尝试创建 1、2、3、4,最多 10 个线程。每种情况下的运行时间变化太大,尽管在线程更多的情况下平均值会更小。我没有运行任何其他应用程序。
谁能解释为什么运行时间变化太大?
如何准确计算运行时间?讲师给出了他的计算机的运行时间,这似乎是一致的。而且他也在使用双核处理器。