请教下面这段代码的差异

dypang

大家好，我的程序里面有这么一段代码：

[Fortran] 纯文本查看 复制代码

     do j=1,nt
      ! csum=G_0t(j)
      do ic0=1,ndims
       ! CP0(j,ic0)=CP0(j,ic0)*csum
       CP0(j,ic0)=CP0(j,ic0)*G_0t(j)
      enddo
     enddo

其中 nt 是比较大的数，比如，640000。以上的代码可以 run。但是如果把上面的注释去掉，像下面这样：

[Fortran] 纯文本查看 复制代码

     do j=1,nt
      csum=G_0t(j)
      do ic0=1,ndims
       CP0(j,ic0)=CP0(j,ic0)*csum
       ! CP0(j,ic0)=CP0(j,ic0)*G_0t(j)
      enddo
     enddo

想来这样会减少对 G_0t 这个数组的内存的访问。但是这样会出现运行时错误，说 segmentation fault。

看起来这两种方式对内存的处理是不一样的。有没有高手帮我解释一下？

谢谢！

dypang

pasuka 发表于 2018-1-8 11:51
csum没有定义吧？改成associate也出错？
do j=1,nt
       associate(csum=>G_0t(j))then

谢谢。不过， csum 是声明成复数类型了的，因此应该不是这个原因

dypang

pasuka 发表于 2018-1-8 13:17
无法提供能够重现错误的最小代码，实在爱莫能助了
非涉密项目，可以考虑加入qq群，然后加群主好友，再远程 ...

谢谢。我问这个问题只是想了解一下fortran的更细致的规则。这两段代码，一般情况下应该是一样的。可是，在我这里，第二种方式确实就会出现段错误。我怀疑是这两种方式对内存的操作是不一样的。但是这种错误大概只是在数组非常大的情况下才会出现。

dypang

dypang 发表于 2018-1-9 13:48
谢谢。我问这个问题只是想了解一下fortran的更细致的规则。这两段代码，一般情况下应该是一样的。可是， ...

实际上，同样的程序，在 Linux 下面用第二种方式也是可以的，但是在 Windows 10 下面就不行。所以，应该跟操作系统也有关系。

dypang

是的，大概跟编译器有关系，都是 pgfortran，但是一个是 windows 下面的版本，一个是 Linux 下面的版本。