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

forcal

1、简单循环

0 t! \/ Y4 {( }" SMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all+ u4 |' D1 O\" w+ [0 t2 l
   
2. tic
   1 X! a) M( a! L+ {7 f) m3 m
3. s=0;2 x6 A$ X1 f  l' j% ~8 f3 y; w
   
4. for k=1:10000
   - W% }( M; L\" b
5.   for i=1:10008 h) `+ e4 \' f/ L/ F' U, j
   
6.     s=s+1;
     A4 f  ?$ O7 L% c+ O$ B9 W
7.   end
   % O2 _' N8 R) ], s) K! b) ~
8. end
   * y9 ?3 t% ?2 i
9. s- A. h2 V3 @7 j; k4 c* [5 n( D
   
10. toc
    : G1 I. }. H& }\" }; y\" c. V& Q. c
11. 
    ( \0 M6 J) e- g2 \1 ^$ w
12. s =
    4 [- X' B\" }\" n; x
13. , ], p+ _5 X+ |2 k- y& l6 R& s
    
14.     10000000' v+ @! v2 `, c4 o) C
    
15. 
    3 f% K( I4 r8 k' m\" n5 x% \
16. Elapsed time is 0.092808 seconds.

复制代码

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

3 @& W& a( S; YForcal没有JIT加速器，代码：

2. !using["sys"];: _8 t- h% z\\" T4 o; J5 \
   
4. mvar:
   
5. ; E- ?9 b) `& f& R
6. t0=clock(),: J. S+ y! K' J
   
8. s=0,
   
9. $ u) f6 y' _& r$ w! ?# |& O4 }% U, R
10. k=0, while{k<10000,
    
11. 7 c; Q4 `' V* r% u4 E; x. Q
12.   i=0, while{i<1000,s=s+1, i++},: r$ J; z2 G; J5 k\\" D
    
14.   k++& z: D3 T1 {2 w7 V. j/ l
    
16. },
    
17. ) M! V# ]0 j5 e: d2 e# R
18. s;5 m6 Z) }5 H# f1 O# c( U
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">: _8 t- h% z\\" T4 o; J5 \</font><br /> mvar:<br /> <span style=\\"display:none\\">; E- ?9 b) `& f& R</span> t0=clock(),<font class=\\"jammer\\">: J. S+ y! K' J</font><br /> s=0,<br /> <span style=\\"display:none\\">$ u) f6 y' _& r$ w! ?# |& O4 }% U, R</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">7 c; Q4 `' V* r% u4 E; x. Q</span>   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">: r$ J; z2 G; J5 k\\" D</font><br />   k++<font class=\\"jammer\\">& z: D3 T1 {2 w7 V. j/ l</font><br /> },<br /> <span style=\\"display:none\\">) M! V# ]0 j5 e: d2 e# R</span> s;<font class=\\"jammer\\">5 m6 Z) }5 H# f1 O# c( U</font><br /> [clock()-t0]/1000;

结果：
10000000.
- u+ g* D  d. P% j9 d# {0.704  秒

二者运行效率相差7倍多。

$ }/ P& O) V. a' {1 v5 F% K7 C2、循环中稍微加点东西

/ X7 u* }" q0 @% j) BMatlab代码：

1. clear all5 C1 e: Z$ Q6 l% H
   
2. tic# v* T1 U+ W  _
   
3. s=0;. |& f/ M7 m* G9 R/ v$ m
   
4. for k=1:1000
   ! S9 c$ {2 H& |) r! x1 U1 c7 z0 G
5.   for i=1:1000
   * |, ]( y% O0 I/ ~4 j9 D: 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);% i  g- V1 x# _; Z
   
7.   end
   \" I8 k+ e( N0 Y8 E
8. end
   9 l0 T/ _  I5 K# \+ S. z. |
9. s
   $ s: h. y* p' u9 X( f0 L
10. toc' I0 v# O/ D' ~& f! d
    
11. 
    % p8 ~4 q% T' b3 [
12. s =6 x8 E0 `8 c& d# z4 O
    
13. 
    . q& d$ A. M; r8 m
14.   1.1009e+0072 a5 O$ o. W! L6 f! y\" h
    
15. 8 N: i8 U4 ~% U: X0 R3 E% B% g
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 2 J2 i% h) ]\\" l
4. mvar:
   
5. 7 X( O5 q1 m, C
6. t0=clock(),5 k; M7 _& R9 R! y, a  a8 Y- g
   
8. s=0,
   
9. ; u8 T) G% u; |; A. e
10. k=1, while{k<=1000,+ P& e9 B$ }8 d- c
    
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++},3 v9 d1 o9 R5 R4 s8 P# f1 X
    
14.   k++
    
15. , O7 I5 A2 w$ s$ K5 E. f
16. },: }/ N3 k9 A2 J$ `+ f% h
    
18. s;* T/ V$ V) H3 O6 s4 _7 y/ t0 O4 p1 [
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">2 J2 i% h) ]\\" l</span> mvar:<br /> <span style=\\"display:none\\">7 X( O5 q1 m, C</span> t0=clock(),<font class=\\"jammer\\">5 k; M7 _& R9 R! y, a  a8 Y- g</font><br /> s=0,<br /> <span style=\\"display:none\\">; u8 T) G% u; |; A. e</span> k=1, while{k<=1000,<font class=\\"jammer\\">+ P& e9 B$ }8 d- c</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\\">3 v9 d1 o9 R5 R4 s8 P# f1 X</font><br />   k++<br /> <span style=\\"display:none\\">, O7 I5 A2 w$ s$ K5 E. f</span> },<font class=\\"jammer\\">: }/ N3 k9 A2 J$ `+ f% h</font><br /> s;<font class=\\"jammer\\">* T/ V$ V) H3 O6 s4 _7 y/ t0 O4 p1 [</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
2 }5 g( A0 q) e, A4 y$ Z3 j1.031

二者效率相差已比较小。
, ^& w5 t9 _: p( _) `. o: E$ d( C
3、再增加一个函数调用

# \* R8 h7 N5 k- t* P, pMatlab代码：

1. clear all4 w0 t0 W' ?3 n5 c/ }! w\" X
   
2. f=@(x)x+1;
   3 l3 I2 v1 z, O8 `. `% I
3. tic, l$ t5 z( c$ G. d# i2 w
   
4. s=0;& f9 D# y\" F% l0 E7 S) c
   
5. for k=1:1000
   ( w, Z3 d' Y6 t- [5 q# o6 n
6.   for i=1:10000 z6 y2 b: j4 I4 }/ e
   
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  P0 ^6 O! M  m5 y$ d- j7 o2 S
   
8.   end
   ) S\" I8 z+ U, n; C
9. end
   ) r' n: i9 i1 f4 O\" x
10. s4 U( l) N) H# Y8 f
    
11. toc6 N9 N, m5 R3 j. ^: y6 h
    
12. . Z9 ~; H5 x, p( P
    
13. s =# b3 B2 N9 \/ ?
    
14. , T+ g4 ~7 v\" c8 h( j# S$ b
    
15.   5.1251e+008
    * H, ?+ a8 T9 j% O$ h7 d\" z
16. ' L$ d+ ~9 r2 w9 [! G2 U% @, B; E2 g+ K
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];, J6 t& R8 E% L5 \4 N
   
4. f(x)=x+1;
   
5. + L1 T, w5 _6 E1 b, h9 y
6. mvar:
   
7.   t8 n  t3 _5 j6 A6 P$ h
8. t0=clock(),
   
9. 5 P5 w7 b) U\\" I% a6 i) c; M' I. D( s
10. s=0,4 J1 {% z\\" G1 [9 A& W
    
12. k=1, while{k<=1000,, p. X- U7 l  @( A. I( 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++},. D$ x9 {, ^- S$ P$ A
    
16.   k++2 \9 f; j9 d4 z, k
    
18. },
    
19.   z$ `( D7 `* \\\" C5 q
20. s;
    
21. 2 V9 I8 w* i; _( ^- r) S! U, g1 f$ o
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">, J6 t& R8 E% L5 \4 N</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">+ L1 T, w5 _6 E1 b, h9 y</span> mvar:<br /> <span style=\\"display:none\\">  t8 n  t3 _5 j6 A6 P$ h</span> t0=clock(),<br /> <span style=\\"display:none\\">5 P5 w7 b) U\\" I% a6 i) c; M' I. D( s</span> s=0,<font class=\\"jammer\\">4 J1 {% z\\" G1 [9 A& W</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">, p. X- U7 l  @( A. I( v</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++},<font class=\\"jammer\\">. D$ x9 {, ^- S$ P$ A</font><br />   k++<font class=\\"jammer\\">2 \9 f; j9 d4 z, k</font><br /> },<br /> <span style=\\"display:none\\">  z$ `( D7 `* \\\" C5 q</span> s;<br /> <span style=\\"display:none\\">2 V9 I8 w* i; _( ^- r) S! U, g1 f$ o</span> [clock()-t0]/1000;

结果：
  P+ ^, C- P* V- _2 {! J8 K512508953.6658993
. d0 F( x! U5 S/ _; I1.078

可以看出，Matlab效率已不及Forcal。
1 z8 K/ F/ f6 l9 u; k$ Q故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

% C" ^8 M2 M) A  bMatlab代码：

1. >> clear all$ U! m6 d; ?1 ~0 ?# b/ r\" O0 i
   
2. a = linspace(1,1,1000);# ~0 `. v/ C5 F2 o1 |# k
   
3. tic8 P  s8 K, P* ?\" \. s: T
   
4. for k=1:10000
   5 v7 E+ t1 U  k! T9 f0 `) {! ?
5.   for i=1:1000\" c9 g* _6 n! @/ v, a
   
6.     a(i)=2;
   7 R& b1 A7 h\" X+ u3 X0 o: ?
7.   end
   , i1 M. _: w6 F- z  d* ~4 @
8. end
   - b, q& V4 h$ r! ~4 [. O
9. toc. j1 b+ V; m# a  ~+ D# B
   
10. 2 L- u3 z; m- F. k1 W7 L3 j
    
11. tic
    ) u% M5 T0 k0 ]. m7 n7 d3 E
12. s=0;\" n+ }- U5 Z5 V  h2 l\" o/ h& o
    
13. for k=1:10000) t5 f' d+ \( ]4 g4 w: T+ P) [
    
14.   for i=1:1000; t+ N1 R  q1 C+ S
    
15.     s=s+a(i);7 B* }2 ]! |) D/ Y2 x# a
    
16.   end
    / G. U2 x3 O7 S- L; l
17. end
    , k! U5 x2 n, ^
18. s
    + M( Y; v1 {- W\" r3 H% P! S6 K; f
19. toc
    + s' E+ b- `& l/ l# L1 i& Z7 p
20. Elapsed time is 0.140658 seconds.
    $ }+ R! L2 ]3 _7 `
21. 
    + z# x3 i: _+ X( R! @% X  Z
22. s =
    ' x% ]. ]# h1 q- ?\" F8 t
23. 
    ! m3 M0 ?/ H3 C+ U  K
24.     20000000
    2 h% F$ ^/ `( A2 N/ @2 e# ~9 E
25. 8 k. R/ {6 E) _' T; Q! r9 n5 W: Y
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

) |$ S! O* x" f( {% JForcal代码：

2. !using["math","sys"];. [\\" U7 H7 Z. G7 O9 [! ~  T! O
   
4. mvar:- r\\" ]! R0 j7 a! \
   
6. main0(::a)=. R; G  b4 P\\" X0 \2 H' K
   
8. a = linspace(1,1,1000),
   
9. 6 U; p3 K# v1 j
10. t0=clock(),
    
11. 9 m# X$ P. Z# x- w# H( A: V& z0 w
12. k=0, while{k<10000,
    
13. 7 f  O! S0 @8 u0 N. u
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. ! S- C8 x: A' l! z
16.   k++
    
17. : F, [4 ~: r! y1 w1 B% U! P
18. },
    
19. 0 x0 |: z4 C3 g( C/ H
20. [clock()-t0]/1000;' |: Z( T* S5 o3 A- M# K- B
    
22. 2 s8 y1 n6 a7 F( w9 ~
    
24. main1(::a)=
    
25. ; d7 K4 [$ h7 m; o6 w5 t% `
26. t0=clock(), s=0,
    
27. 3 U& s# B' a+ u* S, r
28. k=0, while{k<10000,
    
29.   z( {. p- i& C8 g- [* H
30.   i=0, while{i<1000, s=s+a(i), i++},* C+ q: C! D$ F
    
32.   k+++ m1 u& q- j* |$ V$ o
    
34. },
    
35. ! _% N. D( v6 T0 l1 w+ o
36. s;
    
37. ' ^% O' X/ P# q/ L) J; w
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">. [\\" U7 H7 Z. G7 O9 [! ~  T! O</font><br /> mvar:<font class=\\"jammer\\">- r\\" ]! R0 j7 a! \</font><br /> main0(::a)=<font class=\\"jammer\\">. R; G  b4 P\\" X0 \2 H' K</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">6 U; p3 K# v1 j</span> t0=clock(),<br /> <span style=\\"display:none\\">9 m# X$ P. Z# x- w# H( A: V& z0 w</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">7 f  O! S0 @8 u0 N. u</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">! S- C8 x: A' l! z</span>   k++<br /> <span style=\\"display:none\\">: F, [4 ~: r! y1 w1 B% U! P</span> },<br /> <span style=\\"display:none\\">0 x0 |: z4 C3 g( C/ H</span> [clock()-t0]/1000;<font class=\\"jammer\\">' |: Z( T* S5 o3 A- M# K- B</font><br /> <font class=\\"jammer\\">2 s8 y1 n6 a7 F( w9 ~</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">; d7 K4 [$ h7 m; o6 w5 t% `</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">3 U& s# B' a+ u* S, r</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">  z( {. p- i& C8 g- [* H</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">* C+ q: C! D$ F</font><br />   k++<font class=\\"jammer\\">+ m1 u& q- j* |$ V$ o</font><br /> },<br /> <span style=\\"display:none\\">! _% N. D( v6 T0 l1 w+ o</span> s;<br /> <span style=\\"display:none\\">' ^% O' X/ P# q/ L) J; w</span> [clock()-t0]/1000;

结果：
( K7 `2 h- K- X2 ^! Q15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒

6 S2 B9 r; h  c6 d/ d4 ?' U/ i－－－－－－－－－－
( Z- T% o# l. m$ O9 _8 y* x4 ]/ ~
; \4 X% H7 O. L  X; V( S用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. - R+ V  E: c6 w/ n\\" P' s2 A
4. mvar:+ L7 Z) T0 J. g4 D
   
6. main0(::a)=
   
7. ; P9 @; Z% x: i# L6 G+ u7 }3 A0 `
8. a = linspace(1,1,1000),  e( x$ {: c/ c8 z, q2 r6 I; k5 D! I
   
10. t0=clock(),, d$ D8 h3 B5 g& c+ J. `& |; {
    
12. k=0, while{k<10000,3 @8 i6 [* l3 G1 N
    
14.   i=0, while{i<1000, A(a,i)=2, i++},# Y3 o\\" {8 ^- y7 d$ A
    
16.   k++
    
17. + r2 d+ |' p9 e7 `/ ?( L9 z. b
18. },
    
19. & q( [0 L2 s5 j- A\\" `: n/ x, q8 |
20. [clock()-t0]/1000;2 F: E8 ^2 R' I' t: @+ D
    
22. : n3 y# g/ ~+ m7 d+ ]$ T! |* ~
    
24. main1(::a)=% T0 {/ V/ V- i1 `# D$ s; z; A% w  k
    
26. t0=clock(), s=0,$ P+ \/ Z% ^5 ?! v9 W* L0 B  e# v
    
28. k=0, while{k<10000,. {; Y+ `  n+ C
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},( e2 _. n# w: R) e6 y. O& ^7 @
    
32.   k++
    
33. / v# w' Y; U% u& l! O5 \
34. },/ a6 \9 a7 \% k, D
    
36. s;
    
37. ( h8 U4 _# N; ?* v- l
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">- R+ V  E: c6 w/ n\\" P' s2 A</span> mvar:<font class=\\"jammer\\">+ L7 Z) T0 J. g4 D</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">; P9 @; Z% x: i# L6 G+ u7 }3 A0 `</span> a = linspace(1,1,1000),<font class=\\"jammer\\">  e( x$ {: c/ c8 z, q2 r6 I; k5 D! I</font><br /> t0=clock(),<font class=\\"jammer\\">, d$ D8 h3 B5 g& c+ J. `& |; {</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">3 @8 i6 [* l3 G1 N</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\"># Y3 o\\" {8 ^- y7 d$ A</font><br />   k++<br /> <span style=\\"display:none\\">+ r2 d+ |' p9 e7 `/ ?( L9 z. b</span> },<br /> <span style=\\"display:none\\">& q( [0 L2 s5 j- A\\" `: n/ x, q8 |</span> [clock()-t0]/1000;<font class=\\"jammer\\">2 F: E8 ^2 R' I' t: @+ D</font><br /> <font class=\\"jammer\\">: n3 y# g/ ~+ m7 d+ ]$ T! |* ~</font><br /> main1(::a)=<font class=\\"jammer\\">% T0 {/ V/ V- i1 `# D$ s; z; A% w  k</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">$ P+ \/ Z% ^5 ?! v9 W* L0 B  e# v</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">. {; Y+ `  n+ C</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">( e2 _. n# w: R) e6 y. O& ^7 @</font><br />   k++<br /> <span style=\\"display:none\\">/ v# w' Y; U% u& l! O5 \</span> },<font class=\\"jammer\\">/ a6 \9 a7 \% k, D</font><br /> s;<br /> <span style=\\"display:none\\">( h8 U4 _# N; ?* v- l</span> [clock()-t0]/1000;

结果：
/ E& j. a5 E: f9 g/ L10.469
20000000.
5 @! v( P* e5 c3.781
/ J7 _3 Y6 K( K: v7 L+ \
0 Q! j( b$ q- C7 [－－－－－－－

5 U! y; T* t3 ?3 o" E# ]在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
8 U+ U8 N' d/ p! Q$ @在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

; P! f# L7 ^& E" {3 @5 D5、在循环中加点东西
( v1 o: j7 [+ E  K3 u
matlab代码：

1. >> clear all1 r5 N5 g+ O& b1 t9 {
   
2. f=@(x)x+1;, T: A9 h3 a  |5 B$ R
   
3. a = linspace(2,2,1000);6 K! c\" b# k$ X9 D- N, Y' v
   
4. tic
   1 }- b3 M2 e\" C$ U0 e0 S
5. s=0;
   5 o$ N7 W; l2 B( R% g2 m) P
6. for k=1:1000, ?0 n2 K9 w7 H# ^
   
7.   for i=1:1000
   3 \/ i0 B5 T4 Q
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));\" n  t+ ?9 H8 u5 ?- @
   
9.   end( e  \% A+ @# {$ v% o' o\" H7 V
   
10. end
    * D5 z/ G* ?2 @- n( z
11. s- t+ a. o* T0 B: K5 k( s
    
12. toc0 ]! J8 g5 O1 e; A5 t
    
13. 
    : c* Q) u) z+ M/ z0 d
14. s =0 ~+ S, c. Q$ p
    
15. 9 O% R. p& Z  X) }! o
    
16.   2.4002e+007\" g\" B. X6 M! p$ G1 @
    
17. % \/ e& H- q/ e' ]! W* X
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. ! |2 L7 T3 r$ l4 S: V& a4 U
4. mvar:) C. v2 }$ D* I3 I% K
   
6. f(x)=x+1;\\" t: g8 {; @9 c8 n* Y6 o! E
   
8. main0(::a)=
   
9. ! {; x. ?$ m, y0 }  \
10. a = linspace(2,2,1000),
    
11. ; |( L) D/ y3 Q\\" [4 c
12. t0=clock(), s=0,
    
13. * \\\" `0 R5 M0 v$ Q8 p
14. k=0, while{k<1000,6 U8 z% J\\" v4 \* F5 A
    
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 ^- d7 k6 S  I7 C6 J( A
18.   k++
    
19. 4 q\\" H$ N: r6 D1 A& _
20. },
    
21. 6 y+ v2 E' u' z9 a
22. s;
    
23. ; |/ a* f1 ]( M- d+ `
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">! |2 L7 T3 r$ l4 S: V& a4 U</span> mvar:<font class=\\"jammer\\">) C. v2 }$ D* I3 I% K</font><br /> f(x)=x+1;<font class=\\"jammer\\">\\" t: g8 {; @9 c8 n* Y6 o! E</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">! {; x. ?$ m, y0 }  \</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">; |( L) D/ y3 Q\\" [4 c</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">* \\\" `0 R5 M0 v$ Q8 p</span> k=0, while{k<1000,<font class=\\"jammer\\">6 U8 z% J\\" v4 \* F5 A</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\\">9 ^- d7 k6 S  I7 C6 J( A</span>   k++<br /> <span style=\\"display:none\\">4 q\\" H$ N: r6 D1 A& _</span> },<br /> <span style=\\"display:none\\">6 y+ v2 E' u' z9 a</span> s;<br /> <span style=\\"display:none\\">; |/ a* f1 ]( M- d+ `</span> [clock()-t0]/1000;

结果：
3 l% R. g$ ^+ ~! M9 [9 c24002271.69189384
8 G+ X, K" [6 H  t  M" A5.172
# B: X: N# \4 W3 M$ y0 F: [
4 o8 `$ [. ]) j4 e# _( {$ T二者效率现在只有一倍差距。
2 l/ M( }2 U3 B- c" t
) w0 Z7 r- N7 E( ^- A－－－－－－－－－－

3 d/ ~9 F: E1 ]0 w% x) P不过，显然matlab和Forcal的程序应该这样写：
. U/ e! I8 T% ?% }& H# Q
2 A9 r) o$ p7 ?' g( Rmatlab程序：

1. clear all1 U) X1 I! X/ S
   
2. f=@(x)x+1;6 d- I% g\" Z, d# i0 e+ a! {# Q
   
3. a = linspace(2,2,1000);
   + o/ q0 F\" J2 F0 v3 P
4. tic
   $ N+ n* [# }$ B7 T* S5 ~
5. s=0;  L6 F- D* \/ |6 |
   
6. for k=1:10005 }( }5 u, j, X7 U  S7 {
   
7.   for i=1:1000
   / h9 q5 f) q& Q; \  [3 [' @
8.     t=a(i);. n1 ?( C& U9 O: X\" T
   
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);, ]. r; a3 N$ f) D3 ^
   
10.   end$ p' [/ Y  ^! v4 B- D/ }
    
11. end
    ) ?4 r% z/ Q& d  F8 O! ^
12. s( U& g5 E5 F2 t0 |, O1 r. w! i  R
    
13. toc. S/ X7 W  D+ W) Z; h
    
14. 
    ( f) J) P( a3 F0 M3 T( w, E
15. s =
    9 G5 S# N, R\" Y# A/ |) ^
16. $ ]- N1 p* L' g$ h+ f+ _6 g
    
17.   2.4002e+0077 Z5 _% P1 J$ H) l0 }/ Q
    
18. ; l* F+ h8 s! Y% o) m  n* ^
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];2 i; s: d\\" I. T3 U. j, D6 m& G- I
   
4. mvar:# g3 p: K- m$ H0 Q
   
6. f(x)=x+1;; D$ j7 n8 _$ w\\" ]3 i% d) j
   
8. main0(::a)=
   
9. , H% h$ B/ j. q3 S, o# e
10. a = linspace(2,2,1000),/ v\\" M; Z8 f5 ~# h\\" l
    
12. t0=clock(), s=0,: W\\" q7 ^' `  [2 X
    
14. k=0, while{k<1000,, N! k. i5 C) c0 n& @1 R
    
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. 7 @3 ?+ h\\" Z3 M. R) M  C
18.   k++$ D  S8 D3 G\\" [5 D* d  e4 Z1 w. k
    
20. },8 z3 B7 f% u% |8 ^( \
    
22. s;
    
23. ; j' ?$ @# N- U( K# A
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">2 i; s: d\\" I. T3 U. j, D6 m& G- I</font><br /> mvar:<font class=\\"jammer\\"># g3 p: K- m$ H0 Q</font><br /> f(x)=x+1;<font class=\\"jammer\\">; D$ j7 n8 _$ w\\" ]3 i% d) j</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">, H% h$ B/ j. q3 S, o# e</span> a = linspace(2,2,1000),<font class=\\"jammer\\">/ v\\" M; Z8 f5 ~# h\\" l</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">: W\\" q7 ^' `  [2 X</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">, N! k. i5 C) c0 n& @1 R</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++},<br /> <span style=\\"display:none\\">7 @3 ?+ h\\" Z3 M. R) M  C</span>   k++<font class=\\"jammer\\">$ D  S8 D3 G\\" [5 D* d  e4 Z1 w. k</font><br /> },<font class=\\"jammer\\">8 z3 B7 f% u% |8 ^( \</font><br /> s;<br /> <span style=\\"display:none\\">; j' ?$ @# N- U( K# A</span> [clock()-t0]/1000;

结果：
5 {: Q1 B0 p* N# }; Q24002271.69189384
! ~( j( ]4 x4 R8 ]+ m4 G. q8 `, n1.625
8 o/ {5 `3 U& D" J9 }6 T# ?
2 N; i# z4 ^7 X3 xmatlab效率已不及Forcal。

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

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
6 W; w, S; h! i, a. g
Forcal代码：

2. !using["sys"];' c5 t' e) a0 ~! r5 z
   
4. mvar:
   
5. ; A% k- p, Q; \0 ^* V
6. a="\&\[1000]"/4,    //用字符串模拟静态数组) T2 b( @0 v) h+ H
   
8. t0=clock(),
   
9.   B6 {3 w  S& ?% W+ H5 X+ V: U' @\\" z
10. k=0, while{k<10000,
    
11. 6 }/ w\\" Y* Q! k: ~
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. $ E' _+ q\\" J3 i7 e2 j8 {( f
14.   k++& p& o& B& u/ s# j5 `
    
16. },
    
17. + H) {  H' U1 O2 R8 {
18. [clock()-t0]/1000;7 R) f7 }+ Z2 J8 S9 y
    
20. 
    
21. # \* u! ]. f. t& H4 \1 F5 t8 q# m  r
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. $ L* I4 W4 H. k
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值* e( i1 W8 |$ n
    
26. t0=clock(),
    
27. ) N+ M3 i8 j+ |6 X; J8 n% U
28. s=0,
    
29. - j1 |) h: ?3 q) N
30. k=0, while{k<10000,
    
31. . i$ Y- q3 e+ z/ |) ^
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值, _* Q. K# L( l' l* ]% i
    
34.   k++/ C6 \' _; K6 y# }  _# F% g
    
36. },
    
37. ) f& D. k, b, \  z9 j% ^
38. s;
    
39. 0 [# O9 ?3 Q0 e) E- {
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">' c5 t' e) a0 ~! r5 z</font><br /> mvar:<br /> <span style=\\"display:none\\">; A% k- p, Q; \0 ^* V</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">) T2 b( @0 v) h+ H</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">  B6 {3 w  S& ?% W+ H5 X+ V: U' @\\" z</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">6 }/ w\\" Y* Q! k: ~</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">$ E' _+ q\\" J3 i7 e2 j8 {( f</span>   k++<font class=\\"jammer\\">& p& o& B& u/ s# j5 `</font><br /> },<br /> <span style=\\"display:none\\">+ H) {  H' U1 O2 R8 {</span> [clock()-t0]/1000;<font class=\\"jammer\\">7 R) f7 }+ Z2 J8 S9 y</font><br /> <br /> <span style=\\"display:none\\"># \* u! ]. f. t& H4 \1 F5 t8 q# m  r</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">$ L* I4 W4 H. k</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">* e( i1 W8 |$ n</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">) N+ M3 i8 j+ |6 X; J8 n% U</span> s=0,<br /> <span style=\\"display:none\\">- j1 |) h: ?3 q) N</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">. i$ Y- q3 e+ z/ |) ^</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">, _* Q. K# L( l' l* ]% i</font><br />   k++<font class=\\"jammer\\">/ C6 \' _; K6 y# }  _# F% g</font><br /> },<br /> <span style=\\"display:none\\">) f& D. k, b, \  z9 j% ^</span> s;<br /> <span style=\\"display:none\\">0 [# O9 ?3 Q0 e) E- {</span> [clock()-t0]/1000;

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