在Fortran中如何正确将real(kind=8)数组转换为real(kind=4)数组

在gfortran中，real(kind=8)对应IEEE 754双精度浮点数，可保留约15~16位十进制有效数字；real(kind=4)对应IEEE 754单精度浮点数，可保留约6~7位十进制有效数字。双精度转单精度必然丢失超出单精度承载能力的精度，但可以通过正确实现做到仅损失必要信息、无额外精度损失。

核心实现方法 #

直接使用Fortran内置的real()类型转换函数配合整数组操作即可，这是效率最高、精度损失最小的方案，不需要手动写循环遍历数组元素。
代码示例：

program demo_convert
    implicit none
    ! 定义可移植的精度参数，gfortran下与kind=4/8完全等价，避免硬编码kind的兼容性问题
    integer, parameter :: sp = selected_real_kind(p=6, r=37)
    integer, parameter :: dp = selected_real_kind(p=15, r=307)
    
    ! 示例双精度数组
    real(dp), dimension(5) :: arr_dp = [1.234567890123456_dp, 2.345678901234567_dp, &
                                         3.456789012345678_dp, 4.567890123456789_dp, &
                                         5.678901234567890_dp]
    ! 目标单精度数组
    real(sp), dimension(size(arr_dp)) :: arr_sp
    
    ! 核心转换语句：整数组转换，无额外精度损失
    arr_sp = real(arr_dp, kind=sp)
    
    ! 验证输出
    print *, "原双精度数组："
    print '(F22.16)', arr_dp
    print *, "转换后单精度数组："
    print '(F22.16)', arr_sp
end program demo_convert

关键注意事项 #

• 必须显式指定real()函数的第二个kind参数：如果省略该参数，转换结果会跟随gfortran的默认实型设置（编译时添加-fdefault-real-8参数会把默认实型改成kind=8，导致转换不生效），显式指定可以避免编译选项带来的非预期行为
• 不要使用手动乘除10的幂次、手动截断小数位这类自定义转换逻辑：这类操作反而会引入额外的精度损失，gfortran内置的real()转换严格遵循IEEE 754标准的就近舍入规则（ties to even），是当前场景下信息损失最少的转换方式
• 若有特殊舍入需求，可先对双精度数组做舍入处理再转换：比如需要向零舍入、保留固定小数位时，可以配合int()、nint()、floor()等函数实现，示例：

  ! 示例：保留6位小数后转换
  arr_sp = real(nint(arr_dp * 1.0e6_dp) / 1.0e6_dp, kind=sp)
  

• 转换前可做范围校验：单精度浮点数的最大可表示正值约为3.4028235e+38，最小正正规数约为1.17549435e-38，如果原双精度数组存在超出该范围的数值，转换会产生inf/NaN，可提前做截断或缩放处理避免异常。

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