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

forcal

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

* g/ l4 X' g* \: o/ a( gmatlab 2009a代码：

1. clear all
   0 V2 z2 I0 @0 }7 {- G
2. tic) z- Z& |  c+ s
   
3. s=0;  k5 B% e0 B! r
   
4. for k=1:1000$ ^$ U! K% u- W) C8 @
   
5.   for i=1:1000
   9 p0 y9 `& y2 L; l. {9 p
6.     a=linspace(2,2,1000);
   7 T) d# K6 E& {1 j( g
7.     s=s+a(i);
   : X4 H6 b- N8 m4 l! B; `4 _1 e8 Q5 V
8.   end
   $ a# `3 O8 k8 Y$ M4 w3 r$ U
9. end
   2 F; E9 i, Y: s; ?# I
10. s/ k& z7 A2 ?\" w% w! ]\" f, B. u% V
    
11. toc
    ' t2 R- A: N' z  B8 F
12. ! \' p/ l8 L6 O# J
    
13. s =4 r- e% D* P, h* ^9 v\" k
    
14. 
    - N5 [: b3 G% q+ U6 t- A# Q
15.      2000000
    1 i& r& Q9 P+ A8 f) V8 }9 P9 W2 R$ X
16.   @5 ^- m% e$ L3 i$ s( a# i, u( Q
    
17. Elapsed time is 18.325900 seconds.

复制代码

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

－－－－－－－－
' C1 T/ i/ ^* n$ i- C
. m8 _8 t3 ^- x- WForcal代码：

2. !using["math","sys"];, j* ^# V/ S- t
   
4. mvar:
   
5. , X1 p( Q3 N3 Z8 i\\" g  |/ a
6. main(:a)=
   
7. 4 a8 z7 r0 ]* S1 G' Q
8. t0=clock(),
   
9. & B# \4 ^) T0 V; w$ `5 w5 ~8 s1 Y
10. s=0,
    
11. 9 Y0 _4 c- c5 a. \
12. k=0, while{k<1000,3 b4 ?& k4 Y' J6 u2 w
    
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)$ Z' p0 z1 q4 t% ]* X+ x8 l
    
16.   k++1 ^' A. A- P4 E+ @. z2 N
    
18. },* n% G3 Z# y3 D
    
20. s;/ D; z! y$ \' S3 U( C  w
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">, j* ^# V/ S- t</font><br /> mvar:<br /> <span style=\\"display:none\\">, X1 p( Q3 N3 Z8 i\\" g  |/ a</span> main(:a)=<br /> <span style=\\"display:none\\">4 a8 z7 r0 ]* S1 G' Q</span> t0=clock(),<br /> <span style=\\"display:none\\">& B# \4 ^) T0 V; w$ `5 w5 ~8 s1 Y</span> s=0,<br /> <span style=\\"display:none\\">9 Y0 _4 c- c5 a. \</span> k=0, while{k<1000,<font class=\\"jammer\\">3 b4 ?& k4 Y' J6 u2 w</font><br />   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\">$ Z' p0 z1 q4 t% ]* X+ x8 l</font><br />   k++<font class=\\"jammer\\">1 ^' A. A- P4 E+ @. z2 N</font><br /> },<font class=\\"jammer\\">* n% G3 Z# y3 D</font><br /> s;<font class=\\"jammer\\">/ D; z! y$ \' S3 U( C  w</font><br /> [clock()-t0]/1000;

结果：
2000000.
40.766 秒
" f  M) }+ q% Y" z# T0 s4 ?* e
2 G+ L5 B) j" D: Q6 J9 |7 {( S% B3 V程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。
- r* n( N5 ~* x3 I" T2 o1 K& V
－－－－－－－－－－－－

! K% M  A. F1 U以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. ! P* H: ]( l; M+ p
4. mvar:9 S2 j* Q) Q$ |; v
   
6. main(:a)=
   
7. 4 u5 U/ U$ U8 G
8. t0=clock(),+ W' B, @8 V1 G
   
10. s=0,\\" ]1 I' P$ c6 Z# N
    
12. k=0, while{k<1000,
    
13.   u( m1 F; m# x$ ~0 r
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. 7 u% X6 d/ y7 |9 X4 s4 r
16.   k++
    
17. , w4 ^1 h8 m& K- \! e- `; l
18. },
    
19. & G, v, E8 N3 _7 i
20. s;3 e0 T7 W2 o/ _4 U8 U3 O
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">! P* H: ]( l; M+ p</span> mvar:<font class=\\"jammer\\">9 S2 j* Q) Q$ |; v</font><br /> main(:a)=<br /> <span style=\\"display:none\\">4 u5 U/ U$ U8 G</span> t0=clock(),<font class=\\"jammer\\">+ W' B, @8 V1 G</font><br /> s=0,<font class=\\"jammer\\">\\" ]1 I' P$ c6 Z# N</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">  u( m1 F; m# x$ ~0 r</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\\">7 u% X6 d/ y7 |9 X4 s4 r</span>   k++<br /> <span style=\\"display:none\\">, w4 ^1 h8 m& K- \! e- `; l</span> },<br /> <span style=\\"display:none\\">& G, v, E8 N3 _7 i</span> s;<font class=\\"jammer\\">3 e0 T7 W2 o/ _4 U8 U3 O</font><br /> [clock()-t0]/1000;

结果：
2000000.
' Y% X8 e; f5 u( v  |: l# K% U4.609 秒
+ y, P* `  O' J- s9 Q8 h. Y: L. d
+ z9 Y1 Q) x/ C( t/ ^可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
  Y5 U, f1 m+ `+ l& U在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

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

" X2 X- I/ M# ^8 c7 D函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forcal如何处理垃圾对象和进行内存管理。

我就在你背后

forcal 发表于 2011-8-2 18:18
, _/ `' }  n% o8 |/ O0 e. s函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forc ...

恩 我在研究研究，谢谢了

sxjm567

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