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

forcal

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

9 x) y  r, ^) _2 Q& Mmatlab 2009a代码：

1. clear all. ^& n2 r' U5 V$ q9 b$ _1 s
   
2. tic  f6 X, N, i, K: O1 |
   
3. s=0;6 o$ O\" s  c/ _, Y
   
4. for k=1:1000- z; v6 I- G8 k# X2 o
   
5.   for i=1:1000. d5 Y  |% g6 {8 R$ O
   
6.     a=linspace(2,2,1000);
   , l; @) B/ a9 T; h7 Q& W
7.     s=s+a(i);
   2 _$ D! z! A( k& {8 f& A
8.   end
   & ~7 ?! x4 w9 J$ ^7 _7 b# [) Q' v
9. end
   * Y: F+ d: N9 Y2 \
10. s
    4 P- b5 @5 b8 c\" I2 h4 W  g
11. toc! I( _/ ]' ^% P+ s* h8 P
    
12. 
    6 H& P& p3 E+ W7 S
13. s =5 ?' \- M- W6 q3 ^1 A1 T
    
14. 8 ^* e& _* h6 e9 ?7 F8 Q3 j
    
15.      20000003 v+ E; E- T1 A' D* V  ]\" }
    
16. 
    5 f9 U4 d  G2 m
17. Elapsed time is 18.325900 seconds.

复制代码

程序运行时内存比较稳定。
" n$ M; w! o  ^+ ~& V6 h6 r
－－－－－－－－
8 ^/ g7 ^% m% m0 K% h" r
Forcal代码：

2. !using["math","sys"];
   
3. 1 E$ U5 V2 U( A) S9 f- _6 E# w
4. mvar:
   
5. . _3 @* V  n2 J0 F
6. main(:a)=
   
7. , q! U# K/ y% Z; h; n6 i\\" o
8. t0=clock(),4 O% Q4 Q6 L\\" E& r
   
10. s=0,+ X3 S8 W( C8 x  C) T
    
12. k=0, while{k<1000,
    
13. * |3 k/ I6 E1 M  a
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. 4 f  C  P+ o/ q5 J! p: B
16.   k+++ H( L3 `) P$ _% ~% k. B2 g$ K
    
18. },% E  p! q8 p- z( ^/ u5 z\\" z
    
20. s;$ ^0 o. b  U, k5 k1 x
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">1 E$ U5 V2 U( A) S9 f- _6 E# w</span> mvar:<br /> <span style=\\"display:none\\">. _3 @* V  n2 J0 F</span> main(:a)=<br /> <span style=\\"display:none\\">, q! U# K/ y% Z; h; n6 i\\" o</span> t0=clock(),<font class=\\"jammer\\">4 O% Q4 Q6 L\\" E& r</font><br /> s=0,<font class=\\"jammer\\">+ X3 S8 W( C8 x  C) T</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">* |3 k/ I6 E1 M  a</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\\">4 f  C  P+ o/ q5 J! p: B</span>   k++<font class=\\"jammer\\">+ H( L3 `) P$ _% ~% k. B2 g$ K</font><br /> },<font class=\\"jammer\\">% E  p! q8 p- z( ^/ u5 z\\" z</font><br /> s;<font class=\\"jammer\\">$ ^0 o. b  U, k5 k1 x</font><br /> [clock()-t0]/1000;

结果：
! {& I) e) ?# o) I- h, K3 U, U1 d1 v2000000.
40.766 秒
) U* u' Z' A) _: l* ]
& q( G+ z  `# s3 `* P: D. X程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
8 l5 Q# o  z! k+ w: d" R因垃圾收集器多次启动，故耗时较长。
5 `& s" j) b6 I/ U
& S8 ^0 O; Y9 _－－－－－－－－－－－－

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

2. !using["math","sys"];
   
3. 9 R' K# q8 k; G) D0 a8 c) M: J
4. mvar:
   
5. ; h  n4 ~\\" y1 s9 y, r3 }
6. main(:a)=; b3 M; \9 R5 X, F' y\\" p1 m+ J\\" o
   
8. t0=clock(),0 J0 ?6 D# @/ o# G
   
10. s=0,
    
11. 3 ^/ ?\\" E* ~' v: w
12. k=0, while{k<1000,& {+ ?0 W( _. y& v\\" c
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. 5 C; l) v+ w( I0 C0 B
16.   k++6 z' p! ~: w# M
    
18. },2 Y+ \; O7 x- A' @. `+ A
    
20. s;
    
21. / O3 B, V% A0 }' _- G/ f+ {
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">9 R' K# q8 k; G) D0 a8 c) M: J</span> mvar:<br /> <span style=\\"display:none\\">; h  n4 ~\\" y1 s9 y, r3 }</span> main(:a)=<font class=\\"jammer\\">; b3 M; \9 R5 X, F' y\\" p1 m+ J\\" o</font><br /> t0=clock(),<font class=\\"jammer\\">0 J0 ?6 D# @/ o# G</font><br /> s=0,<br /> <span style=\\"display:none\\">3 ^/ ?\\" E* ~' v: w</span> k=0, while{k<1000,<font class=\\"jammer\\">& {+ ?0 W( _. y& v\\" c</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 C; l) v+ w( I0 C0 B</span>   k++<font class=\\"jammer\\">6 z' p! ~: w# M</font><br /> },<font class=\\"jammer\\">2 Y+ \; O7 x- A' @. `+ A</font><br /> s;<br /> <span style=\\"display:none\\">/ O3 B, V% A0 }' _- G/ f+ {</span> [clock()-t0]/1000;

结果：
2000000.
4.609 秒

' |! {" a" l4 R9 g可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
1 V* d* i- X3 O6 W$ d4 r( K在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
1 N$ ^+ G3 Q4 T6 B6 V- H太强大了 看不懂啊

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

我就在你背后

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

恩 我在研究研究，谢谢了

sxjm567

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