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

forcal

1、简单循环
- b+ {3 b* I2 x: w( _6 p0 [9 q
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   : ?% H! E- u  O\" ^3 y/ k% b
2. tic
   \" j\" V9 e5 `% R
3. s=0;5 w: K* C1 w( \$ t
   
4. for k=1:100004 h$ r! o& J. E1 h& I$ U
   
5.   for i=1:1000' P+ F* p% }9 ^* z* v* f8 N, t
   
6.     s=s+1;; J( m( u\" _) ?0 ~& ]
   
7.   end2 v9 M2 O0 Y8 h# o# y  o
   
8. end$ G7 I2 a5 B, M. a; Y& l$ \8 R
   
9. s) A* a% S8 p0 K1 b+ X
   
10. toc
    . R- M6 |5 F\" b/ j
11. 
    ' F& F1 Z+ [  v  s8 C
12. s =* Z8 \' v9 A\" y6 x4 [$ W
    
13. + h& Q% T& Q1 l, e# |+ e
    
14.     10000000
    2 w. b4 b, {$ P+ r
15. 6 O6 I6 y0 L9 g. T9 s  b% |+ A
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

3 v2 f* b0 k3 ]6 [. SForcal没有JIT加速器，代码：

2. !using["sys"];
   
3. 1 W  C  P! o9 x( ^8 i\\" l7 a
4. mvar:. x% i* X6 j) h
   
6. t0=clock(),
   
7. ; e4 P% ?# ^( `. {
8. s=0,
   
9. 3 c6 |\\" o$ s: g  k6 N4 p
10. k=0, while{k<10000,
    
11. ( ^\\" R, v% N4 l# w/ k
12.   i=0, while{i<1000,s=s+1, i++},
    
13. $ i0 l/ e# n7 l1 k( h7 G( R) g
14.   k++
    
15. ( S/ S- |% J$ T6 W) ?
16. },
    
17. ' g% m* J# {2 v7 o; i
18. s;
    
19. \\" z7 O% X) H+ L
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">1 W  C  P! o9 x( ^8 i\\" l7 a</span> mvar:<font class=\\"jammer\\">. x% i* X6 j) h</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">; e4 P% ?# ^( `. {</span> s=0,<br /> <span style=\\"display:none\\">3 c6 |\\" o$ s: g  k6 N4 p</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">( ^\\" R, v% N4 l# w/ k</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">$ i0 l/ e# n7 l1 k( h7 G( R) g</span>   k++<br /> <span style=\\"display:none\\">( S/ S- |% J$ T6 W) ?</span> },<br /> <span style=\\"display:none\\">' g% m* J# {2 v7 o; i</span> s;<br /> <span style=\\"display:none\\">\\" z7 O% X) H+ L</span> [clock()-t0]/1000;

结果：
10000000.
0.704  秒

& f; P' ^% j' _4 h9 r4 `' C4 a二者运行效率相差7倍多。

2、循环中稍微加点东西
0 Y# I$ o9 {; e, T( `; F4 ?
Matlab代码：

1. clear all
   9 j7 n1 d, h, y; Q5 N
2. tic
   5 E! W\" K7 v/ s+ F/ m6 K2 i
3. s=0;
   ' [6 k\" @( Q$ ?6 F9 T+ X! h3 D
4. for k=1:1000
   ( E9 Y4 f: L5 O8 T1 S0 g
5.   for i=1:1000
   - T4 m, X8 i' e/ C
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);: a$ K: r\" t8 y5 U
   
7.   end
   0 V) Z% @9 o4 L
8. end
   6 X+ ?\" P9 R; V: S! A6 A$ R' p# A
9. s# o0 t, V0 A% Q+ O2 k; Z* a, e
   
10. toc
    # j8 J9 b+ ?' D. G/ l- X+ Q
11. 1 l$ D* i5 Y% a1 @\" `+ H
    
12. s =
    ( C  z2 ^3 m( ?7 F, c1 R7 v
13. 
    ' V0 v+ U, g$ i6 x* e1 g
14.   1.1009e+007
    3 {1 r0 j4 Q$ @
15. 8 r# ?( u! z$ p4 V1 b& g* P\" `\" e# V
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. ) c. I- t' ?, L* M5 T) b& G
4. mvar:) v9 m# T. J8 R$ M
   
6. t0=clock(),6 ?2 t7 @, L8 `! ~7 p
   
8. s=0,+ B: \4 a; a  R; k
   
10. k=1, while{k<=1000,
    
11. 4 {0 m& I( y- d$ 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++},
    
13. , O: ?5 C9 Y+ X0 P
14.   k++
    
15. ! }2 }' e! z( N: f8 Q+ k
16. },
    
17. 1 A\\" ~$ Y% [: o3 g+ f/ Q# }
18. s;4 }, `9 K+ D. X\\" D, e/ f7 z- p
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">) c. I- t' ?, L* M5 T) b& G</span> mvar:<font class=\\"jammer\\">) v9 m# T. J8 R$ M</font><br /> t0=clock(),<font class=\\"jammer\\">6 ?2 t7 @, L8 `! ~7 p</font><br /> s=0,<font class=\\"jammer\\">+ B: \4 a; a  R; k</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">4 {0 m& I( y- d$ 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++},<br /> <span style=\\"display:none\\">, O: ?5 C9 Y+ X0 P</span>   k++<br /> <span style=\\"display:none\\">! }2 }' e! z( N: f8 Q+ k</span> },<br /> <span style=\\"display:none\\">1 A\\" ~$ Y% [: o3 g+ f/ Q# }</span> s;<font class=\\"jammer\\">4 }, `9 K+ D. X\\" D, e/ f7 z- p</font><br /> [clock()-t0]/1000;

结果：
! ]2 f; U/ z6 v, j$ {* c" d, H11008953.66602346
1.031
, D1 i0 g# b# W3 l  k2 r* J
二者效率相差已比较小。
4 Q3 E1 s5 U& ]( m9 Z/ ~9 \* f8 a+ _
3、再增加一个函数调用
) ~: {# h2 }' |2 l1 x
; g/ G7 E, T7 A: HMatlab代码：

1. clear all
   3 r: C4 o+ E$ g# f; g2 |
2. f=@(x)x+1;
   7 ^3 X) M! ?9 @* x9 q
3. tic
   . F8 S- f5 m$ C
4. s=0;6 c: K: C6 r5 j
   
5. for k=1:1000# {3 X0 g\" w9 [' X; u& E
   
6.   for i=1:1000- y1 d$ Y1 E) W) S: Q/ W
   
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);7 A$ d# }- B0 y, h* U3 t
   
8.   end( d# r5 I9 j! C+ O3 o# S, X
   
9. end9 v% I5 |  O\" i; n3 W\" k
   
10. s* x5 K' {8 \- Z7 G$ r1 N
    
11. toc
    % O$ S8 F1 A) h& g  d
12. 
    4 |- e+ w  E. z; N5 X: W% W; B3 }
13. s =0 @. {1 ?1 g/ y2 s: n2 Y
    
14. 
    % Y/ ^9 l9 Q% o/ C
15.   5.1251e+008
    # a  z/ S3 Y' J' i) ^6 P
16. 
    ) d5 }7 }\" h' ^1 F% f
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. - x7 G- q' I: H+ `! ?9 O
4. f(x)=x+1;7 B5 T7 b. A8 c
   
6. mvar:+ u' T/ X0 @' y# S8 Q8 `: V6 T2 h
   
8. t0=clock(),, m2 y& t0 c( J. Y) s0 d# P3 X
   
10. s=0,
    
11. 8 m- j6 G4 c7 T- A  K$ _' [  K
12. k=1, while{k<=1000,
    
13. 0 R' O\\" `& t& f
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++},# e6 L' @; p4 p  M
    
16.   k++
    
17. 8 J- _# K6 _\\" [
18. },5 k0 J9 O) I6 m. Z! W9 b
    
20. s;
    
21. 0 `) e  G1 v! U\\" b\\" a: M. d
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">- x7 G- q' I: H+ `! ?9 O</span> f(x)=x+1;<font class=\\"jammer\\">7 B5 T7 b. A8 c</font><br /> mvar:<font class=\\"jammer\\">+ u' T/ X0 @' y# S8 Q8 `: V6 T2 h</font><br /> t0=clock(),<font class=\\"jammer\\">, m2 y& t0 c( J. Y) s0 d# P3 X</font><br /> s=0,<br /> <span style=\\"display:none\\">8 m- j6 G4 c7 T- A  K$ _' [  K</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">0 R' O\\" `& t& f</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\\"># e6 L' @; p4 p  M</font><br />   k++<br /> <span style=\\"display:none\\">8 J- _# K6 _\\" [</span> },<font class=\\"jammer\\">5 k0 J9 O) I6 m. Z! W9 b</font><br /> s;<br /> <span style=\\"display:none\\">0 `) e  G1 v! U\\" b\\" a: M. d</span> [clock()-t0]/1000;

结果：
512508953.6658993
1.078

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

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
- }4 ~$ [6 Y- D. D8 _
Matlab代码：

1. >> clear all
   , O) }; s( d3 [% a& @
2. a = linspace(1,1,1000);' p, w  y# n; x% N# R
   
3. tic
   5 |5 l' Z& P# |  X9 H3 T2 l: M\" K6 B# d
4. for k=1:100003 m3 n* ?, y5 ?4 m5 \, N; ]7 A
   
5.   for i=1:1000& F\" `! j5 T! Y4 v4 x
   
6.     a(i)=2;. ^# z# k) ^9 J1 b; G. V( U: u
   
7.   end
     p1 V3 Z+ m/ @& i5 _\" R
8. end6 Z: k8 z+ i9 u0 J3 C7 V\" k
   
9. toc
   & t- S( p- F\" s6 E! g' q+ w
10. : M: j7 _$ H* I0 b: h
    
11. tic
    $ \8 t0 _9 y  d/ z- D
12. s=0;( c) ~8 c6 q: r* \) W7 Z, _
    
13. for k=1:100007 ?6 s, H. o9 x9 M! s
    
14.   for i=1:1000
    1 o8 `  v3 x* X. [\" s\" R
15.     s=s+a(i);
    + S, s' p5 a: P% N/ s- m9 U
16.   end* O+ i0 g: x  P0 Q6 z: X
    
17. end- R8 S' q* Q' C% b2 H* s* U$ U- N
    
18. s
    * U7 D+ h/ o0 N- Q1 ^: E- m4 v
19. toc, o. A. n( C8 t5 i
    
20. Elapsed time is 0.140658 seconds.\" h9 w3 L, t# F, s- G
    
21. 
    # k8 f. ?: Z2 ~6 P% {
22. s =: _( U' L/ d\" S8 J4 n1 {% B: Z
    
23. 
    - o1 u4 W2 G& d
24.     20000000
    + f. g0 l* z8 a& {% E# Y
25. 
    0 m6 T, z! G, t
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
4 d5 ^/ R9 O% o: F0 s% W- d
- a- i- C3 \( i; G/ qForcal代码：

2. !using["math","sys"];
   
3. ' ~1 t7 W' c6 m2 r% R3 D( k% O8 J
4. mvar:
   
5. ' T$ V3 M8 M. o7 K( Q
6. main0(::a)=
   
7. 0 I* o! D: A0 H9 v\\" j( |
8. a = linspace(1,1,1000),
   
9. & q( j9 x4 T# F
10. t0=clock(),- x$ n* u0 M9 M# f/ Q( H9 z3 |
    
12. k=0, while{k<10000,. Y. P2 J9 N6 F+ o% O3 ~\\" G- g
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. 3 A! V8 k0 C$ M( K
16.   k++8 H: r7 F. u2 J' J& y& X* Z3 c4 o
    
18. },/ f4 O* ]7 i/ G7 G2 {
    
20. [clock()-t0]/1000;0 o& _+ E\\" V\\" N% L; M
    
22. 
    
23. & m: x8 A\\" b; L\\" K6 |$ i
24. main1(::a)=
    
25. 7 W5 f$ b2 n! O$ c3 P
26. t0=clock(), s=0,
    
27. # A) {. R1 J: c) c* z2 |
28. k=0, while{k<10000,' U4 y# Y! G4 u% Y+ {
    
30.   i=0, while{i<1000, s=s+a(i), i++},\\" K9 p7 v$ B* K& S
    
32.   k++
    
33. 8 ~) q+ d- ^. b+ J# }2 V( n
34. },& F8 j9 P7 o/ V3 I
    
36. s;
    
37. 3 x! F- Q+ Z' ~7 t$ l) i
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">' ~1 t7 W' c6 m2 r% R3 D( k% O8 J</span> mvar:<br /> <span style=\\"display:none\\">' T$ V3 M8 M. o7 K( Q</span> main0(::a)=<br /> <span style=\\"display:none\\">0 I* o! D: A0 H9 v\\" j( |</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">& q( j9 x4 T# F</span> t0=clock(),<font class=\\"jammer\\">- x$ n* u0 M9 M# f/ Q( H9 z3 |</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">. Y. P2 J9 N6 F+ o% O3 ~\\" G- g</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">3 A! V8 k0 C$ M( K</span>   k++<font class=\\"jammer\\">8 H: r7 F. u2 J' J& y& X* Z3 c4 o</font><br /> },<font class=\\"jammer\\">/ f4 O* ]7 i/ G7 G2 {</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">0 o& _+ E\\" V\\" N% L; M</font><br /> <br /> <span style=\\"display:none\\">& m: x8 A\\" b; L\\" K6 |$ i</span> main1(::a)=<br /> <span style=\\"display:none\\">7 W5 f$ b2 n! O$ c3 P</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\"># A) {. R1 J: c) c* z2 |</span> k=0, while{k<10000,<font class=\\"jammer\\">' U4 y# Y! G4 u% Y+ {</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">\\" K9 p7 v$ B* K& S</font><br />   k++<br /> <span style=\\"display:none\\">8 ~) q+ d- ^. b+ J# }2 V( n</span> },<font class=\\"jammer\\">& F8 j9 P7 o/ V3 I</font><br /> s;<br /> <span style=\\"display:none\\">3 x! F- Q+ Z' ~7 t$ l) i</span> [clock()-t0]/1000;

结果：
* ]' J; N' B7 N6 s15.484  //存数据所用时间，秒
1 }0 T0 u6 W9 @+ b( r1 e20000000.
5.531   //取数据所用时间，秒
0 H( d7 e8 a9 V, v
. N5 g; p, I9 g9 q－－－－－－－－－－

' j  s5 l2 g9 x用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];- B4 B- I8 d/ V( ~9 |\\" |
   
4. mvar:
   
5. ' U3 [6 u5 z5 W3 O
6. main0(::a)=
   
7. # ]& p; o2 ]' h1 e) E% M, m  H
8. a = linspace(1,1,1000),
   
9. - l9 _& c0 Z. k# A  M
10. t0=clock(),
    
11. , e- `. g: u7 Y2 D; Q) i; v. C
12. k=0, while{k<10000,7 m. }! t\\" V  d* `
    
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. 3 `0 L5 P5 g* |6 a2 Y: Y
16.   k++
    
17. 5 p3 S\\" n9 J( n' {8 I, ]0 V
18. },0 W+ p4 r& k' D
    
20. [clock()-t0]/1000;
    
21. & B2 r7 X& E& E, r) g& y
22. 8 S7 M$ o8 i$ ]: N3 P: H1 K
    
24. main1(::a)=
    
25. % X* |: i2 ^\\" G# K8 w
26. t0=clock(), s=0,5 Y, i3 G  C' L5 I( h% }
    
28. k=0, while{k<10000,4 S9 Q\\" [; q& O) M2 a0 w  n- W
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. 7 `  T1 z) c  ~
32.   k++* q9 i+ f4 m. G5 _( f
    
34. },+ S5 W8 k\\" P3 j$ ^9 S! A4 \
    
36. s;
    
37. 0 h! z0 Z4 z\\" e\\" r6 x\\" Z# b5 d
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">- B4 B- I8 d/ V( ~9 |\\" |</font><br /> mvar:<br /> <span style=\\"display:none\\">' U3 [6 u5 z5 W3 O</span> main0(::a)=<br /> <span style=\\"display:none\\"># ]& p; o2 ]' h1 e) E% M, m  H</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">- l9 _& c0 Z. k# A  M</span> t0=clock(),<br /> <span style=\\"display:none\\">, e- `. g: u7 Y2 D; Q) i; v. C</span> k=0, while{k<10000,<font class=\\"jammer\\">7 m. }! t\\" V  d* `</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">3 `0 L5 P5 g* |6 a2 Y: Y</span>   k++<br /> <span style=\\"display:none\\">5 p3 S\\" n9 J( n' {8 I, ]0 V</span> },<font class=\\"jammer\\">0 W+ p4 r& k' D</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">& B2 r7 X& E& E, r) g& y</span> <font class=\\"jammer\\">8 S7 M$ o8 i$ ]: N3 P: H1 K</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">% X* |: i2 ^\\" G# K8 w</span> t0=clock(), s=0,<font class=\\"jammer\\">5 Y, i3 G  C' L5 I( h% }</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">4 S9 Q\\" [; q& O) M2 a0 w  n- W</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">7 `  T1 z) c  ~</span>   k++<font class=\\"jammer\\">* q9 i+ f4 m. G5 _( f</font><br /> },<font class=\\"jammer\\">+ S5 W8 k\\" P3 j$ ^9 S! A4 \</font><br /> s;<br /> <span style=\\"display:none\\">0 h! z0 Z4 z\\" e\\" r6 x\\" Z# b5 d</span> [clock()-t0]/1000;

结果：
9 q0 ]3 F/ C1 b; w8 K# [10.469
20000000.
1 m# s+ N% z; b3.781
9 \$ @$ o  [1 O. a, O( z; Z
. c& G9 C1 ]) S5 `' u2 r; I－－－－－－－
4 W& w. _3 I+ \+ d! E
在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
# G: u% p* z  B+ B1 `6 C
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

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

5、在循环中加点东西

matlab代码：

1. >> clear all
     N7 t4 U: e: {9 `( B0 ^- }
2. f=@(x)x+1;
   & W2 w\" {2 }9 Y) V3 q0 N
3. a = linspace(2,2,1000);9 n) E, s3 e8 d5 J) R5 I# J5 R4 Z
   
4. tic( g. G7 q( t) m7 X# t% L+ W! y
   
5. s=0;
   ' c; E* k# b- o- Y1 l- v: t
6. for k=1:1000
   * Y4 D) |\" p' ]
7.   for i=1:1000& K3 L$ n: w& I  l6 }. w
   
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));/ J2 @! Q! n/ f5 Y, [
   
9.   end
   . v4 D8 H7 n% E3 E
10. end7 [( M$ k0 q$ H7 D
    
11. s\" j) l6 G4 d; m- d
    
12. toc) E! K# o8 G\" c8 O' ?8 o% s
    
13. 
    - \7 [\" Y8 V, N# P
14. s =
    ! M; {6 \8 C* N9 U
15. 0 W/ E4 |. o' C7 {
    
16.   2.4002e+0079 K9 P9 q  [/ k! N' c5 O
    
17. 
    : m% q! |6 U( A
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. / t2 y8 t2 Z8 i% E, P3 {* l! e
4. mvar:
   
5. ( y1 K/ ?0 c\\" U
6. f(x)=x+1;0 S2 I; |2 l9 R6 h\\" {& R
   
8. main0(::a)=
   
9. $ y5 P2 T- K8 o' s# }
10. a = linspace(2,2,1000),
    
11. 6 C  C6 J3 k5 E# |2 T$ ~( G
12. t0=clock(), s=0,
    
13. 5 ?7 P( q9 i  f/ z0 M! w- E
14. k=0, while{k<1000,
    
15. 6 U6 e% K+ U; y5 x) S\\" [
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++},4 B  A& t' f0 F, H# I
    
18.   k++: |/ b& A  T1 q4 K% [. x
    
20. },
    
21. \\" L! P% `( ~! p/ x
22. s;
    
23. ' r8 j' y) g  V% |$ S+ {& Z
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">/ t2 y8 t2 Z8 i% E, P3 {* l! e</span> mvar:<br /> <span style=\\"display:none\\">( y1 K/ ?0 c\\" U</span> f(x)=x+1;<font class=\\"jammer\\">0 S2 I; |2 l9 R6 h\\" {& R</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">$ y5 P2 T- K8 o' s# }</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">6 C  C6 J3 k5 E# |2 T$ ~( G</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">5 ?7 P( q9 i  f/ z0 M! w- E</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">6 U6 e% K+ U; y5 x) S\\" [</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\\">4 B  A& t' f0 F, H# I</font><br />   k++<font class=\\"jammer\\">: |/ b& A  T1 q4 K% [. x</font><br /> },<br /> <span style=\\"display:none\\">\\" L! P% `( ~! p/ x</span> s;<br /> <span style=\\"display:none\\">' r8 j' y) g  V% |$ S+ {& Z</span> [clock()-t0]/1000;

结果：
24002271.69189384
5.172
9 s& c. B; A6 X# B) s
二者效率现在只有一倍差距。

－－－－－－－－－－

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

. ?( z& B6 O' N5 q* |! pmatlab程序：

1. clear all
   + A/ y9 r( K* v
2. f=@(x)x+1;
   4 s! J- m; G3 x1 m- t7 C
3. a = linspace(2,2,1000);
   ; t. R; U4 o4 C3 D! }
4. tic
   7 P4 _, B' q3 N7 `( t2 X  m: I9 ^5 o
5. s=0;
   8 i- b. A) A/ A) |. r
6. for k=1:10008 C3 ]2 z* f  K4 `2 C6 `+ c/ I- E
   
7.   for i=1:1000
   ! d4 x5 w; _- [* m/ T# O9 R5 _
8.     t=a(i);+ O0 Q5 j( ?5 B5 g/ h
   
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);4 L) y\" _+ i8 W  g
   
10.   end4 W. {. `6 H) a, U, z
    
11. end\" D( Z$ t& p8 ~+ M, A+ _
    
12. s
    . i. ~9 H7 G; g* p( y2 @  l5 b
13. toc
    \" `2 Q  _3 g1 C- b
14. 5 M) Q; ^1 E1 W\" {. v/ m
    
15. s =0 V; m' D  }- M! U6 p( I' H
    
16. 
    ; R/ H  i$ Z+ ~( E, _. t0 c
17.   2.4002e+007
    ' _+ Z0 [. c2 v2 Q! ?( q2 }
18. : X9 w; z/ S0 d
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];$ s- H- t( j8 f4 z  Q
   
4. mvar:+ m! ~2 Q0 W4 b\\" p( J
   
6. f(x)=x+1;
   
7. ! y$ A; i) t; u* E
8. main0(::a)=
   
9. 4 \7 o* I5 \6 v! n- k- u( |
10. a = linspace(2,2,1000),
    
11. 8 v, L$ G# E; P, Z# b2 D7 t; g- }
12. t0=clock(), s=0,
    
13. $ j2 M. y. E9 `5 M* J
14. k=0, while{k<1000,) ~\\" O3 ~( P1 c' G, k\\" L! ?
    
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 x* @- V2 {
18.   k++
    
19. % q' h8 Y( T) ^( [9 `$ ~
20. },
    
21. # d& O6 L! e1 R: q
22. s;. V  g% l8 t5 P8 D& O# b7 T
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">$ s- H- t( j8 f4 z  Q</font><br /> mvar:<font class=\\"jammer\\">+ m! ~2 Q0 W4 b\\" p( J</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">! y$ A; i) t; u* E</span> main0(::a)=<br /> <span style=\\"display:none\\">4 \7 o* I5 \6 v! n- k- u( |</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">8 v, L$ G# E; P, Z# b2 D7 t; g- }</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">$ j2 M. y. E9 `5 M* J</span> k=0, while{k<1000,<font class=\\"jammer\\">) ~\\" O3 ~( P1 c' G, k\\" L! ?</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\\">/ {4 x* @- V2 {</span>   k++<br /> <span style=\\"display:none\\">% q' h8 Y( T) ^( [9 `$ ~</span> },<br /> <span style=\\"display:none\\"># d& O6 L! e1 R: q</span> s;<font class=\\"jammer\\">. V  g% l8 t5 P8 D& O# b7 T</font><br /> [clock()-t0]/1000;

结果：
+ O$ }3 E) E, E3 n1 _24002271.69189384
1.625

matlab效率已不及Forcal。
( w# I9 E, w* H& E* L" `" H5 }; p" ~8 C
/ Y0 q- O7 Z/ t3 Y提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

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

  H2 B2 [2 Y& y# K6 n5 [2 f5 p  c! |Forcal代码：

2. !using["sys"];
   
3. \\" L* D$ I: F2 `/ {5 f
4. mvar:9 A$ N7 m3 D9 n; ~- ^. ?+ h
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. ! d: X# D8 X) f& f1 c
8. t0=clock(),
   
9. , W3 }8 E* i  W8 h# J1 G# m1 P* b
10. k=0, while{k<10000,
    
11. , T: y  u+ \: x6 Z& v
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. , _$ [7 ?2 Y3 Z) I
14.   k++
    
15. - M) b' y: g% N/ F8 u! I/ |
16. },
    
17. ( L% i! O- P9 ]; {
18. [clock()-t0]/1000;- ?. j# W5 ]% c2 Z7 Z2 ]& |
    
20. . K% w3 h) ]) z+ r
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组' o/ J2 t& d. r4 M: a0 a) Y
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值+ r4 c1 K$ j9 F
    
26. t0=clock(),
    
27. + N# X; h; r\\" E$ L/ y
28. s=0,: }) E; E. }8 A) J5 @; h7 O$ ^
    
30. k=0, while{k<10000,
    
31. % F4 y1 J3 P4 @% I
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. & ?: |0 k' A4 ?5 }( {
34.   k++
    
35. $ [  R; M) B9 A5 @# I3 s. w
36. },- j  ?) @1 U% L; f# o! I
    
38. s;
    
39. 0 D5 H2 o- Y* x$ K
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">\\" L* D$ I: F2 `/ {5 f</span> mvar:<font class=\\"jammer\\">9 A$ N7 m3 D9 n; ~- ^. ?+ h</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">! d: X# D8 X) f& f1 c</span> t0=clock(),<br /> <span style=\\"display:none\\">, W3 }8 E* i  W8 h# J1 G# m1 P* b</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">, T: y  u+ \: x6 Z& v</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">, _$ [7 ?2 Y3 Z) I</span>   k++<br /> <span style=\\"display:none\\">- M) b' y: g% N/ F8 u! I/ |</span> },<br /> <span style=\\"display:none\\">( L% i! O- P9 ]; {</span> [clock()-t0]/1000;<font class=\\"jammer\\">- ?. j# W5 ]% c2 Z7 Z2 ]& |</font><br /> <font class=\\"jammer\\">. K% w3 h) ]) z+ r</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">' o/ J2 t& d. r4 M: a0 a) Y</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">+ r4 c1 K$ j9 F</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">+ N# X; h; r\\" E$ L/ y</span> s=0,<font class=\\"jammer\\">: }) E; E. }8 A) J5 @; h7 O$ ^</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">% F4 y1 J3 P4 @% I</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">& ?: |0 k' A4 ?5 }( {</span>   k++<br /> <span style=\\"display:none\\">$ [  R; M) B9 A5 @# I3 s. w</span> },<font class=\\"jammer\\">- j  ?) @1 U% L; f# o! I</font><br /> s;<br /> <span style=\\"display:none\\">0 D5 H2 o- Y* x$ K</span> [clock()-t0]/1000;

结果：
2 x1 s% Y! i* i. i0.5620000000000001  //赋值用时间，秒
20000000.
0.579      //取值用时间，秒