Matlab 就地操作 A = A + (B*C)

Question

当对以下测试进行内存分析时，感兴趣的行上只有一个（相对）小的分配（可能是缓冲结果以添加完整的列），因此我认为这是一个就地操作，其中每个结果组件的乘法直接加到目的矩阵A的对应元素上。

A = rand(4000);
B = rand(4000,40);
C = rand(40,4000);

A = A + (B*C); % Line of interest

eigs(A,1) % to ensure the computation of A is actually done

然而，在更复杂的情况下（例如在一个循环中，前面有额外的操作），这行代码创建了一个完整的临时（大小为 A）

B*C

我想知道为什么会这样。在大规模设置中，这可能决定 Matlab 是否出现内存不足。可能是我高估了 JIT，它甚至没有在测试用例中进行就地操作；或者如我所愿，也有一些技巧可以推动 JIT 就地处理更复杂的代码。

编辑：

简单的测试用例：

复杂的测试用例：

（如果你想知道，这里：R = N = 29008）

复杂代码：

尽管有内存分析器输出，我可以在 htop 中看到，在第 244 行期间，内存使用量激增了 W 的大小，约为 6GB。

htop 在第 244 行：

第 244 行之后的htop：

Answer 1

简单地使用TaskManager （我正在使用Windows 7），无论有没有for 循环，我都会得到B*C 大小的内存峰值：

function [] = toto()
%[
    A = rand(6000);
    B = rand(6000);
    C = rand(6000);

    pause(10);
    %for k = 1:10,
        A = A + (B*C); % Line of interest
    %end
    pause(10);

    eigs(A,1) % to ensure the computation of A is actually done
%]

您确定使用htop 进行的内存分析足够准确还是测试相关？