我正在解决一个问题,其中一个进程对列进行一些计算,我将该列保存在一个临时数组中,并使用 MPI_Bcast 将其广播到所有其他进程,但最终总是在接收中得到数组中的所有零结尾。
以下是我正在使用的代码:
//nProc = no of processors;
//helperA = 2D Array i.e. helperA[size][size]
// colK = 1D Array i.e. colK[size]
for (int k = 0; k < size; ++k) {
if (k % nProc == rank) {
// One of the process will do this calculation
int temp = 0;
for (int j = k + 1; j < size; ++j) {
helperA[j][k] = helperA[j][k]/helperA[k][k];
colK[temp++] = helperA[j][k];
}
}
MPI_Bcast(colK, size - k - 1, MPI_DOUBLE, rank, MPI_COMM_WORLD);
// After this other process should get the colK updated with the calculation
for (int i = k + 1; i < size; ++i) {
if (i % nProc == rank) {
int temp = 0;
for (int j = k + 1; j < size; ++j) {
// Here colK is always zero
printf("%d %f \n", rank, colK[temp]);
helperA[j][i] = helperA[j][i] - (colK[temp++] * helperA[k][i]);
}
}
}
}
我不确定我在这里做错了什么。请提供任何帮助/建议。