调用MPI_Bcast时传入自定义精度dp报错，如何解决？

你犯的核心错误是把Fortran的**精度参数dp**直接传给了MPI函数——MPI要求的是它自己定义的数据类型常量（比如MPI_DOUBLE_PRECISION），而dp只是个整数，MPI根本不认识它，所以才会报“Invalid datatype”。

下面给三种可行的解决方法：

1. 直接用MPI预定义类型（最简单） #

你用selected_real_kind(15,307)定义的dp对应的就是双精度，直接把MPI_Bcast里的dp换成MPI_DOUBLE_PRECISION就行：

call MPI_Bcast(buffername, count, MPI_DOUBLE_PRECISION, 0, MPI_COMM_WORLD, ierr)

如果是单精度（比如selected_real_kind(6,37)），就用MPI_REAL。

2. 动态匹配精度与MPI类型（适配多种精度） #

如果你的代码可能切换不同精度，可以写个简单的判断逻辑，根据dp的值自动选对应的MPI类型：

integer :: mpi_dp_type
if (dp == selected_real_kind(15,307)) then
    mpi_dp_type = MPI_DOUBLE_PRECISION
else if (dp == selected_real_kind(6,37)) then
    mpi_dp_type = MPI_REAL
else
    error stop "当前精度不支持MPI通信"
end if

call MPI_Bcast(buffername, count, mpi_dp_type, 0, MPI_COMM_WORLD, ierr)

3. 创建自定义MPI类型（最严谨跨平台） #

如果要完全保证Fortran自定义实型和MPI类型的匹配（比如某些编译器的特殊精度），可以用MPI_Type_create_f90_real生成专属的MPI数据类型：

integer :: mpi_dp_type, ierr
! 根据Fortran的kind值生成MPI类型
call MPI_Type_create_f90_real(dp, 0, mpi_dp_type, ierr)
! 提交类型让MPI可用
call MPI_Type_commit(mpi_dp_type, ierr)

! 用自定义类型做广播
call MPI_Bcast(buffername, count, mpi_dp_type, 0, MPI_COMM_WORLD, ierr)

! 用完后释放类型（可选，但建议养成习惯）
call MPI_Type_free(mpi_dp_type, ierr)

这里第二个参数0表示用默认对齐方式，这个方法兼容性拉满，适合跨编译器或平台的场景。

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