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

forcal

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

; I3 l/ \2 B+ Y1 M. p7 ~matlab 2009a代码：

1. clear all
   2 O) Z3 k( B% |5 i\" a  M7 p% e
2. tic4 T4 j6 p' H! x) f- m2 [5 K* x
   
3. s=0;1 l! _6 H7 h; q( f: O! m( E; q3 {
   
4. for k=1:1000
   \" x* X; P5 L% ]+ i0 i\" f/ i! \
5.   for i=1:1000) d2 E2 x+ @( o
   
6.     a=linspace(2,2,1000);
   : U' W7 w/ u\" w; D5 U
7.     s=s+a(i);- _  e+ i6 f& k/ Y0 a. V' J
   
8.   end
   ' D, N* T6 u0 s! s: E
9. end
   0 T; Y0 R! G- q0 Y- R3 C. a6 F4 u) ~
10. s/ c; l% ^# i\" g. l
    
11. toc# e4 N7 G2 k5 B# V7 d2 w) j. C\" t
    
12. 
    - B  R0 {. Q4 e7 t. J; T
13. s =! H4 o8 m- e, [: u5 v
    
14. , _\" l! Q4 q/ U, J4 n. M
    
15.      2000000- V% L6 [5 a: [6 B
    
16. \" O4 q* u7 i3 q, D
    
17. Elapsed time is 18.325900 seconds.

复制代码

程序运行时内存比较稳定。
. N2 `" o. I! e9 N
3 c3 S- N( U$ U& M" t9 V) @－－－－－－－－

4 `8 C( q, [! L2 VForcal代码：

2. !using["math","sys"];, W! |6 F' Y$ r- F0 X
   
4. mvar:
   
5. ; `- U/ \2 d1 v
6. main(:a)=
   
7. ' K* ?0 u1 C: s' B. `4 K: y: O1 X
8. t0=clock(),
   
9. & W- }* Z! b2 _5 v. Q
10. s=0,
    
11. . w$ C9 F; x* k\\" d! v2 M; i( e
12. k=0, while{k<1000,
    
13. \\" M% a3 O# N1 o. w. S% e
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)# M$ v: l+ ?; ?5 Z
    
16.   k++' y1 n5 f; Y* G\\" ~& m& J9 p
    
18. },. H- E* i$ @! D\\" Z* |8 D
    
20. s;9 o; H& z) o! s4 }1 R8 G2 K5 m2 f
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">, W! |6 F' Y$ r- F0 X</font><br /> mvar:<br /> <span style=\\"display:none\\">; `- U/ \2 d1 v</span> main(:a)=<br /> <span style=\\"display:none\\">' K* ?0 u1 C: s' B. `4 K: y: O1 X</span> t0=clock(),<br /> <span style=\\"display:none\\">& W- }* Z! b2 _5 v. Q</span> s=0,<br /> <span style=\\"display:none\\">. w$ C9 F; x* k\\" d! v2 M; i( e</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">\\" M% a3 O# N1 o. w. S% e</span>   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\"># M$ v: l+ ?; ?5 Z</font><br />   k++<font class=\\"jammer\\">' y1 n5 f; Y* G\\" ~& m& J9 p</font><br /> },<font class=\\"jammer\\">. H- E* i$ @! D\\" Z* |8 D</font><br /> s;<font class=\\"jammer\\">9 o; H& z) o! s4 }1 R8 G2 K5 m2 f</font><br /> [clock()-t0]/1000;

结果：
; X9 g  A% l  d) X% V6 V2000000.
40.766 秒

程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。

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

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

2. !using["math","sys"];
   
3. ' D; W+ n! n2 W  N& F* _  A
4. mvar:
   
5. 9 x  H& X' T6 ]* y  _9 k5 m
6. main(:a)=
   
7. * U8 |9 C8 _\\" ^6 D
8. t0=clock(),
   
9. 4 n; z6 `/ h6 D: O; Y
10. s=0,* L( \2 ?) J8 g. \# s
    
12. k=0, while{k<1000,
    
13. 1 a6 G$ v. c1 [  l8 f7 W5 Y, y
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. # ?5 z7 j; F3 \4 w, U& w  M
16.   k++
    
17. 5 \7 B( A  g/ H$ ~2 f  g
18. },
    
19. & x/ j' N& U. w; }3 U
20. s;
    
21. 2 [5 C( f- j( d, d( R
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">' D; W+ n! n2 W  N& F* _  A</span> mvar:<br /> <span style=\\"display:none\\">9 x  H& X' T6 ]* y  _9 k5 m</span> main(:a)=<br /> <span style=\\"display:none\\">* U8 |9 C8 _\\" ^6 D</span> t0=clock(),<br /> <span style=\\"display:none\\">4 n; z6 `/ h6 D: O; Y</span> s=0,<font class=\\"jammer\\">* L( \2 ?) J8 g. \# s</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">1 a6 G$ v. c1 [  l8 f7 W5 Y, y</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\\"># ?5 z7 j; F3 \4 w, U& w  M</span>   k++<br /> <span style=\\"display:none\\">5 \7 B( A  g/ H$ ~2 f  g</span> },<br /> <span style=\\"display:none\\">& x/ j' N& U. w; }3 U</span> s;<br /> <span style=\\"display:none\\">2 [5 C( f- j( d, d( R</span> [clock()-t0]/1000;

结果：
2000000.
3 y( ~! K# |. L: J6 F) i8 e2 `4.609 秒
0 s4 m( E" ^( p9 S( K# \/ Q
8 `# V9 H* J+ H* f可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
  u6 G7 I6 r3 I在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

sxjm567

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

我就在你背后

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

7 _* q- P6 ^* _3 v恩 我在研究研究，谢谢了

forcal

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

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

我就在你背后

太强大了 看不懂啊