三维数组在内存的中的存储结构

    在FORTRAN中，数组在内存中是按列连续存储的；一个数组在内存中也是连续的。对三维数组array(x, y, z)来说，x ,y, z分别表示的数组的行、列、页标，在内存中按列连续存储，每一个页面上的元素在内存中的地址也是连续的。以A(2,2,2)为例，它在内存中的存储结构如下图

    目前的计算机硬件体系结构决定了在读取大批量数据时，如果这一批数据都位于临近的内存中，读取操作会执行的较快。因此，在特定的情况下，我们可能通过程序的执行时间，来验证数组的在内存的中的存储结构。

module variable integer, parameter::imax = 800, jmax = 800, kmax = 800 integer a(imax, jmax, kmax), b(imax,jmax, kmax) end module variable program main use variable integer i,j,k real(kind=8) start, end call cpu_time(start) do 10 k = 1, kmax do 10 i = 1, imax do 10 j = 1, jmax a(i,j,k) = i + j + k 10 continue call cpu_time(end) print *, end - start call cpu_time(start) do 20 k = 1, kmax do 20 j = 1, jmax do 20 i = 1, imax b(i,j,k) = i + j + k 20 continue call cpu_time(end) print *, end-start call cpu_time(start) do 30 j = 1, jmax do 30 i = 1, imax do 30 k = 1, kmax b(i,j,k) = i + j + k 30 continue call cpu_time(end) print *, end-start call cpu_time(start) do 40 i = 1, imax do 40 j = 1, jmax do 40 k = 1, kmax b(i,j,k) = i + j + k 40 continue call cpu_time(end) print *, end-start end  

 

结果：

[root@c0109 zlt]# gfortran hello.F -o hello [root@c0109 zlt]# ./hello 7.25245300000000 6.89643000000000 12.9048070000000 66.9721850000000 [root@c0109 zlt]#  

 

 

 

 

附：验证数组的按页存储

program pack integer, parameter::imax=2,jmax=2,kmax = 4 integer arr(imax,jmax,kmax),i, j, k do 10 k = 1, kmax do 10 j = 1, jmax do 10 i = 1, imax arr(i,j,k) = k 10 continue print *,arr stop end  

result:

[root@c0109 zlt]# gfortran abc.F -o abc [root@c0109 zlt]# ./abc 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 [root@c0109 zlt]#