Question

我正在制作一个小代码，它将读取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上。提前谢谢。

Answer 1

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

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

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

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

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

gfortran和ifort给出了不同的结果

1 个答案: