极限测试之Matlab与Forcal有无JIT加速器的比较

forcal

1、简单循环
: y3 f% s% h" I7 R' Q
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   ' A; y, c7 p! e
2. tic
   2 K: l  i2 H1 l/ Y. V6 c
3. s=0;
   / j  z. {- Q. o6 {( ?' V3 e5 k
4. for k=1:100005 S; w' \# {* {; s3 Q9 c' m
   
5.   for i=1:1000
   & }! g, ~, u' j7 E0 R; \
6.     s=s+1;7 g& j- ]% ^$ j
   
7.   end2 z) l( c0 i2 f% P  {# x
   
8. end3 t) K# b+ S$ T. T
   
9. s
   3 U' O' b: n) f& k2 j$ Z' j; C- \
10. toc
    \" P\" B: C$ b! ~1 |3 z
11. 
    + M8 [6 ~\" p5 Q
12. s =
    7 ^( \2 i% K! x( v. b( y9 K
13. 6 d7 n9 @4 X\" K9 G5 T( B
    
14.     10000000
    # h& m% X* d6 `* B0 R% o
15. + j1 j- ~% u\" C9 p7 H
    
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝

Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. : ~7 s3 ^. G4 {5 ]9 N
4. mvar:, K; a3 _, p$ B9 C\\" ^8 b4 C. ?
   
6. t0=clock(),0 m( V7 O* N7 [6 j/ Q0 a9 }/ W
   
8. s=0,6 E; s) G( Q9 _3 h5 S  _
   
10. k=0, while{k<10000,) m\\" b) g% q6 f) `, U8 _. n
    
12.   i=0, while{i<1000,s=s+1, i++},
    
13. 1 g4 H5 u9 z( M1 l
14.   k++
    
15. 4 K) _- K* P# M) D5 P) ^4 C
16. },
    
17. 5 U, S4 k& i( z1 G5 R
18. s;
    
19. ( B! i2 |/ s' \1 g
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">: ~7 s3 ^. G4 {5 ]9 N</span> mvar:<font class=\\"jammer\\">, K; a3 _, p$ B9 C\\" ^8 b4 C. ?</font><br /> t0=clock(),<font class=\\"jammer\\">0 m( V7 O* N7 [6 j/ Q0 a9 }/ W</font><br /> s=0,<font class=\\"jammer\\">6 E; s) G( Q9 _3 h5 S  _</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">) m\\" b) g% q6 f) `, U8 _. n</font><br />   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">1 g4 H5 u9 z( M1 l</span>   k++<br /> <span style=\\"display:none\\">4 K) _- K* P# M) D5 P) ^4 C</span> },<br /> <span style=\\"display:none\\">5 U, S4 k& i( z1 G5 R</span> s;<br /> <span style=\\"display:none\\">( B! i2 |/ s' \1 g</span> [clock()-t0]/1000;

结果：
3 u: p: w, ?5 v7 h# r6 D10000000.
0.704  秒
% O7 k; X7 y, [8 [0 N  C
, Q  k0 o7 J0 e3 u9 Q2 t# }5 b二者运行效率相差7倍多。
% s" v, n- Y/ J3 m2 K" C0 q
% s2 K$ M/ f$ E/ @9 Q7 ~2、循环中稍微加点东西
! i2 o/ |2 C3 `, [) B( T
# A; n  G( y' D* }; n( ]) _Matlab代码：

1. clear all
   ( x% e. _4 A! g3 d
2. tic
   8 V7 @7 B- W9 J! |, o
3. s=0;+ R9 a5 F$ u- d7 x& B\" y% w
   
4. for k=1:1000! X+ j' P) F9 d4 @& d
   
5.   for i=1:1000
   \" g& p) o  }! d; q2 R0 z9 A
6.     s=s+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i);
   + K8 k: ]$ Y0 ]2 j$ k
7.   end
   & `1 Y  v, i! o) `( U
8. end. O/ w( V6 i( t$ M
   
9. s
   * Q7 Y+ ~  w# s& ]' O) ?  W! F! f8 B
10. toc0 p  Z  L/ o# j+ h- C+ h
    
11. 
    7 c& g5 }9 p* M2 ]6 F) a- i
12. s =8 G$ I5 Q& ~) i3 Q
    
13. # a. w\" o# t3 h  e& s! h* l
    
14.   1.1009e+007
    & c; x  c+ w4 U5 `7 V- Z& M
15. 9 E9 }2 ^  R; x
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 5 p2 c* C# J, Z# ]$ p( k( V
4. mvar:1 S: S. \4 V3 z9 l) @0 u0 t9 Q$ q
   
6. t0=clock(),% Z0 e- @) ^; n- T; j1 V
   
8. s=0,  z6 @+ ~5 j% n5 ^
   
10. k=1, while{k<=1000,' @# ^\\" a3 @- n* p8 H
    
12.   i=1, while{i<=1000, s=s+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i), i++},& ?# h3 G$ b% \. c) y
    
14.   k++
    
15. . {7 ]' j% H# S8 S2 _- M9 x
16. },! v* D1 p8 k9 _$ S& ^, i
    
18. s;0 P* E) {  x' S9 X
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">5 p2 c* C# J, Z# ]$ p( k( V</span> mvar:<font class=\\"jammer\\">1 S: S. \4 V3 z9 l) @0 u0 t9 Q$ q</font><br /> t0=clock(),<font class=\\"jammer\\">% Z0 e- @) ^; n- T; j1 V</font><br /> s=0,<font class=\\"jammer\\">  z6 @+ ~5 j% n5 ^</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">' @# ^\\" a3 @- n* p8 H</font><br />   i=1, while{i<=1000, s=s+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i), i++},<font class=\\"jammer\\">& ?# h3 G$ b% \. c) y</font><br />   k++<br /> <span style=\\"display:none\\">. {7 ]' j% H# S8 S2 _- M9 x</span> },<font class=\\"jammer\\">! v* D1 p8 k9 _$ S& ^, i</font><br /> s;<font class=\\"jammer\\">0 P* E) {  x' S9 X</font><br /> [clock()-t0]/1000;

结果：
/ D! ~$ {! y% Y, l+ F11008953.66602346
1.031
! _  k  H& [! w" I! h
二者效率相差已比较小。
' R  U4 `, ^2 z3 p; \& r
* }! H, B% h4 H3 m4 O3、再增加一个函数调用

Matlab代码：

1. clear all
     U* p. `, t% G\" u, F$ Q- r
2. f=@(x)x+1;, p  D  Q& t4 |/ ?' i
   
3. tic
   0 m* z0 W( }/ C( c
4. s=0;7 O0 g9 `0 K' V, ^0 V2 P* y. K5 s
   
5. for k=1:1000* s7 v( R! v1 j% k
   
6.   for i=1:1000* l1 i0 G( D+ w, X1 B\" a2 w
   
7.     s=s+f(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i);# S7 [, B1 {$ `; e9 B+ T' _+ e
   
8.   end
   , X* Q' p\" }) Q  r; A. k
9. end
     ?* F: `, s5 X/ @7 H8 V# X
10. s/ n. S$ n! X  V+ Y( T+ w9 C
    
11. toc
    % V& {2 g0 C\" m6 _& b5 A4 R9 a
12. 
      P% p7 G8 G6 o) G; h2 ]! ^1 P4 o3 j
13. s =
      B% ~! N! |- |7 _6 C6 V+ f+ D+ U
14. 2 ~* i% _1 E  x* m2 I- N\" ?, R
    
15.   5.1251e+0089 |2 y) W% {3 K2 h
    
16. # t\" c  o, r: F# W
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 8 t3 c0 ?( f, F+ O: }: q' t
4. f(x)=x+1;\\" F) F6 N' r& J' t& a
   
6. mvar:7 P1 _: A: s; n+ u  C
   
8. t0=clock(),3 X% J% @/ y) B: E8 E/ t
   
10. s=0,2 H, \9 T4 h6 G$ J& I
    
12. k=1, while{k<=1000,( P9 O; f4 A4 ?1 u/ m: Y: J$ L4 v$ r
    
14.   i=1, while{i<=1000, s=s+f(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i), i++},
    
15. 4 {+ o  X! L9 l\\" z
16.   k++
    
17.   C  M1 X/ e; g
18. },
    
19. - ?- w, e0 M& F  s9 _# F' A
20. s;+ c& V( B* u  `9 q. c- d3 p& H
    
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">8 t3 c0 ?( f, F+ O: }: q' t</span> f(x)=x+1;<font class=\\"jammer\\">\\" F) F6 N' r& J' t& a</font><br /> mvar:<font class=\\"jammer\\">7 P1 _: A: s; n+ u  C</font><br /> t0=clock(),<font class=\\"jammer\\">3 X% J% @/ y) B: E8 E/ t</font><br /> s=0,<font class=\\"jammer\\">2 H, \9 T4 h6 G$ J& I</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">( P9 O; f4 A4 ?1 u/ m: Y: J$ L4 v$ r</font><br />   i=1, while{i<=1000, s=s+f(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i), i++},<br /> <span style=\\"display:none\\">4 {+ o  X! L9 l\\" z</span>   k++<br /> <span style=\\"display:none\\">  C  M1 X/ e; g</span> },<br /> <span style=\\"display:none\\">- ?- w, e0 M& F  s9 _# F' A</span> s;<font class=\\"jammer\\">+ c& V( B* u  `9 q. c- d3 p& H</font><br /> [clock()-t0]/1000;

结果：
, @% E0 ]5 n/ j' M! g9 o& K% s5 v9 [512508953.6658993
1.078
/ a8 W! t3 j; x, N9 q; R, C) a
可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
8 ^+ D% b% T0 Y0 L$ W; R* q+ k8 a3 x* {8 T
Matlab代码：

1. >> clear all. ]( w! X4 k: A. ]' @% b
   
2. a = linspace(1,1,1000);6 _$ n  K) C. M/ t. @
   
3. tic7 F/ R) F' c0 L' s( b4 e1 ~
   
4. for k=1:10000
   6 g  \7 m; d3 p. |
5.   for i=1:1000
   & ~. H% F. v) l: l: k2 M
6.     a(i)=2;
   # W# s8 E: h' `! ^' _# X1 Y7 y
7.   end- M9 H0 R' L$ C4 j  g, s; V
   
8. end' W. J+ c9 J; N
   
9. toc
   & p7 y# I* u  q- x* ]0 n
10. 
    , w& ?8 f$ e% l/ E% F$ z
11. tic
    ; V# i  w8 \6 r\" L: l- v3 m  l8 [. a
12. s=0;+ l6 Q7 r. y1 [% l1 X/ T# K
    
13. for k=1:10000
    , L\" ?2 B! i* i! F5 j* T
14.   for i=1:1000. R5 @( ^' ~! f1 }0 }
    
15.     s=s+a(i);! L3 m3 G' b- X) ~9 T
    
16.   end, U4 m7 I7 L5 ?, i1 ]) o. v
    
17. end
    / g- O* o6 i5 C6 s7 Z6 ]; ?- v\" Q
18. s
    \" I* Q. {% S+ F4 I3 e: E6 ^6 y
19. toc
    3 g! T4 W! `1 v: Z: {
20. Elapsed time is 0.140658 seconds.9 O/ V\" a; v% A& @4 \; g# q\" b
    
21. 
    : W8 d( _) h, y- l8 _! `
22. s =
      R5 L0 H. g  E3 a7 b6 U
23. $ ~$ r! Y6 l6 C
    
24.     20000000
    9 A2 p/ ~' ^\" U
25. 
    3 Q3 o' M/ y- Y4 ^! G6 F( I4 Z, \
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝

, {4 n% J3 Q4 H! \% WForcal代码：

2. !using["math","sys"];' [% `; G  w! A& ~
   
4. mvar:
   
5. $ w3 x# \9 J+ ]3 G' n& \* {5 p' |
6. main0(::a)=
   
7. % ?3 T+ n8 x+ F\\" {6 l
8. a = linspace(1,1,1000),8 {0 ]5 B: U6 u% m
   
10. t0=clock(),
    
11. - n2 o$ w, ?4 q* I0 [
12. k=0, while{k<10000,
    
13. 2 h+ T- m* t* g- i2 |1 A1 r9 _
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. + p5 x8 j: m+ y\\" h2 l+ p
16.   k++
    
17. & G; A5 _! e3 I3 F
18. },
    
19. / c% L7 \8 c* \# P* v+ a3 L! ~) `
20. [clock()-t0]/1000;% n( t1 F& l- u
    
22. 3 q9 c& {/ B# g% X3 b
    
24. main1(::a)=
    
25. , s- N' ^\\" w3 c
26. t0=clock(), s=0,
    
27. % g/ g* r; I' u7 J
28. k=0, while{k<10000,
    
29. * t2 v4 x* K( D) o. g
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. 7 w2 b0 l5 L+ e* c$ b( b
32.   k++% L3 T. B/ j: s( Z  |+ a7 D; h% j1 h
    
34. },
    
35. 1 _; J\\" q2 h4 k\\" a8 s. d  A* E
36. s;9 s2 {- w/ s- l# G8 @% {
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">' [% `; G  w! A& ~</font><br /> mvar:<br /> <span style=\\"display:none\\">$ w3 x# \9 J+ ]3 G' n& \* {5 p' |</span> main0(::a)=<br /> <span style=\\"display:none\\">% ?3 T+ n8 x+ F\\" {6 l</span> a = linspace(1,1,1000),<font class=\\"jammer\\">8 {0 ]5 B: U6 u% m</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">- n2 o$ w, ?4 q* I0 [</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">2 h+ T- m* t* g- i2 |1 A1 r9 _</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">+ p5 x8 j: m+ y\\" h2 l+ p</span>   k++<br /> <span style=\\"display:none\\">& G; A5 _! e3 I3 F</span> },<br /> <span style=\\"display:none\\">/ c% L7 \8 c* \# P* v+ a3 L! ~) `</span> [clock()-t0]/1000;<font class=\\"jammer\\">% n( t1 F& l- u</font><br /> <font class=\\"jammer\\">3 q9 c& {/ B# g% X3 b</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">, s- N' ^\\" w3 c</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">% g/ g* r; I' u7 J</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">* t2 v4 x* K( D) o. g</span>   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">7 w2 b0 l5 L+ e* c$ b( b</span>   k++<font class=\\"jammer\\">% L3 T. B/ j: s( Z  |+ a7 D; h% j1 h</font><br /> },<br /> <span style=\\"display:none\\">1 _; J\\" q2 h4 k\\" a8 s. d  A* E</span> s;<font class=\\"jammer\\">9 s2 {- w/ s- l# G8 @% {</font><br /> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒

－－－－－－－－－－
! v- x- K, Z* ?% z+ B  S* G/ H
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3.   _- c  o  G% d
4. mvar:8 K0 D* v1 p: }4 _/ C
   
6. main0(::a)=
   
7. 6 ?( x) I3 j$ Y
8. a = linspace(1,1,1000),
   
9. - z* H! j$ J2 F/ t8 D: H
10. t0=clock(),0 p3 p$ V% L9 R4 I! h+ p
    
12. k=0, while{k<10000,
    
13. , A- T, m* o- A9 a4 d
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. 6 }/ {9 H* a1 h( _3 E$ L% R
16.   k++( L. M, r& ?/ }! i( I
    
18. },2 p( h\\" A  u! j. M4 G+ h
    
20. [clock()-t0]/1000;$ J/ x% s\\" X3 I
    
22. 
    
23. 2 W6 V. S; k- I3 S7 \
24. main1(::a)=
    
25. 9 l( F, n. m- p
26. t0=clock(), s=0,. |) H; v: B% ~( [
    
28. k=0, while{k<10000,
    
29. 6 {. H+ k  C8 Q( A2 h1 I& Q; u5 r' x5 S
30.   i=0, while{i<1000, s=s+A(a,i), i++},1 Q- S$ w% O8 ~, E7 `
    
32.   k++
    
33. 0 k' g( ^  g6 {1 }1 k) j, U& s
34. },
    
35. ( o\\" D8 B# k. `% M
36. s;% R. ^& g0 g: {: l0 B* W9 t
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">  _- c  o  G% d</span> mvar:<font class=\\"jammer\\">8 K0 D* v1 p: }4 _/ C</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">6 ?( x) I3 j$ Y</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">- z* H! j$ J2 F/ t8 D: H</span> t0=clock(),<font class=\\"jammer\\">0 p3 p$ V% L9 R4 I! h+ p</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">, A- T, m* o- A9 a4 d</span>   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">6 }/ {9 H* a1 h( _3 E$ L% R</span>   k++<font class=\\"jammer\\">( L. M, r& ?/ }! i( I</font><br /> },<font class=\\"jammer\\">2 p( h\\" A  u! j. M4 G+ h</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">$ J/ x% s\\" X3 I</font><br /> <br /> <span style=\\"display:none\\">2 W6 V. S; k- I3 S7 \</span> main1(::a)=<br /> <span style=\\"display:none\\">9 l( F, n. m- p</span> t0=clock(), s=0,<font class=\\"jammer\\">. |) H; v: B% ~( [</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">6 {. H+ k  C8 Q( A2 h1 I& Q; u5 r' x5 S</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">1 Q- S$ w% O8 ~, E7 `</font><br />   k++<br /> <span style=\\"display:none\\">0 k' g( ^  g6 {1 }1 k) j, U& s</span> },<br /> <span style=\\"display:none\\">( o\\" D8 B# k. `% M</span> s;<font class=\\"jammer\\">% R. ^& g0 g: {: l0 B* W9 t</font><br /> [clock()-t0]/1000;

结果：
10.469
! b% W, \* K- @: q5 C/ \" q20000000.
, K4 T  Z) p: ^8 x7 ^: _; C3.781
; ^! |; C$ V) m
+ V$ f+ q" C$ x0 m－－－－－－－

0 \, J1 t  ]! t! G  @% Z在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
1 S% N$ w* Q% G! c
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
, g5 |1 [& y, o+ j# x; j
" g4 k) }, _  H+ tForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
4 y* ^0 I  h+ Z3 ]1 s
5、在循环中加点东西

matlab代码：

1. >> clear all2 C% u4 i0 v9 M
   
2. f=@(x)x+1;8 O; p\" Q2 A8 ~8 W& v/ g1 J; p1 P
   
3. a = linspace(2,2,1000);
   % F/ R2 M\" p' o
4. tic/ a1 E6 A: r7 Z2 Z( `4 R+ m
   
5. s=0;
   ( U6 q\" j/ @8 b
6. for k=1:1000
   . S8 f' L* [\" I( v# ~/ @1 @0 B) V
7.   for i=1:1000
   7 H0 _2 N9 [( J7 [. `\" s6 I5 t
8.     s=s+f(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i));
   % @, Q\" N! ^2 Q\" W- G& r* H. b
9.   end; x2 z: W; h\" u! u7 H  r0 F1 _\" ]* q
   
10. end& k- ]* Q: S2 F) |) X' E
    
11. s+ a4 d  s. x. H
    
12. toc5 c2 k; \% i9 e* {  q, b5 @$ S
    
13. 
    # u( C5 D' X( k- F. A
14. s =& m. C% {: e+ ^7 ^. c8 |* m
    
15. 
    + n' w9 ?9 q) [' t\" o# o- R% ]
16.   2.4002e+007
    7 T4 R4 K) D9 R9 L8 c, D4 O6 P& Q
17. 2 H/ K9 R2 K4 a9 R9 U* Z
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];# _/ I) I7 ^9 h8 o! j
   
4. mvar:4 P3 [! M- k3 i8 N% `
   
6. f(x)=x+1;7 e: Q; f# B$ x, S& \
   
8. main0(::a)=) l  _% C/ I) s# w+ P- a
   
10. a = linspace(2,2,1000),' Z9 t% ^1 B5 E1 }5 h1 B; o\\" i
    
12. t0=clock(), s=0,
    
13. $ J/ [3 \$ `* z9 P( o8 s
14. k=0, while{k<1000,
    
15. \\" e5 ?) n0 D; k% l  V1 b
16.   i=0, while{i<1000, s=s+f(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i)), i++},
    
17. + \) Z\\" \5 z8 x. {4 T7 Z1 D
18.   k++# U0 @' {1 p3 i. {$ F* U7 r
    
20. },\\" H! T# D; n  m
    
22. s;
    
23. + Q! h; P- H# ?8 w+ _; i5 P( L( ]6 N
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\"># _/ I) I7 ^9 h8 o! j</font><br /> mvar:<font class=\\"jammer\\">4 P3 [! M- k3 i8 N% `</font><br /> f(x)=x+1;<font class=\\"jammer\\">7 e: Q; f# B$ x, S& \</font><br /> main0(::a)=<font class=\\"jammer\\">) l  _% C/ I) s# w+ P- a</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">' Z9 t% ^1 B5 E1 }5 h1 B; o\\" i</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">$ J/ [3 \$ `* z9 P( o8 s</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">\\" e5 ?) n0 D; k% l  V1 b</span>   i=0, while{i<1000, s=s+f(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i)), i++},<br /> <span style=\\"display:none\\">+ \) Z\\" \5 z8 x. {4 T7 Z1 D</span>   k++<font class=\\"jammer\\"># U0 @' {1 p3 i. {$ F* U7 r</font><br /> },<font class=\\"jammer\\">\\" H! T# D; n  m</font><br /> s;<br /> <span style=\\"display:none\\">+ Q! h; P- H# ?8 w+ _; i5 P( L( ]6 N</span> [clock()-t0]/1000;

结果：
24002271.69189384
$ }8 f' T6 Q1 r5.172
; c. Q$ H4 \+ [4 f
9 T5 W. I# A( g, R) p二者效率现在只有一倍差距。
, v" d, R- [& Q! h/ x
; h4 {  V( q4 x5 j& X- Y# Q－－－－－－－－－－
6 r! L9 ~) ]4 `5 q2 ?$ i4 ?
/ r$ c/ D: w8 i8 K' Z/ a# B  m  ?: X不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

1. clear all
   8 E) y0 r+ ^% u* z* Q: `$ O% k
2. f=@(x)x+1;6 L& F# O% i6 L2 e
   
3. a = linspace(2,2,1000);4 F* ~- P  w0 K) i
   
4. tic
   - a# |% Y7 A6 m1 j) h\" j
5. s=0;
   * ~9 h\" u' x3 K8 \/ r! x
6. for k=1:1000
   9 `8 r3 [8 P3 n0 U& F
7.   for i=1:1000
   5 O\" _* Y' F) c0 m$ ?+ ]3 t0 Z
8.     t=a(i);- @  p6 B/ W3 Z, f; }  n
   
9.     s=s+f(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t);
   ! [& ]8 s+ H( i
10.   end
    + S# D0 K; W\" k7 n- l6 {
11. end
    \" J2 E8 w9 a) s; Z) C+ o
12. s( u! L' N. o. \$ b1 |
    
13. toc
    + x; D7 `, B; e  C% E, x- s
14. 
    5 F, K\" S% W* _9 c( R& o
15. s =\" R% e+ g* l1 I6 \
    
16. 0 r4 d3 ~( ]2 {% c5 t7 j2 I9 L
    
17.   2.4002e+0079 R3 D  x/ j0 f0 k, |. V
    
18. 
    1 Y  f1 s/ l1 x
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3.   r3 X, r: C2 f& B  V
4. mvar:
   
5. # }! d( q' X. M
6. f(x)=x+1;
   
7. $ ~- U' C2 b! O, @/ [
8. main0(::a)=
   
9. , n! }$ d% R2 ]  m# x
10. a = linspace(2,2,1000),  D9 t7 o- x8 u! B\\" w4 X3 Q
    
12. t0=clock(), s=0,. h0 Y: F: U; H
    
14. k=0, while{k<1000,- k* q\\" t0 M. O) B# L
    
16.   i=0, while{i<1000, t=a(i), s=s+f(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t), i++},. H3 }* J( e  I( E
    
18.   k++
    
19. * U& w, {& H1 t2 Z* Y
20. },; F& \, {: v3 t4 D' @
    
22. s;6 M* I5 n' L% @\\" o8 m9 p8 a
    
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">  r3 X, r: C2 f& B  V</span> mvar:<br /> <span style=\\"display:none\\"># }! d( q' X. M</span> f(x)=x+1;<br /> <span style=\\"display:none\\">$ ~- U' C2 b! O, @/ [</span> main0(::a)=<br /> <span style=\\"display:none\\">, n! }$ d% R2 ]  m# x</span> a = linspace(2,2,1000),<font class=\\"jammer\\">  D9 t7 o- x8 u! B\\" w4 X3 Q</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">. h0 Y: F: U; H</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">- k* q\\" t0 M. O) B# L</font><br />   i=0, while{i<1000, t=a(i), s=s+f(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t), i++},<font class=\\"jammer\\">. H3 }* J( e  I( E</font><br />   k++<br /> <span style=\\"display:none\\">* U& w, {& H1 t2 Z* Y</span> },<font class=\\"jammer\\">; F& \, {: v3 t4 D' @</font><br /> s;<font class=\\"jammer\\">6 M* I5 n' L% @\\" o8 m9 p8 a</font><br /> [clock()-t0]/1000;

结果：
3 P$ o! b0 T$ Z* F- Q. ]24002271.69189384
2 R- R/ b- T9 S8 H  a8 {; S& L1.625
3 J/ Q3 U6 Q4 X) T# A: ~
& Z' }  ?" d; x  T/ Ematlab效率已不及Forcal。

提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
- C0 R) d% O  R& @; ^5 v2 v6 U
1 A6 t$ T3 m+ y$ Q1 |/ ^4 B& qForcal代码：

2. !using["sys"];* h2 ^# O& a% q: {0 {5 p' @
   
4. mvar:
   
5. 1 v6 I! [& f) ?: `. N8 m+ X\\" ?
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. 2 n+ J$ k) n; b# s9 p9 Z8 z\\" Y
8. t0=clock(),4 w3 @& C2 B8 S
   
10. k=0, while{k<10000,
    
11. - R9 ?2 p# x2 T. B
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值: C; `9 L8 f, A& Y2 j! G
    
14.   k+++ {* l# V0 u, c3 u3 ?7 u
    
16. },# ~, o7 G0 B6 O\\" |& R0 ~\\" P- i$ O
    
18. [clock()-t0]/1000;2 @* \7 f; k$ p/ D8 o! T
    
20. 
    
21. + c/ L% S4 ?+ e; g& f1 R5 L
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. ! H6 a  R: _% r, r4 M
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值3 @/ O2 {' b6 y) \& `- x
    
26. t0=clock(),
    
27. 6 k8 K' y3 t; H' C. G\\" {
28. s=0,& H* O5 w! r7 _
    
30. k=0, while{k<10000,9 B; t# D) s2 I: ~3 I
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. \\" m! P# B$ y1 t6 q  n& C, O8 l
34.   k++
    
35. \\" w' k\\" B- m6 c  o! y1 O
36. },1 c1 i% M# N. U1 b% i$ h
    
38. s;9 `\\" a9 N! U+ R. m* E9 l+ x
    
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">* h2 ^# O& a% q: {0 {5 p' @</font><br /> mvar:<br /> <span style=\\"display:none\\">1 v6 I! [& f) ?: `. N8 m+ X\\" ?</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">2 n+ J$ k) n; b# s9 p9 Z8 z\\" Y</span> t0=clock(),<font class=\\"jammer\\">4 w3 @& C2 B8 S</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">- R9 ?2 p# x2 T. B</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">: C; `9 L8 f, A& Y2 j! G</font><br />   k++<font class=\\"jammer\\">+ {* l# V0 u, c3 u3 ?7 u</font><br /> },<font class=\\"jammer\\"># ~, o7 G0 B6 O\\" |& R0 ~\\" P- i$ O</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">2 @* \7 f; k$ p/ D8 o! T</font><br /> <br /> <span style=\\"display:none\\">+ c/ L% S4 ?+ e; g& f1 R5 L</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">! H6 a  R: _% r, r4 M</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">3 @/ O2 {' b6 y) \& `- x</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">6 k8 K' y3 t; H' C. G\\" {</span> s=0,<font class=\\"jammer\\">& H* O5 w! r7 _</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">9 B; t# D) s2 I: ~3 I</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">\\" m! P# B$ y1 t6 q  n& C, O8 l</span>   k++<br /> <span style=\\"display:none\\">\\" w' k\\" B- m6 c  o! y1 O</span> },<font class=\\"jammer\\">1 c1 i% M# N. U1 b% i$ h</font><br /> s;<font class=\\"jammer\\">9 `\\" a9 N! U+ R. m* E9 l+ x</font><br /> [clock()-t0]/1000;

结果：
( X0 W. m) [: U0.5620000000000001  //赋值用时间，秒
! S5 P* t( @/ G4 \* P20000000.
1 C% i) }9 h3 B: Z7 {7 v0.579      //取值用时间，秒