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

forcal

这里用一个不太常用的小例子说明一下。

! P  g& S3 T  i6 ~, p  M. W' umatlab 2009a代码：

1. clear all
   7 b; \$ H( D; d& c. f1 t+ R- ?+ D
2. tic
   5 _( p  q8 B- V5 @\" l
3. s=0;& ]. s- Z) A3 \8 I/ D* q2 m
   
4. for k=1:1000
   \" f7 V: M7 o& h: F
5.   for i=1:1000! ~! H2 M: ]0 V, H' g9 r
   
6.     a=linspace(2,2,1000);' Y\" d/ M. [; e2 L* s
   
7.     s=s+a(i);* }- W/ |+ I. O9 q
   
8.   end
   : B: J0 c1 S0 M: j# q0 s6 i3 g\" L
9. end
   ( _: C+ K\" \) J. t1 _
10. s
    , x8 F: t' Q! b7 n, R3 Q
11. toc
    * [4 s, z0 k5 H\" N# j& g2 z
12. 
    . S% z' `0 j0 G
13. s =
    $ X' ]  i2 ~# j2 T; r! o! ?\" M
14.   a! c6 V5 F1 c! b4 R1 E$ ~5 D
    
15.      2000000
    3 O& w% W4 d' N) `
16. 
    ; h\" ?, _\" ~& ]1 i
17. Elapsed time is 18.325900 seconds.

复制代码

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

－－－－－－－－

; i; p5 s( z4 ZForcal代码：

2. !using["math","sys"];
   
3. 1 I1 _+ v) O7 C+ t9 C) A
4. mvar:
   
5. 8 |/ M9 x: r! K3 s; ]6 a
6. main(:a)=; z8 C- E3 k, t+ o! p\\" J
   
8. t0=clock(),1 M8 K# b/ d. r% }' r& N
   
10. s=0,* }+ ]% X* ?/ f/ |; ]2 E1 y( a) w
    
12. k=0, while{k<1000,
    
13. & j3 t( ^. i\\" ^: y
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. - v0 Z2 K0 f. y  n8 x( Z; j: P
16.   k++; e0 ^* U* [4 o/ |3 L. b
    
18. },+ b4 e! {$ Q  u3 M) W) U' ]/ C! E
    
20. s;
    
21. 0 B  X; g9 h0 W( U+ K: j
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">1 I1 _+ v) O7 C+ t9 C) A</span> mvar:<br /> <span style=\\"display:none\\">8 |/ M9 x: r! K3 s; ]6 a</span> main(:a)=<font class=\\"jammer\\">; z8 C- E3 k, t+ o! p\\" J</font><br /> t0=clock(),<font class=\\"jammer\\">1 M8 K# b/ d. r% }' r& N</font><br /> s=0,<font class=\\"jammer\\">* }+ ]% X* ?/ f/ |; ]2 E1 y( a) w</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">& j3 t( ^. i\\" ^: y</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\\">- v0 Z2 K0 f. y  n8 x( Z; j: P</span>   k++<font class=\\"jammer\\">; e0 ^* U* [4 o/ |3 L. b</font><br /> },<font class=\\"jammer\\">+ b4 e! {$ Q  u3 M) W) U' ]/ C! E</font><br /> s;<br /> <span style=\\"display:none\\">0 B  X; g9 h0 W( U+ K: j</span> [clock()-t0]/1000;

结果：
! F/ |0 q. P0 I, u2000000.
40.766 秒

5 j* |" `8 c1 x$ j程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
8 k% j3 `0 L9 G因垃圾收集器多次启动，故耗时较长。

6 V$ w! P4 |. w- K  Z* _; O－－－－－－－－－－－－
" O$ k1 u& K: _- `2 E0 t  |* P
以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. \\" f8 I8 b4 F# n8 V: S6 H! ^- C
4. mvar:& h* s) c' V' r  J7 f
   
6. main(:a)=. f4 j0 Q$ ^) X9 H: V, D: M
   
8. t0=clock(),
   
9. + e. i6 O, v4 C3 H2 G7 _
10. s=0,
    
11. 8 n& ?5 I5 A. }\\" B5 e* o# S
12. k=0, while{k<1000,. x  W) V& e0 Z1 p
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. . y' L! G3 L5 e8 a3 a
16.   k++
    
17. . w: Y: \- L- }) P  X) l) \5 {
18. },
    
19. 4 x7 ^5 d2 {% v1 z
20. s;
    
21. . k! \- G- t; g5 t6 Y
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">\\" f8 I8 b4 F# n8 V: S6 H! ^- C</span> mvar:<font class=\\"jammer\\">& h* s) c' V' r  J7 f</font><br /> main(:a)=<font class=\\"jammer\\">. f4 j0 Q$ ^) X9 H: V, D: M</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">+ e. i6 O, v4 C3 H2 G7 _</span> s=0,<br /> <span style=\\"display:none\\">8 n& ?5 I5 A. }\\" B5 e* o# S</span> k=0, while{k<1000,<font class=\\"jammer\\">. x  W) V& e0 Z1 p</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\\">. y' L! G3 L5 e8 a3 a</span>   k++<br /> <span style=\\"display:none\\">. w: Y: \- L- }) P  X) l) \5 {</span> },<br /> <span style=\\"display:none\\">4 x7 ^5 d2 {% v1 z</span> s;<br /> <span style=\\"display:none\\">. k! \- G- t; g5 t6 Y</span> [clock()-t0]/1000;

结果：
2000000.
3 c  h- p  o. C( `0 d7 e4.609 秒
0 S+ C7 J. f' l6 i$ G
4 x; O0 y3 j& ?  E: H可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

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

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

我就在你背后

forcal 发表于 2011-8-2 18:18
; |3 f% R! y+ C3 h* m$ Y函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forc ...

恩 我在研究研究，谢谢了

sxjm567

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