我在这行 MATLAB 代码中遇到内存不足错误:
result = (A(1:xmax,1:ymax,1:zmax) .* B(2:xmax+1,2:ymax+1,2:zmax+1) +
A(2:xmax+1,2:ymax+1,2:zmax+1) .* B(1:xmax,1:ymax,1:zmax)) ./ C
其中 C 是另一个数组。这是在 32 位 MATLAB 上(我目前似乎无法获得 64 位版本,这将暂时解决我的问题)。
数组结果 A、B 和 C 是预先初始化的并且永远不会改变大小。然后我猜这个计算不是在恒定空间中执行的。
这个对吗?有没有办法让它运行或检查它是否在恒定空间中运行?
这些数组的大小约为(250、250、250)。
如果 MATLAB 没有以恒定大小运行它,是否有人对 Octave 或 Julia 或(插入类似语言)是否有任何经验?
编辑1:
我消除了多余的数组。有 10 个 258 x 258 x 338 的阵列,对应 1.67 GB。还有很多其他变量,但它们要小得多。给出的计算被简化,计算形式为:
R = (A(3Drange) .* B(3Drange) + A(new_3Drange) .* D(new_3Drange) + . . . ) ./ C
其中范围通常仅相差正负 1 或 2。
内存命令的输出:
Maximum possible array: 669 MB (7.013e+08 bytes) *
Memory available for all arrays: 1541 MB (1.616e+09 bytes) **
Memory used by MATLAB: 2209 MB (2.316e+09 bytes)
Physical Memory (RAM): 8154 MB (8.550e+09 bytes)
* Limited by contiguous virtual address space available.
** Limited by virtual address space available.
显然我应该违反第二条线。但是,代码运行良好,直到我对数组实际执行的第一个操作。也许 MATLAB 在我输入时很懒惰并且没有分配:
A=zeros(xmax+2,ymax+2,zmax+2);
但仍然在工作区告诉我变量已分配。
此代码以前曾用于较小的数组。(编辑:但似乎实际内存大小是问题,而不是每个单独数组的大小)。
对我来说非常奇怪的是为什么它在分配过程中不会出错,而是在第一次计算时出错。
编辑2:
我已经确认循环在空间中运行不恒定。在计算过程中分配了大约 0.8 GB 的内存。这是循环执行命令时资源使用情况的图像:
但是,我尝试将计算分解为多行。我在每次添加时拆分计算,并在新命令中添加每个部分,将R
其视为累加器。结果是一次分配的内存更少,但可能更频繁。这是图片:
我仍然很好奇为什么 MATLAB 不想在恒定空间中执行它。我认为这可能与索引被转移有关 - 我计划稍后对其进行更多调查,然后将所有这些放在一起作为答案,但有人可能会击败我,这也很棒。不过,现在我可以运行我正在寻找的数组大小并可以完成我的项目。