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

forcal

1、简单循环
) [( J; @6 f# P! `* Q
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   4 Y, y9 i! b9 R
2. tic
   4 b' G  u6 b* A0 S; c
3. s=0;
   3 R, [% C+ l& p+ U- E
4. for k=1:10000' q1 c) Y! }9 l1 D. _5 t( ?
   
5.   for i=1:1000
   5 A* o7 [  I7 z, U
6.     s=s+1;\" M7 t. N; B% U# ^' H
   
7.   end
   0 I2 u8 D+ k6 o' P/ ~
8. end' W8 c; @/ H# I( u5 H( k6 ^
   
9. s+ n; I7 r6 |4 L; T) [, N! O
   
10. toc\" t; p1 i& q1 l; ~6 W* W( Z
    
11. 9 |8 |  ?  g, l) w4 {
    
12. s =  d% o9 y/ Q* ^$ E5 w  ^7 E
    
13. 
    ! ?# O' ~) i; ~( S5 M, Q
14.     10000000# Q5 o% ~! }  L
    
15. 
    . H4 M( @2 E& x: G* x  l
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
  V3 O* i4 w5 m% |
- I* [( [5 L0 b2 M. V0 Z' IForcal没有JIT加速器，代码：

2. !using["sys"];
   
3. * x+ ?, A$ B* D1 {& n7 x* F
4. mvar:! R\\" N; f4 J/ o. r+ Q
   
6. t0=clock(),
   
7. & y! k* v6 A8 h+ ~0 S) |
8. s=0,
   
9. - G  ]\\" J& \1 M4 H' z) n8 G0 m
10. k=0, while{k<10000,8 e) _2 z7 R2 ^\\" M- h, A/ a
    
12.   i=0, while{i<1000,s=s+1, i++},
    
13. ( d5 W  n5 w* S3 ]; W! o) X. `
14.   k++* W( C8 T& f0 i% z1 ]
    
16. },: L9 X& S5 x5 {) G( P( O
    
18. s;! A- B; j5 M1 C: k3 {
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">* x+ ?, A$ B* D1 {& n7 x* F</span> mvar:<font class=\\"jammer\\">! R\\" N; f4 J/ o. r+ Q</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">& y! k* v6 A8 h+ ~0 S) |</span> s=0,<br /> <span style=\\"display:none\\">- G  ]\\" J& \1 M4 H' z) n8 G0 m</span> k=0, while{k<10000,<font class=\\"jammer\\">8 e) _2 z7 R2 ^\\" M- h, A/ a</font><br />   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">( d5 W  n5 w* S3 ]; W! o) X. `</span>   k++<font class=\\"jammer\\">* W( C8 T& f0 i% z1 ]</font><br /> },<font class=\\"jammer\\">: L9 X& S5 x5 {) G( P( O</font><br /> s;<font class=\\"jammer\\">! A- B; j5 M1 C: k3 {</font><br /> [clock()-t0]/1000;

结果：
0 m0 _* _) e8 K10000000.
0.704  秒

" e: ?" ^9 t7 }% \# Z& h5 L二者运行效率相差7倍多。
$ r* F: q; Y& w  o8 ]
2、循环中稍微加点东西

6 d0 e2 i. t4 o( KMatlab代码：

1. clear all
   ) t; \0 |9 u' W* |8 r# Q5 M% I  D  r& X
2. tic
   8 e& H3 L& o- x8 P1 h. {
3. s=0;
   . G* e5 c4 c. S
4. for k=1:1000
   ( y* u. z\" N, w' Q$ r
5.   for i=1:1000/ `5 T0 Q( B\" v  K) r$ f# w  w
   
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);9 r) m- N1 c5 D# W
   
7.   end
   4 e. G: g) N2 r  [
8. end
   . C9 Q9 r- C! S6 T4 d' Q
9. s4 i  J* |) @* l% w$ z$ @; R
   
10. toc4 j' K7 E; Z0 G, V
    
11. 
      g! s5 x2 X) O( D9 u2 a
12. s =
    3 v( A  `& P  `8 M( @
13. 
    6 p9 t$ W0 U/ O
14.   1.1009e+007
      A1 k4 h5 ?% Q5 I7 Y. |
15. % Y5 o5 A& |' m. O% z4 \
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];: P) F& ]6 n. I
   
4. mvar:6 a* L, G9 ]2 s3 w
   
6. t0=clock(),
   
7. 7 y& q& P8 X0 t; Z! S8 W4 g, F
8. s=0,9 W; t  `, l: ]) v1 Q9 V9 `
   
10. k=1, while{k<=1000,- y7 ]  ^. e* a\\" t1 ^
    
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++},
    
13. : W8 w\\" R4 @$ P7 Y$ [$ v7 q8 _) ~# m1 M
14.   k++
    
15. 9 i; V8 G2 @1 |: j5 |
16. },
    
17. ' P+ t( Y& h; y1 l# c# E* L
18. s;; L: g7 q/ z8 Q) H3 y& c2 j% s' V
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">: P) F& ]6 n. I</font><br /> mvar:<font class=\\"jammer\\">6 a* L, G9 ]2 s3 w</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">7 y& q& P8 X0 t; Z! S8 W4 g, F</span> s=0,<font class=\\"jammer\\">9 W; t  `, l: ]) v1 Q9 V9 `</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">- y7 ]  ^. e* a\\" t1 ^</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++},<br /> <span style=\\"display:none\\">: W8 w\\" R4 @$ P7 Y$ [$ v7 q8 _) ~# m1 M</span>   k++<br /> <span style=\\"display:none\\">9 i; V8 G2 @1 |: j5 |</span> },<br /> <span style=\\"display:none\\">' P+ t( Y& h; y1 l# c# E* L</span> s;<font class=\\"jammer\\">; L: g7 q/ z8 Q) H3 y& c2 j% s' V</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
1.031
) b  o, l$ g  Z5 a  V% `
5 |; k' Q0 ^6 ]: p4 G& O二者效率相差已比较小。

3、再增加一个函数调用
7 d. G* t4 J* i( c9 G9 w6 j! F
Matlab代码：

1. clear all
   5 t* s: S( [% ~+ t* W
2. f=@(x)x+1;
   2 x. u8 @3 |/ A! q6 `5 [
3. tic- N& k5 S; f2 H7 s1 t
   
4. s=0;
   1 c- a* q* M7 c6 |* J. {8 C0 _7 n
5. for k=1:1000
   # z8 R7 q6 q3 Y; g& u6 p
6.   for i=1:1000' e' t* I8 G/ W& P% D# T% C
   
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);( k\" V( ?! c5 n) r8 K) B
   
8.   end
   , ]2 _3 C- S  ~9 V
9. end
   * v6 Y; u& u. o3 i: u
10. s6 [; s+ K9 s2 q4 B/ P0 {- @
    
11. toc
    5 E3 A\" S; o8 w1 e. \( l% ^' f
12. : P# N' j' }& \  ~$ V7 Z
    
13. s =5 D\" s9 _: h- s3 @, r0 R
    
14. 
    2 F2 k& D! U# Z& {1 [
15.   5.1251e+008! T5 q' t) H) v/ A1 u- u
    
16. 4 a: ^\" a: r+ K8 z9 {
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. ) L, q! b! ?) a( j
4. f(x)=x+1;
   
5. # s7 W+ p8 z# G% W1 a/ O, |
6. mvar:, i4 b& Q5 d1 ]
   
8. t0=clock(),8 }; Z  B, X* }
   
10. s=0,
    
11. $ {: i! a( H5 P; K( r) a' ~4 S0 e8 E
12. k=1, while{k<=1000,
    
13. ' w7 @' f% ?; p- B) v
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. ! S: b! s3 B: i' x\\" n6 {
16.   k++* I' G; Q/ _7 o
    
18. },
    
19. 5 q' g# q0 l5 S7 D\\" T
20. s;
    
21. 4 J0 L3 M3 N# d) S2 O9 [
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">) L, q! b! ?) a( j</span> f(x)=x+1;<br /> <span style=\\"display:none\\"># s7 W+ p8 z# G% W1 a/ O, |</span> mvar:<font class=\\"jammer\\">, i4 b& Q5 d1 ]</font><br /> t0=clock(),<font class=\\"jammer\\">8 }; Z  B, X* }</font><br /> s=0,<br /> <span style=\\"display:none\\">$ {: i! a( H5 P; K( r) a' ~4 S0 e8 E</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">' w7 @' f% ?; p- B) v</span>   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\\">! S: b! s3 B: i' x\\" n6 {</span>   k++<font class=\\"jammer\\">* I' G; Q/ _7 o</font><br /> },<br /> <span style=\\"display:none\\">5 q' g# q0 l5 S7 D\\" T</span> s;<br /> <span style=\\"display:none\\">4 J0 L3 M3 N# d) S2 O9 [</span> [clock()-t0]/1000;

结果：
' ?' K1 [2 K1 w% h$ M+ Z% C512508953.6658993
1.078
4 f) U$ i, ?9 S
8 S0 y  g1 r1 [9 O+ r+ r可以看出，Matlab效率已不及Forcal。
9 \- M" |5 z% J2 }故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高

Matlab代码：

1. >> clear all
   . Y6 K! U% K) r( j. Q  i- I
2. a = linspace(1,1,1000);, P' g0 d% X- r
   
3. tic
   / Y3 \! n; t8 y4 ?: e+ W
4. for k=1:10000
   % {5 q4 w  m; V& V! H, |
5.   for i=1:1000
   ) S& r1 K0 ~) v6 F' j3 k; n- `0 c\" R
6.     a(i)=2;
   & n6 I8 a1 ^1 C- D
7.   end
   3 u% c+ T0 d; k8 z$ u- A' T, Y
8. end
   0 y! v9 }* L9 q2 v: F, @
9. toc
   ( W! r* e' k9 {  G. s0 U/ z, L8 f
10. 
    \" U1 f4 a7 H% L- ~8 Q2 [2 |
11. tic6 B; z* p  o8 \5 t/ W! k) ?+ ^  A1 o) U
    
12. s=0;- p3 Y7 E# c  B  s
    
13. for k=1:10000\" `/ n- h$ _: I
    
14.   for i=1:1000
    3 E0 ~& m! V8 t8 j
15.     s=s+a(i);
    % {, O* r1 k5 g2 S# I$ @% _
16.   end) A; d( \; y9 n, u  ?: c( f; Q
    
17. end
    $ R\" I  T$ P7 J
18. s5 l& u) A  E7 @* c8 O% X! g5 ]
    
19. toc
    3 e. N\" h& c8 c' I
20. Elapsed time is 0.140658 seconds.+ N! J5 a2 J1 L% W  }- W* e! V* d* ]7 q
    
21. 
    7 q5 _/ x6 Z, e2 N
22. s =
    ! `: ]' H4 L9 U7 T
23. 
    5 g4 ]  |% A\" \# e+ o/ X5 w$ j
24.     200000006 H4 P6 s- t) N0 G* ]2 B
    
25. 3 g& g5 s( i  b  |. r$ a) F
    
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
+ J. [$ G* F1 N: b
4 Y& }/ B9 d; F% k- w% t: {Forcal代码：

2. !using["math","sys"];  E4 U+ L2 ]( a  R3 ]+ M
   
4. mvar:
   
5. # a\\" w, c/ @4 O. a( Z! @3 }; F
6. main0(::a)=
   
7. 7 j5 \* d# l! X2 K- ^  N
8. a = linspace(1,1,1000),) K/ p7 p5 V; U
   
10. t0=clock(),. n7 x, ?$ E5 j) I
    
12. k=0, while{k<10000,
    
13. \\" e  \  K% |) k( p- M7 k8 {+ T
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. , H( E6 x; |2 D9 M. j7 _& G
16.   k++5 N7 s8 e0 K\\" O' S
    
18. },; t0 W* j7 T. Y! x  p% R! d
    
20. [clock()-t0]/1000;
    
21. ( I( \9 K' y! F7 h1 R, G
22. 3 u$ {. Z' V5 V+ k6 I
    
24. main1(::a)=
    
25. + t# Y7 b; N# c4 ^* C6 J# p4 `
26. t0=clock(), s=0,
    
27. & P, ?( E. W7 G
28. k=0, while{k<10000,, G# R9 E9 t4 _$ e
    
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. , }/ B* S6 L8 e# p6 i% R
32.   k++. B6 \( {7 j; f# E$ o
    
34. },
    
35. % A4 X+ ~0 \% P/ q
36. s;8 F- I; g  k, H5 T  h\\" G. _
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">  E4 U+ L2 ]( a  R3 ]+ M</font><br /> mvar:<br /> <span style=\\"display:none\\"># a\\" w, c/ @4 O. a( Z! @3 }; F</span> main0(::a)=<br /> <span style=\\"display:none\\">7 j5 \* d# l! X2 K- ^  N</span> a = linspace(1,1,1000),<font class=\\"jammer\\">) K/ p7 p5 V; U</font><br /> t0=clock(),<font class=\\"jammer\\">. n7 x, ?$ E5 j) I</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">\\" e  \  K% |) k( p- M7 k8 {+ T</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">, H( E6 x; |2 D9 M. j7 _& G</span>   k++<font class=\\"jammer\\">5 N7 s8 e0 K\\" O' S</font><br /> },<font class=\\"jammer\\">; t0 W* j7 T. Y! x  p% R! d</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">( I( \9 K' y! F7 h1 R, G</span> <font class=\\"jammer\\">3 u$ {. Z' V5 V+ k6 I</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">+ t# Y7 b; N# c4 ^* C6 J# p4 `</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">& P, ?( E. W7 G</span> k=0, while{k<10000,<font class=\\"jammer\\">, G# R9 E9 t4 _$ e</font><br />   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">, }/ B* S6 L8 e# p6 i% R</span>   k++<font class=\\"jammer\\">. B6 \( {7 j; f# E$ o</font><br /> },<br /> <span style=\\"display:none\\">% A4 X+ ~0 \% P/ q</span> s;<font class=\\"jammer\\">8 F- I; g  k, H5 T  h\\" G. _</font><br /> [clock()-t0]/1000;

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

6 V8 k: C" p- L4 G－－－－－－－－－－
9 F; N0 y: G, K( Z* ?9 V
0 T+ i2 e: q& ?5 q9 y2 l用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. % V' P! o6 k' G# \
4. mvar:
   
5. 3 Z\\" ^3 n4 R  Q' s
6. main0(::a)=
   
7. : T* @; e5 \' p3 p$ A8 i# ~/ {
8. a = linspace(1,1,1000),( b; M! B$ U# B0 P5 K
   
10. t0=clock(),
    
11. . p5 `3 L. A$ ?+ N  K- j6 ]
12. k=0, while{k<10000,
    
13. : L5 A. R0 L3 `9 i. y! _4 N* K
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. - J/ `% o( E2 b. k
16.   k++% G5 |1 ~  S$ [6 x+ A3 q5 Y
    
18. },
    
19. ( v1 c\\" X' O# _/ w' b+ e
20. [clock()-t0]/1000;
    
21. 1 N9 N  z, m2 Y2 _* f4 N
22. 
    
23. ) i9 H+ g2 L/ T: G0 d) u7 U3 e' \( q  G
24. main1(::a)=
    
25. ( x% u( E, P. k7 ?
26. t0=clock(), s=0,  G8 t+ l% s% v5 @) j: L
    
28. k=0, while{k<10000,7 C6 B' V: \! p\\" W% C2 k/ r6 d
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. 7 F8 f& U) ~) v4 _  ^5 F
32.   k++
    
33. : K% i3 x3 u9 l# N6 s7 v
34. },7 z4 V! b3 v+ C
    
36. s;' T6 i\\" H- b% v6 k  `, \
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">% V' P! o6 k' G# \</span> mvar:<br /> <span style=\\"display:none\\">3 Z\\" ^3 n4 R  Q' s</span> main0(::a)=<br /> <span style=\\"display:none\\">: T* @; e5 \' p3 p$ A8 i# ~/ {</span> a = linspace(1,1,1000),<font class=\\"jammer\\">( b; M! B$ U# B0 P5 K</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">. p5 `3 L. A$ ?+ N  K- j6 ]</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">: L5 A. R0 L3 `9 i. y! _4 N* K</span>   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">- J/ `% o( E2 b. k</span>   k++<font class=\\"jammer\\">% G5 |1 ~  S$ [6 x+ A3 q5 Y</font><br /> },<br /> <span style=\\"display:none\\">( v1 c\\" X' O# _/ w' b+ e</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">1 N9 N  z, m2 Y2 _* f4 N</span> <br /> <span style=\\"display:none\\">) i9 H+ g2 L/ T: G0 d) u7 U3 e' \( q  G</span> main1(::a)=<br /> <span style=\\"display:none\\">( x% u( E, P. k7 ?</span> t0=clock(), s=0,<font class=\\"jammer\\">  G8 t+ l% s% v5 @) j: L</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">7 C6 B' V: \! p\\" W% C2 k/ r6 d</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">7 F8 f& U) ~) v4 _  ^5 F</span>   k++<br /> <span style=\\"display:none\\">: K% i3 x3 u9 l# N6 s7 v</span> },<font class=\\"jammer\\">7 z4 V! b3 v+ C</font><br /> s;<font class=\\"jammer\\">' T6 i\\" H- b% v6 k  `, \</font><br /> [clock()-t0]/1000;

结果：
. i- L( O$ B& Y! W- d' e$ q1 i10.469
% O9 W: C2 u/ x; Z& }' V20000000.
8 E6 @0 e- d- ~% M3.781

" [4 b6 |' [7 T- F# Z－－－－－－－

/ B* O' L* W, f- T  H* G: v" N& z在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
( o4 ^' S& ]% A5 B6 I7 |  R在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
( P4 {% S8 K! u  y6 O' `9 a
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
. o' G( L  O* w
Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
0 w2 G( M  G4 b) |2 g5 ]- h7 G1 z+ a
7 k$ G8 k1 `; S! p5、在循环中加点东西

matlab代码：

1. >> clear all
   5 O0 i$ [+ A: l( x/ X0 t/ y4 K$ z
2. f=@(x)x+1;
   6 e: m/ l; W1 _( q% J8 f
3. a = linspace(2,2,1000);
   : n+ Y# Z; z. E; O  {* R2 O$ {
4. tic
   8 o: @! c+ t& A
5. s=0;1 ]' C  g* q- t; i  `( G( e
   
6. for k=1:1000
   $ m3 R# O0 @. @3 ^! x* f' m( @3 F
7.   for i=1:1000( j% ^$ {; N# q3 T7 |5 F! Z5 P
   
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));
   ! D1 _5 A& u2 e% n, M
9.   end' M% {1 I/ w9 D
   
10. end8 |  O\" F9 ]/ n7 d1 ~
    
11. s. V1 p8 {# Q. `& Z& V# w+ \
    
12. toc/ U& p3 y8 Y7 v. m2 L5 D' c
    
13. 
    # @+ ]% |6 [' H' P\" [4 ?6 {; d
14. s =) q, |5 Z& |/ y8 M6 r/ H
    
15. 
    + S( t! G; M# Y. e! s
16.   2.4002e+007+ m' U! |) _8 b, a7 Z. {
    
17. 4 z/ ~/ J5 d+ W$ s& |+ d
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];7 ~, [  q3 N5 Y
   
4. mvar:' }! u5 C7 Z8 I' `* y
   
6. f(x)=x+1;! {  I+ L1 U  X
   
8. main0(::a)=- q2 [5 b' m; v
   
10. a = linspace(2,2,1000),7 t: z. E# b6 S& K
    
12. t0=clock(), s=0,1 n) y9 P! a5 j
    
14. k=0, while{k<1000,
    
15. ! f7 D% N) ^! `. d! C# n/ ~) y
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. 9 X0 n3 r/ s9 A  i) u
18.   k++
    
19. % b8 e; D' n+ A( J: I5 `- y
20. },
    
21. - J! N; \8 o) y0 T7 v
22. s;
    
23. 9 }; }4 N6 p7 r. P4 v0 }  P\\" @
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">7 ~, [  q3 N5 Y</font><br /> mvar:<font class=\\"jammer\\">' }! u5 C7 Z8 I' `* y</font><br /> f(x)=x+1;<font class=\\"jammer\\">! {  I+ L1 U  X</font><br /> main0(::a)=<font class=\\"jammer\\">- q2 [5 b' m; v</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">7 t: z. E# b6 S& K</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">1 n) y9 P! a5 j</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">! f7 D% N) ^! `. d! C# n/ ~) y</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\\">9 X0 n3 r/ s9 A  i) u</span>   k++<br /> <span style=\\"display:none\\">% b8 e; D' n+ A( J: I5 `- y</span> },<br /> <span style=\\"display:none\\">- J! N; \8 o) y0 T7 v</span> s;<br /> <span style=\\"display:none\\">9 }; }4 N6 p7 r. P4 v0 }  P\\" @</span> [clock()-t0]/1000;

结果：
24002271.69189384
5.172

二者效率现在只有一倍差距。

, p1 Q7 d; B8 w9 R' ?－－－－－－－－－－
, @: ^. F' d1 c7 C' p, O6 M
$ f8 l# ^/ l# O: g) ^, r9 \- r) f不过，显然matlab和Forcal的程序应该这样写：
) \' ]  G: _% j- T+ i. ^1 P
8 I8 ?2 ^% m/ ^$ T, {matlab程序：

1. clear all
   % ?, L& f& Q; {: L; e\" I. N
2. f=@(x)x+1;4 s4 ?! v9 x' U4 I4 _2 Z3 f* v+ s; p
   
3. a = linspace(2,2,1000);# Y3 }% J- H4 u9 i& s3 ?: |
   
4. tic- z4 w5 {6 D; l' r6 L* P7 T, j
   
5. s=0;2 t3 o+ M3 A) p# j2 E% k! V
   
6. for k=1:1000$ U: y' {  C/ N, p
   
7.   for i=1:1000
   $ t) k3 }& o! g) K
8.     t=a(i);* x4 N, H) i! t$ V  t0 A\" a  }
   
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);
   . E; i6 n0 ~$ e% A- Q1 b$ y- n. L+ t
10.   end
    6 C9 P3 M$ {4 ~, c1 V: N' ]
11. end) z; H0 }- X$ i% \- r, k3 }
    
12. s' G% Z$ `% h; ~' U! @5 P) ?
    
13. toc$ |: C4 V$ i2 X6 B
    
14. 
    9 B9 R0 J\" A7 _# ?/ N3 n5 @
15. s =
    7 ]3 z\" X3 U, x! ]8 l) d3 ^. U9 B
16. : x2 x. N; e% @# m7 {3 E
    
17.   2.4002e+007% A. Q- A2 s( I: ^9 Y% ?- p! Z' n
    
18. # I# L% ?; K1 R5 _% D3 k\" A
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];* e+ V: i& E, S8 Z' H! S% Q* g+ U6 a
   
4. mvar:5 \' j\\" |  J. v' P5 _3 G9 G: C6 A
   
6. f(x)=x+1;
   
7. 0 D# b% G3 E4 K# G1 d. D\\" l
8. main0(::a)=$ q7 U: l+ |: x+ h9 G
   
10. a = linspace(2,2,1000),
    
11. 2 R\\" n% T+ }# r+ \% n0 ]3 N
12. t0=clock(), s=0,8 T% r% S( ?* t; \
    
14. k=0, while{k<1000,
    
15. 0 _4 q; C) m, w7 j; n
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++},
    
17. 4 s- \7 `0 e) O* @- y; [4 S  N
18.   k++. g- h9 L% b4 t; m7 j( ]
    
20. },( a+ Z: P& ]+ Q  A) t4 [
    
22. s;6 ^7 p  k4 q/ J
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">* e+ V: i& E, S8 Z' H! S% Q* g+ U6 a</font><br /> mvar:<font class=\\"jammer\\">5 \' j\\" |  J. v' P5 _3 G9 G: C6 A</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">0 D# b% G3 E4 K# G1 d. D\\" l</span> main0(::a)=<font class=\\"jammer\\">$ q7 U: l+ |: x+ h9 G</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">2 R\\" n% T+ }# r+ \% n0 ]3 N</span> t0=clock(), s=0,<font class=\\"jammer\\">8 T% r% S( ?* t; \</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">0 _4 q; C) m, w7 j; n</span>   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++},<br /> <span style=\\"display:none\\">4 s- \7 `0 e) O* @- y; [4 S  N</span>   k++<font class=\\"jammer\\">. g- h9 L% b4 t; m7 j( ]</font><br /> },<font class=\\"jammer\\">( a+ Z: P& ]+ Q  A) t4 [</font><br /> s;<font class=\\"jammer\\">6 ^7 p  k4 q/ J</font><br /> [clock()-t0]/1000;

结果：
% m6 C9 O/ l9 e5 L* Q24002271.69189384
+ B* p4 ]6 V/ W$ z' H( r1.625

matlab效率已不及Forcal。

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

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
3 ~5 i$ P. |6 |
Forcal代码：

2. !using["sys"];
   
3. 0 W/ s* M- l( v% s6 @& F( V9 T
4. mvar:
   
5. ! w9 \# {5 Q' B
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. , X3 T9 m- b) `5 G, g, [
8. t0=clock(),
   
9. . v5 Y) f& T* C$ v+ U( V7 }1 g; `4 l
10. k=0, while{k<10000,
    
11. 1 r- Z& O$ F5 d6 u! T4 a
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. : x) }* h& O1 K0 p
14.   k++- Z\\" _! r2 T5 X; y1 P
    
16. },# [& _( u' u  @
    
18. [clock()-t0]/1000;
    
19. 4 q\\" Z# k. [( ^+ L! G1 ]
20. 
    
21. ; A- Q2 X* X: i- }4 d5 {$ B
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23.   P. p( k' K\\" ^1 K: y8 M
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值/ w0 R! }! w. w+ l\\" t) b' t( e3 ?7 s
    
26. t0=clock(),& O! T+ D+ v4 ]  _2 C3 r: C
    
28. s=0,
    
29. % g' f* D0 a4 J0 \. R
30. k=0, while{k<10000,
    
31. % z9 H0 J3 g& o5 y. s
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. 4 L0 w! T% J\\" ]1 E. {& y
34.   k++( M  |3 m* ]$ ?- y+ m\\" a+ b8 I8 k
    
36. },
    
37.   P+ Y7 R/ \% J) V4 v: k0 j
38. s;
    
39. ) D8 ^0 R# x% Y! G0 h) g; Q
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">0 W/ s* M- l( v% s6 @& F( V9 T</span> mvar:<br /> <span style=\\"display:none\\">! w9 \# {5 Q' B</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">, X3 T9 m- b) `5 G, g, [</span> t0=clock(),<br /> <span style=\\"display:none\\">. v5 Y) f& T* C$ v+ U( V7 }1 g; `4 l</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">1 r- Z& O$ F5 d6 u! T4 a</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">: x) }* h& O1 K0 p</span>   k++<font class=\\"jammer\\">- Z\\" _! r2 T5 X; y1 P</font><br /> },<font class=\\"jammer\\"># [& _( u' u  @</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">4 q\\" Z# k. [( ^+ L! G1 ]</span> <br /> <span style=\\"display:none\\">; A- Q2 X* X: i- }4 d5 {$ B</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">  P. p( k' K\\" ^1 K: y8 M</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">/ w0 R! }! w. w+ l\\" t) b' t( e3 ?7 s</font><br /> t0=clock(),<font class=\\"jammer\\">& O! T+ D+ v4 ]  _2 C3 r: C</font><br /> s=0,<br /> <span style=\\"display:none\\">% g' f* D0 a4 J0 \. R</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">% z9 H0 J3 g& o5 y. s</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">4 L0 w! T% J\\" ]1 E. {& y</span>   k++<font class=\\"jammer\\">( M  |3 m* ]$ ?- y+ m\\" a+ b8 I8 k</font><br /> },<br /> <span style=\\"display:none\\">  P+ Y7 R/ \% J) V4 v: k0 j</span> s;<br /> <span style=\\"display:none\\">) D8 ^0 R# x% Y! G0 h) g; Q</span> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
20000000.
0.579      //取值用时间，秒