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

forcal

1、简单循环
, d' I  Q" y8 r
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   4 w  ^/ q5 [1 m  A1 \: Q2 y
2. tic$ t+ n* q' K) P! C7 \4 r5 C% }
   
3. s=0;
   9 q' F; t! m3 B4 l, e! Z( b* X
4. for k=1:10000
   9 u* P0 R+ |) G* n
5.   for i=1:1000
   6 i3 \* ?; @& Q8 c  h- R0 U: P% p
6.     s=s+1;% |, B+ B) E2 o  x' V# R
   
7.   end& E( B; H% ~1 G\" D4 k* K+ s. X
   
8. end8 T8 N7 G0 n6 f, @! {8 V$ N# H
   
9. s% r# q+ S9 x! o5 _3 [* \
   
10. toc
    . `4 g: g6 Z) y' v6 D
11. 
    6 u6 [: _2 x( _& l% u$ M
12. s =2 _0 p& B7 b: N8 F( D( i
    
13. + Z4 n9 Z' v0 y9 B3 n+ o
    
14.     10000000  C% H' D- o1 x1 K
    
15. / Z# R) Q& i0 s& P, [
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

2. !using["sys"];5 o  n+ j& `1 p; {6 y) h
   
4. mvar:4 k6 s8 R\\" j. ]
   
6. t0=clock(),9 E/ M/ N7 O) Q8 E; d
   
8. s=0,
   
9. ( g7 H; Q, w3 ~9 ~. Z! P8 N( a
10. k=0, while{k<10000,
    
11. 9 E* v' |% t% I  S. g; F
12.   i=0, while{i<1000,s=s+1, i++},; K4 ^\\" H  m/ b1 J$ g; x
    
14.   k++
    
15. : n! m& V7 [6 J  L0 V
16. },3 h- j9 k% @8 J. i\\" j
    
18. s;7 |; P* N4 `8 {7 f) p- D
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">5 o  n+ j& `1 p; {6 y) h</font><br /> mvar:<font class=\\"jammer\\">4 k6 s8 R\\" j. ]</font><br /> t0=clock(),<font class=\\"jammer\\">9 E/ M/ N7 O) Q8 E; d</font><br /> s=0,<br /> <span style=\\"display:none\\">( g7 H; Q, w3 ~9 ~. Z! P8 N( a</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">9 E* v' |% t% I  S. g; F</span>   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">; K4 ^\\" H  m/ b1 J$ g; x</font><br />   k++<br /> <span style=\\"display:none\\">: n! m& V7 [6 J  L0 V</span> },<font class=\\"jammer\\">3 h- j9 k% @8 J. i\\" j</font><br /> s;<font class=\\"jammer\\">7 |; P* N4 `8 {7 f) p- D</font><br /> [clock()-t0]/1000;

结果：
: m8 G& \7 a1 n! M3 m/ l10000000.
0.704  秒

二者运行效率相差7倍多。
, M  L8 D, Y( e& [! }
2、循环中稍微加点东西
: e+ i! N" _- c8 n$ P
Matlab代码：

1. clear all
   \" J. [, I) w5 `
2. tic$ ~1 O7 Q# E\" b: Z4 |, j& B4 p( N
   
3. s=0;
   ) f+ X) e4 S+ w+ M
4. for k=1:1000
   ( L, |# v( `* R* c
5.   for i=1:1000
   / g! H2 v: v$ A8 e1 f5 T% o
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);
   * j# p8 a. U$ |8 G5 C/ x
7.   end
   . w5 S; X! p& i- v5 ~
8. end+ p6 N% ?+ h/ G4 H8 A- d6 c* @
   
9. s
   1 ]\" |- z! }. K5 Z
10. toc6 Y- M) D3 s# W) u6 y5 o; n
    
11. * ]6 [2 @0 F1 O. I
    
12. s =
    2 E! y8 ^& U7 C8 P( a/ x
13. / B& H( J* b, E7 B' B$ u. r' R
    
14.   1.1009e+007
    4 c* E# \- x8 b6 O, |8 a7 j3 M
15. 4 n- Y9 D1 M% Z: T. S- N
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. ' q! u) U+ L. M  O2 O3 N
4. mvar:
   
5. 0 b' T& A6 J6 {# A/ o
6. t0=clock(),8 S! y# |. R( y2 x! W1 ?- w
   
8. s=0,( K  {6 a3 X& q' m\\" V
   
10. k=1, while{k<=1000,* ?% x0 M9 o5 m9 S\\" ?. w$ t
    
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++},: @, j6 p\\" h! M! f+ k\\" b4 v
    
14.   k++
    
15. 1 o3 c% D. B! t  f
16. },
    
17. . g, w+ S# Y5 [# L4 Q\\" j
18. s;
    
19. ! m! A2 ^; I$ I# O+ O
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">' q! u) U+ L. M  O2 O3 N</span> mvar:<br /> <span style=\\"display:none\\">0 b' T& A6 J6 {# A/ o</span> t0=clock(),<font class=\\"jammer\\">8 S! y# |. R( y2 x! W1 ?- w</font><br /> s=0,<font class=\\"jammer\\">( K  {6 a3 X& q' m\\" V</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">* ?% x0 M9 o5 m9 S\\" ?. w$ t</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\\">: @, j6 p\\" h! M! f+ k\\" b4 v</font><br />   k++<br /> <span style=\\"display:none\\">1 o3 c% D. B! t  f</span> },<br /> <span style=\\"display:none\\">. g, w+ S# Y5 [# L4 Q\\" j</span> s;<br /> <span style=\\"display:none\\">! m! A2 ^; I$ I# O+ O</span> [clock()-t0]/1000;

结果：
& L% v+ [  h0 W# L: b0 Y11008953.66602346
& a; ]/ o# i6 X1.031

二者效率相差已比较小。

3、再增加一个函数调用

Matlab代码：

1. clear all
   7 e7 b' i! g0 U/ [9 r
2. f=@(x)x+1;% x2 t* @$ \& ?# w\" g3 Z. Q+ g
   
3. tic8 S' C! T( ^1 H# f9 \9 W
   
4. s=0;
   & I: @8 h9 m. E: n5 L) m$ D$ T  `5 S
5. for k=1:1000
   - V+ S9 e8 Y+ A$ ?8 D  G& ^
6.   for i=1:1000# }( ?6 B& Q% e9 H  @# n
   
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);) y. e7 C! Q+ I$ f) n1 P: R, x# {
   
8.   end
   : }; [, ?6 n3 v0 C, k; n0 R
9. end' D5 J0 }+ L) R5 T0 J* G( N8 x
   
10. s
    , ^% ~' h9 G  N- c3 D\" a. ?  O
11. toc  K7 m3 K\" C! Q) ]3 J4 B
    
12. 
    ; |. M# \3 r' k* U
13. s =+ o* b+ [! i4 D7 `  u
    
14. 
    : x) d- U  X: K% @$ g
15.   5.1251e+0083 E; }2 [$ G* i6 d* s( o9 w
    
16. 
    7 ^  B% ~3 s! p3 O
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];4 a) X1 m9 Y7 a0 G% @. y\\" \
   
4. f(x)=x+1;: D; V7 }' w$ y4 `/ Z
   
6. mvar:+ b! J+ ~5 |5 J9 B2 a3 D\\" c
   
8. t0=clock(),
   
9. ) e9 e0 U* J! m/ ?2 o
10. s=0,3 ?; }' U% J& O0 }
    
12. k=1, while{k<=1000,
    
13. \\" h- e; w! ?  j5 Z2 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++},4 N/ v! `9 {% j! c7 @
    
16.   k++
    
17. 1 J0 i' ?0 t  Z
18. },( X' d* S* ~  `2 A5 ]
    
20. s;
    
21. 7 L+ `' A\\" ^1 C: C' C( L
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">4 a) X1 m9 Y7 a0 G% @. y\\" \</font><br /> f(x)=x+1;<font class=\\"jammer\\">: D; V7 }' w$ y4 `/ Z</font><br /> mvar:<font class=\\"jammer\\">+ b! J+ ~5 |5 J9 B2 a3 D\\" c</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">) e9 e0 U* J! m/ ?2 o</span> s=0,<font class=\\"jammer\\">3 ?; }' U% J& O0 }</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">\\" h- e; w! ?  j5 Z2 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++},<font class=\\"jammer\\">4 N/ v! `9 {% j! c7 @</font><br />   k++<br /> <span style=\\"display:none\\">1 J0 i' ?0 t  Z</span> },<font class=\\"jammer\\">( X' d* S* ~  `2 A5 ]</font><br /> s;<br /> <span style=\\"display:none\\">7 L+ `' A\\" ^1 C: C' C( L</span> [clock()-t0]/1000;

结果：
512508953.6658993
1.078

可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

' R) B) b; k" |2 C/ b+ G" zMatlab代码：

1. >> clear all
   ( w' l% C( c5 q& ?2 z/ ?
2. a = linspace(1,1,1000);& x1 V5 a* u, G6 j6 m
   
3. tic  \- J5 T; i) z/ m
   
4. for k=1:10000
   , l, v\" i6 M  K& [
5.   for i=1:1000, N( M; |/ ]! x2 x: F7 R\" L
   
6.     a(i)=2;/ R* a) D9 B# P. W; s' u1 q0 [5 s
   
7.   end
     J$ E! J! ~8 u
8. end
   + q! }% {- F& [9 M. a$ [
9. toc
   8 y2 O5 E6 x. D8 Z! q# H
10. 
    + F. v2 ~) a/ r; `+ \
11. tic
    & x) U1 V4 z0 x% I9 }% ]/ y
12. s=0;5 Z# c* E. P) R3 D% `1 d
    
13. for k=1:10000: j: `9 B1 P2 @
    
14.   for i=1:1000  j7 {! N) s: _6 e
    
15.     s=s+a(i);
    2 w. ]2 i1 ~5 q\" U9 E4 e9 h6 O$ @: L5 A
16.   end
    : R8 Q& [\" s/ H8 @- z5 |
17. end
    * c/ t0 I9 }  Z! r& L2 ]
18. s- p\" i, F3 g4 G7 _' w
    
19. toc3 t  l% d& n& z- {/ J; }
    
20. Elapsed time is 0.140658 seconds.
    , Y9 P( F! t$ J& Z8 X* q' ?; v7 b0 e
21. ; Y- Y9 m' x/ a# }* m
    
22. s =5 y0 t' _$ Y) c
    
23. 
    \" W- w! y5 m$ r* y
24.     20000000; z7 }\" x7 }6 X5 o* |
    
25. : l# ]1 S4 B1 K3 G5 M
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];+ l7 |3 }2 Z' p6 f! F; [7 P( T
   
4. mvar:+ E% B/ l\\" j; s8 g& G
   
6. main0(::a)=  Q3 z  |. r3 F8 w% z
   
8. a = linspace(1,1,1000),
   
9. # S: f$ E0 v$ R: \5 }
10. t0=clock(),
    
11. # R; \! V- s  Q\\" c$ c9 [\\" n
12. k=0, while{k<10000,
    
13. # G3 G- E3 {* K2 L$ x
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. 8 ~9 q7 ]# L% i* m4 L
16.   k++
    
17. # y, u1 T8 l: f1 G7 G
18. },- t' f8 c. p7 E- {' s
    
20. [clock()-t0]/1000;
    
21. ! z9 J4 n: s) M$ S1 n: N
22. 
    
23. 4 Z( `' d% @- U! j
24. main1(::a)=
    
25. 8 l3 v3 v: f8 U/ a) J' u7 \4 d
26. t0=clock(), s=0,8 x2 Z6 L0 v/ K! x! @
    
28. k=0, while{k<10000,% ^' p( y  _5 f3 M7 b$ N, F
    
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. 0 T; B1 a8 O; F! m7 F; i- ^\\" ^
32.   k++
    
33. ( p( T4 [% `) B\\" Z# e
34. },
    
35. - _( d3 M8 q6 B, M& l$ L7 C
36. s;
    
37. $ Q6 M( B% F7 a& w+ l( ]6 a
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ l7 |3 }2 Z' p6 f! F; [7 P( T</font><br /> mvar:<font class=\\"jammer\\">+ E% B/ l\\" j; s8 g& G</font><br /> main0(::a)=<font class=\\"jammer\\">  Q3 z  |. r3 F8 w% z</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\"># S: f$ E0 v$ R: \5 }</span> t0=clock(),<br /> <span style=\\"display:none\\"># R; \! V- s  Q\\" c$ c9 [\\" n</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\"># G3 G- E3 {* K2 L$ x</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">8 ~9 q7 ]# L% i* m4 L</span>   k++<br /> <span style=\\"display:none\\"># y, u1 T8 l: f1 G7 G</span> },<font class=\\"jammer\\">- t' f8 c. p7 E- {' s</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">! z9 J4 n: s) M$ S1 n: N</span> <br /> <span style=\\"display:none\\">4 Z( `' d% @- U! j</span> main1(::a)=<br /> <span style=\\"display:none\\">8 l3 v3 v: f8 U/ a) J' u7 \4 d</span> t0=clock(), s=0,<font class=\\"jammer\\">8 x2 Z6 L0 v/ K! x! @</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">% ^' p( y  _5 f3 M7 b$ N, F</font><br />   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">0 T; B1 a8 O; F! m7 F; i- ^\\" ^</span>   k++<br /> <span style=\\"display:none\\">( p( T4 [% `) B\\" Z# e</span> },<br /> <span style=\\"display:none\\">- _( d3 M8 q6 B, M& l$ L7 C</span> s;<br /> <span style=\\"display:none\\">$ Q6 M( B% F7 a& w+ l( ]6 a</span> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
1 I8 q  r- S3 p/ l/ ?+ \20000000.
5.531   //取数据所用时间，秒
, i$ t. p8 F! ^7 ~
5 ]- `/ _" C. y5 H. [% d－－－－－－－－－－
5 w* `3 d: f( S% ~+ Y9 z
1 O2 m1 g4 Y- R% D4 f* ]9 t7 ]用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. # ]2 w3 d# v5 }8 l$ b/ r\\" n' K
4. mvar:( T6 p$ I8 M3 B
   
6. main0(::a)=- W) ^5 `5 B  c! _/ S/ K& k# B/ ^+ r
   
8. a = linspace(1,1,1000),
   
9. + A4 \+ E3 `% q0 o; k4 y4 u4 V) v) _  E4 o
10. t0=clock(),
    
11. \\" A; K8 Q- u9 U2 M. d\\" E8 u
12. k=0, while{k<10000,8 ]4 B: U: |) m* ^* B$ H
    
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. % S! w7 Y7 V1 r8 B6 j+ w6 a$ @8 b- F
16.   k++
    
17. / Q+ {4 c$ b$ P! B. n, x. o
18. },
    
19. 3 I8 g+ Y. t* v5 c0 I
20. [clock()-t0]/1000;
    
21. 3 H* d% I# ~1 g4 _4 I
22.   T* W, ]+ {) S
    
24. main1(::a)=\\" ~2 f, b( z5 b% M
    
26. t0=clock(), s=0,. b9 G4 ?( O3 P
    
28. k=0, while{k<10000,
    
29. . B3 v: M1 O% ?
30.   i=0, while{i<1000, s=s+A(a,i), i++},( q! T  f1 X* s! m
    
32.   k++) j0 q% @; R! U
    
34. },& I6 c+ E* E! f% w
    
36. s;  k& Q- J3 y& y: Z( U
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\"># ]2 w3 d# v5 }8 l$ b/ r\\" n' K</span> mvar:<font class=\\"jammer\\">( T6 p$ I8 M3 B</font><br /> main0(::a)=<font class=\\"jammer\\">- W) ^5 `5 B  c! _/ S/ K& k# B/ ^+ r</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">+ A4 \+ E3 `% q0 o; k4 y4 u4 V) v) _  E4 o</span> t0=clock(),<br /> <span style=\\"display:none\\">\\" A; K8 Q- u9 U2 M. d\\" E8 u</span> k=0, while{k<10000,<font class=\\"jammer\\">8 ]4 B: U: |) m* ^* B$ H</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">% S! w7 Y7 V1 r8 B6 j+ w6 a$ @8 b- F</span>   k++<br /> <span style=\\"display:none\\">/ Q+ {4 c$ b$ P! B. n, x. o</span> },<br /> <span style=\\"display:none\\">3 I8 g+ Y. t* v5 c0 I</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">3 H* d% I# ~1 g4 _4 I</span> <font class=\\"jammer\\">  T* W, ]+ {) S</font><br /> main1(::a)=<font class=\\"jammer\\">\\" ~2 f, b( z5 b% M</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">. b9 G4 ?( O3 P</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">. B3 v: M1 O% ?</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">( q! T  f1 X* s! m</font><br />   k++<font class=\\"jammer\\">) j0 q% @; R! U</font><br /> },<font class=\\"jammer\\">& I6 c+ E* E! f% w</font><br /> s;<font class=\\"jammer\\">  k& Q- J3 y& y: Z( U</font><br /> [clock()-t0]/1000;

结果：
) _5 ^' x  ^8 Z! |/ q9 U! ^10.469
20000000.
3.781

/ s& K& K5 i  h3 \; e2 T－－－－－－－

在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
3 E+ x" U1 U5 Y& Q9 N! z3 j在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

7 D& ~- o; f; ]& lmatlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
) e: V/ |, |6 A+ ]3 r2 b% @
+ f' R1 C$ R% d: L+ oForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

; f; Y% ~& |- X8 t' V5、在循环中加点东西

, \( N7 W5 J; u4 w% E9 n. ?matlab代码：

1. >> clear all
   ; P0 b8 f$ u9 _( r  k\" e& b# d
2. f=@(x)x+1;
   ) }- o- D& K2 S! k
3. a = linspace(2,2,1000);
   5 c* e( A! e( o* R, ^
4. tic
   ; M7 ?! @' X! u8 b. w$ ?
5. s=0;0 h6 L, A7 l, z  j* n
   
6. for k=1:1000
   8 G1 X* @9 \) }5 M: g' c* ]
7.   for i=1:10001 ~! s8 g/ c% u  e) X% S
   
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));6 y# A1 i: c# X) R$ G8 _' r! y. f
   
9.   end4 R% [  o6 U: k
   
10. end
    . _2 l8 I8 H$ T' Y' ?
11. s
    \" ?9 E5 x+ L4 Q& q
12. toc( G$ j$ c% T  d2 X' h
    
13. $ R; z0 W: y7 r7 @' o\" ~
    
14. s =
    ) C# V/ ~/ V4 M) W7 @
15. 
    \" R* Z) I' p\" c+ m% n( B. A4 |
16.   2.4002e+0071 F; c) W/ f. M; R; I( `8 q( ?* p
    
17. 
    ! _1 w\" @& G8 e# C: k+ D
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];/ d0 Z5 ?* B' U, G/ C7 H* U% `& y) h
   
4. mvar:
   
5. , o0 }2 J/ c: V+ I* Q
6. f(x)=x+1;2 S, `7 Y$ @/ g; j+ z/ c6 X
   
8. main0(::a)=
   
9. \\" a+ V$ A0 H/ a
10. a = linspace(2,2,1000),
    
11. $ o7 O& Z\\" O* v: x/ d
12. t0=clock(), s=0,# @5 `4 A. ~! b; `! L0 u
    
14. k=0, while{k<1000,
    
15. 6 C. s2 L6 k6 T7 F1 A% r
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++},# W1 \+ F1 O: K' A# V
    
18.   k++\\" s/ B+ [- Y: o' ^
    
20. },\\" N) Y/ Y5 x* f; ?5 r
    
22. s;/ z9 Z1 r& ]\\" t0 ]- E: }5 \( u
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">/ d0 Z5 ?* B' U, G/ C7 H* U% `& y) h</font><br /> mvar:<br /> <span style=\\"display:none\\">, o0 }2 J/ c: V+ I* Q</span> f(x)=x+1;<font class=\\"jammer\\">2 S, `7 Y$ @/ g; j+ z/ c6 X</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">\\" a+ V$ A0 H/ a</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">$ o7 O& Z\\" O* v: x/ d</span> t0=clock(), s=0,<font class=\\"jammer\\"># @5 `4 A. ~! b; `! L0 u</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">6 C. s2 L6 k6 T7 F1 A% r</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++},<font class=\\"jammer\\"># W1 \+ F1 O: K' A# V</font><br />   k++<font class=\\"jammer\\">\\" s/ B+ [- Y: o' ^</font><br /> },<font class=\\"jammer\\">\\" N) Y/ Y5 x* f; ?5 r</font><br /> s;<font class=\\"jammer\\">/ z9 Z1 r& ]\\" t0 ]- E: }5 \( u</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
4 x5 _/ a# @9 h" E! ]6 v( T5.172
* Z2 `! y+ ?6 e9 o) x; l
+ a( i* g4 R) W( T, V5 e+ t二者效率现在只有一倍差距。

8 X: C# e/ Y  ~2 e－－－－－－－－－－
6 F" E, B/ K6 f- v- J( ^: f; h( [
不过，显然matlab和Forcal的程序应该这样写：
% \4 _7 a( ]# h4 W( J& v5 O, @0 h" G% t' j
# U/ R, R) c; A& mmatlab程序：

1. clear all
   9 d1 ^, E: B2 f  X
2. f=@(x)x+1;4 [+ K- K% z  N/ k9 c8 \
   
3. a = linspace(2,2,1000);7 }- s2 j  K2 M- a  _; s\" U
   
4. tic+ r( F2 K+ Y& W2 Q$ H, e7 E
   
5. s=0;( _4 u9 _9 M4 s# u
   
6. for k=1:1000
   , m5 m! K7 o4 g' u' T
7.   for i=1:1000( m4 q+ P3 O  P( `, S
   
8.     t=a(i);
   / R5 A: f# H% C. l1 M9 N6 @
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);
   2 _: t; B6 J0 h, T3 z\" L0 Z
10.   end
    ) x6 j+ K# C! O: A+ I5 f  \9 N
11. end% h5 c, o( d0 Y& S
    
12. s
      S0 ^2 E; {8 `* r- `
13. toc# V0 R; z2 o) G/ r
    
14. 
    1 D8 d, K) d9 B' R
15. s =
    2 q; e) `4 O/ ?, w\" l
16. 
    9 T7 k3 L6 ^\" g# J
17.   2.4002e+007
    # f( F: x, B4 k( D
18. . H4 n* h1 A9 R
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];4 K' W, p* y* i2 d/ j
   
4. mvar:4 c& I0 u! N% ?+ {( T4 U  B
   
6. f(x)=x+1;
   
7. & q7 X* Q; L) L6 ]( t$ j
8. main0(::a)=7 U7 r- S4 e! K) b
   
10. a = linspace(2,2,1000),% X8 ]# M( x4 E: |\\" H
    
12. t0=clock(), s=0,
    
13. ) b. N2 M* T6 V1 @! L- H\\" t- a& U
14. k=0, while{k<1000,
    
15. \\" S7 R6 ?- _0 `1 Y\\" U! b
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. + t3 Y\\" s( r1 C1 R0 l/ a
18.   k++- m\\" n5 z6 R\\" i8 {; [- F3 [
    
20. },
    
21. % O- Q* C. t; e# _; ^2 D
22. s;8 u6 A/ \! l\\" b. W: i4 W! x
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">4 K' W, p* y* i2 d/ j</font><br /> mvar:<font class=\\"jammer\\">4 c& I0 u! N% ?+ {( T4 U  B</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">& q7 X* Q; L) L6 ]( t$ j</span> main0(::a)=<font class=\\"jammer\\">7 U7 r- S4 e! K) b</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">% X8 ]# M( x4 E: |\\" H</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">) b. N2 M* T6 V1 @! L- H\\" t- a& U</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">\\" S7 R6 ?- _0 `1 Y\\" U! b</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\\">+ t3 Y\\" s( r1 C1 R0 l/ a</span>   k++<font class=\\"jammer\\">- m\\" n5 z6 R\\" i8 {; [- F3 [</font><br /> },<br /> <span style=\\"display:none\\">% O- Q* C. t; e# _; ^2 D</span> s;<font class=\\"jammer\\">8 u6 A/ \! l\\" b. W: i4 W! x</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
: A3 c: l, O" n" H1.625
4 I0 ?% O$ l+ d2 b  M
matlab效率已不及Forcal。
7 ?) Q& c# y, v9 t9 _
6 P+ ~) K: j  n2 S" r. p8 Y; \提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。

Forcal代码：

2. !using["sys"];% n: _& V7 r$ K, ~; ~7 l5 A
   
4. mvar:$ W! E) y5 H; ?+ T4 D  g
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. / `! |$ R( j5 e\\" U$ j: \7 G
8. t0=clock(),- U! N- t) o# \: J4 ?- h5 z4 e
   
10. k=0, while{k<10000,& I, v( A4 Z# ]- v8 K. w+ Y0 t3 u\\" l. n
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值0 S7 _( r1 R) H5 k2 B, d9 X9 V
    
14.   k++
    
15. ! r; \- s9 S1 K8 u7 @' w
16. },- l7 q. [# i5 e2 l' ]
    
18. [clock()-t0]/1000;' w8 V9 n6 c6 g8 R, h2 K
    
20. 
    
21. . ^( r1 D  A* Y+ Q- D8 u0 v
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. % ?\\" K# {- L  p# U
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值  ]- c4 {, ]. V% n9 e
    
26. t0=clock(),
    
27. - ^# g  `) r( u  L' K8 q
28. s=0,3 x8 F3 Z% a4 s. @( ^- E$ ~' c3 }8 O
    
30. k=0, while{k<10000,
    
31. 3 d5 t! `. L9 f( V' s+ \5 e% u
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值. l6 O% }8 w8 h& N' |/ \( s! \% F
    
34.   k++& l\\" m( Z+ y1 m; p
    
36. },3 J& J) B7 ]. A4 [
    
38. s;' ~: v, [\\" r7 W* U4 n. l, k
    
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">% n: _& V7 r$ K, ~; ~7 l5 A</font><br /> mvar:<font class=\\"jammer\\">$ W! E) y5 H; ?+ T4 D  g</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">/ `! |$ R( j5 e\\" U$ j: \7 G</span> t0=clock(),<font class=\\"jammer\\">- U! N- t) o# \: J4 ?- h5 z4 e</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">& I, v( A4 Z# ]- v8 K. w+ Y0 t3 u\\" l. n</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">0 S7 _( r1 R) H5 k2 B, d9 X9 V</font><br />   k++<br /> <span style=\\"display:none\\">! r; \- s9 S1 K8 u7 @' w</span> },<font class=\\"jammer\\">- l7 q. [# i5 e2 l' ]</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">' w8 V9 n6 c6 g8 R, h2 K</font><br /> <br /> <span style=\\"display:none\\">. ^( r1 D  A* Y+ Q- D8 u0 v</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">% ?\\" K# {- L  p# U</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">  ]- c4 {, ]. V% n9 e</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">- ^# g  `) r( u  L' K8 q</span> s=0,<font class=\\"jammer\\">3 x8 F3 Z% a4 s. @( ^- E$ ~' c3 }8 O</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">3 d5 t! `. L9 f( V' s+ \5 e% u</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">. l6 O% }8 w8 h& N' |/ \( s! \% F</font><br />   k++<font class=\\"jammer\\">& l\\" m( Z+ y1 m; p</font><br /> },<font class=\\"jammer\\">3 J& J) B7 ]. A4 [</font><br /> s;<font class=\\"jammer\\">' ~: v, [\\" r7 W* U4 n. l, k</font><br /> [clock()-t0]/1000;

结果：
. v. Y6 |. T0 g# z0.5620000000000001  //赋值用时间，秒
* e4 C; Q4 O& H8 ^20000000.
0.579      //取值用时间，秒