[已解决]gfortran常数问题

bian222

一段fortran 代码：

1. 
   
2. program exam
   
3. implicit none
   
4. real*8,parameter::pi=3.1415926
   
5. write(*,*) pi
   
6. end

复制代码

编译完了，给个结果是3.1415925025939941，怎么回事？在win下用visual fortran 6.5  编译结果是正确的。

fmying

你用fortran 77试试看，也一样的。

2.        program exam
3.        parameter (pi=3.1415926)
4.        write(*,*)pi
5.        end

复制代码

2. program exam
3. implicit none
4. real,parameter :: pi=3.1415926
5. write(*,*) pi
6. end

复制代码

两段代码的结果都是 3.1415925。上面的是77，下面的是90。
我测试了一下，不知道下面的结果是不是能让lz满意：

2. program exam
3. implicit none
4. real*8,parameter :: pi=3.1415926d0
5. write(*,*) pi
6. end

复制代码

输出

1. 3.1415926000000001

复制代码

bian222

又装了intel fortran，编译结果也是3.1415925。我用的是Debian,不知道在其他linux版本下结果如何，在win下结果是3.1415926。想搞清楚究竟是什么原因造成的。

tjpm

real*8的：
program exam
implicit none
real*8,parameter :: pi=3.1415926D0
write(*,*) pi
end
intel的编译器和gfortran都能正确输出

real*16的，仅限intel编译器 for linux 支持
program exam
implicit none
real*16,parameter :: pi=3.1415926Q0
write(*,*) pi
end
也能正确输出

常数赋值的时候需要给编译器说明这个数的精度

Post by bian222;1967807
又装了intel fortran，编译结果也是3.1415925。我用的是Debian,不知道在其他linux版本下结果如何，在win下结果是3.1415926。想搞清楚究竟是什么原因造成的。

bian222

明白了，谢谢两位。