fortran - Fortran 中的无穷大

Question

在 Fortran 中将变量设置为 +Infinity 的最安全方法是什么？目前我正在使用:

program test
  implicit none
  print *,infinity()
contains
  real function infinity()
    implicit none
    real :: x
    x = huge(1.)
    infinity = x + x
  end function infinity
end program test

但我想知道是否有更好的方法？

Answer 1

如果你的编译器 supports ISO TR 15580 IEEE 算术是所谓的 Fortran 2003 标准的一部分，您可以使用 ieee_* 模块中的程序。

PROGRAM main

  USE ieee_arithmetic

  IMPLICIT NONE

  REAL :: r

  IF (ieee_support_inf(r)) THEN
    r = ieee_value(r,  ieee_negative_inf)
  END IF

  PRINT *, r

END PROGRAM main