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

forcal

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

matlab 2009a代码：

1. clear all6 c+ b. Q7 U  F5 v) A* P
   
2. tic
   , F9 r+ }3 ~( P9 z
3. s=0;
     ?( W# E; ?1 x% ]3 h* H, Y
4. for k=1:1000
   * P5 _- {1 J* b$ J$ w% y* E
5.   for i=1:10007 d: o6 N9 h) C3 m
   
6.     a=linspace(2,2,1000);; \7 e8 k3 A' M+ }7 S9 T* b
   
7.     s=s+a(i);
   6 K9 C  O. \, S
8.   end- z2 H  M% d) u+ R& P7 Q
   
9. end# h. g- [: J4 Z4 ]5 L: D
   
10. s
    3 t& e. y! S- N+ d& I& C7 G
11. toc
    7 u+ s+ W! \\" ^4 a- Y7 M! J$ Z; E
12. 4 D; C) F1 F0 Z( I: A, Q4 [# F  ^
    
13. s =8 ^, d. b8 Q' D, [$ G1 v
    
14. . z2 L1 S: n0 p3 {. s4 J0 f! Y
    
15.      2000000/ }: y$ f, P, E; Q1 ?) [
    
16. 
    # m+ I/ A. v8 D8 K7 l, {- r
17. Elapsed time is 18.325900 seconds.

复制代码

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

* N; N# a3 a8 I* h, m－－－－－－－－
5 O0 T( p! P) @/ u$ T
5 h+ x$ m4 ?4 H$ d/ z2 G9 ZForcal代码：

2. !using["math","sys"];
   
3. . F8 B) Z& Z& x$ ^. v\\" [$ S
4. mvar:# {' N' V! e. }\\" V6 `
   
6. main(:a)=
   
7. / H% A& \; a9 }# d$ v: l! j
8. t0=clock(),
   
9. ( `/ n, M9 \$ ?) T; v6 ?. v
10. s=0,
    
11. ) W4 l( q4 f4 h\\" l# [8 u& V
12. k=0, while{k<1000,1 `; s! g, F+ @1 h& s6 M0 G: Z* V9 l
    
14.   i=0, while{i<1000, a=linspace(2,2,1000), s=s+a.i, i++},  // 嗯，a.i 就是 a(i)
    
15. 1 }. b\\" W/ E8 D
16.   k++3 w- \9 }3 O/ k7 \, T: J
    
18. },
    
19. + u. ~5 ~  c) F6 G6 E+ o* s7 a, _
20. s;) M; s8 X! S1 T/ U
    
22. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">. F8 B) Z& Z& x$ ^. v\\" [$ S</span> mvar:<font class=\\"jammer\\"># {' N' V! e. }\\" V6 `</font><br /> main(:a)=<br /> <span style=\\"display:none\\">/ H% A& \; a9 }# d$ v: l! j</span> t0=clock(),<br /> <span style=\\"display:none\\">( `/ n, M9 \$ ?) T; v6 ?. v</span> s=0,<br /> <span style=\\"display:none\\">) W4 l( q4 f4 h\\" l# [8 u& V</span> k=0, while{k<1000,<font class=\\"jammer\\">1 `; s! g, F+ @1 h& s6 M0 G: Z* V9 l</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\\">1 }. b\\" W/ E8 D</span>   k++<font class=\\"jammer\\">3 w- \9 }3 O/ k7 \, T: J</font><br /> },<br /> <span style=\\"display:none\\">+ u. ~5 ~  c) F6 G6 E+ o* s7 a, _</span> s;<font class=\\"jammer\\">) M; s8 X! S1 T/ U</font><br /> [clock()-t0]/1000;

结果：
& o2 c6 Q; T/ b1 J8 Z2000000.
) g/ B$ O, _" m  T, M40.766 秒
2 D5 ~9 R. j4 _. ~/ O, h
程序运行时内存不稳定，时大时小。该程序导致Forcal垃圾收集器启动约10次，垃圾收集前内存消耗大，垃圾收集后内存小。
因垃圾收集器多次启动，故耗时较长。
* A. o6 o* w' a6 |& c" f) `1 e
－－－－－－－－－－－－
# |: D" Z  ?, u9 j# J; A8 l
/ E2 d5 w; |' l8 }. \以下将a=linspace(2,2,1000)放到了oo{  }函数中：

2. !using["math","sys"];( u5 p. I; K; X6 h3 Q
   
4. mvar:
   
5. 0 J8 x! m1 s  U+ ?0 K8 O
6. main(:a)=
   
7. ( K' b2 t2 ^) i! C) f! f# ?
8. t0=clock(),5 {# k+ l: t. ~9 r& @: @$ B
   
10. s=0,
    
11. 5 o- P1 P- H% t& ~/ R
12. k=0, while{k<1000,+ X( ?9 u\\" ^* }$ n
    
14.   i=0, while{i<1000, oo{a=linspace(2,2,1000)}, s=s+a.i, i++},  // 嗯，a.i就是a(i)
    
15. % b\\" O- [4 N2 K8 l\\" J
16.   k++' T' r) s# y& B2 l1 l$ r9 `. `& S
    
18. },$ ^  y  ?9 ^\\" N' D0 L0 K
    
20. s;, M3 w% Z: ^* b# k\\" ?- D  }: L
    
22. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">( u5 p. I; K; X6 h3 Q</font><br /> mvar:<br /> <span style=\\"display:none\\">0 J8 x! m1 s  U+ ?0 K8 O</span> main(:a)=<br /> <span style=\\"display:none\\">( K' b2 t2 ^) i! C) f! f# ?</span> t0=clock(),<font class=\\"jammer\\">5 {# k+ l: t. ~9 r& @: @$ B</font><br /> s=0,<br /> <span style=\\"display:none\\">5 o- P1 P- H% t& ~/ R</span> k=0, while{k<1000,<font class=\\"jammer\\">+ X( ?9 u\\" ^* }$ n</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\\">% b\\" O- [4 N2 K8 l\\" J</span>   k++<font class=\\"jammer\\">' T' r) s# y& B2 l1 l$ r9 `. `& S</font><br /> },<font class=\\"jammer\\">$ ^  y  ?9 ^\\" N' D0 L0 K</font><br /> s;<font class=\\"jammer\\">, M3 w% Z: ^* b# k\\" ?- D  }: L</font><br /> [clock()-t0]/1000;

结果：
# \7 m. N9 R: T1 `) z: w% B2000000.
4.609 秒

! N% G* @+ C5 j* P可以发现，Forcal耗时大大减少了，远小于matlab耗时，这要归功于oo{ }函数。
: g( ]7 d5 @- g  A8 ?' o4 O在Forcal中，函数oo{ }不仅提供了运算符重载功能，而且会自动高效地回收垃圾。通常，所有涉及到对象的代码都要放到函数oo{ }中。

我就在你背后

太强大了 看不懂啊

forcal

我就在你背后 发表于 2011-8-2 16:08
2 q$ E, e+ f8 b! p6 ]& u太强大了 看不懂啊

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

我就在你背后

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

( x3 ?5 W1 w4 q8 F  D恩 我在研究研究，谢谢了

sxjm567

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