Apple M2芯片编译双精度Fortran程序（AVL等）报错求助

问题概述 #

在2023款搭载M2 Pro芯片、运行macOS Sonoma系统的MacBook Pro上，使用gcc/gfortran编译Mark Drela开发的AVL程序时，单精度编译版本可正常运行，但双精度编译版本出现异常：程序能正常启动、加载测试用例hershey.avl并查看几何模型，但进入oper菜单执行x命令时，出现配平失败错误，伴随浮点运算异常提示。

复现步骤 #

执行以下命令序列可复现问题：

./avl
load ../runs/hershey.avl
oper
x

完整报错信息 #

.OPER (case 1/1)   c>  x
Building normalwash AIC matrix...
Factoring normalwash AIC matrix...
Building source+doublet strength AIC matrix...
Building source+doublet velocity AIC matrix...
Building bound-vortex velocity matrix...

iter d(alpha)   d(beta)    d(pb/2V)   d(qc/2V)   d(rb/2V)   LEflap
1   Infinity   Infinity -0.225+308 -0.000E+00  0.000E+00 -0.000E+00
Cannot trim.  Alpha too large.  a =                  Infinity

.OPER (case 1/1)   c>

AVL   c>  quit
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_OVERFLOW_FLAG IEEE_UNDERFLOW_FLAG

当前编译参数 #

使用的gfortran编译参数为：
-m64 -fdefault-real-8 -O -fallow-argument-mismatch
可提供完整Makefile用于进一步排查。

排查方向与解决方案 #

1. 统一依赖库的双精度编译配置 #

确保plotlib、eispack等所有依赖库均使用与AVL相同的双精度编译参数（-fdefault-real-8）编译。若依赖库为单精度编译，与双精度AVL链接会导致类型不匹配，引发数值计算错误。

2. 调整浮点运算行为适配Apple Silicon #

Apple Silicon的ARM架构浮点处理与x86存在差异，尝试添加以下编译参数调整浮点行为：

• -fno-signaling-nans：关闭信号NaN处理，避免因IEEE标准严格性触发异常
• -fno-trapping-math：禁用浮点陷阱，防止程序因溢出/除零异常终止
• -ffpe-summary=none：关闭浮点异常的汇总提示，仅保留核心错误信息

3. 降低优化等级排查编译优化问题 #

当前使用-O优化等级，尝试降低至-O0或-O1编译，排查是否为ARM优化器对旧Fortran代码的双精度处理存在bug。优化等级过高可能导致数值计算精度丢失或逻辑错误。

4. 适配Apple Silicon的编译参数 #

检查Makefile中是否存在针对x86的特定参数（如-march=native），替换为ARM架构适配的-mcpu=native，确保编译参数与M2 Pro芯片匹配。

5. 调试定位数值异常点 #

添加-g -Wall编译参数开启调试模式，使用lldb或gdb启动AVL，在配平计算环节设置断点，查看导致Infinity的具体变量与计算步骤。重点检查AIC矩阵的条件数、求逆运算，确认是否存在病态矩阵或除数为0的情况。

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