macOS编译pyradex/myradex遇Fortran类型不匹配错误求助

错误原因 #

这个类型不匹配错误本质是老Fortran77代码的不规范写法与现代gfortran严格类型检查的冲突，结合Rosetta2转译环境下的编译器行为差异导致：

• 早期Fortran77依赖隐式类型规则（I-N开头变量默认整数，其他默认实数），且允许通过内存地址直接传递不同类型的参数（比如用实数数组存储整数数据，直接传给要求整数数组的子例程），这类写法在旧版编译器中被容忍。
• 新版gfortran（如15.x版本）默认启用了更严格的类型检查，在x86架构通过Rosetta2转译的环境下，不会再忽略这种明显的类型不匹配；而ARM架构下可能编译时隐含了宽松检查选项，或代码内存布局巧合绕过了检查，因此无报错。

解决方法 #

1. 调整编译选项（最快临时方案） #

给gfortran添加-fallow-argument-mismatch选项，允许参数类型不匹配，兼容老代码的写法：

gfortran -fallow-argument-mismatch [其他编译选项] opkda1.f ...

该选项是GNU Fortran的扩展，专门用于处理这类老代码的兼容性问题，无需修改源码。

2. 修改源码（长期规范方案） #

若要彻底解决类型不匹配问题，需调整调用DPREP/DPREPI的代码：

• 确认RWORK(LWM)对应的内存区域实际存储整数数据，通过EQUIVALENCE声明将整数数组与该区域绑定：

  INTEGER IWK(...)
  REAL*8 RWORK(...)
  EQUIVALENCE (IWK, RWORK(LWM))
  

• 调用子例程时传递IWK而非RWORK(LWM)，确保参数类型与子例程声明匹配。

3. 不建议降级编译器 #

降级编译器会丢失新版gfortran的性能优化、bug修复和新特性，优先采用前两种方法解决问题。

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