Fortran MPI_Scatter处理大数组卡住、小数组正常运行问题求助

1. 核心错误：MPI_Scatter 根参数配置错误 #

MPI_Scatter 的第7个参数是发送根进程的rank，必须固定为持有发送缓冲区的进程rank（你的场景下就是0），但你填了当前进程的rank mpirank，相当于所有进程都认为自己是发送根，各自等待其他进程接收数据，直接触发死锁。
小数据量时可能因为MPI实现的内部缓冲机制凑巧跑完，大数据量缓冲占满后就会永久卡住，和你遇到的现象完全吻合。

2. 次要错误：MPI_Scatter 不支持不均匀分块 #

你给rank 0和其他rank分配的ni、nj长度不一致（rank 0多了mod余数值），但标准MPI_Scatter要求给所有接收进程发送相同长度的数据。如果你确实需要不均匀分块，需要改用MPI_Scatterv接口，手动定义每个进程的接收长度和偏移量。
如果你的总长度刚好可以被进程数整除（比如你测试用的np=10，ni_t=50000刚好是10的倍数），可以把所有进程的ni、nj都设为int(ni_t/mpisize)，不用加mod，就可以正常用MPI_Scatter。

修复后的核心代码片段 #

把你两个MPI_Scatter调用里的根参数从mpirank改成0即可，示例：

! 错误写法
call mpi_scatter(xi_t_x, ni, MPI_DOUBLE_PRECISION, temp_xi, ni, &
    MPI_DOUBLE_PRECISION, mpirank, mpi_comm_world, IERR)

! 正确写法
call mpi_scatter(xi_t_x, ni, MPI_DOUBLE_PRECISION, temp_xi, ni, &
    MPI_DOUBLE_PRECISION, 0, mpi_comm_world, IERR)

xj_t_x的scatter调用也要做同样修改。

额外优化建议 #

• xi_t_x、xj_t_x的分配和初始化只需要在rank 0执行即可，其他进程不需要持有这两个数组，可以减少内存占用。
• 调用MPI接口后可以判断ierr返回值，方便排查接口调用错误。

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