MPI_COMM_SPAWN父进程逐行发矩阵给子进程出现零值输出问题

让我们一步步拆解你的问题，主要有几个关键错误导致了子进程输出全0且不符合“每行对应一个子进程”的预期：

1. 矩阵生成函数的内存分配错误 #

你的generate_matrix函数里，分配行指针数组时犯了一个低级错误：

x = (float **) malloc(m * sizeof(float));

这里你分配的是m个float的空间，但实际上x是float**类型，每个元素是指向float数组的指针，应该分配m个float*的空间。这个错误会导致内存越界，矩阵的行指针指向无效内存，最终生成的矩阵数据都是垃圾值（表现为全0）。

修复后的函数：

float** generate_matrix(int n, int m) {
    int i, j;
    float **x;
    // 修正：分配m个float指针的空间
    x = (float **) malloc(m * sizeof(float*));
    for (i = 0; i < m; i++) {
        x[i] = (float *) malloc(n * sizeof(float));
    }
    for (i = 0; i < m; i++) {
        for (j = 0; j < n; j++) {
            x[i][j] = random(100);
        }
    }
    return x;
}

2. 父进程错误使用MPI_Bcast #

MPI_Bcast的作用是把根进程的数据广播给通信器内的所有进程，你的循环执行m次广播，意味着每个子进程会收到m行数据（最后一次覆盖之前的），这完全不符合“每个子进程对应一行”的需求。

正确的做法是给每个子进程单独发送对应的行，用MPI_Send即可：

// 替换原来的MPI_Bcast循环
for (int i = 0; i < m; i++) {
    // 给远程组rank为i的子进程发送第i行
    MPI_Send(matrix[i], n, MPI_FLOAT, i, tag, workercomm);
}

3. 子进程的两处核心错误 #

错误1：获取子进程rank的方式不对 #

你用MPI_Comm_rank(MPI_COMM_WORLD, &myid);，但每个通过MPI_Comm_spawn启动的子进程，它的MPI_COMM_WORLD只包含自身，所以myid永远是0（你的输出里出现其他数字是环境异常导致的）。正确的方式是从parent通信器中获取自己在子进程组的rank：

// 替换原来的MPI_Comm_rank调用
MPI_Comm_rank(parent, &myid);

错误2：用MPI_Bcast接收数据 #

子进程不需要接收广播，应该单独接收父进程发送给自己的那一行，改用MPI_Recv：

if (parent != MPI_COMM_NULL) {
    // 接收父进程（parent通信器中rank为0的进程）发送的数据
    MPI_Recv(vector, n, MPI_FLOAT, 0, tag, parent, &status);
}

修复后的完整子进程代码：

int tag = 0;
MPI_Status status;

int main(int argc, char** argv) {

    MPI_Init(&argc, &argv);
    MPI_Comm parent;
    MPI_Comm_get_parent(&parent);
    int myid;
    MPI_Comm_rank(parent, &myid);

    int n = 4;
    float* vector = (float *) malloc(n * sizeof(float));

    if (parent != MPI_COMM_NULL) {
        MPI_Recv(vector, n, MPI_FLOAT, 0, tag, parent, &status);
    }
    printf("%d ->", myid);
    for (int i = 0; i < n; i++) {
        printf("%.2f ", vector[i]);
    }
    printf("\n");

    MPI_Comm_free(&parent);
    free(vector);
    MPI_Finalize();
    return 0;
}

额外优化建议 #

你用rand()生成随机数，但没有设置种子，每次运行生成的矩阵都是一样的。可以在父进程main函数开头添加：

srand(time(NULL));

这样每次运行都会生成不同的随机矩阵。

修复完成后，每个子进程会正确打印矩阵对应的行，输出类似：

0 ->45.00 78.00 23.00 91.00 
1 ->12.00 67.00 34.00 89.00 
2 ->56.00 21.00 78.00 43.00 
3 ->90.00 14.00 52.00 68.00 
4 ->37.00 82.00 19.00 55.00 

内容的提问来源于stack exchange，提问作者Buculei Razvan