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

forcal

这里用一个不太常用的小例子说明一下。
/ V! w( ^* @5 b& v, u
matlab 2009a代码：

1. clear all
   ! z6 t/ O( i0 ?. H
2. tic9 h- E- Y- G\" h. x1 x0 c
   
3. s=0;
   \" V) Y8 P( p4 q, G
4. for k=1:1000
     @7 Y\" M- {: b! B# y: T
5.   for i=1:1000
   , [# U$ h2 S8 @5 k, d6 i
6.     a=linspace(2,2,1000);9 f3 G5 n! v, n6 ^' ~
   
7.     s=s+a(i);; Z# q5 _$ C2 U6 j
   
8.   end) v' }3 C! [- x6 n/ S! |5 i' m- f
   
9. end8 B# z3 G+ V0 j  w6 H2 O1 f) X
   
10. s
    ) `5 ]+ O1 u; ]' q
11. toc$ t+ }0 i- u5 V$ |& j; [0 r, g
    
12. 
    0 D8 T4 s/ T/ {6 u7 ]) `3 Z' G
13. s =% i! U' ~5 I% g+ D& l
    
14. + T0 c) z: q! ?* P+ T- \0 r3 M
    
15.      2000000. h9 M2 P) U& @) A\" r4 D
    
16. 
    \" T& o1 Z/ Q6 c$ ^
17. Elapsed time is 18.325900 seconds.

复制代码

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

/ P' j. V7 G, I( }7 e. b3 ?－－－－－－－－

Forcal代码：

2. !using["math","sys"];' Q% a3 b4 z+ q1 R7 d; p
   
4. mvar:
   
5. * S6 D! ?3 e: g3 p! [
6. main(:a)=( H. ?* p% x% q: V
   
8. t0=clock(),
   
9. % C- A0 l\\" |- z1 S' j& O* A6 w3 V
10. s=0,+ g* D: |' l8 N2 {- n$ p
    
12. k=0, while{k<1000,3 p# u2 X& o/ T7 R3 i' T& O
    
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)6 `+ t- L2 c% z/ l! y: O6 u/ r7 s
    
16.   k++* J/ d9 n, U6 F( ^& Z
    
18. },6 `2 K8 @' t  i+ ~4 P& c& r
    
20. s;\\" |7 B; z: o. h9 K\\" O7 f
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">' Q% a3 b4 z+ q1 R7 d; p</font><br /> mvar:<br /> <span style=\\"display:none\\">* S6 D! ?3 e: g3 p! [</span> main(:a)=<font class=\\"jammer\\">( H. ?* p% x% q: V</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">% C- A0 l\\" |- z1 S' j& O* A6 w3 V</span> s=0,<font class=\\"jammer\\">+ g* D: |' l8 N2 {- n$ p</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">3 p# u2 X& o/ T7 R3 i' T& O</font><br />   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\">6 `+ t- L2 c% z/ l! y: O6 u/ r7 s</font><br />   k++<font class=\\"jammer\\">* J/ d9 n, U6 F( ^& Z</font><br /> },<font class=\\"jammer\\">6 `2 K8 @' t  i+ ~4 P& c& r</font><br /> s;<font class=\\"jammer\\">\\" |7 B; z: o. h9 K\\" O7 f</font><br /> [clock()-t0]/1000;

结果：
  D3 C3 x8 i5 ?$ B- }9 b2000000.
40.766 秒
+ G& z" Z+ T8 n1 o) z
程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。
9 X( u! t$ A3 W3 h6 c0 Z8 {
* j7 Y" \  g* E+ f; t$ B5 T2 r－－－－－－－－－－－－
  z# G1 ?1 W3 W. ~( R+ j5 |
以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. # t* S8 t0 K! I' [8 r
4. mvar:, D, P$ ?7 Y6 l4 F5 S$ J% c2 @! w
   
6. main(:a)=
   
7. / S: D/ }\\" }. F1 C. J. ^. {0 H) i5 y
8. t0=clock(),; q5 o/ a\\" V6 \+ X* o
   
10. s=0,
    
11. ) Z4 k( `# |# B0 V% p) o! x
12. k=0, while{k<1000,
    
13. 6 ?- K8 k\\" S3 L+ x
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. 1 L( ~) g( B# O
16.   k++
    
17. ) Y! S( P8 m' Y$ b; L+ W
18. },/ u! a5 ~  t1 ]( r1 D  N
    
20. s;% s\\" K: W' j1 B1 |1 t* E, z
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\"># t* S8 t0 K! I' [8 r</span> mvar:<font class=\\"jammer\\">, D, P$ ?7 Y6 l4 F5 S$ J% c2 @! w</font><br /> main(:a)=<br /> <span style=\\"display:none\\">/ S: D/ }\\" }. F1 C. J. ^. {0 H) i5 y</span> t0=clock(),<font class=\\"jammer\\">; q5 o/ a\\" V6 \+ X* o</font><br /> s=0,<br /> <span style=\\"display:none\\">) Z4 k( `# |# B0 V% p) o! x</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">6 ?- K8 k\\" S3 L+ x</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\\">1 L( ~) g( B# O</span>   k++<br /> <span style=\\"display:none\\">) Y! S( P8 m' Y$ b; L+ W</span> },<font class=\\"jammer\\">/ u! a5 ~  t1 ]( r1 D  N</font><br /> s;<font class=\\"jammer\\">% s\\" K: W' j1 B1 |1 t* E, z</font><br /> [clock()-t0]/1000;

结果：
1 O% ]5 B+ \$ R, H* B9 C2 Z  T& O2000000.
8 k; z; @8 P0 o+ E  i4.609 秒
6 Y7 ^& R) P2 v4 J) h: l1 z) c
可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
- f$ F" b( E6 Q2 g  d' T太强大了 看不懂啊

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

我就在你背后

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

* i+ `/ @0 [* B, [* _+ F' L恩 我在研究研究，谢谢了

sxjm567

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