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

forcal

1、简单循环
7 a/ w  X2 t2 h: {* j0 g5 m' q
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all3 [% g% X: A( m  L' ~
   
2. tic' v- S' q, D2 `8 [$ `- o3 z
   
3. s=0;
     x) ?9 z4 s* x2 s2 o) q$ P
4. for k=1:10000
   + W8 h% k  I- A5 f
5.   for i=1:1000
   * ]  g& h) G$ {0 x- s& `6 X
6.     s=s+1;, I% Z' B\" U$ C' O  N# ^4 Y
   
7.   end9 I6 g  Y9 }/ ?
   
8. end- ?; K9 G. Q: c# i7 ~4 j1 f2 e
   
9. s2 E( i  n8 q7 T9 c( x1 H
   
10. toc, j4 L, A  V7 p' V4 D
    
11. 
    & i8 q# _* i3 ]3 X) ~
12. s =
    - \8 ]3 |. m\" j: d  F
13. * ?$ b2 S- l+ o; \' P
    
14.     100000001 [( z. x\" B  h5 a
    
15. 
    ( A3 d' B9 e0 _' B# z( S  i, i
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
$ m5 c! C7 b% @4 o+ L
% C  d8 P! R- V! G+ DForcal没有JIT加速器，代码：

2. !using["sys"];5 S# H9 f3 e) N2 _; D
   
4. mvar:
   
5. 3 f+ R7 A' c8 W0 C
6. t0=clock(),( P6 Y* C7 U+ f) d3 U
   
8. s=0,
   
9. 4 U9 ?\\" d4 [# f9 f! c0 @; S8 C7 P+ K
10. k=0, while{k<10000,, {# e5 N% g& T# s
    
12.   i=0, while{i<1000,s=s+1, i++},
    
13. 7 L$ F- N% l. j% Y\\" p
14.   k++
    
15. $ X# j1 S9 I  s- A$ h2 \! e
16. },
    
17. ' K# G& r# x+ f* P; h: V
18. s;\\" z  {0 g2 ]& p! G* B5 e9 c
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">5 S# H9 f3 e) N2 _; D</font><br /> mvar:<br /> <span style=\\"display:none\\">3 f+ R7 A' c8 W0 C</span> t0=clock(),<font class=\\"jammer\\">( P6 Y* C7 U+ f) d3 U</font><br /> s=0,<br /> <span style=\\"display:none\\">4 U9 ?\\" d4 [# f9 f! c0 @; S8 C7 P+ K</span> k=0, while{k<10000,<font class=\\"jammer\\">, {# e5 N% g& T# s</font><br />   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">7 L$ F- N% l. j% Y\\" p</span>   k++<br /> <span style=\\"display:none\\">$ X# j1 S9 I  s- A$ h2 \! e</span> },<br /> <span style=\\"display:none\\">' K# G& r# x+ f* P; h: V</span> s;<font class=\\"jammer\\">\\" z  {0 g2 ]& p! G* B5 e9 c</font><br /> [clock()-t0]/1000;

结果：
10000000.
  v: i2 b% e3 T0 l0 _! Z, U0.704  秒
- m; _& }) m. G/ L" y1 B% S- i; y2 w( x9 ?
: o' R2 e% h4 {) T% d二者运行效率相差7倍多。
! U& C- `7 F$ Z5 i# A. s
2、循环中稍微加点东西

Matlab代码：

1. clear all, {2 @1 a: h2 W% B
   
2. tic+ H4 q, k3 n) u1 Y\" s' ?8 `4 Z1 i
   
3. s=0;
   ' g\" s( e- F! I: S4 ]
4. for k=1:1000* Z0 g: c  y9 P1 S% O2 s0 j3 ~
   
5.   for i=1:1000
   % W; r6 M9 I  L* b
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);8 o) c: [% }& [4 t
   
7.   end# x' _; [/ K5 X. W* y* j& W
   
8. end
   4 u4 f( S, Q! g
9. s4 g\" k& C, x* [; r% f\" [$ B' N* L
   
10. toc% e2 g& g* C) h8 ?7 P
    
11. 
    % }8 A# L7 {! J! }, ^
12. s =
    9 h; n$ Z/ R/ k5 [; Q
13. 
      c8 R9 t# K; T+ p
14.   1.1009e+007$ a8 k& {$ l9 w$ P7 Z
    
15. & A; V. Q6 c! d  k9 ^
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. - R) n( n' \9 o; U; |+ E9 U% P0 d
4. mvar:. P, p: `* B7 a4 H) v# l! ]\\" W
   
6. t0=clock(),
   
7. # ~1 A5 r8 _; C$ u- q4 M. R
8. s=0,0 a3 K0 i# r; k0 I
   
10. k=1, while{k<=1000,
    
11. ' K% o% ^+ ]4 P# D1 K  {: [
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++},4 Z8 b  M/ y  d2 W6 `
    
14.   k++
    
15. % g  K; o6 U5 A% z+ O+ P
16. },1 `/ q: W$ w4 e1 [( n\\" }
    
18. s;
    
19. ; z. q- i9 F6 |$ A+ ~, \4 u
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">- R) n( n' \9 o; U; |+ E9 U% P0 d</span> mvar:<font class=\\"jammer\\">. P, p: `* B7 a4 H) v# l! ]\\" W</font><br /> t0=clock(),<br /> <span style=\\"display:none\\"># ~1 A5 r8 _; C$ u- q4 M. R</span> s=0,<font class=\\"jammer\\">0 a3 K0 i# r; k0 I</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">' K% o% ^+ ]4 P# D1 K  {: [</span>   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\\">4 Z8 b  M/ y  d2 W6 `</font><br />   k++<br /> <span style=\\"display:none\\">% g  K; o6 U5 A% z+ O+ P</span> },<font class=\\"jammer\\">1 `/ q: W$ w4 e1 [( n\\" }</font><br /> s;<br /> <span style=\\"display:none\\">; z. q- i9 F6 |$ A+ ~, \4 u</span> [clock()-t0]/1000;

结果：
  N4 W8 P0 Q5 l* ]' U- x8 Q+ s" [11008953.66602346
0 m3 W2 J" f9 K! P1.031
1 j8 l" H: A3 X7 X5 N3 ?/ c( c
二者效率相差已比较小。
- e1 @5 Y% }! w8 m6 d% l1 o" l
3、再增加一个函数调用

Matlab代码：

1. clear all
   2 @; w' d8 R0 d6 |
2. f=@(x)x+1;
   ' E: y$ m* j+ Y# C- T3 F  E
3. tic
   - I2 M( Z\" m& Z; h
4. s=0;
   \" M& j1 L  v' m- ]6 l8 R* e
5. for k=1:10002 ~6 r. ^: U$ s+ ?0 k+ g
   
6.   for i=1:1000! {) a4 l* p) `& ]+ P: {
   
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);
   % Q* m\" [6 a/ Q- R7 q1 Z
8.   end/ m\" z$ ~7 Q( V/ M# S2 a: R% _# A
   
9. end  O+ c( e4 j+ z& W' m
   
10. s2 e% j* [5 P4 A; V. I: Q; E& n
    
11. toc
    \" [  [0 R) B/ U+ n8 x9 R
12. 
    ; u  H+ Z$ ?  x1 Q+ o
13. s =4 x/ x6 @2 `+ P; s: m
    
14. 
    1 U+ s( N# g# Z+ _
15.   5.1251e+008
    , \  }6 z7 ~9 \8 v9 A. C
16. ( w9 |7 W3 m3 b, d2 S
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];# `0 \7 J( K8 t; k4 F6 e
   
4. f(x)=x+1;* S5 C: Y% i$ q: c' u# g. y
   
6. mvar:2 Y5 `9 }  x- F, }& u2 J) J
   
8. t0=clock(),
   
9. , _( [: K- I: _) ]* P
10. s=0,$ L2 \) X! ~9 U; O9 M2 P8 C
    
12. k=1, while{k<=1000,
    
13. , B9 l- w# X7 z9 s2 ~7 |
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++},) i+ ]# I3 C4 Q' C
    
16.   k++
    
17. 9 e# ]: Y; v1 A2 g. M
18. },6 G5 N3 G5 y/ d& Q4 w
    
20. s;
    
21. ) E1 O% x5 E, N/ ]- w
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\"># `0 \7 J( K8 t; k4 F6 e</font><br /> f(x)=x+1;<font class=\\"jammer\\">* S5 C: Y% i$ q: c' u# g. y</font><br /> mvar:<font class=\\"jammer\\">2 Y5 `9 }  x- F, }& u2 J) J</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">, _( [: K- I: _) ]* P</span> s=0,<font class=\\"jammer\\">$ L2 \) X! ~9 U; O9 M2 P8 C</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">, B9 l- w# X7 z9 s2 ~7 |</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\\">) i+ ]# I3 C4 Q' C</font><br />   k++<br /> <span style=\\"display:none\\">9 e# ]: Y; v1 A2 g. M</span> },<font class=\\"jammer\\">6 G5 N3 G5 y/ d& Q4 w</font><br /> s;<br /> <span style=\\"display:none\\">) E1 O% x5 E, N/ ]- w</span> [clock()-t0]/1000;

结果：
512508953.6658993
  v( C; @0 p. U1 n1.078

( q7 J4 c2 `( j: t& J可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

Matlab代码：

1. >> clear all. f\" Q. D$ G7 H  b* V) }, [
   
2. a = linspace(1,1,1000);# h/ K& x' Y) M& z5 f
   
3. tic8 `1 f  T) H# I3 u
   
4. for k=1:10000
   : S% L\" c7 g' _
5.   for i=1:1000\" ~- ?\" t6 h( e( d, P7 ~) r
   
6.     a(i)=2;
   & c1 w  p$ b; V
7.   end
   ' v, D6 d# q; ]7 O/ s( t4 E% D( x
8. end
   % g* a. q. E4 X; C: U- q
9. toc
   % H4 U+ O+ e! Z; F1 G& b, [( v
10. 
    % i0 B/ A, G( N, R' d3 M; e. F1 @
11. tic
    8 x6 _6 [- s  R* v+ s4 M( `6 c
12. s=0;
    3 N! S8 @/ J# J$ D2 R6 L1 K
13. for k=1:10000
    * H8 S) s* U+ {! u/ g# v. Q& Q
14.   for i=1:1000
    ) f% w$ ~/ q\" `
15.     s=s+a(i);. }: x. T- j) W\" J$ f
    
16.   end\" m; H- U2 N- I
    
17. end
    5 Q$ m5 T( |: U
18. s2 j% U( n+ I& _% }' X5 C
    
19. toc( m& u, }8 b4 \/ V; c
    
20. Elapsed time is 0.140658 seconds.
    $ O+ `' A1 `6 U. Q( Z& q
21. 0 I7 i3 M+ T2 @' U' j
    
22. s =. `- H. e/ U: L6 F# f9 y
    
23. 8 g) u* D4 y- A6 r  F0 Y- H: i
    
24.     20000000
    7 ^; t% y! R\" n3 [7 y# Q: q\" q
25. 
    6 C- k% {; ]6 T6 n. Q: |
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
. B) \5 j& _# m% ]: _
Forcal代码：

2. !using["math","sys"];
   
3.   c9 `0 f: f- k6 |2 _5 N) l
4. mvar:
   
5. 1 }1 j! U7 b, R6 n7 P# \4 q
6. main0(::a)=
   
7. : ~& [. w& q8 o* P$ N
8. a = linspace(1,1,1000),* @; S- K# T\\" W3 J& T- g+ t
   
10. t0=clock(),8 R( _& d, \  t2 {- ~
    
12. k=0, while{k<10000,
    
13. 4 J8 d4 L% B) E9 I2 D
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. , ~2 z1 F( b) j
16.   k++
    
17. ' b6 c  n% l) d! X) O8 s6 t
18. },; E$ T8 o/ M% J6 z
    
20. [clock()-t0]/1000;
    
21. ) [  L! c/ M' J. Q\\" e, N. t7 |1 k. V
22. 8 G$ l! W\\" O: u\\" {' z* r\\" b
    
24. main1(::a)=
    
25. ' A- r# p\\" `+ J8 ~- A1 u  f
26. t0=clock(), s=0,
    
27. 2 {, \( e$ C6 J; H& r\\" N# @; }
28. k=0, while{k<10000,# x4 E5 N- E  r! Z
    
30.   i=0, while{i<1000, s=s+a(i), i++},  b4 M6 Y5 ~/ Q  S5 k
    
32.   k++- q- u# v+ Q1 t( h+ R
    
34. },
    
35. ( E! ~/ d9 z7 D, Y1 R7 y
36. s;
    
37. 9 c% |, z- h' l( S2 N. ]- A
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">  c9 `0 f: f- k6 |2 _5 N) l</span> mvar:<br /> <span style=\\"display:none\\">1 }1 j! U7 b, R6 n7 P# \4 q</span> main0(::a)=<br /> <span style=\\"display:none\\">: ~& [. w& q8 o* P$ N</span> a = linspace(1,1,1000),<font class=\\"jammer\\">* @; S- K# T\\" W3 J& T- g+ t</font><br /> t0=clock(),<font class=\\"jammer\\">8 R( _& d, \  t2 {- ~</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">4 J8 d4 L% B) E9 I2 D</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">, ~2 z1 F( b) j</span>   k++<br /> <span style=\\"display:none\\">' b6 c  n% l) d! X) O8 s6 t</span> },<font class=\\"jammer\\">; E$ T8 o/ M% J6 z</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">) [  L! c/ M' J. Q\\" e, N. t7 |1 k. V</span> <font class=\\"jammer\\">8 G$ l! W\\" O: u\\" {' z* r\\" b</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">' A- r# p\\" `+ J8 ~- A1 u  f</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">2 {, \( e$ C6 J; H& r\\" N# @; }</span> k=0, while{k<10000,<font class=\\"jammer\\"># x4 E5 N- E  r! Z</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">  b4 M6 Y5 ~/ Q  S5 k</font><br />   k++<font class=\\"jammer\\">- q- u# v+ Q1 t( h+ R</font><br /> },<br /> <span style=\\"display:none\\">( E! ~/ d9 z7 D, Y1 R7 y</span> s;<br /> <span style=\\"display:none\\">9 c% |, z- h' l( S2 N. ]- A</span> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
20000000.
& o) k  q- w; a9 G1 X5.531   //取数据所用时间，秒

－－－－－－－－－－
4 `* e* s3 o* [! o/ O
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];+ a0 g8 M# A7 X/ j# M2 I: y4 e
   
4. mvar:
   
5. 2 O; Q+ r+ v$ S, v) a4 E
6. main0(::a)=$ o\\" w1 l) [\\" m  P; K* F5 Z
   
8. a = linspace(1,1,1000),
   
9. 4 Z) P) |* _% x
10. t0=clock(),
    
11. - _' m3 N! z# k* Y/ z+ D
12. k=0, while{k<10000,
    
13. ! r1 _4 ^\\" T5 d0 p4 @$ ^
14.   i=0, while{i<1000, A(a,i)=2, i++},5 ^- k. v* I1 m
    
16.   k++
    
17. & \/ m* }9 w* |$ L
18. },9 O\\" H1 h' {% ^
    
20. [clock()-t0]/1000;
    
21. : N) R6 l0 C* Q5 D; [+ {$ l
22. 
    
23. / a' ?) T# [2 [( u# U& {, c
24. main1(::a)=) b8 w\\" }( E: e
    
26. t0=clock(), s=0,' T& b8 M; m7 g% A
    
28. k=0, while{k<10000,0 D6 O3 V8 W  I: X+ ~( v/ T
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. + b. |\\" v; y. {6 n
32.   k++
    
33. : P9 |7 P1 o# |' b) f& Q\\" l' G
34. },
    
35. ) y4 V; i( n! B\\" O7 \, M, d
36. s;
    
37. % e1 D! z% K4 {8 {6 u1 H- V) f
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ a0 g8 M# A7 X/ j# M2 I: y4 e</font><br /> mvar:<br /> <span style=\\"display:none\\">2 O; Q+ r+ v$ S, v) a4 E</span> main0(::a)=<font class=\\"jammer\\">$ o\\" w1 l) [\\" m  P; K* F5 Z</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">4 Z) P) |* _% x</span> t0=clock(),<br /> <span style=\\"display:none\\">- _' m3 N! z# k* Y/ z+ D</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">! r1 _4 ^\\" T5 d0 p4 @$ ^</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">5 ^- k. v* I1 m</font><br />   k++<br /> <span style=\\"display:none\\">& \/ m* }9 w* |$ L</span> },<font class=\\"jammer\\">9 O\\" H1 h' {% ^</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">: N) R6 l0 C* Q5 D; [+ {$ l</span> <br /> <span style=\\"display:none\\">/ a' ?) T# [2 [( u# U& {, c</span> main1(::a)=<font class=\\"jammer\\">) b8 w\\" }( E: e</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">' T& b8 M; m7 g% A</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">0 D6 O3 V8 W  I: X+ ~( v/ T</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">+ b. |\\" v; y. {6 n</span>   k++<br /> <span style=\\"display:none\\">: P9 |7 P1 o# |' b) f& Q\\" l' G</span> },<br /> <span style=\\"display:none\\">) y4 V; i( n! B\\" O7 \, M, d</span> s;<br /> <span style=\\"display:none\\">% e1 D! z% K4 {8 {6 u1 H- V) f</span> [clock()-t0]/1000;

结果：
10.469
20000000.
* F6 X" B9 ]6 T% L/ M3.781
, F! ?2 X7 I8 K8 j  D( O' ^
$ g1 n. L: U" {7 R' }9 U－－－－－－－
6 A" z& e; `: m6 D
. J& H+ ?: i; O; h6 ]3 Q$ y7 _在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
9 M/ D+ i0 B7 e! j) X在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
$ R7 L% r/ P# ]  @4 d# ]
% O) M) f9 h, ]+ T# }- @1 o0 Wmatlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

& c9 i% L4 T. m! pForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
. L2 N9 m" s9 L
5、在循环中加点东西
/ t( @- u' K& x5 W* C( d
matlab代码：

1. >> clear all
   0 K* z* A/ }. B0 u7 R( Q1 n5 w& l7 K$ `
2. f=@(x)x+1;' y. U& F. o\" v
   
3. a = linspace(2,2,1000);1 \! M+ \/ y6 |
   
4. tic
   % H\" K* R$ b( z
5. s=0;) k& w* i\" g  D& G/ _+ N4 j  _
   
6. for k=1:1000
   9 _6 o& ?: ~9 G' `
7.   for i=1:1000/ }4 t+ b. \4 r: J\" R: u
   
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));! \$ @+ k5 {( I0 Y% T0 G
   
9.   end+ d; _: T& W: B
   
10. end
    8 l( K7 i* Z! t/ @6 b' N+ I$ K
11. s
    / `\" u1 m) {\" G* {% {* f+ `\" o
12. toc
    2 ^1 u\" z. R8 X* ]
13. 
    5 Y4 U0 G2 K/ _1 T& H: `
14. s =
    0 I\" A0 c7 G0 S3 V
15. ! {3 |; ?( G2 d1 s2 Q& l
    
16.   2.4002e+007
    % q# {6 p4 b' u6 U6 b9 r7 I& t- L
17. : Q- D0 h; C6 F/ |. }$ u8 P4 e) O! F
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. : u# |1 U7 H( k1 i/ J, A/ p
4. mvar:- E7 a# B, \6 G5 n, {. t+ r) V+ ?
   
6. f(x)=x+1;& y) ]3 v, T0 R, A
   
8. main0(::a)=
   
9. 9 h\\" O. A8 @! {' O1 V
10. a = linspace(2,2,1000),
    
11.   V# o' \1 d5 e; B
12. t0=clock(), s=0,% I4 C+ W' H9 R
    
14. k=0, while{k<1000,% S. S5 A; h2 H/ Q4 `2 G
    
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. : }5 k! Y\\" F\\" _4 O- R
18.   k++
    
19. # c$ p1 z: C- H- A\\" w\\" T/ v
20. },
    
21. ) c# O! m; u! L9 `! }! e
22. s;
    
23. 8 o9 x9 r/ f& {2 b7 `$ |* i1 t2 T& ^
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">: u# |1 U7 H( k1 i/ J, A/ p</span> mvar:<font class=\\"jammer\\">- E7 a# B, \6 G5 n, {. t+ r) V+ ?</font><br /> f(x)=x+1;<font class=\\"jammer\\">& y) ]3 v, T0 R, A</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">9 h\\" O. A8 @! {' O1 V</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">  V# o' \1 d5 e; B</span> t0=clock(), s=0,<font class=\\"jammer\\">% I4 C+ W' H9 R</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">% S. S5 A; h2 H/ Q4 `2 G</font><br />   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\\">: }5 k! Y\\" F\\" _4 O- R</span>   k++<br /> <span style=\\"display:none\\"># c$ p1 z: C- H- A\\" w\\" T/ v</span> },<br /> <span style=\\"display:none\\">) c# O! m; u! L9 `! }! e</span> s;<br /> <span style=\\"display:none\\">8 o9 x9 r/ f& {2 b7 `$ |* i1 t2 T& ^</span> [clock()-t0]/1000;

结果：
5 R& J; Y: G; J1 z  C24002271.69189384
5.172

二者效率现在只有一倍差距。
7 D' Q: g, l& Y. T# ~7 H
$ A5 R8 {, A: n0 p- Q. @) i－－－－－－－－－－

不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

1. clear all2 H/ i5 w7 G( q' V5 c# R/ N( s3 H7 v
   
2. f=@(x)x+1;\" l, z* v* B% E2 j3 I, B
   
3. a = linspace(2,2,1000);4 g5 T2 }! R* L4 s! }$ z2 h
   
4. tic9 ], c\" ^\" J. K( ]
   
5. s=0;2 L& U. T+ W9 G# C, N
   
6. for k=1:1000  @6 D. D% k6 x\" ^. y0 t1 t/ }/ `
   
7.   for i=1:1000
   8 o( s, g  O: Y$ _
8.     t=a(i);
   6 p9 n- {& \7 s; |, b
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);
   \" k; g9 a* J- n* _$ |6 K
10.   end
    + W/ G* _0 E! k
11. end: _4 ~4 k  D\" A2 G
    
12. s( `2 |/ D* V( S+ @( X8 D6 m& d- }
    
13. toc/ l- ~, s% A- b, e. H9 H' X
    
14. 
    2 X) H/ B7 |+ z7 n
15. s =7 B& h6 w; I( i; b/ F3 Z3 [3 R
    
16. . _' w+ A0 e% E
    
17.   2.4002e+007+ c- Z9 {! S3 X! f( S* l
    
18. 
    6 ^5 U9 o% T  \; A$ ?0 k* I
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];% @) k' o& S, U
   
4. mvar:# ~! o$ R7 u+ X- \5 Q0 r) s
   
6. f(x)=x+1;/ Y5 n, {9 V8 R; g, ?
   
8. main0(::a)=\\" H: J2 L& q( `: f+ C0 Q
   
10. a = linspace(2,2,1000),. K4 ^7 @# J# J9 h& H
    
12. t0=clock(), s=0,
    
13. # K0 M; K% c' s( A! r
14. k=0, while{k<1000,3 ]; J  A9 f; S2 {
    
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++},' M9 h& J# B  c0 I, x' }
    
18.   k++
    
19. : k, }0 o( U0 q$ H  V4 ^
20. },
    
21.   E0 g* @0 P$ u5 g3 ^% h& C
22. s;# a\\" `) r! X& K' _1 t6 a
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">% @) k' o& S, U</font><br /> mvar:<font class=\\"jammer\\"># ~! o$ R7 u+ X- \5 Q0 r) s</font><br /> f(x)=x+1;<font class=\\"jammer\\">/ Y5 n, {9 V8 R; g, ?</font><br /> main0(::a)=<font class=\\"jammer\\">\\" H: J2 L& q( `: f+ C0 Q</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">. K4 ^7 @# J# J9 h& H</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\"># K0 M; K% c' s( A! r</span> k=0, while{k<1000,<font class=\\"jammer\\">3 ]; J  A9 f; S2 {</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\\">' M9 h& J# B  c0 I, x' }</font><br />   k++<br /> <span style=\\"display:none\\">: k, }0 o( U0 q$ H  V4 ^</span> },<br /> <span style=\\"display:none\\">  E0 g* @0 P$ u5 g3 ^% h& C</span> s;<font class=\\"jammer\\"># a\\" `) r! X& K' _1 t6 a</font><br /> [clock()-t0]/1000;

结果：
. Y9 J( P2 |+ s24002271.69189384
1.625

% `5 X! J' d1 X8 ?matlab效率已不及Forcal。
' T4 }& Y% f4 X
提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
4 E, N7 u4 @, L7 X
" [5 O5 S& E) k5 x0 GForcal代码：

2. !using["sys"];8 O# R/ b+ j0 Y$ D! {- Q: [; i
   
4. mvar:4 F6 h) G$ D- x/ D
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组8 y2 V# T  u5 z, l( d3 I5 c
   
8. t0=clock(),
   
9. % ?: ]0 \  h8 e1 q. @
10. k=0, while{k<10000,
    
11. / i8 r0 }3 k, E: r! d
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值1 K2 `\\" O\\" T5 z
    
14.   k++
    
15. 0 n/ E5 c\\" K/ Z
16. },% O3 P0 \- W5 h; K
    
18. [clock()-t0]/1000;
    
19. , i2 p( B1 t4 n3 S! [' }
20. 
    
21. 4 ^: a: U5 |, Q/ \
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. & v; E$ W2 R9 t% L* ]
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值( X, ^# G0 |( b\\" x( L$ F' j
    
26. t0=clock(),% ]. X) Q& H- ~* D\\" I5 ~
    
28. s=0,9 |4 V4 J9 W\\" I\\" U( V: Y( F; C
    
30. k=0, while{k<10000,
    
31. \\" Q4 f& R/ _/ T3 |& m+ V2 \
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. ! \  _0 d2 X/ H2 x& X1 ?
34.   k++
    
35. + N+ w, q  q' e/ Q1 {; p
36. },
    
37. # I\\" Q5 ~! ?8 R3 I7 r# n. g
38. s;# |) U3 E! M+ U$ Q$ k9 Y
    
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">8 O# R/ b+ j0 Y$ D! {- Q: [; i</font><br /> mvar:<font class=\\"jammer\\">4 F6 h) G$ D- x/ D</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">8 y2 V# T  u5 z, l( d3 I5 c</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">% ?: ]0 \  h8 e1 q. @</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">/ i8 r0 }3 k, E: r! d</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">1 K2 `\\" O\\" T5 z</font><br />   k++<br /> <span style=\\"display:none\\">0 n/ E5 c\\" K/ Z</span> },<font class=\\"jammer\\">% O3 P0 \- W5 h; K</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">, i2 p( B1 t4 n3 S! [' }</span> <br /> <span style=\\"display:none\\">4 ^: a: U5 |, Q/ \</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">& v; E$ W2 R9 t% L* ]</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">( X, ^# G0 |( b\\" x( L$ F' j</font><br /> t0=clock(),<font class=\\"jammer\\">% ]. X) Q& H- ~* D\\" I5 ~</font><br /> s=0,<font class=\\"jammer\\">9 |4 V4 J9 W\\" I\\" U( V: Y( F; C</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">\\" Q4 f& R/ _/ T3 |& m+ V2 \</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">! \  _0 d2 X/ H2 x& X1 ?</span>   k++<br /> <span style=\\"display:none\\">+ N+ w, q  q' e/ Q1 {; p</span> },<br /> <span style=\\"display:none\\"># I\\" Q5 ~! ?8 R3 I7 r# n. g</span> s;<font class=\\"jammer\\"># |) U3 E! M+ U$ Q$ k9 Y</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
$ T( H/ n4 g: i/ h20000000.
: I' T$ \7 u2 k& E" a0.579      //取值用时间，秒