Question

在<>Fortran上,对阵列的初始化和分配价值如下:

<>Initializing:

real(kind=8):: r(3,4)
...
r(:,:) = 0.0_8

如果我们只使用

real(kind=8):: r(3,4)
...
r = 0.0_8

and what if we do as:

real(kind=8):: r(3,4)
...
r = 0

此外,情况如下:

real(kind=8):: r(3,4), q(3,4), p(30,40)
...
q = 0
r = q
r = p(1:3,21:24)

我们更愿意这样做:

real(kind=8):: r(3,4), q(3,4), p(30,40)
...
q = 0.0_8
r(:,:) = q(:,:)
r(:,:) = p(1:3,21:24)

我们确信,hope你为我们提供了你选择的每一个理由。

Answer 1

你们是否认为“实物=8”是指8个。所有汇编者都是如此。 1. 最大可转移性是利用所选择的---------------------------------------------------------------------------------------------------------

integer, parameter :: DRK = selected_real_kind (14)

Of some of the options that you list, I prefer r = 0.0_8, or better r=0.0_DRK. This defines the entire array. There is no need to designate array sections in this case, since you are calling out the entire array: r (:, :). Steve Lionel has a discussion of why trying to make arrays obvious with ":" isn t a good idea -- there are differences between array 以及 array (:). As an argument, the first has it declared dimensions, while the second always begins at 1 -- see http://software.intel.com/en-us/blogs/2008/03/31/doctor-it-hurts-when-i-do-this/

由于r是真实的,r=0意味着转换。对零值而言,这几乎肯定是偶然的,但它可能与其他价值发生重要变化。例如,

r = 3.1415927654

以及

r = 3.141592654_8

将为<代码>r提供不同的数值,因为固定数值不同,第一个数值在转让前转换为单一精确度。

Answer 2

关于有效使用Fortran阵列的一般考虑(包括声明和指派),我建议改为:。

更确切地回答您的问题:。几个月前,我进行了一些可能感兴趣的检测。这里是结果。这是对我手提电脑进行测试的,该笔手提电脑在Lino Archlinux x86-64上,GNU Fortran (GCC) 4.6.1 20110819(预发)没有优化选择。

do i = 1 , 100
  do j = 1 , 100
    do k = 1 , 100 ! innest is fastest
      array ( i , j , k ) = 0.0d0
    end do
  end do
end do
! reference time : 1.00

纽约总部

do i = 1 , 100
  do j = 1 , 100
    do k = 1 , 100
      array ( k , j , i ) = 0.0d0
    end do
  end do
end do
! time : 0.499

纽约总部

array = 0.0d0
! time : 0.250

纽约总部

array ( : , : , : ) = 0.0d0
! time : 0.250

友情链接