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

forcal

这里用一个不太常用的小例子说明一下。
) z. K2 w  f; g+ i, O$ V( X
1 }8 @  |, D  ?) x6 M' _matlab 2009a代码：

1. clear all
   % Z' Z% @; i6 G6 n
2. tic! Z. ^4 F9 q2 w4 T+ c
   
3. s=0;
   + V# s  _1 V0 q4 F, F
4. for k=1:1000
   ! ^# A% o7 M! _\" ^6 G
5.   for i=1:1000
   / I4 j, k6 K* e9 X
6.     a=linspace(2,2,1000);
   ! g7 F  I% R9 ?& y  r& C! t
7.     s=s+a(i);) |\" R- Z5 _) U. u: ~. @
   
8.   end
   ( k\" Q  L- D( G1 N3 f3 R* p
9. end- @' o5 e! k4 c) x
   
10. s
    + a1 S2 o8 u4 h0 _
11. toc
    2 E! U' b1 b7 A+ \) J/ K9 N
12. 
    % i2 ~7 u; b: a/ X; L) ~
13. s =
    \" ?% ^6 K) e! K
14. 
    # m\" {3 R) t( ]
15.      20000005 d5 o: H* }# S: s
    
16. 1 Q1 K; D\" X, |
    
17. Elapsed time is 18.325900 seconds.

复制代码

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

－－－－－－－－

Forcal代码：

2. !using["math","sys"];+ O; L) y5 H* ~  q
   
4. mvar:
   
5. - A  l6 H2 z+ Y( r' Y1 w6 B% e
6. main(:a)=
   
7. % ~% Z+ p1 l7 h' I, X4 D1 P2 O
8. t0=clock(),
   
9. : N2 e7 ]. i$ _  Q8 T# N4 e4 \
10. s=0,
    
11. % P9 x( x. Z, Y* V4 x( k
12. k=0, while{k<1000,
    
13. / S& i. b3 i) D1 d2 F
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. : @% {( N+ k# N
16.   k++
    
17. % J! ~3 E! j- S% T
18. },/ z9 B( Y- B2 k1 m( ~$ h0 I
    
20. s;  ]. Y\\" R3 c+ O0 d; \4 A0 s% d& V
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ O; L) y5 H* ~  q</font><br /> mvar:<br /> <span style=\\"display:none\\">- A  l6 H2 z+ Y( r' Y1 w6 B% e</span> main(:a)=<br /> <span style=\\"display:none\\">% ~% Z+ p1 l7 h' I, X4 D1 P2 O</span> t0=clock(),<br /> <span style=\\"display:none\\">: N2 e7 ]. i$ _  Q8 T# N4 e4 \</span> s=0,<br /> <span style=\\"display:none\\">% P9 x( x. Z, Y* V4 x( k</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">/ S& i. b3 i) D1 d2 F</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\\">: @% {( N+ k# N</span>   k++<br /> <span style=\\"display:none\\">% J! ~3 E! j- S% T</span> },<font class=\\"jammer\\">/ z9 B( Y- B2 k1 m( ~$ h0 I</font><br /> s;<font class=\\"jammer\\">  ]. Y\\" R3 c+ O0 d; \4 A0 s% d& V</font><br /> [clock()-t0]/1000;

结果：
2000000.
40.766 秒

+ [& P- {* _5 y' h程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。
( |% o: J; w% K4 J, a: G
9 b1 ^4 U( v% U. |2 h. [－－－－－－－－－－－－
# i- d$ x8 }( I% h# e1 |9 @! `
1 j, t- C# @' \9 @, _. E% [7 A以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];+ P6 H( z% K; h# F3 K, \# k
   
4. mvar:$ n/ ~8 B- m8 g, h0 W
   
6. main(:a)=
   
7. 8 [- }. W9 S0 f! b4 N
8. t0=clock(),
   
9. 3 Z1 V6 X; m+ R6 b9 L4 `\\" G
10. s=0,: @: j* J* r) i% f5 f
    
12. k=0, while{k<1000,
    
13. 6 U3 V7 p\\" i& w\\" w& K9 S5 M% O
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. 6 \( ?% I- n2 j2 e: f: j. l7 O/ M, x
16.   k++
    
17. ! e+ O1 a/ c) O, n# V% L, c
18. },) p1 B8 O# D; v
    
20. s;
    
21. # N* u. X) a- f+ M\\" E; G9 W
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ P6 H( z% K; h# F3 K, \# k</font><br /> mvar:<font class=\\"jammer\\">$ n/ ~8 B- m8 g, h0 W</font><br /> main(:a)=<br /> <span style=\\"display:none\\">8 [- }. W9 S0 f! b4 N</span> t0=clock(),<br /> <span style=\\"display:none\\">3 Z1 V6 X; m+ R6 b9 L4 `\\" G</span> s=0,<font class=\\"jammer\\">: @: j* J* r) i% f5 f</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">6 U3 V7 p\\" i& w\\" w& K9 S5 M% O</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\\">6 \( ?% I- n2 j2 e: f: j. l7 O/ M, x</span>   k++<br /> <span style=\\"display:none\\">! e+ O1 a/ c) O, n# V% L, c</span> },<font class=\\"jammer\\">) p1 B8 O# D; v</font><br /> s;<br /> <span style=\\"display:none\\"># N* u. X) a- f+ M\\" E; G9 W</span> [clock()-t0]/1000;

结果：
2000000.
! }8 |0 ]6 ^3 i4.609 秒

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

我就在你背后

太强大了 看不懂啊

forcal

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

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

我就在你背后

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

恩 我在研究研究，谢谢了

sxjm567

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