gpt4 book ai didi

c++ - 为什么并行版本更慢?

转载 作者:搜寻专家 更新时间:2023-10-31 00:06:28 24 4
gpt4 key购买 nike

<分区>

我想在矩阵上应用特定的过滤器。 (从[0][0]顺序到结束)

A[i][j] = 0.2 * (A[i][j] + A[i+1][j] + A[i-1][j] + A[i][ j+1] + A[i][j-1])

例如,如果 [i],[j] 是 [0][0](矩阵中的第一个值),我使用零作为左侧和上侧的值。

我试图理解为什么我的代码的并行版本比顺序版本慢。

当我使用多线程进行计算时,我使用的是沿对角线有独立工作的事实。我故意将矩阵扩展两行和两列(用零填充)以简化过滤器的计算。

我还尝试了各种尺寸的矩阵(最大 7000x7000)。

我的问题: http://15418.courses.cs.cmu.edu/fall2017/lecture/progbasics/slide_032

后续版本:

for (int i = 1; i < r-1; i++) {
for (int j = 1; j < c-1; j++) {
arr[i][j] = 0.2f * (arr[i][j] + arr[i][j - 1] + arr[i - 1][j] + arr[i][j + 1] + arr[i + 1][j]);
}
}

并行版本:

int n = r - 2;  
for (int slice = 0; slice < 2 * n - 1; ++slice) { //along the diagonals
int z = (slice < n) ? 0 : slice - n + 1;
#pragma omp parallel for schedule(static) //spawns threads
for (int j = z; j <= slice - z; ++j) {
pl_arr[j + 1][slice - j + 1] = 0.2f * (pl_arr[j + 1][slice - j + 1] + pl_arr[j + 1][slice - j] + pl_arr[j][slice - j + 1] + pl_arr[j + 1][slice - j + 1 + 1] + pl_arr[j + 1 + 1][slice - j + 1]);
}
}

其余代码:

int r = 7000, c = 7000;

r = r + 2;
c = c + 2;

/* initialize random seed: */
srand(time(NULL));

float **arr = (float **)malloc(r * sizeof(float *));
for (int i = 0; i < r; i++)
arr[i] = (float *)malloc(c * sizeof(float));

float **pl_arr = (float **)malloc(r * sizeof(float *));
for (int i = 0; i < r; i++)
pl_arr[i] = (float *)malloc(c * sizeof(float));


for (int i = 0; i < r; i++) {
for (int j = 0; j < c; j++) {
if ((i == 0) || (i == (r - 1)) || (j == 0) || (j == (c - 1)) ){
arr[i][j] = 0;
pl_arr[i][j] = 0;
}
else {
arr[i][j] = rand() % 99 + 1;
pl_arr[i][j] = arr[i][j];
}
}
}

#pragma omp parallel for schedule(static) - for 构造拆分 for 循环,以便当前组中的每个线程处理循环的不同部分。

结果:并行版本总是比顺序版本慢

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com