linux fortran 大数组,数组 – Fortran 90中双精度的最大数组大小是多少？

好的,我将解释为什么字节数有限,而不仅仅是元素数.在数组索引期间,必须计算元素的地址.当然它必须适合intptr_t C变量.此外,以字节为单位的数组大小必须符合size_t C变量.这些在现代机器上的32位和64位程序上都是32位或64位.同样也适用于程序可寻址的虚拟内存！还有OS和CPU可寻址的内存,即使程序是32位也可以是64位.

这是32位程序和操作系统无法处理超过4 GB内存的根本原因.即使您可以使用比所选CPU字大小更宽的Fortran变量以某种方式计算地址,CPU也无法访问它.

最后,我在32位模式的英特尔Fortran中进行了一项实验,其中包含32字节元素的数组：

complex(16), allocatable :: a(:)

do i=1,100

allocate(a(2**i))

a(size(a)) = 1

deallocate(a)

write(*,*) i

end do

end

ifort arraysize.f90 -m32 -check -traceback -g

输出如预期：

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

forrtl: severe (179): Cannot allocate array - overflow on array size calculation.

正如预期的那样,数组的大小以字节为单位溢出,程序在索引变量溢出之前很久就崩溃了.这不是编译器特定的功能,但有一个根本原因.