Question

我正在尝试编写一个简单的Fortran库来进行区间算术练习。我想显式设置舍入模式，做一些工作，然后将舍入模式恢复到原来的状态。但是，在使用gfortran，gcc的Fortran前端进行编译时，我无法确定需要将哪个库链接到生成的可执行文件。

! get_rounding_mode.f03
! print the rounding mode

program get_rounding_mode
  f = IEEE_GET_ROUNDING_MODE()
  print *,f
end program get_rounding_mode

尝试可能可行的最简单方法会给我

gfortran get_rounding_mode.f03 
/usr/bin/ld: /tmp/ccTLaxeN.o: in function `MAIN__':
get_rounding_mode.f03:(.text+0x20): undefined reference to `ieee_get_rounding_mode_'
collect2: error: ld returned 1 exit status
Exit 1

我到处寻找ieee_get_rounding时发现了它，但我不知道如何引导gfortran链接它，因为它似乎已经在libgfortran中。

find /usr/ -exec nm --print-file-name '{}' '+' 2>&1 | grep 'ieee_get_rounding'
/usr/lib/libgfortran.so.5:000000000023edc0 T __ieee_arithmetic_MOD_ieee_get_rounding_mode
/usr/lib/libgfortran.so:000000000023edc0 T __ieee_arithmetic_MOD_ieee_get_rounding_mode

Answer 1

IEEE_GET_ROUNDING_MODE不是函数。这是一个子例程。您需要做类似的事情

program get_rounding_mode
   use ieee_arithmetic
   implicit none
   ieee_rounding_type mode
   real x
   if (ieee_support_rounding(x)) then
      call ieee_get_rounding_mode(mode)       ! Get current rounding mode
      call ieee_set_rounding_mode(IEEE_TO_UP) ! Set rounding up
      !
      ! Do your work here!
      !
      call ieee_set_rounding_mode(mode)       ! Reset rounding mode
   end if
 end program get_rounding_mode

哇，忘记了implicit none和声明x

包含IEEE_GET_ROUNDING_MODE

1 个答案: