错误原因分析 #

你遇到的There is no specific subroutine for the generic ‘mpi_file_write’错误，本质是编译器无法找到与你传入参数类型完全匹配的MPI_FILE_WRITE具体实现。通常由以下几个原因导致：

1. MPI数据类型与变量类型不匹配 #

Fortran的自定义类型别名（比如r8、i4）和MPI预定义数据类型的映射可能存在差异：

• 你的代码中用MPI_REAL8对应REAL(r8)变量，但部分MPI库（如OpenMPI）中，双精度浮点数对应的标准MPI类型是MPI_DOUBLE_PRECISION，而非MPI_REAL8（MPI_REAL8是某些实现的非标准别名）。
• 同理，INTEGER(i4)对应的标准MPI类型通常是MPI_INT或MPI_INTEGER，而非MPI_INTEGER4。

2. MPI接口模块使用问题 #

传统的USE mpi模块可能没有正确导出所有泛型子程序的接口，部分MPI库推荐使用现代Fortran接口USE mpi_f08来获得更完整的类型检查支持。

3. Cygwin环境的MPI库兼容性 #

在Windows的Cygwin环境下编译Linux目标代码，可能存在MPI库（如OpenMPI）的Fortran接口实现不完整，或者编译选项未正确指定目标平台的问题。

具体修复步骤 #

步骤1：统一MPI数据类型与变量类型的映射 #

修改代码中MPI_FILE_WRITE调用的MPI数据类型参数，替换为标准MPI类型：

• 将MPI_REAL8替换为MPI_DOUBLE_PRECISION（对应REAL(r8)双精度）
• 将MPI_INTEGER4替换为MPI_INT（对应INTEGER(i4)4字节整数）
• 将MPI_INTEGER8替换为MPI_INTEGER8（如果你的MPI库支持，否则用MPI_LONG_LONG_INT）

修改后的示例代码片段：

! 原调用
CALL MPI_FILE_WRITE(h%filehandle, b%dim_mults, ndims, MPI_REAL8, &
 MPI_STATUS_IGNORE, errcode)
! 修改后
CALL MPI_FILE_WRITE(h%filehandle, b%dim_mults, ndims, MPI_DOUBLE_PRECISION, &
 MPI_STATUS_IGNORE, errcode)

! 原调用
CALL MPI_FILE_WRITE(h%filehandle, b%geometry, 1, MPI_INTEGER4, &
 MPI_STATUS_IGNORE, errcode)
! 修改后
CALL MPI_FILE_WRITE(h%filehandle, b%geometry, 1, MPI_INT, &
 MPI_STATUS_IGNORE, errcode)

! 原调用
CALL MPI_FILE_WRITE(h%filehandle, b%npoints, 1, MPI_INTEGER8, &
 MPI_STATUS_IGNORE, errcode)
! 修改后（如果MPI_INTEGER8不可用，替换为MPI_LONG_LONG_INT）
CALL MPI_FILE_WRITE(h%filehandle, b%npoints, 1, MPI_INTEGER8, &
 MPI_STATUS_IGNORE, errcode)

步骤2：检查并更新MPI模块引用 #

尝试将USE mpi替换为USE mpi_f08（如果你的代码兼容现代Fortran MPI接口）：

! 原代码
USE mpi
! 修改为
USE mpi_f08

注意：如果代码中使用了大量旧版MPI的语法（比如MPI_COMM_WORLD的旧定义），可能需要逐步调整，或者继续使用USE mpi但确保MPI库的接口文件正确安装。

步骤3：调整Cygwin编译环境的MPI配置 #

1. 确保Cygwin中已安装完整的OpenMPI Fortran支持包：

   apt-cyg install openmpi openmpi-devel openmpi-fortran
   

2. 编译时指定正确的MPI库链接选项，在Makefile中确保mpif90编译器正确调用，并且没有遗漏任何MPI相关的编译参数。
3. 如果目标是Linux平台，考虑使用交叉编译工具链，而非直接在Cygwin中编译，避免环境差异导致的兼容性问题。

步骤4：验证变量类型定义 #

检查代码中r8、i4的定义是否正确：

! 确保r8是双精度
INTEGER, PARAMETER :: r8 = SELECTED_REAL_KIND(15, 307)
! 确保i4是4字节整数
INTEGER, PARAMETER :: i4 = SELECTED_INT_KIND(9)

如果类型定义不匹配，也会导致MPI数据类型无法正确映射。

测试修复 #

修改完成后，重新执行make命令，检查错误是否消失。如果仍有问题，可以增加编译的调试信息，比如：

mpif90 -c -O3 -I/usr/include -Iobj -Jobj -o obj/sdf_output_point_ru.o -Wall -Wextra src/io/sdf_output_point_ru.f90

通过-Wall -Wextra获取更多类型不匹配的警告信息，进一步定位问题。

内容的提问来源于stack exchange，提问作者J.Zhao