Question

我在用openMP并行化程序时遇到问题。第一个for循环大约需要10毫秒，但是第二个大约需要45秒。我不确定是否在浪费时间的循环中做错了什么。

float A[M][M];
float B[M][M];
float C[M][M];

main(int argc, char** argv) {
float temp;
float real;
float error = 0;
int i,j,k;
double time_start;
double time_end;
double time_mid;
int n  = 12;

omp_set_num_threads(n);
time_start = omp_get_wtime();


#pragma omp parallel default(shared) private(i,j,k,temp,real) reduction(+:error)
#pragma omp for
for (i=0; i<M; i++) {
        for (j=0; j<M; j++) {
                A[i][j] = ((i+1)*(j+1))/(float)M;
                B[i][j] = (j+1)/(float)(i+1);
        }
}

time_mid = omp_get_wtime();
#pragma omp for
for (i=0; i<M; i++) {
        for (j=0; j<M; j++) {
                temp = 0;
                for (k=0; k<M; k++) {
                        temp += A[i][k]*B[k][j];
                }
            C[i][j] = temp;
            real =(float) (i+1)*(j+1);
            error = error + (float) fabs(temp-real)/real;

}
}


time_end = omp_get_wtime();
error = (100/(float)(M*M))*error;

printf("Percent error for C[][] is: %f\n", error);
printf("Time is: %f\n%f\n%f\n%f\n", time_end-time_start, time_start, time_mid, time_end);

return 0;
}

Answer 1

根据OpenMP规范（第35页，2.1指令格式C / C ++）
https://www.openmp.org/wp-content/uploads/openmp-4.5.pdf

OpenMP可执行指令最多可应用于一个后续指令语句，它必须是结构化的块。

C ++中块的定义为stmt.block

因此，#pragma omp parallel default(shared) private(i,j,k,temp,real) reduction(+:error)仅适用于第一个块（您的第一个for循环）

其他循环不在'#pragma omp parallel'上下文中。

使用#pragma omp parallel{}包围第二个循环。

在openMP程序中嵌套for循环花费的时间太长

1 个答案: