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

forcal

这里用一个不太常用的小例子说明一下。
# d5 f- ?* @3 K- n* x7 {0 p- P
! [! Y/ k7 B, ~7 n% H/ `matlab 2009a代码：

1. clear all8 x5 P/ O( i, v/ N: @7 e4 ?. N5 o5 T
   
2. tic  M5 D; i\" ]1 C$ e
   
3. s=0;
   & R6 U3 h) u/ J
4. for k=1:1000
   8 M& W0 F) O. G; l3 _
5.   for i=1:10002 m6 }\" B  X# K# @1 x! R/ h6 L+ z
   
6.     a=linspace(2,2,1000);; J. l$ @' x' T, o  X
   
7.     s=s+a(i);, x; X  S& c, s. z# }1 B4 t5 Y
   
8.   end/ ]0 p! ?. ?) @, {
   
9. end
   . a0 l* Y1 {1 s% P' N( `+ H
10. s
    / ?! k8 U& k+ \) x
11. toc2 y  T* T1 H& q& h& O6 b
    
12. 6 x. Y+ r, u  a
    
13. s =
    , d3 v; W3 ?/ p9 x2 B/ m4 W
14. 
    ( y\" c\" w4 v! a, t
15.      2000000
    0 Q# ^4 a, G+ ]; v' J: Y9 A
16. 
    8 j. x- a! v# F/ x: {  X
17. Elapsed time is 18.325900 seconds.

复制代码

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

# [& N' T( g! I* C; n# M- }2 t8 C' E8 r－－－－－－－－

Forcal代码：

2. !using["math","sys"];
   
3. 6 P3 B. ~! c( ?- m# }6 A
4. mvar:; I; X6 O6 Y+ n7 Q
   
6. main(:a)=
   
7. 0 I# n: n4 S) Y/ B0 E
8. t0=clock(),! y  }/ ]- K% l4 f\\" V$ x. v1 p! I
   
10. s=0,
    
11. 4 U7 t: b! q4 M* Y7 f
12. k=0, while{k<1000,
    
13. ) V( F# q4 r6 Q4 M& a8 M
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. 1 y& K  z# Z' R* f\\" E
16.   k+++ A3 S. r0 k  b, d
    
18. },- I: r$ J0 {# t+ R$ c$ e- C! e- o
    
20. s;
    
21. : m8 Y% }* Q. d/ z
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">6 P3 B. ~! c( ?- m# }6 A</span> mvar:<font class=\\"jammer\\">; I; X6 O6 Y+ n7 Q</font><br /> main(:a)=<br /> <span style=\\"display:none\\">0 I# n: n4 S) Y/ B0 E</span> t0=clock(),<font class=\\"jammer\\">! y  }/ ]- K% l4 f\\" V$ x. v1 p! I</font><br /> s=0,<br /> <span style=\\"display:none\\">4 U7 t: b! q4 M* Y7 f</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">) V( F# q4 r6 Q4 M& a8 M</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\\">1 y& K  z# Z' R* f\\" E</span>   k++<font class=\\"jammer\\">+ A3 S. r0 k  b, d</font><br /> },<font class=\\"jammer\\">- I: r$ J0 {# t+ R$ c$ e- C! e- o</font><br /> s;<br /> <span style=\\"display:none\\">: m8 Y% }* Q. d/ z</span> [clock()-t0]/1000;

结果：
2000000.
4 h2 h# Z  C3 o6 o# Z5 c  {! J9 P40.766 秒
3 A0 F3 ^2 C- W, k; k3 J1 N
' m1 Y9 y% j% k; f程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
1 }/ q6 J3 o- {  b  R' t6 Y) \因垃圾收集器多次启动，故耗时较长。

/ a' n) n! {& T( z5 R－－－－－－－－－－－－

" U! w4 f5 B% R, Q, z7 T以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];
   
3. # I2 S, N1 F  M% }$ g# D
4. mvar:! M9 r9 ~. k/ J; G1 |% u\\" ?4 `% t7 D
   
6. main(:a)=
   
7. 0 k. l- `8 h; ]; c6 I
8. t0=clock(),
   
9.   b+ {! e+ Y; y& X
10. s=0,
    
11.   I/ B8 a5 D& i) Y
12. k=0, while{k<1000,
    
13. 0 N3 N, r# c! f& ^* s
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. % J9 z\\" e, h6 V
16.   k++
    
17. % s* G* C9 a9 l
18. },
    
19. ! v% }' v  Z9 Z( ]
20. s;
    
21. 9 q& p& C& w: s2 z( x
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\"># I2 S, N1 F  M% }$ g# D</span> mvar:<font class=\\"jammer\\">! M9 r9 ~. k/ J; G1 |% u\\" ?4 `% t7 D</font><br /> main(:a)=<br /> <span style=\\"display:none\\">0 k. l- `8 h; ]; c6 I</span> t0=clock(),<br /> <span style=\\"display:none\\">  b+ {! e+ Y; y& X</span> s=0,<br /> <span style=\\"display:none\\">  I/ B8 a5 D& i) Y</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">0 N3 N, r# c! f& ^* s</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\\">% J9 z\\" e, h6 V</span>   k++<br /> <span style=\\"display:none\\">% s* G* C9 a9 l</span> },<br /> <span style=\\"display:none\\">! v% }' v  Z9 Z( ]</span> s;<br /> <span style=\\"display:none\\">9 q& p& C& w: s2 z( x</span> [clock()-t0]/1000;

结果：
2000000.
4.609 秒
+ n, h1 S! V( u2 d1 z
* _6 S, D8 T, G! J* q1 A3 {可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
+ l7 K& O7 e, |7 T在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

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

. G' j! i+ d& U! R% \- _  [0 y+ D函数linspace(2,2,1000)放到循环里，每循环一次就会产生一个垃圾对象（一维数组），主要对比Matlab和Forcal如何处理垃圾对象和进行内存管理。

我就在你背后

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

/ C7 c; N4 N! {* {8 A- ?恩 我在研究研究，谢谢了

sxjm567

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