Fortran G格式输出星号异常原因及解决方法问询

问题复现 #

以下代码在Intel Fortran Compiler（w_oneAPI_2022.0.0.161 + Visual Studio 2019）环境下，输出结果为星号（*），而非预期的-10.00：

WRITE(*,'(G10.4)') -9.99961061241612  
WRITE(*,'(G10.4)') -9.9996

原因分析 #

这是该版本Intel Fortran编译器G格式处理逻辑的特定bug：

• 当负数的小数部分进位导致整数位数增加时（如-9.9996进位后变为-10.00），编译器未正确重新计算输出所需的字段宽度。
• 内部仍按进位前的整数位数（1位）分配小数位数，导致格式匹配逻辑误判为指定的10位总宽度不足以容纳输出内容，最终用星号表示字段宽度不足。
• 正数无此问题是因为编译器对正数的进位宽度计算逻辑正常；极端值触发G格式的科学计数法分支（Ew.d），绕过了该bug。

解决办法 #

• 加宽格式总宽度：将G格式的总宽度增加1位，改用G11.4即可容纳进位后的负数输出：

  WRITE(*,'(G11.4)') -9.99961061241612  ! 输出：-10.00
  WRITE(*,'(G11.4)') -9.9996            ! 输出：-10.00
  

• 替换为F格式：若数值范围明确在0.1~10000之间（符合G格式使用F分支的场景），直接用F10.4格式可规避该bug：

  WRITE(*,'(F10.4)') -9.99961061241612  ! 输出：-10.0000
  WRITE(*,'(F10.4)') -9.9996            ! 输出：-10.0000
  

• 升级编译器版本：Intel oneAPI 2022.1及后续版本已修复该G格式负数进位处理bug，升级至最新稳定版可彻底解决问题。

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