极限测试之Matlab与Forcal动态内存管理效率

forcal

这里用一个不太常用的小例子说明一下。

! L7 O5 M2 {. wmatlab 2009a代码：

1. clear all
   0 `5 J6 r. G4 a) t  l; X
2. tic
   . R' j; \) q) B5 I3 g
3. s=0;
   / O3 e+ Q( A/ S1 |
4. for k=1:1000
   , ?$ T, O0 s0 k, x# t* q
5.   for i=1:1000( ^2 p\" d  X0 ?! Y
   
6.     a=linspace(2,2,1000);  `$ b* B0 }1 f( g8 y
   
7.     s=s+a(i);
   \" N+ L9 `+ O8 b, H  P
8.   end) h5 G9 f6 X- U) R: Q. S2 ?, K2 \3 L! q
   
9. end, o+ g  \7 A! x
   
10. s+ h5 F, S: q) _+ y# ]
    
11. toc6 T9 g, z1 h; h* U8 Z2 A/ k8 D! d
    
12. 
    & j( V4 ^  ]3 S
13. s =
    3 Q4 c7 G0 t\" ~. y7 d' j( G+ ~5 m) ]
14. 
    ' f. e. N. T\" I: a( b5 J- H
15.      2000000
    / l. C- D& a- K& I
16. . a* u) {\" \6 {\" H' L
    
17. Elapsed time is 18.325900 seconds.

复制代码

程序运行时内存比较稳定。
) z' R; d: A* K! ?: H0 Q1 o
* r* \7 W5 D9 e+ `" J－－－－－－－－

Forcal代码：

2. !using["math","sys"];. P3 r8 I, d\\" |( \5 i* Y( R6 J: D( ^
   
4. mvar:# G  B! }. I, m; j3 Z0 I
   
6. main(:a)=; `\\" D+ Y3 R# \! q1 z; Z, t
   
8. t0=clock(),+ y9 ^& U: r- Q4 U
   
10. s=0,
    
11.   L5 ^9 e' m9 Z+ I\\" v# e) _3 [
12. k=0, while{k<1000,
    
13. : n2 b  p, r6 C( p9 p
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. - m' g1 V% V9 T. T
16.   k++
    
17. . [! l9 Y: P2 U# [
18. },0 m; Q% D8 A$ O$ b3 X
    
20. s;
    
21. % x! A8 K( g6 ~3 S: U5 Z! Q6 A' z
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">. P3 r8 I, d\\" |( \5 i* Y( R6 J: D( ^</font><br /> mvar:<font class=\\"jammer\\"># G  B! }. I, m; j3 Z0 I</font><br /> main(:a)=<font class=\\"jammer\\">; `\\" D+ Y3 R# \! q1 z; Z, t</font><br /> t0=clock(),<font class=\\"jammer\\">+ y9 ^& U: r- Q4 U</font><br /> s=0,<br /> <span style=\\"display:none\\">  L5 ^9 e' m9 Z+ I\\" v# e) _3 [</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">: n2 b  p, r6 C( p9 p</span>   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<br /> <span style=\\"display:none\\">- m' g1 V% V9 T. T</span>   k++<br /> <span style=\\"display:none\\">. [! l9 Y: P2 U# [</span> },<font class=\\"jammer\\">0 m; Q% D8 A$ O$ b3 X</font><br /> s;<br /> <span style=\\"display:none\\">% x! A8 K( g6 ~3 S: U5 Z! Q6 A' z</span> [clock()-t0]/1000;

结果：
2000000.
+ K1 L' s/ P+ D( g40.766 秒

程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。

/ d( O& d2 J0 k, [2 g* @7 T3 b－－－－－－－－－－－－
% @/ m1 j* a( ~1 R) T& \. y
以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. 0 k6 r9 M+ c, D. [
4. mvar:
   
5. 1 u0 K$ h\\" V6 ~4 K  J) D
6. main(:a)=# y8 A+ m- m0 \$ k
   
8. t0=clock(),* B\\" h! o% Q8 w/ N
   
10. s=0,  X: F# v2 M  \: a4 i1 W
    
12. k=0, while{k<1000,5 s9 D5 }, r$ y0 ~$ W3 J1 l
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. + x$ o5 N4 f! c& X
16.   k++
    
17. : @' h1 ^7 Q8 T2 q9 k
18. },, H, V: D9 ~! l\\" @, _
    
20. s;6 ]% C* @3 |  K
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">0 k6 r9 M+ c, D. [</span> mvar:<br /> <span style=\\"display:none\\">1 u0 K$ h\\" V6 ~4 K  J) D</span> main(:a)=<font class=\\"jammer\\"># y8 A+ m- m0 \$ k</font><br /> t0=clock(),<font class=\\"jammer\\">* B\\" h! o% Q8 w/ N</font><br /> s=0,<font class=\\"jammer\\">  X: F# v2 M  \: a4 i1 W</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">5 s9 D5 }, r$ y0 ~$ W3 J1 l</font><br />   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)<br /> <span style=\\"display:none\\">+ x$ o5 N4 f! c& X</span>   k++<br /> <span style=\\"display:none\\">: @' h1 ^7 Q8 T2 q9 k</span> },<font class=\\"jammer\\">, H, V: D9 ~! l\\" @, _</font><br /> s;<font class=\\"jammer\\">6 ]% C* @3 |  K</font><br /> [clock()-t0]/1000;

结果：
4 ]( z, `# F" Q2 X5 [2 b* y5 B1 x/ I2000000.
4.609 秒

可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
2 h1 \* U1 T. U1 h$ ^3 g在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
太强大了 看不懂啊

函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forcal如何处理垃圾对象和进行内存管理。

我就在你背后

forcal 发表于 2011-8-2 18:18
函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forc ...

恩 我在研究研究，谢谢了

sxjm567

我反复看了多遍，好帖，得支持