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

forcal

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

& J' m; k/ `5 B& O3 Qmatlab 2009a代码：

1. clear all
   8 e( _5 Q' }% c9 x, k. m
2. tic
   3 @0 \) }5 x; @6 m1 @: |1 V
3. s=0;
   ( `6 @5 _, P9 W% ?
4. for k=1:1000* N( d, y+ F- c( h& y+ d
   
5.   for i=1:1000
   ' M& j  ]3 K9 u5 w/ m
6.     a=linspace(2,2,1000);
   7 u( @9 X' B\" L* K$ ~
7.     s=s+a(i);+ P- T9 l9 h. ^
   
8.   end% f3 B, q4 }; c0 p3 [
   
9. end
   0 t* u& j# S+ i7 l' |2 j
10. s
    % C, ^% G( X5 v6 E% g* Y
11. toc
    # r* u. O9 M- y2 J2 J' j1 ~
12. 
    , n1 G1 I& L# y( m! J6 N# w: h
13. s =+ b3 u9 j4 g5 ^
    
14. 
    6 l# z7 d, u( D1 c
15.      2000000
    ( U+ Z5 l8 B7 I6 O7 N2 S( c) M
16. 
    8 `& d2 c- }; m, D0 d4 i0 |. e! S
17. Elapsed time is 18.325900 seconds.

复制代码

程序运行时内存比较稳定。

－－－－－－－－

+ H; B8 F/ ?! J/ q; R4 v0 o& iForcal代码：

2. !using["math","sys"];
   
3. 6 |, O0 V2 @& l( Q
4. mvar:
   
5. 2 \9 K# T: C* \0 ]
6. main(:a)=
   
7. ! @  o: N$ v3 y# j$ {
8. t0=clock(),% F( U0 E3 }- {  B
   
10. s=0,
    
11. $ v& ]3 }! l( e5 _/ |) M6 M+ `* y4 N
12. k=0, while{k<1000,. n( d( \; P. y3 R; _
    
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)* w5 H9 G* p3 v+ m) I( G# _. a
    
16.   k++
    
17. 0 w$ b7 F& N: |0 F3 _/ }+ v
18. },: o) q/ ^- {* [- E& E; I  C
    
20. s;$ z3 l+ R9 B: B
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">6 |, O0 V2 @& l( Q</span> mvar:<br /> <span style=\\"display:none\\">2 \9 K# T: C* \0 ]</span> main(:a)=<br /> <span style=\\"display:none\\">! @  o: N$ v3 y# j$ {</span> t0=clock(),<font class=\\"jammer\\">% F( U0 E3 }- {  B</font><br /> s=0,<br /> <span style=\\"display:none\\">$ v& ]3 }! l( e5 _/ |) M6 M+ `* y4 N</span> k=0, while{k<1000,<font class=\\"jammer\\">. n( d( \; P. y3 R; _</font><br />   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\">* w5 H9 G* p3 v+ m) I( G# _. a</font><br />   k++<br /> <span style=\\"display:none\\">0 w$ b7 F& N: |0 F3 _/ }+ v</span> },<font class=\\"jammer\\">: o) q/ ^- {* [- E& E; I  C</font><br /> s;<font class=\\"jammer\\">$ z3 l+ R9 B: B</font><br /> [clock()-t0]/1000;

结果：
! j+ ?* q- H1 Q2 g+ }6 R2000000.
1 s$ h0 M) P% p8 y! s" q40.766 秒
: {1 ^3 F" B; s) N% I' t& j+ v
8 g4 ?: C, L+ Q/ f' T' g程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。

－－－－－－－－－－－－
5 Q# U( }  B0 U* w+ w. X, J
以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. : Z) W\\" r0 t1 }. `# o6 r5 R
4. mvar:# l; `. C7 F& g$ o\\" ^: ]  C9 Y. l
   
6. main(:a)=( k* B% I. P( W0 z  X  u
   
8. t0=clock(),
   
9. 0 P6 ^9 l/ m/ \  ?; {. \7 V* L
10. s=0,# P4 A1 y* d6 h
    
12. k=0, while{k<1000,
    
13. 0 m2 f) F% C) B: ~1 E: x0 [4 j4 y6 r
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. ' t4 d: Z- K1 O# T4 o+ S\\" V: T
16.   k++
    
17. 6 ]\\" }\\" e+ _* m\\" U) W; P4 h8 \+ F
18. },
    
19. 9 ]- p# X/ W' ^! d+ h& n
20. s;
    
21. ( E: F5 w) c' _6 R  b# [0 ^
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">: Z) W\\" r0 t1 }. `# o6 r5 R</span> mvar:<font class=\\"jammer\\"># l; `. C7 F& g$ o\\" ^: ]  C9 Y. l</font><br /> main(:a)=<font class=\\"jammer\\">( k* B% I. P( W0 z  X  u</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">0 P6 ^9 l/ m/ \  ?; {. \7 V* L</span> s=0,<font class=\\"jammer\\"># P4 A1 y* d6 h</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">0 m2 f) F% C) B: ~1 E: x0 [4 j4 y6 r</span>   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\\">' t4 d: Z- K1 O# T4 o+ S\\" V: T</span>   k++<br /> <span style=\\"display:none\\">6 ]\\" }\\" e+ _* m\\" U) W; P4 h8 \+ F</span> },<br /> <span style=\\"display:none\\">9 ]- p# X/ W' ^! d+ h& n</span> s;<br /> <span style=\\"display:none\\">( E: F5 w) c' _6 R  b# [0 ^</span> [clock()-t0]/1000;

结果：
4 s1 f# y7 e4 t2000000.
+ u- n2 B* I5 F; |* ^* Z% j9 O  ]4 |4.609 秒

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

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
' m9 f+ m& F# V6 g$ [太强大了 看不懂啊

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

我就在你背后

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

恩 我在研究研究，谢谢了

sxjm567

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