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

forcal

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

matlab 2009a代码：

1. clear all9 ~4 ]8 v) o% t
   
2. tic% l, K+ ^. P+ \* q6 P\" y
   
3. s=0;4 l# @; M' C) U. j- H' p- b( L/ {
   
4. for k=1:1000( ?0 ~3 \5 C2 l. ]8 r/ t( V
   
5.   for i=1:1000$ b0 V1 S) ]) |- x% D2 z
   
6.     a=linspace(2,2,1000);% Z- W& R( ?4 i& F! l\" H/ W8 r
   
7.     s=s+a(i);1 b) Z% x( `# g4 Q5 J9 X# f9 d
   
8.   end
   : |/ [% w' x( m
9. end\" A  _\" @( Q9 n9 |  a
   
10. s
    # B& N% v\" ?; W: |( v2 [4 M% L
11. toc
    + @% q+ y( p1 Y, \5 W$ J2 B2 |% V' |/ ^
12.   ~4 i8 v2 F2 ~. k& ~3 J
    
13. s =  F  r  f% U( Q: c( h' I4 s
    
14. ' n( _4 Q6 w4 K
    
15.      2000000
    4 e* r* \( g6 B7 q; t( J\" l
16. 2 S\" ?8 h/ f: {: _4 l: E* t6 L
    
17. Elapsed time is 18.325900 seconds.

复制代码

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

5 R2 k+ X- @/ d9 y0 z/ z: Y! U) u－－－－－－－－

Forcal代码：

2. !using["math","sys"];5 B* I. ^1 q' [! w
   
4. mvar:
   
5. . c$ S# }- D9 K9 b/ |1 `
6. main(:a)=- L7 r; Y. p+ H. p9 N
   
8. t0=clock(),# {, r: ]\\" O& Y) j! A2 A# M
   
10. s=0,, }4 F3 G0 s8 B
    
12. k=0, while{k<1000,
    
13. 2 X5 b4 U+ m% F; p7 c' O3 Z
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)3 J- {$ K0 b, X8 _  F3 p
    
16.   k++
    
17. / t, b; R( d$ E/ F8 B- h0 U
18. },
    
19. & R5 p4 Y6 k! Z% V0 o9 H\\" S5 a
20. s;
    
21. 6 k- w: s; {& K+ g1 C
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">5 B* I. ^1 q' [! w</font><br /> mvar:<br /> <span style=\\"display:none\\">. c$ S# }- D9 K9 b/ |1 `</span> main(:a)=<font class=\\"jammer\\">- L7 r; Y. p+ H. p9 N</font><br /> t0=clock(),<font class=\\"jammer\\"># {, r: ]\\" O& Y) j! A2 A# M</font><br /> s=0,<font class=\\"jammer\\">, }4 F3 G0 s8 B</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">2 X5 b4 U+ m% F; p7 c' O3 Z</span>   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<font class=\\"jammer\\">3 J- {$ K0 b, X8 _  F3 p</font><br />   k++<br /> <span style=\\"display:none\\">/ t, b; R( d$ E/ F8 B- h0 U</span> },<br /> <span style=\\"display:none\\">& R5 p4 Y6 k! Z% V0 o9 H\\" S5 a</span> s;<br /> <span style=\\"display:none\\">6 k- w: s; {& K+ g1 C</span> [clock()-t0]/1000;

结果：
2000000.
40.766 秒

程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
$ C. o+ C$ N/ n# e, Q" ]8 [5 D* P因垃圾收集器多次启动，故耗时较长。

5 m% s8 ^: @- q7 l4 l/ y+ `2 ~－－－－－－－－－－－－
# R3 N8 E$ L3 c8 }1 B2 k
; g/ x0 i' j' e1 Q' h# p以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];; E+ q+ r  I6 B5 j7 n) U
   
4. mvar:
   
5. 9 ^0 {1 w7 H/ m' A6 G
6. main(:a)=
   
7.   g3 X# k& A* u; U5 U
8. t0=clock(),. _8 K% @; x' \5 d# e# c* U
   
10. s=0,
    
11. 8 [+ \7 ~9 x! g* j  Y0 Q$ A
12. k=0, while{k<1000,
    
13.   w& Q2 f7 `# j. v% _/ q, t
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. 4 U( B, s- x/ x6 }: J& [6 F8 l( R
16.   k++
    
17. 0 R0 q) \# ~% {1 G; h
18. },
    
19. . ?: P1 Y8 g9 _  r; R( q# G5 r) r& B
20. s;
    
21. \\" Q' C! H  A# I7 U\\" }
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">; E+ q+ r  I6 B5 j7 n) U</font><br /> mvar:<br /> <span style=\\"display:none\\">9 ^0 {1 w7 H/ m' A6 G</span> main(:a)=<br /> <span style=\\"display:none\\">  g3 X# k& A* u; U5 U</span> t0=clock(),<font class=\\"jammer\\">. _8 K% @; x' \5 d# e# c* U</font><br /> s=0,<br /> <span style=\\"display:none\\">8 [+ \7 ~9 x! g* j  Y0 Q$ A</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">  w& Q2 f7 `# j. v% _/ q, t</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\\">4 U( B, s- x/ x6 }: J& [6 F8 l( R</span>   k++<br /> <span style=\\"display:none\\">0 R0 q) \# ~% {1 G; h</span> },<br /> <span style=\\"display:none\\">. ?: P1 Y8 g9 _  r; R( q# G5 r) r& B</span> s;<br /> <span style=\\"display:none\\">\\" Q' C! H  A# I7 U\\" }</span> [clock()-t0]/1000;

结果：
2000000.
0 ]0 m8 C5 v& }* j4.609 秒

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

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
0 q% m% x9 R! J* l4 e太强大了 看不懂啊

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

我就在你背后

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

8 ]+ P7 T6 _  x" b恩 我在研究研究，谢谢了

sxjm567

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