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

forcal

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

, u: ^2 n/ t% n* H( B: q+ bmatlab 2009a代码：

1. clear all0 @0 s: C& W' c; c7 G$ h% F
   
2. tic
   8 }: @6 l. d) n/ k4 J8 O
3. s=0;
   # t$ n, @! k) U  r' W: q% D0 E
4. for k=1:1000
   6 X& l/ l9 E* ~! g+ E1 W$ C* v
5.   for i=1:1000$ W7 }: E, A& L
   
6.     a=linspace(2,2,1000);3 M\" q4 @' L4 i& A
   
7.     s=s+a(i);
     R( n# m, n\" s\" r% B
8.   end/ i5 k' J% M& O9 L5 v- h
   
9. end* X$ O/ |' i# c/ a
   
10. s
    - q  C4 C, |) p. J/ c4 _
11. toc
    : [  f0 S2 f  K- y2 \5 ]7 a
12. / G7 j: ?9 r\" @- a
    
13. s =
    1 O1 X5 _/ e- M3 P+ W9 ]7 K
14. * ^\" `/ ^: {% q
    
15.      2000000
    . O9 m9 D/ ?9 ]
16. 
    $ l0 j/ {# o4 H  u; H
17. Elapsed time is 18.325900 seconds.

复制代码

程序运行时内存比较稳定。
- b. _: v( ~, ^+ B& k
+ w% F' S2 d* z/ F3 M% J& R－－－－－－－－
' a' i& [+ q9 Z' F) q6 f
Forcal代码：

2. !using["math","sys"];
   
3. ! n1 e- Q  \  x
4. mvar:
   
5. 5 Y; w- }1 {: l1 y\\" t. h8 v
6. main(:a)=
   
7. $ y9 ]# J- M, L' g7 E$ l& S
8. t0=clock(),% {, g) q& B' J( s1 X1 v- }2 Y) B
   
10. s=0,
    
11. 3 U' d4 b+ `# c0 i
12. k=0, while{k<1000,8 u3 j\\" T9 g% z
    
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. . N9 r$ x' [2 _$ W$ ~9 n
16.   k++! V( g) \5 i% a8 A/ O' _! g: L( O: G
    
18. },
    
19. ' f\\" d% X: S( j9 e0 B
20. s;
    
21.   N) t5 H4 q' J+ R% E) E\\" G; W
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">! n1 e- Q  \  x</span> mvar:<br /> <span style=\\"display:none\\">5 Y; w- }1 {: l1 y\\" t. h8 v</span> main(:a)=<br /> <span style=\\"display:none\\">$ y9 ]# J- M, L' g7 E$ l& S</span> t0=clock(),<font class=\\"jammer\\">% {, g) q& B' J( s1 X1 v- }2 Y) B</font><br /> s=0,<br /> <span style=\\"display:none\\">3 U' d4 b+ `# c0 i</span> k=0, while{k<1000,<font class=\\"jammer\\">8 u3 j\\" T9 g% z</font><br />   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)<br /> <span style=\\"display:none\\">. N9 r$ x' [2 _$ W$ ~9 n</span>   k++<font class=\\"jammer\\">! V( g) \5 i% a8 A/ O' _! g: L( O: G</font><br /> },<br /> <span style=\\"display:none\\">' f\\" d% X: S( j9 e0 B</span> s;<br /> <span style=\\"display:none\\">  N) t5 H4 q' J+ R% E) E\\" G; W</span> [clock()-t0]/1000;

结果：
2 c; Q5 ^1 g. d8 r9 m2000000.
40.766 秒
# `! D$ j# }  e" w
# Q( k4 Q' i: |! z2 ~1 ~- u程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
# j' ~. T* z7 X2 P因垃圾收集器多次启动，故耗时较长。

( Q) @; U) V% J" X& p－－－－－－－－－－－－

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

2. !using["math","sys"];& {. y- P3 Y5 @2 T
   
4. mvar:4 j% G7 i. E5 n, V+ [
   
6. main(:a)=. D. K8 ^- I3 b& Y' X% B+ Q; D
   
8. t0=clock(),
   
9.   j; s. c- H/ X) j, _) j
10. s=0,
    
11. / f0 r. X6 V% q: a
12. k=0, while{k<1000,# N8 F8 [0 R. ?6 Q7 c
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. ! {$ t- z, L7 j# g1 ~
16.   k++
    
17. ! @2 k/ |4 M7 n: U& T
18. },
    
19. ; T- d- u, u! `/ ~
20. s;
    
21. . I+ E6 ]- @8 G8 F1 h. G7 }
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">& {. y- P3 Y5 @2 T</font><br /> mvar:<font class=\\"jammer\\">4 j% G7 i. E5 n, V+ [</font><br /> main(:a)=<font class=\\"jammer\\">. D. K8 ^- I3 b& Y' X% B+ Q; D</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">  j; s. c- H/ X) j, _) j</span> s=0,<br /> <span style=\\"display:none\\">/ f0 r. X6 V% q: a</span> k=0, while{k<1000,<font class=\\"jammer\\"># N8 F8 [0 R. ?6 Q7 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\\">! {$ t- z, L7 j# g1 ~</span>   k++<br /> <span style=\\"display:none\\">! @2 k/ |4 M7 n: U& T</span> },<br /> <span style=\\"display:none\\">; T- d- u, u! `/ ~</span> s;<br /> <span style=\\"display:none\\">. I+ E6 ]- @8 G8 F1 h. G7 }</span> [clock()-t0]/1000;

结果：
2000000.
4.609 秒
' n* ~, Z  @8 O  k. I
可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
$ K8 m: |2 k! G在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
# d0 B* ~3 N. N$ h; \9 F- w; h8 Q太强大了 看不懂啊

+ \0 ]  _, `6 h. i' r函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forcal如何处理垃圾对象和进行内存管理。

我就在你背后

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

7 l7 k9 r% ^0 z. Y5 U- y恩 我在研究研究，谢谢了

sxjm567

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