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

forcal

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

8 z) V! a0 O0 \' ^9 Smatlab 2009a代码：

1. clear all
   ' s  {! U; i. m8 @( B4 t. z2 F7 S+ d
2. tic1 f) Y/ p. `/ [# a$ L
   
3. s=0;
   5 N, V  g( a+ C) o& I! u& P
4. for k=1:10006 w$ a/ m- l/ m
   
5.   for i=1:10001 R4 k7 U! H3 T/ a% f
   
6.     a=linspace(2,2,1000);/ U& ~& _( M' J0 \
   
7.     s=s+a(i);0 h' U0 n+ d2 H7 z: _+ t\" e$ U. i/ f
   
8.   end
   / k4 l0 `2 Z- p7 T6 v
9. end1 `6 l) Y6 @( n! L  ]
   
10. s
    & F$ w) V* M! H& |: d0 y' C1 l2 `# p7 V
11. toc! \2 ~5 B% a8 m+ s  m
    
12. ' E1 P' f4 G\" x1 Z1 y) D
    
13. s =
    0 R& v+ I- `7 S0 B! ?2 V9 b
14. 
    - w5 T9 _& m) t: {3 O
15.      2000000
    , X: O2 \1 y/ s& V6 b! C& O
16. 7 a! k. F, \7 {+ N5 p$ n' B
    
17. Elapsed time is 18.325900 seconds.

复制代码

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

8 Z/ z8 T: ]% K# Y－－－－－－－－
9 q2 d: Q' B4 d/ W) ?
; Z6 s  R1 i- b% N- DForcal代码：

2. !using["math","sys"];
   
3. 2 @3 G0 }8 A4 q& ?4 p' m3 B( l
4. mvar:
   
5. - U5 F) o0 Q! T) C5 E0 b
6. main(:a)=4 m* t$ a( o+ T; L' G
   
8. t0=clock(),% V. }5 C1 X* U
   
10. s=0,
    
11. 9 @  v\\" Q4 a  d1 C3 w
12. k=0, while{k<1000,
    
13. 5 x: l! C% X5 B4 n' N
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)5 ]  N4 g0 H1 H) Q1 p5 v
    
16.   k++
    
17. : D. W5 m0 Y: T( Z3 p& O
18. },
    
19. 0 A4 q$ ~, k/ _$ H3 W
20. s;
    
21. # y& ^2 I6 c+ \+ g$ \5 L5 a2 M
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">2 @3 G0 }8 A4 q& ?4 p' m3 B( l</span> mvar:<br /> <span style=\\"display:none\\">- U5 F) o0 Q! T) C5 E0 b</span> main(:a)=<font class=\\"jammer\\">4 m* t$ a( o+ T; L' G</font><br /> t0=clock(),<font class=\\"jammer\\">% V. }5 C1 X* U</font><br /> s=0,<br /> <span style=\\"display:none\\">9 @  v\\" Q4 a  d1 C3 w</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">5 x: l! C% X5 B4 n' N</span>   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\">5 ]  N4 g0 H1 H) Q1 p5 v</font><br />   k++<br /> <span style=\\"display:none\\">: D. W5 m0 Y: T( Z3 p& O</span> },<br /> <span style=\\"display:none\\">0 A4 q$ ~, k/ _$ H3 W</span> s;<br /> <span style=\\"display:none\\"># y& ^2 I6 c+ \+ g$ \5 L5 a2 M</span> [clock()-t0]/1000;

结果：
( k  Q1 _8 A% g2000000.
40.766 秒

程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。
3 G" p# D! i6 |( E% O
3 |2 P7 K2 V$ P1 ?  P( V  l0 A－－－－－－－－－－－－

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

2. !using["math","sys"];1 t$ x  v& I) }: i, N! v) ~- k! O' Y7 u
   
4. mvar:
   
5. # U) @2 Y% B  M/ S: h1 ?
6. main(:a)=
   
7. 1 E) L; m  D% f
8. t0=clock(),
   
9. 4 ]$ j, \* F: u5 `\\" K! R' ?/ t% m# c
10. s=0,
    
11. ' e. P& ?( v- G( w7 ^6 j% T
12. k=0, while{k<1000,1 C# r7 W' O+ m5 K) [
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15.   {0 o+ ^. p; u# x: c- L
16.   k++
    
17. , T5 i8 f\\" r+ o% T& H
18. },
    
19. + X# o\\" ^+ J$ r0 ]
20. s;' w2 M7 X5 q: B\\" U) H
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">1 t$ x  v& I) }: i, N! v) ~- k! O' Y7 u</font><br /> mvar:<br /> <span style=\\"display:none\\"># U) @2 Y% B  M/ S: h1 ?</span> main(:a)=<br /> <span style=\\"display:none\\">1 E) L; m  D% f</span> t0=clock(),<br /> <span style=\\"display:none\\">4 ]$ j, \* F: u5 `\\" K! R' ?/ t% m# c</span> s=0,<br /> <span style=\\"display:none\\">' e. P& ?( v- G( w7 ^6 j% T</span> k=0, while{k<1000,<font class=\\"jammer\\">1 C# r7 W' O+ m5 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\\">  {0 o+ ^. p; u# x: c- L</span>   k++<br /> <span style=\\"display:none\\">, T5 i8 f\\" r+ o% T& H</span> },<br /> <span style=\\"display:none\\">+ X# o\\" ^+ J$ r0 ]</span> s;<font class=\\"jammer\\">' w2 M7 X5 q: B\\" U) H</font><br /> [clock()-t0]/1000;

结果：
! O! A5 b% Y1 |, \& q2 x. ]/ ^2000000.
8 o6 n& l7 p5 w1 v( m4.609 秒

9 `. |, N. r- @* j/ L* X; t可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
, m' C: I; b6 z太强大了 看不懂啊

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

我就在你背后

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

: u# m$ |  ~' [  ?恩 我在研究研究，谢谢了

sxjm567

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