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

forcal

这里用一个不太常用的小例子说明一下。
2 e" ^% t$ ?7 |0 B9 p4 C3 i( K, W; x$ m2 o2 p
- {/ O2 P- |" \. O3 z5 J: nmatlab 2009a代码：

1. clear all
   7 n& {+ Q0 ]4 p: g# M8 _- ]8 a- R; z
2. tic
   9 N# q2 b7 Q3 i$ f: p
3. s=0;- V\" ?& D4 E& g; [, z$ w
   
4. for k=1:1000
   , y+ k2 q' {/ }% ?
5.   for i=1:10006 u) _# c7 f; T, ~& q
   
6.     a=linspace(2,2,1000);
   % S8 B  Y$ e9 o6 r
7.     s=s+a(i);* ]! V2 e% U! w0 f+ H
   
8.   end
   ( {% P( j+ f, v
9. end) X& ]  j6 J6 d! u
   
10. s8 y- X. q, I6 H/ F: @% f
    
11. toc
    - v# k, O9 r% y\" b; \
12. 
    5 D0 h7 u  K5 e8 a$ T$ m
13. s =
    ! H8 g  }6 o  M2 n* N: I
14. - t. V5 \9 q' S
    
15.      2000000, A\" t, i$ V' @: f\" r+ b
    
16. 
    1 `\" e% X) p+ m8 n
17. Elapsed time is 18.325900 seconds.

复制代码

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

, |$ p5 T$ k9 s" w－－－－－－－－
4 u6 I$ n- u! N! l) G0 y% M
Forcal代码：

2. !using["math","sys"];
   
3. 9 V# [$ g$ u  m4 b+ q$ U
4. mvar:; r; _+ z\\" @; @% ?' w, k
   
6. main(:a)=
   
7. ; k) B/ a* @; M) ^) Q6 R. ?. N
8. t0=clock(),) f: x6 Y3 H, n) D2 I
   
10. s=0,+ C. q3 c) r% Y8 |( z% f0 d
    
12. k=0, while{k<1000,; ?# H% V7 n# H# h7 b0 Q
    
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. . ]2 ]6 T3 V) X& K5 Y0 f* c
16.   k++' H& G$ H6 p4 H$ U\\" n$ }. U' y
    
18. },. w* T* Z/ P. O1 _- p/ [2 s. Y
    
20. s;% |' e3 U5 j. k+ i* D2 `
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">9 V# [$ g$ u  m4 b+ q$ U</span> mvar:<font class=\\"jammer\\">; r; _+ z\\" @; @% ?' w, k</font><br /> main(:a)=<br /> <span style=\\"display:none\\">; k) B/ a* @; M) ^) Q6 R. ?. N</span> t0=clock(),<font class=\\"jammer\\">) f: x6 Y3 H, n) D2 I</font><br /> s=0,<font class=\\"jammer\\">+ C. q3 c) r% Y8 |( z% f0 d</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">; ?# H% V7 n# H# h7 b0 Q</font><br />   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<br /> <span style=\\"display:none\\">. ]2 ]6 T3 V) X& K5 Y0 f* c</span>   k++<font class=\\"jammer\\">' H& G$ H6 p4 H$ U\\" n$ }. U' y</font><br /> },<font class=\\"jammer\\">. w* T* Z/ P. O1 _- p/ [2 s. Y</font><br /> s;<font class=\\"jammer\\">% |' e3 U5 j. k+ i* D2 `</font><br /> [clock()-t0]/1000;

结果：
2000000.
( d$ E1 J( s4 f& q" ?40.766 秒
7 i8 x; H# l$ Z0 ?
程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
5 V3 g/ {. {$ U0 h8 L因垃圾收集器多次启动，故耗时较长。

－－－－－－－－－－－－

6 O* p1 e9 n. h& L; w2 w以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. ! d\\" a9 o( ^) h( ]
4. mvar:% ]& m7 Y! j4 y( y; j
   
6. main(:a)=
   
7. ( U) ~8 y2 H3 R\\" l! n
8. t0=clock(),
   
9.   Q6 a4 G. ^! v/ w/ J- }& A
10. s=0,
    
11. ' o2 o& s2 ]& b
12. k=0, while{k<1000,- k/ _. C; _+ y& g. m- w( x
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. $ W  p1 N+ H- Y: N
16.   k++
    
17. 0 ]/ E& X2 U7 G0 J9 ]+ h, g
18. },) X2 K6 u( F4 E# _) y5 B( W
    
20. s;3 X: F% Z2 @4 H6 t9 k7 c8 s
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">! d\\" a9 o( ^) h( ]</span> mvar:<font class=\\"jammer\\">% ]& m7 Y! j4 y( y; j</font><br /> main(:a)=<br /> <span style=\\"display:none\\">( U) ~8 y2 H3 R\\" l! n</span> t0=clock(),<br /> <span style=\\"display:none\\">  Q6 a4 G. ^! v/ w/ J- }& A</span> s=0,<br /> <span style=\\"display:none\\">' o2 o& s2 ]& b</span> k=0, while{k<1000,<font class=\\"jammer\\">- k/ _. C; _+ y& g. m- w( x</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\\">$ W  p1 N+ H- Y: N</span>   k++<br /> <span style=\\"display:none\\">0 ]/ E& X2 U7 G0 J9 ]+ h, g</span> },<font class=\\"jammer\\">) X2 K6 u( F4 E# _) y5 B( W</font><br /> s;<font class=\\"jammer\\">3 X: F% Z2 @4 H6 t9 k7 c8 s</font><br /> [clock()-t0]/1000;

结果：
2000000.
4.609 秒

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

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
. r3 s. M1 V8 U- ~% `/ v太强大了 看不懂啊

, k9 O# i: |  T. ?7 C7 z5 e函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forcal如何处理垃圾对象和进行内存管理。

我就在你背后

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

恩 我在研究研究，谢谢了

sxjm567

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