gfortran和ifort给出了不同的结果

时间:2018-06-12 09:13:33

标签: fortran gfortran intel-fortran fortran95

我正在制作一个小代码,它将读取2.txt中的数字,这是一堆随机实数,将其保存在二维数组中并将其传递给子程序。在子例程中,我将此数组的特定维度分配给另一个矩阵,以便将其保存以供将来计算。

   implicit none
   double precision, allocatable :: x1(:,:)
   integer :: nstep, nheatoms, k, i, j
   open(unit=999,file="2.txt",status="old")

   nstep = 5
   nheatoms = 2

   allocate(x1(3,nheatoms))
   do k = 1, nstep
   do i = 1, 3
   do j = 1, nheatoms
   read(999,*) x1(i,j)
   end do
   read(999,*)
   end do
   call diffusion(nstep,nheatoms,x1)
   end do

   stop
   end

  subroutine diffusion(nstep,nheatoms,qhecent)
  implicit none
  integer, intent (in) :: nheatoms, nstep
  double precision :: qhecent(3,nheatoms)
  integer :: j
  double precision, allocatable  :: qhestep(:,:)
  integer :: nstepdiff, ncross
  integer, save :: l = 0
  double precision :: diff

  l = l + 1
   allocate(qhestep(nstep,nheatoms))
  do j = 1, nheatoms
  qhestep(l,j) = qhecent(3,j)
  end do

  if (l .gt. 1) then
  do j = 1, nheatoms
  write(*,*)qhestep(l,j), qhestep(l-1,j)
  end do
  end if
  end subroutine

2.txt文件如下:

1.0
2.1

3.2
-1.1

-2.2
-3.3

5.0
3.5

4.4
1.9

2.1
1.5

6.0
3.5

4.4
1.9

2.8
2.5

6.0
3.5

4.4
1.9

2.1
3.2

6.0
3.5

4.4
1.9

-4.3
7.9

现在,如果我使用gfortran编译,大多数时候我获得输出为:

   2.1000000000000001       -2.2000000000000002     
   1.5000000000000000       -3.2999999999999998     
   2.7999999999999998        2.1000000000000001     
   2.5000000000000000        1.5000000000000000     
   2.1000000000000001        2.7999999999999998     
   3.2000000000000002        2.5000000000000000     
  -4.2999999999999998        2.1000000000000001     
   7.9000000000000004        3.2000000000000002  

这是预料之中的。但是如果我多次运行代码,输出中会出现一个未知数字,例如:

   2.1000000000000001       -2.2000000000000002     
   1.5000000000000000       -3.2999999999999998     
   2.7999999999999998        1.2882297539194267E-231

   2.5000000000000000        0.0000000000000000     
   2.1000000000000001        2.7999999999999998     
   3.2000000000000002        2.5000000000000000     
  -4.2999999999999998        2.1000000000000001     
   7.9000000000000004        3.2000000000000002 

如果我用ifort编译,那么我得到的输出有两个零,这是错误的。

   2.10000000000000       -2.20000000000000     
   1.50000000000000       -3.30000000000000     
   2.80000000000000       0.000000000000000E+000

   2.50000000000000       0.000000000000000E+000

   2.10000000000000        2.80000000000000     
   3.20000000000000        2.50000000000000     
  -4.30000000000000        2.10000000000000     
   7.90000000000000        3.20000000000000 

我必须强调,在大多数情况下对ifort进行编译会给出正确的结果,就像gfortran一样。我在Mac OS High Sierra上。 提前谢谢。

1 个答案:

答案 0 :(得分:0)

两个编译器都应该给出垃圾结果,因为你在扩散中为每个DO k循环分配qhestep。在退出扩散时,qhestep被解除分配,因此应该丢失累积的信息。 在每次调用时返回qhestep的相同堆地址是误导性的。

要解决这个问题,你应该在进入DO k循环之前使用allocate(qhestep(nstep,nheatoms)),然后在调用扩散中传输累积数组。

ifort给出了你的预期结果是不幸的,因为它伪装了一个重大的逻辑错误。每个ALLOCATE基本上为数组提供了一个新的内存位置。

读取文件时也不提供错误测试。我建议使用iostat =进行读取,并在读取时报告数据。

总之,您的方法给出了错误的结果,这两个编译器都提供了这些结果。您提供的ifort结果在某些情况下也是错误的,但在其他情况下是正确的,这会产生误导。