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

forcal

这里用一个不太常用的小例子说明一下。
! P& N& e2 n# _% ~' O
* w, ?; g. H8 _( A' N8 e# imatlab 2009a代码：

1. clear all4 v4 _: ~2 p* |
   
2. tic8 f& t& v  F7 N( z9 o\" a' x
   
3. s=0;
   + w* N( C, {8 k$ `0 F
4. for k=1:1000
   6 R4 g8 k* N\" a/ s
5.   for i=1:1000
   . ]( g4 Z\" c# r, a  w
6.     a=linspace(2,2,1000);
   5 H  g8 {; W: T4 C, A5 D8 ?  r
7.     s=s+a(i);
   : I% B- d2 Y1 |0 _3 [. S: B
8.   end3 e2 ?0 [/ ~; Z1 J% E6 w
   
9. end7 V) @7 O+ h- r. C; S& g6 T
   
10. s
    $ I- r/ D6 T\" _$ e! g0 d9 o
11. toc8 x) o/ M8 `) P
    
12. ) t, l\" I\" K+ z8 P  Y+ ]+ `% a5 p( ~
    
13. s =
    \" K# t/ @0 p4 \/ w7 D
14. ! \1 T' }2 D4 \( \% h  X) l
    
15.      2000000
      y8 D; `3 g\" n4 \% `2 c; d
16. 
    6 f( D9 i) A# j- g
17. Elapsed time is 18.325900 seconds.

复制代码

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

! ^# t0 x. ]+ g5 t( @% y. k－－－－－－－－

Forcal代码：

2. !using["math","sys"];4 S\\" W8 y1 Q( v
   
4. mvar:  D$ @/ B8 y6 t& t
   
6. main(:a)=
   
7. + B8 l. h, R3 v( [, A
8. t0=clock(),+ d6 W' h, x: h7 o
   
10. s=0,\\" Z' F% x8 H! x8 O$ {+ C5 g5 [
    
12. k=0, while{k<1000,
    
13. 4 T! v/ t+ N: _# t8 M  Z% C
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)/ E' X* s9 N; j8 V- b\\" o# h4 z
    
16.   k++1 t. G4 i* b; A5 {
    
18. },/ ]\\" r9 o: i/ k: H3 y+ c: L
    
20. s;9 ?0 ]3 ?9 o  R! n4 _5 Z
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">4 S\\" W8 y1 Q( v</font><br /> mvar:<font class=\\"jammer\\">  D$ @/ B8 y6 t& t</font><br /> main(:a)=<br /> <span style=\\"display:none\\">+ B8 l. h, R3 v( [, A</span> t0=clock(),<font class=\\"jammer\\">+ d6 W' h, x: h7 o</font><br /> s=0,<font class=\\"jammer\\">\\" Z' F% x8 H! x8 O$ {+ C5 g5 [</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">4 T! v/ t+ N: _# t8 M  Z% C</span>   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\">/ E' X* s9 N; j8 V- b\\" o# h4 z</font><br />   k++<font class=\\"jammer\\">1 t. G4 i* b; A5 {</font><br /> },<font class=\\"jammer\\">/ ]\\" r9 o: i/ k: H3 y+ c: L</font><br /> s;<font class=\\"jammer\\">9 ?0 ]3 ?9 o  R! n4 _5 Z</font><br /> [clock()-t0]/1000;

结果：
* M! O% B2 b9 X( Y5 h( D2000000.
40.766 秒

7 ^# p5 d0 u% T: [$ s2 ]% U: Q: U程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
5 H7 K; C0 s: i因垃圾收集器多次启动，故耗时较长。

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

以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. . b0 P# p+ M8 E9 |\\" z
4. mvar:: v\\" |5 M) ]+ z4 O2 T
   
6. main(:a)=
   
7. ! F% |! c& @5 t! c
8. t0=clock(),7 |! B9 G8 Y( q: l  Q0 C+ s7 Y
   
10. s=0,
    
11. , E  D1 O) h) X1 s
12. k=0, while{k<1000,
    
13. 6 a$ u: G6 t: |& e
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. 3 {+ I0 o. P  u5 d3 o
16.   k++\\" _- D0 N\\" a7 U: l' I( j$ h
    
18. },
    
19. & ]0 c' X/ k- u
20. s;
    
21. * t/ q8 M# V\\" S; V
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">. b0 P# p+ M8 E9 |\\" z</span> mvar:<font class=\\"jammer\\">: v\\" |5 M) ]+ z4 O2 T</font><br /> main(:a)=<br /> <span style=\\"display:none\\">! F% |! c& @5 t! c</span> t0=clock(),<font class=\\"jammer\\">7 |! B9 G8 Y( q: l  Q0 C+ s7 Y</font><br /> s=0,<br /> <span style=\\"display:none\\">, E  D1 O) h) X1 s</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">6 a$ u: G6 t: |& e</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\\">3 {+ I0 o. P  u5 d3 o</span>   k++<font class=\\"jammer\\">\\" _- D0 N\\" a7 U: l' I( j$ h</font><br /> },<br /> <span style=\\"display:none\\">& ]0 c' X/ k- u</span> s;<br /> <span style=\\"display:none\\">* t/ q8 M# V\\" S; V</span> [clock()-t0]/1000;

结果：
5 d; w: t8 o1 |! Z" g6 E- o2000000.
1 Z  Y! g! D/ p4.609 秒
3 v' R# i) o. [. w" a" ^
可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
& [* i6 ~! e. p太强大了 看不懂啊

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

我就在你背后

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

恩 我在研究研究，谢谢了

sxjm567

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