MPI分布式数组视图读取输入值异常问题排查

核心问题：写入时的MPI_File_write参数错误 #

你遇到的最后一个元素值错误（1而非8），根源完全在写入环节的代码错误：
在循环写入矩阵元素时，MPI_File_write的第三个参数（要写入的元素个数）被错误设置为sizeof(&matrix[i][j])：

// 错误代码
MPI_File_write(file, &matrix[i][j], sizeof (&matrix[i][j]), MPI_FLOAT, MPI_STATUS_IGNORE);

&matrix[i][j]是指向float的指针，在64位系统上sizeof返回8，这意味着每次调用会尝试写入8个float元素。你的矩阵只有9个元素，循环执行9次后，会写入远超实际数据量的内容，直接破坏了文件末尾的数据，导致读取时最后一位出现乱值。

修复步骤 #

1. 修正写入代码 #

将MPI_File_write的count参数改为1，因为每次只需要写入1个float元素：

// 修正后的写入循环
if (rank == 0) {        
    for (int i = 0; i < 3; i++) {
        for (int j = 0; j < 3; j++) {
            MPI_File_write(file, &matrix[i][j], 1, MPI_FLOAT, MPI_STATUS_IGNORE);
        }
    }
}

2. 验证修复效果 #

修复后重新运行程序，读取输出会恢复正常：

Proc 0 received:
0 1 2 
3 4 5 
6 7 8

额外优化建议（单进程场景） #

单进程测试时，MPI_Type_create_darray的使用属于冗余操作，你可以直接用更简单的MPI_File_write_at一次性写入整个矩阵，避免循环出错：

if (rank == 0) {
    MPI_File_write_at(file, 0, matrix, 9, MPI_FLOAT, MPI_STATUS_IGNORE);
}

读取时也可以直接用MPI_File_read_at一次性读取，代码更简洁且不易出错。

内容的提问来源于stack exchange，提问作者G. Ianni