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

forcal

这里用一个不太常用的小例子说明一下。
# y+ N! N% z  m) L- ^
matlab 2009a代码：

1. clear all
   / `) b2 R2 n2 i  V- d7 x7 L2 y
2. tic  `2 z& J: P3 R5 h; s
   
3. s=0;2 l* O2 {; _$ K5 o+ F. Y
   
4. for k=1:1000
   $ o6 Z8 e' X2 C  @
5.   for i=1:1000
   % N- g) e: B+ X0 V9 D. X
6.     a=linspace(2,2,1000);; \4 s7 S- A- X- g# E
   
7.     s=s+a(i);
   8 b) Q+ p\" i% I9 r$ M) r! _
8.   end0 o0 n) t8 }* G1 \- }9 F; [\" {
   
9. end
   0 X  t& A- w4 v, g. u' H
10. s( L+ J. }% @! T5 n5 z$ h) I
    
11. toc7 q: _1 V5 E9 k$ N/ W- U
    
12. 0 l, P+ d- N2 I* Q& P
    
13. s =4 F  h! F; P! q$ L% X
    
14. * I& L; R1 s9 i6 z7 ]1 s8 ~
    
15.      2000000
    $ z1 `2 x6 {5 T+ z( u* e5 Q
16.   L) a! L) a' o$ [& K3 i
    
17. Elapsed time is 18.325900 seconds.

复制代码

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

' Z( [0 J8 q! ]# h7 F9 H－－－－－－－－

. T6 @+ Z) q, x* {Forcal代码：

2. !using["math","sys"];
   
3. . @; y2 ]! `# |! A# M6 Z
4. mvar:- I, l' K) D) v; u' L$ [2 _
   
6. main(:a)=1 A1 S8 {( S- u% U$ ^/ _4 o
   
8. t0=clock(),' Q+ @: [5 }  X# G$ @1 W! F9 [. Y# k
   
10. s=0,+ [: U; I6 ]2 {+ r( G' C4 @
    
12. k=0, while{k<1000,5 k( j: `) ]\\" H7 y  Z
    
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)# j\\" Y% D5 e  Q3 ]4 s
    
16.   k++/ }# n\\" s! L( {, c! E# v
    
18. },: R* C- A) W$ @+ `+ m9 Y: c$ Z
    
20. s;9 Q' k6 C, D9 Z7 C
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">. @; y2 ]! `# |! A# M6 Z</span> mvar:<font class=\\"jammer\\">- I, l' K) D) v; u' L$ [2 _</font><br /> main(:a)=<font class=\\"jammer\\">1 A1 S8 {( S- u% U$ ^/ _4 o</font><br /> t0=clock(),<font class=\\"jammer\\">' Q+ @: [5 }  X# G$ @1 W! F9 [. Y# k</font><br /> s=0,<font class=\\"jammer\\">+ [: U; I6 ]2 {+ r( G' C4 @</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">5 k( j: `) ]\\" H7 y  Z</font><br />   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\"># j\\" Y% D5 e  Q3 ]4 s</font><br />   k++<font class=\\"jammer\\">/ }# n\\" s! L( {, c! E# v</font><br /> },<font class=\\"jammer\\">: R* C- A) W$ @+ `+ m9 Y: c$ Z</font><br /> s;<font class=\\"jammer\\">9 Q' k6 C, D9 Z7 C</font><br /> [clock()-t0]/1000;

结果：
2000000.
  O) K" I8 k+ V. p40.766 秒
- @: }+ o: a4 Q
程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
0 q4 D1 o' O. S. f" f因垃圾收集器多次启动，故耗时较长。

4 S4 ?% v, E3 C+ c% U! [  p－－－－－－－－－－－－
! T& f% y3 f4 d' K$ Z
: u" y' q1 ]6 w以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. $ V* N- l3 `, _\\" C3 s0 u; L$ ^4 I
4. mvar:% J$ V\\" u4 y% Q
   
6. main(:a)=- F, |9 ]1 Y; b\\" m6 S+ s
   
8. t0=clock(),+ [6 K, |$ m/ {, l! N
   
10. s=0,
    
11. ! `$ C- z+ @9 ?/ V6 B9 Y& d
12. k=0, while{k<1000,3 N4 |( }+ n5 q3 L, A
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. ' L6 g  z, }9 S9 K7 b+ d+ ~
16.   k++4 W8 O0 o5 q& u) @
    
18. },
    
19. 8 s; m/ u- g; U
20. s;* u; |$ {9 V7 Q( e3 J% x% Q& ~
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">$ V* N- l3 `, _\\" C3 s0 u; L$ ^4 I</span> mvar:<font class=\\"jammer\\">% J$ V\\" u4 y% Q</font><br /> main(:a)=<font class=\\"jammer\\">- F, |9 ]1 Y; b\\" m6 S+ s</font><br /> t0=clock(),<font class=\\"jammer\\">+ [6 K, |$ m/ {, l! N</font><br /> s=0,<br /> <span style=\\"display:none\\">! `$ C- z+ @9 ?/ V6 B9 Y& d</span> k=0, while{k<1000,<font class=\\"jammer\\">3 N4 |( }+ n5 q3 L, A</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\\">' L6 g  z, }9 S9 K7 b+ d+ ~</span>   k++<font class=\\"jammer\\">4 W8 O0 o5 q& u) @</font><br /> },<br /> <span style=\\"display:none\\">8 s; m/ u- g; U</span> s;<font class=\\"jammer\\">* u; |$ {9 V7 Q( e3 J% x% Q& ~</font><br /> [clock()-t0]/1000;

结果：
2000000.
& j' @8 z5 Z+ Q4.609 秒
& y. V/ _" E' E2 H
( D" P+ ?$ |5 z可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
5 J+ }/ `$ J+ W1 H+ d: v5 O在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

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

8 I5 ^+ u; _" n. `. l( M6 P: O函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forcal如何处理垃圾对象和进行内存管理。

我就在你背后

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

) T/ \$ E+ ]' ~4 d恩 我在研究研究，谢谢了

sxjm567

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