Home/ Questions/Q 6188835
In Process

The Archive Base Latest Questions

Editorial Team
  • 0
Asked: 2026-05-24T02:18:17+00:00

This question is inspired by the question on memory leaks in Mathematica due to

  • 0

This question is inspired by the question on memory leaks in Mathematica due to internal caching of results of intermediate computations. All these things are undocumented but are important for anyone who runs memory-intensive computations.

When trying to understand the logic of the internal caching mechanism I found something interesting. Consider the following:

$HistoryLength = 0;
(*dummy command for loading of the Root package*)
Root[# &, 1];
d = 13;
f[z_, i_] := Sum[(2 Mod[Floor[(i - 1)/2^k], 2] - 1) z^(d - k), {k, 0, d}];
(memLog = Flatten[
     Table[Root[f[z, i], j]; {SessionTime[], MemoryInUse[]/1024.^2}, {j, 1, 
       d}, {i, 1, 2^d}], 1];) // Timing
pl1 = ListLinePlot[memLog, 
  FrameLabel -> {"SessionTime, sec", "MemoryInUse, Mb"}, PlotRange -> All, 
  Frame -> True, Axes -> False]
pl2 = ListLinePlot[memLog[[All, 2]], 
  FrameLabel -> {"Point", "MemoryInUse, Mb"}, PlotRange -> All, Frame -> True,
   Axes -> False]

In fresh kernel session the output on my machine (Mathematica 7.0.1 for Windows) is always as follows:

screenshot

Can anyone explain why there is a break of the curve near the point number 8400?

Share

Leave an answer

You must login to add an answer.

Need An Account,

1 Answer

  1. Editorial Team

    The point at which the change in slope occurs is no. 8192. This equals 2^13, which is the upper bound of the fastest of your Table indices, i. When mma arrives at this point it has calculated all possible f[z, i] values. Since this call is independent of j, the next value of j (2) will generate the same range (f[z, 1], ..., f[z, 8192]) again. I assume this is where internal caching kicks in and fewer additional memory is needed.

    • 0