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

forcal

1、简单循环
( _3 p3 e4 m/ d  t+ ?: J
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all! {* u. c8 r  L: Z1 K( w4 n
   
2. tic1 r9 x: J/ m! T0 f
   
3. s=0;
   , c3 a$ T( t) {, M. G
4. for k=1:100009 `; M% `& X  C$ i\" l) d3 d
   
5.   for i=1:1000
   6 h5 M\" X0 T$ G* a3 a1 y\" W
6.     s=s+1;
   6 L3 s/ Z; r: j: l
7.   end6 p5 \* w# M$ M
   
8. end
   / M$ z, ?, Z. ]: u3 z9 x' m
9. s
   * |& k9 ?5 D) I
10. toc& n7 m! W+ ~# K! U5 k; z
    
11. $ @! k) A) V- l. b) _% p+ R2 a% ~7 x
    
12. s =% r% \7 M* a/ ]6 x
    
13. 7 B/ O\" h5 [# v\" @
    
14.     10000000
    7 j, C# m7 ~. W) ~& Q
15. & Y7 O  `* I$ q6 d  z! j
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

  `7 v- T3 i  y% yForcal没有JIT加速器，代码：

2. !using["sys"];
   
3. 7 U& _2 l2 H: Z2 Y. X
4. mvar:
   
5. ! L; n0 |! _9 F0 n5 h5 u% i
6. t0=clock(),- ?/ g: K% Q+ X  R& y
   
8. s=0,
   
9. 8 }9 U! R$ t- q5 L  C8 p4 ?8 O
10. k=0, while{k<10000,\\" j* d8 ?! V* B9 R0 x
    
12.   i=0, while{i<1000,s=s+1, i++},: u) b( e' o, c
    
14.   k++
    
15. % W3 s/ @4 A' |( ?$ m- [
16. },
    
17.   }& C) k+ E  h9 J
18. s;
    
19. 7 t\\" G9 A+ x% n8 S5 c8 c7 x- Q
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">7 U& _2 l2 H: Z2 Y. X</span> mvar:<br /> <span style=\\"display:none\\">! L; n0 |! _9 F0 n5 h5 u% i</span> t0=clock(),<font class=\\"jammer\\">- ?/ g: K% Q+ X  R& y</font><br /> s=0,<br /> <span style=\\"display:none\\">8 }9 U! R$ t- q5 L  C8 p4 ?8 O</span> k=0, while{k<10000,<font class=\\"jammer\\">\\" j* d8 ?! V* B9 R0 x</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">: u) b( e' o, c</font><br />   k++<br /> <span style=\\"display:none\\">% W3 s/ @4 A' |( ?$ m- [</span> },<br /> <span style=\\"display:none\\">  }& C) k+ E  h9 J</span> s;<br /> <span style=\\"display:none\\">7 t\\" G9 A+ x% n8 S5 c8 c7 x- Q</span> [clock()-t0]/1000;

结果：
10000000.
0.704  秒
& |/ G# p' m: ]" `, X+ u' f% k9 v
1 q, Q2 q) f8 Q$ V" S0 `二者运行效率相差7倍多。
9 }; u7 l% Q+ f' Q1 w
2、循环中稍微加点东西
5 d1 Z! B  F8 R8 H. J" P  c, j0 D
Matlab代码：

1. clear all4 o8 c- A- ~6 r) M( i( b. a1 {
   
2. tic; r1 N& f\" b% m
   
3. s=0;
   4 ^. s% U2 g0 S: t) Q- f4 Q* i
4. for k=1:1000
   ; t4 a5 n% ^4 l2 u6 A, Z. u7 R
5.   for i=1:1000
   1 ?  E/ W4 ]$ J
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);
   , e3 ]) Z& l, P7 i- R& Q
7.   end
   ! q' j/ S/ i* ^9 L. P+ Y, u
8. end3 b8 q% ^# q7 X$ @* g
   
9. s/ N! y; g: n  H' b
   
10. toc9 q4 [5 A) @. z  |. ?# S
    
11. 
    6 }; u% o* k3 m4 b% c! X0 U
12. s =
    + B/ z. p. a- B  k2 D* I
13. \" ]2 f& U: L2 d* [: l* y/ d
    
14.   1.1009e+007
    ; H; ^) P2 K( c/ G
15. 
    ' l7 v\" S5 l8 x9 _9 F! N
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];' s+ R, G/ I* ?4 n1 D+ Q: j
   
4. mvar:! s9 I* p, ?3 x6 i# p\\" K5 Z$ Y
   
6. t0=clock(),# Q, X( l9 K8 L; H
   
8. s=0,
   
9. \\" ^3 ^0 k, t, v2 ^: N
10. k=1, while{k<=1000,
    
11. 3 g6 g) w: s# b
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++},9 k! J: F$ Z\\" u
    
14.   k++
    
15. . }2 z. M/ k/ A- A. B
16. },
    
17. 4 z% [! n\\" A4 I3 W& Q
18. s;+ O' I5 F  A5 F
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">' s+ R, G/ I* ?4 n1 D+ Q: j</font><br /> mvar:<font class=\\"jammer\\">! s9 I* p, ?3 x6 i# p\\" K5 Z$ Y</font><br /> t0=clock(),<font class=\\"jammer\\"># Q, X( l9 K8 L; H</font><br /> s=0,<br /> <span style=\\"display:none\\">\\" ^3 ^0 k, t, v2 ^: N</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">3 g6 g) w: s# b</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\\">9 k! J: F$ Z\\" u</font><br />   k++<br /> <span style=\\"display:none\\">. }2 z. M/ k/ A- A. B</span> },<br /> <span style=\\"display:none\\">4 z% [! n\\" A4 I3 W& Q</span> s;<font class=\\"jammer\\">+ O' I5 F  A5 F</font><br /> [clock()-t0]/1000;

结果：
. p0 V  D8 v' u9 F# q% b11008953.66602346
& u# ]+ i* G2 k/ m5 B/ b1.031

- g/ ?! S7 \- R4 m% i# R- }二者效率相差已比较小。

3、再增加一个函数调用

Matlab代码：

1. clear all4 ^  c1 T8 K# \- T# n
   
2. f=@(x)x+1;, K& q: G4 a+ S) D
   
3. tic1 S0 k2 u# T' s\" B3 H\" [9 `
   
4. s=0;% x4 ?4 @; l2 w/ L) }+ D\" o
   
5. for k=1:1000
   1 E5 Q\" |) A& u, b( ?9 k8 b
6.   for i=1:1000
   , f+ b* Q! d( H  q5 \8 Y* L- b
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);' ]: f& ~: @* g+ I! s$ y! q
   
8.   end
   , G9 H5 o( i/ m0 j* F* X\" X/ {1 l& y
9. end
   7 @& K0 p3 P9 G9 t
10. s
    1 t0 g8 q( _7 u5 L. S
11. toc1 s8 |3 q/ @9 Y) z+ u3 L! n
    
12. - y* M1 K0 w( t2 G5 h$ k  j
    
13. s =4 Q! z! G9 N7 R* Q) L( w. P0 l3 c: M
    
14. ; r7 }/ D: r: t3 o1 J; A1 ~
    
15.   5.1251e+008) W7 R: b: c8 W, w: ^, M* [
    
16.   n* J9 ]8 p% Q6 N: D- c9 X
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. / D) d! Y9 k- L) [0 s; q, k
4. f(x)=x+1;* D6 t1 o/ I$ v: ?
   
6. mvar:
   
7.   t. g2 D8 _( u! h
8. t0=clock(),* g5 {8 j, x9 l6 ^/ \# B
   
10. s=0,6 L6 B0 A! ~% [0 Y; ?4 v
    
12. k=1, while{k<=1000,
    
13. - [$ A& \: g1 S4 \
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++},# `/ ~* p7 M& F% B. g
    
16.   k++
    
17. 3 D4 P  S/ k3 R! N5 p
18. },
    
19. ; o# d5 j: u; g: I
20. s;  {) ?- @. w6 E8 @' X
    
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">/ D) d! Y9 k- L) [0 s; q, k</span> f(x)=x+1;<font class=\\"jammer\\">* D6 t1 o/ I$ v: ?</font><br /> mvar:<br /> <span style=\\"display:none\\">  t. g2 D8 _( u! h</span> t0=clock(),<font class=\\"jammer\\">* g5 {8 j, x9 l6 ^/ \# B</font><br /> s=0,<font class=\\"jammer\\">6 L6 B0 A! ~% [0 Y; ?4 v</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">- [$ A& \: g1 S4 \</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\\"># `/ ~* p7 M& F% B. g</font><br />   k++<br /> <span style=\\"display:none\\">3 D4 P  S/ k3 R! N5 p</span> },<br /> <span style=\\"display:none\\">; o# d5 j: u; g: I</span> s;<font class=\\"jammer\\">  {) ?- @. w6 E8 @' X</font><br /> [clock()-t0]/1000;

结果：
512508953.6658993
1.078
9 W. g" i' ^1 S! p- O$ i5 V! |
可以看出，Matlab效率已不及Forcal。
% k5 v  a8 g, ~& b" O6 [2 P故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

Matlab代码：

1. >> clear all
   - l( [1 K) f% Z0 f
2. a = linspace(1,1,1000);  H7 p# m+ n& v5 W0 e* }$ `$ Q
   
3. tic
   , f5 _' B6 o& u1 }
4. for k=1:100007 W9 ?\" s3 P. ^, L8 F\" g
   
5.   for i=1:1000* d+ M+ ^  h8 X# n( i+ G
   
6.     a(i)=2;6 v! }7 T9 w  `6 J' P
   
7.   end\" F# S0 \8 H  |! Z
   
8. end
   ( H' F0 l0 S\" G\" @( P' L
9. toc
   - M* W# x: y; W4 m: u$ Z, E
10. ' k, W$ n% d& Y
    
11. tic
    ) v( P4 W+ c  q& D- Q- X( ^  W
12. s=0;0 I7 \+ E3 B3 D. D( [% d9 r4 o
    
13. for k=1:10000
    # ^: I0 S: B# h* E7 E9 j
14.   for i=1:1000! E% l0 }1 E) p\" n8 \7 Q) z
    
15.     s=s+a(i);
    % r\" t- r$ w5 b9 v9 a, o+ u
16.   end# I3 b) f, y9 Q8 K
    
17. end( h$ s\" L, u( N7 V$ U1 u
    
18. s+ G9 f9 {5 G6 u* M
    
19. toc7 |% |$ [# y8 C  @5 |
    
20. Elapsed time is 0.140658 seconds.
    $ [  F6 N3 N/ i3 |. C: V1 N
21. 
    ; m\" A  S6 x. |+ x9 v5 _
22. s =
    : L+ ?% W7 o* p
23. 6 ?9 ^1 s8 U3 c
    
24.     20000000
    5 X9 f' \- q6 \; E& `9 M5 W
25. 
    - V& O, n) t2 q/ Z# e3 M
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
# O7 f' {2 u1 e/ U1 a
. L) _% b. v2 C% lForcal代码：

2. !using["math","sys"];
   
3. / v# I0 @* U6 _
4. mvar:  r6 w) E5 z\\" ]: l; O, s8 a
   
6. main0(::a)=
   
7. ( _' x$ S+ F7 v. v0 F7 x
8. a = linspace(1,1,1000),2 _\\" z\\" |! @& Y# W* E( t9 B
   
10. t0=clock(),
    
11. / l& Y8 K  F3 J, i& _& z5 o
12. k=0, while{k<10000,$ h: i5 }  M( `) O
    
14.   i=0, while{i<1000, a(i)=2, i++},2 X$ }' t( n- b$ @) v+ P\\" p
    
16.   k++) r  K) `5 Z\\" {( N5 {
    
18. },; B% `\\" u7 _3 u+ @
    
20. [clock()-t0]/1000;
    
21. + M8 P: j7 r. D. x: b
22. ' G1 R* @  N, Y( Z6 o0 p, O  f
    
24. main1(::a)=% Z# ~9 K3 U2 x4 `: L4 Y
    
26. t0=clock(), s=0,2 y3 B\\" N+ d9 ^5 T' z
    
28. k=0, while{k<10000,9 u6 t& Q& Z7 P
    
30.   i=0, while{i<1000, s=s+a(i), i++},# Y, H2 s/ h& A9 J) p
    
32.   k++\\" e9 L/ A# I+ A0 l! q; s4 V
    
34. },
    
35. # P* q) j% M' v
36. s;- Z# x' A\\" [7 D* Z) d. |! v
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">/ v# I0 @* U6 _</span> mvar:<font class=\\"jammer\\">  r6 w) E5 z\\" ]: l; O, s8 a</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">( _' x$ S+ F7 v. v0 F7 x</span> a = linspace(1,1,1000),<font class=\\"jammer\\">2 _\\" z\\" |! @& Y# W* E( t9 B</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">/ l& Y8 K  F3 J, i& _& z5 o</span> k=0, while{k<10000,<font class=\\"jammer\\">$ h: i5 }  M( `) O</font><br />   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">2 X$ }' t( n- b$ @) v+ P\\" p</font><br />   k++<font class=\\"jammer\\">) r  K) `5 Z\\" {( N5 {</font><br /> },<font class=\\"jammer\\">; B% `\\" u7 _3 u+ @</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">+ M8 P: j7 r. D. x: b</span> <font class=\\"jammer\\">' G1 R* @  N, Y( Z6 o0 p, O  f</font><br /> main1(::a)=<font class=\\"jammer\\">% Z# ~9 K3 U2 x4 `: L4 Y</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">2 y3 B\\" N+ d9 ^5 T' z</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">9 u6 t& Q& Z7 P</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\"># Y, H2 s/ h& A9 J) p</font><br />   k++<font class=\\"jammer\\">\\" e9 L/ A# I+ A0 l! q; s4 V</font><br /> },<br /> <span style=\\"display:none\\"># P* q) j% M' v</span> s;<font class=\\"jammer\\">- Z# x' A\\" [7 D* Z) d. |! v</font><br /> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
5 j1 b# K/ k: O; K! `' U3 m20000000.
5.531   //取数据所用时间，秒

－－－－－－－－－－
+ ]5 E2 Y* _. z3 S  K) J- F+ `
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];+ B( R/ q1 D' n) g; G
   
4. mvar:1 s+ S( Q1 p# Y, ?  X$ @8 e- m
   
6. main0(::a)=
   
7. . X3 ^; D- x5 K6 r; N' Y2 l* G0 @
8. a = linspace(1,1,1000),
   
9. # {: N: y\\" y* o0 k9 t6 e3 }
10. t0=clock(),\\" r) }* c% K) Z. t
    
12. k=0, while{k<10000,
    
13. * k% @' V- H, Y8 h8 V& k/ _
14.   i=0, while{i<1000, A(a,i)=2, i++},- \, V; N9 N: Q, s' }* S4 g. F
    
16.   k++
    
17. % @+ e+ A4 n5 ^* J6 F% [6 C4 z, o
18. },6 W7 ?' O; m4 R0 [
    
20. [clock()-t0]/1000;
    
21. 2 v( m0 U; h# W7 Q7 v
22. 3 H* X: }; e7 g7 g
    
24. main1(::a)=& g7 V+ A# r8 y) t, n/ a0 g
    
26. t0=clock(), s=0,; W; j! m6 @% ^- d6 T4 c% q# V8 L! l
    
28. k=0, while{k<10000,
    
29. 2 z\\" b\\" \6 T! w) H, O7 E) l1 x( k
30.   i=0, while{i<1000, s=s+A(a,i), i++},+ P! R5 [) _1 A3 `0 t
    
32.   k++  [! ]6 ^0 d  S3 U3 L) V
    
34. },
    
35. 4 H4 X/ q% z8 a- @/ p& D
36. s;
    
37. 2 J9 V9 Z. d  q& Q5 B1 K% w9 F
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ B( R/ q1 D' n) g; G</font><br /> mvar:<font class=\\"jammer\\">1 s+ S( Q1 p# Y, ?  X$ @8 e- m</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">. X3 ^; D- x5 K6 r; N' Y2 l* G0 @</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\"># {: N: y\\" y* o0 k9 t6 e3 }</span> t0=clock(),<font class=\\"jammer\\">\\" r) }* c% K) Z. t</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">* k% @' V- H, Y8 h8 V& k/ _</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">- \, V; N9 N: Q, s' }* S4 g. F</font><br />   k++<br /> <span style=\\"display:none\\">% @+ e+ A4 n5 ^* J6 F% [6 C4 z, o</span> },<font class=\\"jammer\\">6 W7 ?' O; m4 R0 [</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">2 v( m0 U; h# W7 Q7 v</span> <font class=\\"jammer\\">3 H* X: }; e7 g7 g</font><br /> main1(::a)=<font class=\\"jammer\\">& g7 V+ A# r8 y) t, n/ a0 g</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">; W; j! m6 @% ^- d6 T4 c% q# V8 L! l</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">2 z\\" b\\" \6 T! w) H, O7 E) l1 x( k</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">+ P! R5 [) _1 A3 `0 t</font><br />   k++<font class=\\"jammer\\">  [! ]6 ^0 d  S3 U3 L) V</font><br /> },<br /> <span style=\\"display:none\\">4 H4 X/ q% z8 a- @/ p& D</span> s;<br /> <span style=\\"display:none\\">2 J9 V9 Z. d  q& Q5 B1 K% w9 F</span> [clock()-t0]/1000;

结果：
10.469
20000000.
3.781

! r7 h; \2 G  u" E－－－－－－－
8 ~1 A9 ?6 y( w3 }4 B% d2 [8 I
+ p" K( B' I# [0 M4 _7 b在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
% b2 N% j; p# U7 @  N  V+ V在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

) b" O. Y$ N! P# Bmatlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

2 a- Q& Y- Z# }: dForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
4 Q1 u8 _1 h+ Q0 k. a4 j' j
9 F6 ^3 a' L( @- B% m1 B4 v5、在循环中加点东西
+ T7 s  z5 ^, O0 V
9 A4 ^* g) F$ ~# k# [" N* mmatlab代码：

1. >> clear all$ G0 A7 r# k0 O, W+ I! R( {
   
2. f=@(x)x+1;
   9 Z- k, S9 `9 y/ S
3. a = linspace(2,2,1000);2 @) x$ V. f8 e6 C* W8 b
   
4. tic. Y0 P  r\" \\" B6 I9 N# u
   
5. s=0;
   ' `3 p! o  [& m6 B
6. for k=1:1000; N2 E/ w4 Y4 x3 @. Z# a4 z
   
7.   for i=1:1000+ d# h* V$ {& H1 z( ^
   
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));\" B6 g0 p2 Y1 T2 o$ s' M/ O
   
9.   end
   9 E' W6 M) x2 y0 f: c3 x2 O
10. end2 G\" v0 g% V2 @+ D' r
    
11. s
    $ C$ P* s* V/ x+ L6 Y
12. toc
    - t( N( H; c8 \; }) ?8 f
13. 
    9 M7 S, Z' u1 a! W: T6 t) [, H
14. s =6 w- o$ g! ?' W\" p
    
15. ) [% P4 M' [/ B3 o! @
    
16.   2.4002e+007! _$ ?' q3 R% C2 _( Y
    
17. 
    5 R2 v5 c6 _# ]. |+ d9 F
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. 2 G4 C* H! a$ U! u2 R/ z6 l
4. mvar:
   
5. 2 S6 i* ]8 R0 y: l! _
6. f(x)=x+1;$ M  y2 e# g7 P# r6 k  t
   
8. main0(::a)=\\" W1 x' d) g1 \* s  Z, m6 }% R
   
10. a = linspace(2,2,1000),: ?9 \8 R) X3 ]+ A. o\\" c! r. N  S
    
12. t0=clock(), s=0,+ v4 ]3 d. r# o8 X
    
14. k=0, while{k<1000,4 k. A( i1 M. q; ^5 K
    
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++},1 s, m# J9 y, o% v& K. r/ f4 [7 _8 S& [
    
18.   k++% d4 v, q. L& X! v( A0 Y
    
20. },
    
21. 2 z+ X: ~' }# m2 q
22. s;
    
23. 1 b6 E8 ^5 d- A1 O# E5 }6 u
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">2 G4 C* H! a$ U! u2 R/ z6 l</span> mvar:<br /> <span style=\\"display:none\\">2 S6 i* ]8 R0 y: l! _</span> f(x)=x+1;<font class=\\"jammer\\">$ M  y2 e# g7 P# r6 k  t</font><br /> main0(::a)=<font class=\\"jammer\\">\\" W1 x' d) g1 \* s  Z, m6 }% R</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">: ?9 \8 R) X3 ]+ A. o\\" c! r. N  S</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">+ v4 ]3 d. r# o8 X</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">4 k. A( i1 M. q; ^5 K</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++},<font class=\\"jammer\\">1 s, m# J9 y, o% v& K. r/ f4 [7 _8 S& [</font><br />   k++<font class=\\"jammer\\">% d4 v, q. L& X! v( A0 Y</font><br /> },<br /> <span style=\\"display:none\\">2 z+ X: ~' }# m2 q</span> s;<br /> <span style=\\"display:none\\">1 b6 E8 ^5 d- A1 O# E5 }6 u</span> [clock()-t0]/1000;

结果：
  G, S8 x/ E1 q1 @2 R6 G24002271.69189384
: l, b  c5 F$ z( m* x8 |: B$ x( l5.172
& O4 `& a. s2 c
二者效率现在只有一倍差距。

$ J4 j* P' D6 s－－－－－－－－－－
( @! `( X7 d; M8 h7 S; y5 I
6 }, i' ]2 g  q0 l* X不过，显然matlab和Forcal的程序应该这样写：
7 b/ ?7 C3 H. e6 h. J9 M) F
& E/ z# o3 ~* b& p+ ?. \7 V1 Omatlab程序：

1. clear all4 x3 ]  Y( ?+ ~9 J$ v+ v6 D
   
2. f=@(x)x+1;
   6 L* {4 Y8 g\" G7 s
3. a = linspace(2,2,1000);( ]\" B& M3 y: H/ b
   
4. tic
   6 N6 J6 T1 j+ X5 u
5. s=0;
   - b  R' x* i) w
6. for k=1:1000
   & [9 z& c! D4 n
7.   for i=1:1000$ F. x+ ?  d1 Q
   
8.     t=a(i);
   ! I+ H# w$ g$ i; ?& d7 c( |+ t5 o
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);
   # f7 C4 C, B2 @* z! l
10.   end
    & ~3 x4 q4 K. X8 L8 U% n
11. end, _2 y9 O* q\" q# ~+ J\" H. K
    
12. s
      Q2 \# i' h+ H9 K( Z5 Q
13. toc( D0 u7 T& M/ S5 |! |
    
14. # w% P/ V1 E) S: O' ?' ?
    
15. s =\" y\" F\" W4 c4 ~! v/ t; {, c
    
16. 
    2 v+ Z) m7 k, s1 i3 W& J
17.   2.4002e+007% d) P' n5 B( C5 U! N2 J\" f\" p
    
18. ) P5 g1 H3 v3 _+ @4 _
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];3 r7 A8 C5 p* S  x. v! w
   
4. mvar:
   
5. , f* j8 {6 B4 Q& ^
6. f(x)=x+1;+ i- U3 a6 c: f! p+ t4 y
   
8. main0(::a)=
   
9. ! a. t5 d* L# Y6 H8 L7 ^6 Z9 j
10. a = linspace(2,2,1000),
    
11. . I9 x5 w# K& D
12. t0=clock(), s=0,% M( B2 f& Q3 P! S  G
    
14. k=0, while{k<1000,/ s) m8 I+ X! ~- h. h: r5 J
    
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. 0 E9 N4 w6 ~6 u3 \! _/ x
18.   k++
    
19.   k% j2 ]0 N7 p* f
20. },
    
21. 4 m4 p: I$ y$ ^6 w& ~& I
22. s;/ l  q2 E& C* E: }
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">3 r7 A8 C5 p* S  x. v! w</font><br /> mvar:<br /> <span style=\\"display:none\\">, f* j8 {6 B4 Q& ^</span> f(x)=x+1;<font class=\\"jammer\\">+ i- U3 a6 c: f! p+ t4 y</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">! a. t5 d* L# Y6 H8 L7 ^6 Z9 j</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">. I9 x5 w# K& D</span> t0=clock(), s=0,<font class=\\"jammer\\">% M( B2 f& Q3 P! S  G</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">/ s) m8 I+ X! ~- h. h: r5 J</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\\">0 E9 N4 w6 ~6 u3 \! _/ x</span>   k++<br /> <span style=\\"display:none\\">  k% j2 ]0 N7 p* f</span> },<br /> <span style=\\"display:none\\">4 m4 p: I$ y$ ^6 w& ~& I</span> s;<font class=\\"jammer\\">/ l  q2 E& C* E: }</font><br /> [clock()-t0]/1000;

结果：
& j5 e) K/ u6 V6 A: a24002271.69189384
1.625
% B/ D; j# S( y- Z0 R& L
* J8 L3 y- v& c" n3 cmatlab效率已不及Forcal。
7 g7 F1 u  N& J4 I- |6 H# ]
4 {  T( e0 S( {5 T提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

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

Forcal代码：

2. !using["sys"];6 |6 @1 ^1 D$ }5 f
   
4. mvar:
   
5. * j* _/ j! d; e* o/ }
6. a="\&\[1000]"/4,    //用字符串模拟静态数组# B/ f: B- c8 i1 n# P: I6 O. w
   
8. t0=clock(),
   
9. 2 f6 F5 H7 V3 L% a
10. k=0, while{k<10000,
    
11. ) o1 v( e* `: @- P
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值( u% K) [! V3 W5 U/ N. C$ B8 U
    
14.   k++/ D# s' w3 [: f7 w# |* B
    
16. },
    
17. 0 P# L2 a& g4 H( n# @
18. [clock()-t0]/1000;1 B0 U5 c' w+ M2 g* u0 D6 a. W
    
20. 
    
21. 8 R/ E\\" I: r% l0 W' v+ Z: M0 F
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. 5 w1 c1 r' C4 ]$ x$ Z
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. 0 G0 s7 I& z3 U5 e2 H
26. t0=clock(),7 H  U1 b- P8 e6 N( e0 |, _
    
28. s=0,
    
29. $ ^! S0 v& G: w: H\\" X
30. k=0, while{k<10000,
    
31. / U- Z: S) j\\" }/ F4 @( }2 F; z6 f0 I
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. ; x1 H. A( M7 |- {
34.   k++
    
35. ; O! e3 f( E9 }\\" R) k$ J
36. },
    
37. 5 r4 d0 _2 m1 w5 a* {  G9 `
38. s;
    
39. \\" ]* c) ^3 ^9 c
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">6 |6 @1 ^1 D$ }5 f</font><br /> mvar:<br /> <span style=\\"display:none\\">* j* _/ j! d; e* o/ }</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\"># B/ f: B- c8 i1 n# P: I6 O. w</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">2 f6 F5 H7 V3 L% a</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">) o1 v( e* `: @- P</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">( u% K) [! V3 W5 U/ N. C$ B8 U</font><br />   k++<font class=\\"jammer\\">/ D# s' w3 [: f7 w# |* B</font><br /> },<br /> <span style=\\"display:none\\">0 P# L2 a& g4 H( n# @</span> [clock()-t0]/1000;<font class=\\"jammer\\">1 B0 U5 c' w+ M2 g* u0 D6 a. W</font><br /> <br /> <span style=\\"display:none\\">8 R/ E\\" I: r% l0 W' v+ Z: M0 F</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">5 w1 c1 r' C4 ]$ x$ Z</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">0 G0 s7 I& z3 U5 e2 H</span> t0=clock(),<font class=\\"jammer\\">7 H  U1 b- P8 e6 N( e0 |, _</font><br /> s=0,<br /> <span style=\\"display:none\\">$ ^! S0 v& G: w: H\\" X</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">/ U- Z: S) j\\" }/ F4 @( }2 F; z6 f0 I</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">; x1 H. A( M7 |- {</span>   k++<br /> <span style=\\"display:none\\">; O! e3 f( E9 }\\" R) k$ J</span> },<br /> <span style=\\"display:none\\">5 r4 d0 _2 m1 w5 a* {  G9 `</span> s;<br /> <span style=\\"display:none\\">\\" ]* c) ^3 ^9 c</span> [clock()-t0]/1000;

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