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

forcal

这里用一个不太常用的小例子说明一下。
  J( H9 O6 M* z3 S# ^9 S# T
matlab 2009a代码：

1. clear all( T) d\" u9 X2 [5 n4 z) i6 T( k
   
2. tic
   8 k; \( E' l6 B% I: b, Y  K
3. s=0;
   5 N\" h/ f8 N5 P, ~3 L8 j
4. for k=1:1000
   . p4 n8 {7 {* N2 l# W7 c5 N. K
5.   for i=1:10001 ?( k0 n/ Z3 B4 p, |7 T
   
6.     a=linspace(2,2,1000);
   ) t$ b4 f% }- C$ Z% q
7.     s=s+a(i);  u) b+ O5 q3 K/ {+ O5 b
   
8.   end
   ; _9 m2 v+ ?, N7 z
9. end# T% B# b\" v% G* f9 Q, j
   
10. s
    \" M; e  F& u3 ?
11. toc/ Q9 a; Q$ Q5 }2 R) f5 x
    
12. + D. a: ^* p0 _4 \5 D. y
    
13. s =
    / {9 `: |6 N6 ^1 j+ B: P& @0 k
14. 
    ; F& Y, _( z4 Y
15.      2000000
    3 h  z% c; a: W7 S* }
16. 
    ) T\" B' G) P( R8 p  b! C: ]/ z
17. Elapsed time is 18.325900 seconds.

复制代码

程序运行时内存比较稳定。
% \1 Q/ R4 i: i
# I# s; ~0 l/ E% _－－－－－－－－

Forcal代码：

2. !using["math","sys"];3 H: \! ^; X! E$ _
   
4. mvar:# }  Z- s( h  }  n  x8 C
   
6. main(:a)=
   
7. 3 h1 S/ Z) ^9 p0 E% c
8. t0=clock(),6 J& m/ h/ _& E' ^* g
   
10. s=0,. @7 w\\" \# }, j5 b4 p. y8 I5 h
    
12. k=0, while{k<1000,
    
13. ( e: X0 R9 G; y  F
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)\\" E( s, A' e. M* F$ \: G0 |
    
16.   k++% n7 s. Z5 x& s) ~  x3 ~' B
    
18. },2 d& O; P# v1 J% t8 \* a
    
20. s;2 }5 u9 }/ g7 c% s* l
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">3 H: \! ^; X! E$ _</font><br /> mvar:<font class=\\"jammer\\"># }  Z- s( h  }  n  x8 C</font><br /> main(:a)=<br /> <span style=\\"display:none\\">3 h1 S/ Z) ^9 p0 E% c</span> t0=clock(),<font class=\\"jammer\\">6 J& m/ h/ _& E' ^* g</font><br /> s=0,<font class=\\"jammer\\">. @7 w\\" \# }, j5 b4 p. y8 I5 h</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">( e: X0 R9 G; y  F</span>   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\">\\" E( s, A' e. M* F$ \: G0 |</font><br />   k++<font class=\\"jammer\\">% n7 s. Z5 x& s) ~  x3 ~' B</font><br /> },<font class=\\"jammer\\">2 d& O; P# v1 J% t8 \* a</font><br /> s;<font class=\\"jammer\\">2 }5 u9 }/ g7 c% s* l</font><br /> [clock()-t0]/1000;

结果：
) R: T& [& _# \& h$ s2000000.
4 p3 P! F  S, \0 @" c, s8 Y* n8 p8 r40.766 秒

程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。
% Z$ i( R6 q. e8 |' L
, \' f8 |5 l$ G－－－－－－－－－－－－
& u& f4 X0 t1 h, p. z. n! R2 L
: b4 Z) \" q& y以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];6 u! Y5 N5 I  R
   
4. mvar:
   
5. % R% w8 z, @' v( }+ U
6. main(:a)=
   
7. & ?) ?8 n9 a) [; q\\" G5 ]# y2 o6 F
8. t0=clock(),- K5 ~$ h+ b/ h1 d0 ?+ y. _7 Z
   
10. s=0,6 W% D4 q! N' L0 ^5 {& D! S+ ], c& h1 L
    
12. k=0, while{k<1000,' u5 B0 C  \\\" ~( ]) @* k
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. 5 K6 y) }- v5 F) X7 g) P
16.   k++8 ^# G3 [: x- t: Z# p
    
18. },, N0 w( ]- W* M  A  o
    
20. s;
    
21. \\" f  k7 D5 o+ K+ @6 e2 p$ B\\" ?# i
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">6 u! Y5 N5 I  R</font><br /> mvar:<br /> <span style=\\"display:none\\">% R% w8 z, @' v( }+ U</span> main(:a)=<br /> <span style=\\"display:none\\">& ?) ?8 n9 a) [; q\\" G5 ]# y2 o6 F</span> t0=clock(),<font class=\\"jammer\\">- K5 ~$ h+ b/ h1 d0 ?+ y. _7 Z</font><br /> s=0,<font class=\\"jammer\\">6 W% D4 q! N' L0 ^5 {& D! S+ ], c& h1 L</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">' u5 B0 C  \\\" ~( ]) @* k</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\\">5 K6 y) }- v5 F) X7 g) P</span>   k++<font class=\\"jammer\\">8 ^# G3 [: x- t: Z# p</font><br /> },<font class=\\"jammer\\">, N0 w( ]- W* M  A  o</font><br /> s;<br /> <span style=\\"display:none\\">\\" f  k7 D5 o+ K+ @6 e2 p$ B\\" ?# i</span> [clock()-t0]/1000;

结果：
0 q* r& o0 o4 K9 [3 H2000000.
0 Q9 z8 L7 G) x2 H9 B: P  b4.609 秒

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

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
0 a% }$ T+ f# d3 h( d4 O* `太强大了 看不懂啊

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

我就在你背后

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

, v6 c' E! B: d3 w恩 我在研究研究，谢谢了

sxjm567

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