如何在Fortran中为C_DOUBLE类型实型字面量使用科学计数法

在Fortran里用科学计数法写浮点数确实挺顺手的，比如1.5d-10这种写法，但如果要严格对应C语言里的double类型（也就是Fortran的C_DOUBLE），可不能只靠d后缀凑活——虽然多数编译器里double precision和C_DOUBLE是同一个kind，但标准没强制要求，所以得用更靠谱的方法：

方法1：直接使用_c_double后缀（推荐） #

首先必须引入iso_c_binding模块，这个模块提供了和C类型对应的Fortran类型参数。然后给你的科学计数法字面量加上_c_double后缀，编译器就会直接把它识别为C_DOUBLE类型的字面量：

use, intrinsic :: iso_c_binding
implicit none

real(c_double) :: my_val = 1.5e-10_c_double  ! 直接指定为C_DOUBLE类型

这里用e或者d都可以（比如1.5d-10_c_double），后缀才是决定类型的关键。

方法2：显式类型转换 #

如果已经有一个普通的双精度字面量（比如1.5d-10），可以用real()函数显式转换为C_DOUBLE类型：

use, intrinsic :: iso_c_binding
implicit none

real(c_double) :: my_val
my_val = real(1.5d-10, c_double)  ! 把普通双精度转成C兼容的双精度

这种方式适合你需要处理已有字面量，或者动态转换的场景。

关键提醒 #

千万别默认d后缀的双精度就是C_DOUBLE——虽然大多数情况下两者一致，但Fortran的double precision是语言自身的类型，而C_DOUBLE是严格对齐C语言double的类型，为了代码的可移植性（尤其是跨编译器、跨平台的时候），一定要通过iso_c_binding模块来指定类型。

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