Question

这是一个简短的问题，但我试图提供尽可能详细的信息。

我正在编译一个旧的，但仍然积极开发的fortran代码（f77标准）在Scientific Linux上。这段代码的推荐编译器是ifort和gfortran。使用gfortran我可以编译并运行此代码。但是，如果我使用DEBUG = 1标志进行编码，则代码会编译，但会以SEG FAULT终止。单步执行gdb会导致以下错误来源：

REAL*4 TIME_CURRENT     
CALL CPU_TIME(TIME_CURRENT) 
ISECS = 100*INT(TIME_CURRENT)

程序终止于：

Program received signal SIGFPE, Arithmetic exception. 
timer (init=1, isecs=0) at myprog.f:1818
1818 ISECS = 100*INT(TIME_CURRENT)

如果我在第1818行停止执行并检查ISECS和TIME_CURRENT，我会得到：

(gdb) ptype(TIME_CURRENT)
type = real(kind=4)
(gdb) ptype(ISECS)
type = integer(kind=4)

我尝试过更具体和使用：

ISECS = 100*INT(TIME_CURRENT,4)

但它没有帮助。我没有看到这可能等同于算术错误？

我的（makefile生成的）调试编译标志是：

gfortran -fno-automatic -m32 -O0 -g \
         -ffpe-trap=invalid,zero,overflow,underflow,precision -static-libgfortran

当我编译出调试时，我不再收到SEG FAULT但是我的编译标志是

gfortran -fno-automatic -m32 -O2 -ffast-math -static-libgfortran

我不是一名强大的程序员，所以非常感谢任何帮助。请注意，我正在64位系统上进行编译，但强制进行32位编译，因为这是必需的。

Answer 1

您的代码看起来不像FORTRAN 77，CPU_TIME来自Fortran 95.无论如何，您的调试选项非常严格。 -ffpe-trap=invalid,zero,overflow,underflow,precision表示浮点运算的许多合法用法都会导致异常。我建议只使用-ffpe-trap=invalid,zero,overflow。

具体来自gfortran手册：

Some of the routines in the Fortran runtime library, like CPU_TIME ,
are likely to trigger floating point exceptions when "ffpe-trap=precision"
is used. For this reason, the use of "ffpe-trap=precision" is not recommended.

Fortran算术异常

1 个答案: