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

forcal

1、简单循环
5 }; ?' `+ \* ]9 M/ ]+ B% a
+ t' y: L4 p+ j" LMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all6 z# E. W- u+ ^4 F2 z
   
2. tic
   0 q' H* K/ }$ [
3. s=0;
   ! E  Z, |$ x1 y9 ~
4. for k=1:10000- _2 B) r# ^, E# G2 Y/ u7 X0 G% h
   
5.   for i=1:1000  \) G, A. m+ H
   
6.     s=s+1;
   ) B7 O$ _; X8 \( P+ N; x# N
7.   end% p7 F& O0 ?\" l2 V
   
8. end
   4 P% X9 ?; K+ W8 c& F- ]% n0 M
9. s. b) I; T: h0 w: M* V- Y3 D+ A4 |
   
10. toc  y) s$ m0 m\" E
    
11. 
    . ]% T7 r& K# q& Q. F# j! A
12. s =3 d: S  e! Z( ^- o8 n1 J5 b
    
13. / P# k) n, c3 ^3 V9 s* B1 v+ ~
    
14.     100000004 R( N- K2 w! G! C9 o( j. q
    
15. 
    , G! ~. m- T( D% n# y+ r* ]
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
- r% z% h5 F/ i0 Q7 u
Forcal没有JIT加速器，代码：

2. !using["sys"];5 r% u& Q4 s: N2 U( R
   
4. mvar:
   
5. 6 g8 C) ^# \4 [! `8 T
6. t0=clock(),7 v3 m' @6 f! v! w$ @
   
8. s=0,
   
9. ) y. p\\" ]7 Y6 @; ~) n: `4 g$ e
10. k=0, while{k<10000,
    
11. 9 {* G/ X' U- U2 y* Z
12.   i=0, while{i<1000,s=s+1, i++},1 t0 |8 L+ }8 t* r0 q- T4 k/ y. k
    
14.   k++  r* ~! X  Q; q) e  A\\" b
    
16. },% \  o0 e8 F# G- B
    
18. s;
    
19. & u8 A3 @2 l& h$ [6 Z0 N' I
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">5 r% u& Q4 s: N2 U( R</font><br /> mvar:<br /> <span style=\\"display:none\\">6 g8 C) ^# \4 [! `8 T</span> t0=clock(),<font class=\\"jammer\\">7 v3 m' @6 f! v! w$ @</font><br /> s=0,<br /> <span style=\\"display:none\\">) y. p\\" ]7 Y6 @; ~) n: `4 g$ e</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">9 {* G/ X' U- U2 y* Z</span>   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">1 t0 |8 L+ }8 t* r0 q- T4 k/ y. k</font><br />   k++<font class=\\"jammer\\">  r* ~! X  Q; q) e  A\\" b</font><br /> },<font class=\\"jammer\\">% \  o0 e8 F# G- B</font><br /> s;<br /> <span style=\\"display:none\\">& u8 A3 @2 l& h$ [6 Z0 N' I</span> [clock()-t0]/1000;

结果：
4 ]3 ~, }0 a9 R1 o# I: n' p/ Z% M9 _10000000.
4 ~4 B9 ~: h- |1 Q1 B7 {: o0.704  秒

二者运行效率相差7倍多。
3 z% p/ _3 K% D( [
2、循环中稍微加点东西
" |" o3 @% F: v# j0 ~
" z4 ^% l) R4 i0 t, U1 EMatlab代码：

1. clear all
   6 X  \* j( O: W1 u1 E\" K6 j% q
2. tic
   \" S/ t! Q\" W4 I1 H  [6 z7 c
3. s=0;& o' L# ^$ F& \! e' r2 C+ `  _
   
4. for k=1:1000
   / w  d2 u: h. M6 W/ C( a8 q
5.   for i=1:10000 `& i1 F' D' R! a( S  i! t! V
   
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);
   / Q3 m* t. R1 T+ I. d
7.   end
   + J  R! ~$ G. u; d' @
8. end7 R  j5 j+ E% s  H
   
9. s
   ! q  e1 v$ Q\" t
10. toc
    6 P8 @- R0 J\" k: x+ ]- b
11. 7 a# a9 R' C9 T6 n
    
12. s =
    ; u* E7 S0 ^7 C: ?\" t
13. ( R* P5 d5 e) c: k% k5 ^
    
14.   1.1009e+0073 N  q: O+ m! m8 j3 A2 Q; v' ~
    
15. 
    9 v# q  o3 \$ e; d* y3 i# d: l8 E0 Y
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. ( z6 e+ u$ L6 Y4 a- a, ~
4. mvar:
   
5. , d+ L* n+ H$ T9 ~5 J8 z
6. t0=clock(),# F4 o. |  a+ U, C
   
8. s=0,
   
9. \\" u% B2 T0 v( x# t$ e
10. k=1, while{k<=1000,# s% G+ W2 y, B. h, W! b! u
    
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. % H% D0 v. t7 I7 n0 n6 \
14.   k++# d; [% ~4 W' ?/ J
    
16. },; x/ m/ I# A0 E) m& h
    
18. s;* B  B7 c9 t% J! O5 T
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">( z6 e+ u$ L6 Y4 a- a, ~</span> mvar:<br /> <span style=\\"display:none\\">, d+ L* n+ H$ T9 ~5 J8 z</span> t0=clock(),<font class=\\"jammer\\"># F4 o. |  a+ U, C</font><br /> s=0,<br /> <span style=\\"display:none\\">\\" u% B2 T0 v( x# t$ e</span> k=1, while{k<=1000,<font class=\\"jammer\\"># s% G+ W2 y, B. h, W! b! u</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++},<br /> <span style=\\"display:none\\">% H% D0 v. t7 I7 n0 n6 \</span>   k++<font class=\\"jammer\\"># d; [% ~4 W' ?/ J</font><br /> },<font class=\\"jammer\\">; x/ m/ I# A0 E) m& h</font><br /> s;<font class=\\"jammer\\">* B  B7 c9 t% J! O5 T</font><br /> [clock()-t0]/1000;

结果：
) g  I& D* D& v) E7 f: U) _, {11008953.66602346
1.031

, p1 \+ N; e& _; t二者效率相差已比较小。
4 N; {, z2 b/ @5 N$ W7 ]3 I" D
3、再增加一个函数调用
4 ~% ~6 S. k5 p* E/ |8 Z' r
! ?8 F0 E$ {) C0 }: E- h/ OMatlab代码：

1. clear all  s9 x: a# r6 }+ a, p& v! l
   
2. f=@(x)x+1;' n8 @2 H; _% E\" l% b% r& z0 ^
   
3. tic7 B. F; L( r+ G0 P
   
4. s=0;3 {$ ?9 s- Y9 i7 P3 a
   
5. for k=1:1000
   ) A5 L$ y9 Z' U8 x7 V
6.   for i=1:1000
   . o& ]( Y$ o/ }, ^  r* j4 s  |
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);
   \" |3 g: d/ u+ p7 z! V
8.   end0 D\" W: Z1 Q6 y
   
9. end) }2 a2 o0 f, W
   
10. s1 e- Z7 V+ C9 J% c
    
11. toc
    : l- u0 K3 u\" F  B  A; X. ^! W
12. \" |4 r  T# e( p' \; H) J7 T
    
13. s =2 m: n- [9 h  h' W4 Z  h7 \* o
    
14. 
    : d3 ?3 h1 Z8 _\" U
15.   5.1251e+008
    8 K) h+ ]. J% ^5 D+ p* N$ N
16. 5 n: l- F3 K( s
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. / ?' v2 ]6 K8 E
4. f(x)=x+1;2 ~2 X8 h* h\\" m$ ]7 r
   
6. mvar:
   
7. * I# [7 X) f$ G% ]
8. t0=clock(),
   
9.   X( z. z  H; |$ m2 H
10. s=0,* `9 e# p' d# W& P! k* S
    
12. k=1, while{k<=1000,
    
13. % \; M) B2 a! i% X' ^8 ]: T
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++},
    
15. 2 h$ |5 ?: S, ^# V
16.   k++) H\\" [  _: b\\" y: \
    
18. },
    
19. 8 g1 {# H# s0 P0 ^6 g$ a7 v
20. s;
    
21. & a; ]2 ?3 _/ j/ y9 T
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">/ ?' v2 ]6 K8 E</span> f(x)=x+1;<font class=\\"jammer\\">2 ~2 X8 h* h\\" m$ ]7 r</font><br /> mvar:<br /> <span style=\\"display:none\\">* I# [7 X) f$ G% ]</span> t0=clock(),<br /> <span style=\\"display:none\\">  X( z. z  H; |$ m2 H</span> s=0,<font class=\\"jammer\\">* `9 e# p' d# W& P! k* S</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">% \; M) B2 a! i% X' ^8 ]: T</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++},<br /> <span style=\\"display:none\\">2 h$ |5 ?: S, ^# V</span>   k++<font class=\\"jammer\\">) H\\" [  _: b\\" y: \</font><br /> },<br /> <span style=\\"display:none\\">8 g1 {# H# s0 P0 ^6 g$ a7 v</span> s;<br /> <span style=\\"display:none\\">& a; ]2 ?3 _/ j/ y9 T</span> [clock()-t0]/1000;

结果：
- V. G2 Z) g: U( {1 B/ d512508953.6658993
8 Q5 b2 X( F' |1 M8 l1.078

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

海水

forcal

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

4 ^. }* c" u! e$ M* nMatlab代码：

1. >> clear all
   ' a1 t# }' D5 r6 _$ n/ F8 P9 T/ {
2. a = linspace(1,1,1000);
   , B5 H( f% T3 O& y
3. tic7 I9 z4 J# F1 q! o
   
4. for k=1:100004 w1 Z\" o+ Q; m\" y# k3 ~
   
5.   for i=1:1000
   5 T+ Q% f\" q/ H5 p( Z
6.     a(i)=2;1 p$ p3 S4 @4 T' y* |
   
7.   end
     [$ N4 Y. t4 I( r
8. end* [) Q2 k# |7 [. d2 y+ ?; b1 k& d& c9 o
   
9. toc
   5 J4 k$ u& r# r6 e0 s( C+ P+ l
10. 
    2 @# f/ g$ N4 J+ o4 O! N, X
11. tic
    1 ~' \) i) J5 B( h
12. s=0;! P( z+ B8 l5 G% I0 X
    
13. for k=1:10000
    ( g/ K$ F, Z0 w  J$ N/ f, ]# H8 J
14.   for i=1:1000: k3 @& h- \2 F
    
15.     s=s+a(i);
    , C+ o1 z  L5 Y8 j
16.   end
    7 K. T5 m+ a+ F3 {
17. end
    % p\" ]: i, m# _\" X0 t1 e
18. s
    : a$ l. }6 c2 f5 i& _, ?' B
19. toc
    : l: t) o' ~! I4 m( S
20. Elapsed time is 0.140658 seconds.
    2 C6 {. H2 H0 W1 M; d4 S
21. 
    % M* X. l; T0 Z3 l' I+ S$ P
22. s =
    * @2 y/ F  K8 \7 {$ k$ L% V; R9 a
23. 
    , _1 I$ F2 B5 a\" X
24.     20000000
    6 J- i/ k$ J- B+ q\" F3 j0 W
25. 
    1 H9 ]5 c+ ~- {) l1 D
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];* R4 T; Q- H! O1 r* q
   
4. mvar:4 _. h6 \( d! E\\" }
   
6. main0(::a)=8 Q: y5 }. u0 F5 N& J, V
   
8. a = linspace(1,1,1000),
   
9. * H: J: _) U% a( m3 D+ b
10. t0=clock(),
    
11. : B4 q, K, k; U6 Y) Q
12. k=0, while{k<10000,4 R: S! |- s% I0 c: N9 u3 y/ |+ C\\" D# J
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. - O- G\\" e+ H; [, E$ M2 @7 \
16.   k++
    
17. + ~# ^) ^1 F$ f. {
18. },
    
19. ( ]9 u$ ^' a2 W& t# U
20. [clock()-t0]/1000;
    
21. 1 M1 c6 W  Z4 C' `9 ^: k
22. 
    
23. 8 N4 W, C% ]$ j) N
24. main1(::a)=
    
25. 6 o) H1 I' _5 q\\" U7 n3 s
26. t0=clock(), s=0,/ n2 c' V2 ]0 d7 Y; F3 w: G) P
    
28. k=0, while{k<10000,& J# H+ r. h! A& P: e
    
30.   i=0, while{i<1000, s=s+a(i), i++},6 W1 U( B7 J. A2 t& M+ T1 p
    
32.   k++
    
33. 3 [9 Y- R9 p* q
34. },% t; x5 \  ]# O
    
36. s;
    
37. 0 N\\" ?1 f( x  Y, K4 S: ]
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">* R4 T; Q- H! O1 r* q</font><br /> mvar:<font class=\\"jammer\\">4 _. h6 \( d! E\\" }</font><br /> main0(::a)=<font class=\\"jammer\\">8 Q: y5 }. u0 F5 N& J, V</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">* H: J: _) U% a( m3 D+ b</span> t0=clock(),<br /> <span style=\\"display:none\\">: B4 q, K, k; U6 Y) Q</span> k=0, while{k<10000,<font class=\\"jammer\\">4 R: S! |- s% I0 c: N9 u3 y/ |+ C\\" D# J</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">- O- G\\" e+ H; [, E$ M2 @7 \</span>   k++<br /> <span style=\\"display:none\\">+ ~# ^) ^1 F$ f. {</span> },<br /> <span style=\\"display:none\\">( ]9 u$ ^' a2 W& t# U</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">1 M1 c6 W  Z4 C' `9 ^: k</span> <br /> <span style=\\"display:none\\">8 N4 W, C% ]$ j) N</span> main1(::a)=<br /> <span style=\\"display:none\\">6 o) H1 I' _5 q\\" U7 n3 s</span> t0=clock(), s=0,<font class=\\"jammer\\">/ n2 c' V2 ]0 d7 Y; F3 w: G) P</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">& J# H+ r. h! A& P: e</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">6 W1 U( B7 J. A2 t& M+ T1 p</font><br />   k++<br /> <span style=\\"display:none\\">3 [9 Y- R9 p* q</span> },<font class=\\"jammer\\">% t; x5 \  ]# O</font><br /> s;<br /> <span style=\\"display:none\\">0 N\\" ?1 f( x  Y, K4 S: ]</span> [clock()-t0]/1000;

结果：
, O: o# E$ I6 P$ {4 X15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒
4 L; V- |1 C* |; q$ O9 J$ I
, ?% N+ B& o9 q; u7 Z+ J, W－－－－－－－－－－
$ q; X" a0 l7 M* d5 {" w1 L& ]
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];: ]\\" }+ h\\" _9 j8 \. R
   
4. mvar:- t# Z$ N; y) }# x# G
   
6. main0(::a)=0 ]# F; D$ I1 i+ v
   
8. a = linspace(1,1,1000),* O# e: `, X, K. N0 t6 p$ Q
   
10. t0=clock(),
    
11. ; I6 P( M! V7 ~5 X
12. k=0, while{k<10000,
    
13. - |6 N+ `& N% P3 ]; V
14.   i=0, while{i<1000, A(a,i)=2, i++},+ f, m: c+ R* i8 n0 F4 m\\" N
    
16.   k++
    
17. 3 s  a% s0 e9 C: N$ k+ z; s$ i
18. },$ ?( q: Z! o/ g, A2 Z* o
    
20. [clock()-t0]/1000;
    
21. % u5 @1 i4 _' @) T% B; |
22. # k6 @; S6 H, k8 M
    
24. main1(::a)=
    
25. % b) a, e2 O# h3 v
26. t0=clock(), s=0,3 X\\" q, g+ Z( a! F1 r
    
28. k=0, while{k<10000,
    
29. 1 c7 v: \  T: E; Y& [
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. 6 ]# Z; ]\\" r  b. `& R1 H: g7 E
32.   k++
    
33. 4 C# V  L+ c$ @
34. },& z1 M% U9 h% Y9 A
    
36. s;5 S; z' p  j( S# }; V* H0 o
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">: ]\\" }+ h\\" _9 j8 \. R</font><br /> mvar:<font class=\\"jammer\\">- t# Z$ N; y) }# x# G</font><br /> main0(::a)=<font class=\\"jammer\\">0 ]# F; D$ I1 i+ v</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">* O# e: `, X, K. N0 t6 p$ Q</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">; I6 P( M! V7 ~5 X</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">- |6 N+ `& N% P3 ]; V</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">+ f, m: c+ R* i8 n0 F4 m\\" N</font><br />   k++<br /> <span style=\\"display:none\\">3 s  a% s0 e9 C: N$ k+ z; s$ i</span> },<font class=\\"jammer\\">$ ?( q: Z! o/ g, A2 Z* o</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">% u5 @1 i4 _' @) T% B; |</span> <font class=\\"jammer\\"># k6 @; S6 H, k8 M</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">% b) a, e2 O# h3 v</span> t0=clock(), s=0,<font class=\\"jammer\\">3 X\\" q, g+ Z( a! F1 r</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">1 c7 v: \  T: E; Y& [</span>   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">6 ]# Z; ]\\" r  b. `& R1 H: g7 E</span>   k++<br /> <span style=\\"display:none\\">4 C# V  L+ c$ @</span> },<font class=\\"jammer\\">& z1 M% U9 h% Y9 A</font><br /> s;<font class=\\"jammer\\">5 S; z' p  j( S# }; V* H0 o</font><br /> [clock()-t0]/1000;

结果：
1 h4 a, o% Z6 W: B1 p9 v10.469
& [7 L9 H  d# ^2 A! a/ J20000000.
3.781
& M0 o1 @( W. v0 V4 U
0 R6 x& C  q5 ?/ u) M0 e( F; Q－－－－－－－

, s* W! r4 S' f" ^在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
# c' F3 ?; b0 a* s在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
# v3 x7 {3 F  |
: Q" y2 n7 L$ W9 y. D6 H6 Rmatlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
* S6 `* I( r  V+ x2 T
3 ?9 {6 V' I/ G: r& R% y5、在循环中加点东西
$ Q; j: J4 @1 W! U2 b9 ^
matlab代码：

1. >> clear all
   9 |. i, F' F* D: T6 R: t. T' M
2. f=@(x)x+1;
   7 ~; E\" E$ R5 C$ J% E
3. a = linspace(2,2,1000);! y  W: s$ l& ~! V0 F$ A
   
4. tic
   ' o% P& D! F1 f
5. s=0;/ t+ V) o' X6 ~5 A+ |3 r, c& G
   
6. for k=1:1000
   * @, F\" T, D\" I* y! {: e  z  s
7.   for i=1:1000
   \" _# O' F2 A+ o' _
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));/ a1 r' B3 r$ t* H+ b% M6 s, b
   
9.   end6 ~8 L# T( r4 S- y\" s0 ?\" D& L
   
10. end4 Q, P+ E5 o) i2 m1 a
    
11. s
    4 J; e: A$ U; M' Z$ R; K
12. toc% o% c( [# {! d8 g+ E! \% i
    
13. 3 F5 D- f, R; Z4 ^6 f  |5 G4 n
    
14. s =9 Z- t6 o$ s# T8 p- Q( @. w
    
15. 7 _+ h2 J) D2 h0 k
    
16.   2.4002e+007
    ' c' h3 v% a; g. @  L
17. 
    ' P5 \% m9 x1 x8 F; ]: |6 w
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];5 i( Y: b2 P7 w5 p
   
4. mvar:
   
5. 5 [; C: {2 F' ^6 D& v& g
6. f(x)=x+1;
   
7. 9 H: P& M1 ]/ S' i6 v
8. main0(::a)=
   
9. ; r. x) B3 h. q( U$ v+ i
10. a = linspace(2,2,1000),
    
11. 3 N. Y  U4 |$ _! h& }. g
12. t0=clock(), s=0,
    
13. % y) @1 ^; q+ I
14. k=0, while{k<1000,
    
15. ) B# n% U8 W7 \\\" V5 t# A7 U
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++},% N8 \. Z3 Q- |6 f, o
    
18.   k++
    
19. ! i: q! W; B2 h8 [
20. },
    
21. ; ]# Z- R, n) q/ I2 M- A! D
22. s;
    
23. ( F* U8 ?8 ~+ d  N
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">5 i( Y: b2 P7 w5 p</font><br /> mvar:<br /> <span style=\\"display:none\\">5 [; C: {2 F' ^6 D& v& g</span> f(x)=x+1;<br /> <span style=\\"display:none\\">9 H: P& M1 ]/ S' i6 v</span> main0(::a)=<br /> <span style=\\"display:none\\">; r. x) B3 h. q( U$ v+ i</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">3 N. Y  U4 |$ _! h& }. g</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">% y) @1 ^; q+ I</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">) B# n% U8 W7 \\\" V5 t# A7 U</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\\">% N8 \. Z3 Q- |6 f, o</font><br />   k++<br /> <span style=\\"display:none\\">! i: q! W; B2 h8 [</span> },<br /> <span style=\\"display:none\\">; ]# Z- R, n) q/ I2 M- A! D</span> s;<br /> <span style=\\"display:none\\">( F* U8 ?8 ~+ d  N</span> [clock()-t0]/1000;

结果：
24002271.69189384
5.172

- ~. f4 C- t0 T" J2 `+ u. O二者效率现在只有一倍差距。
# {2 P  ?+ K3 W; n! M/ i
－－－－－－－－－－

8 K& |8 X4 l  R! o0 @不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

1. clear all* ]' n0 h8 J& r' }! i6 d
   
2. f=@(x)x+1;
   6 h( ^9 ?9 B; w8 u
3. a = linspace(2,2,1000);/ G6 P- B4 ?& e. ^; [% E8 v* j
   
4. tic
   9 t# V  w- L/ h& e
5. s=0;
   2 j% D; \  C* v; ]2 [; Z+ I
6. for k=1:1000
   & P7 H# @; W9 Y( m
7.   for i=1:10008 z/ t4 {! U0 R4 e* `\" D\" r6 t
   
8.     t=a(i);
   % W9 b; h& y- p! c3 y4 B8 x. f$ ~
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 a9 A, Y+ D7 @4 J; i. f# u
10.   end( F6 I2 S* {$ E+ t' T  t' ~
    
11. end
    7 m- M: ^3 i/ Z) l, A\" S: p
12. s4 X' D6 t; B# q
    
13. toc
    % w3 L8 Q( z2 p: M. j
14. 
    & [- O, Q; n! y% [) y
15. s =& g8 l$ a* s+ \9 i' p. L
    
16. / m. L' J. K  S$ \; ^  r# m
    
17.   2.4002e+0071 H0 R' d, o4 m& h
    
18. 2 l) G) m. q  b; L( \. g1 F
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. 5 O* X) |) a* Y3 ?- W! @
4. mvar:# p8 R* z' J. B$ d; p& C
   
6. f(x)=x+1;
   
7. 8 E& i& w6 d! H: G+ y# f
8. main0(::a)=
   
9. 2 P0 T0 m, `4 j2 u
10. a = linspace(2,2,1000),  d8 r2 }: e! E
    
12. t0=clock(), s=0,
    
13. 4 ^- O4 p6 D5 f) j0 I+ u
14. k=0, while{k<1000,
    
15. \\" }2 k3 a1 M7 S7 T
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 X% l. {\\" D% h. R8 R
18.   k++
    
19. $ e4 H) m/ \  W* ?
20. },& {) l3 X, B9 C8 T: U6 t# v
    
22. s;
    
23. \\" X. h# s1 |) P; _2 X5 `
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">5 O* X) |) a* Y3 ?- W! @</span> mvar:<font class=\\"jammer\\"># p8 R* z' J. B$ d; p& C</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">8 E& i& w6 d! H: G+ y# f</span> main0(::a)=<br /> <span style=\\"display:none\\">2 P0 T0 m, `4 j2 u</span> a = linspace(2,2,1000),<font class=\\"jammer\\">  d8 r2 }: e! E</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">4 ^- O4 p6 D5 f) j0 I+ u</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">\\" }2 k3 a1 M7 S7 T</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\\">7 X% l. {\\" D% h. R8 R</span>   k++<br /> <span style=\\"display:none\\">$ e4 H) m/ \  W* ?</span> },<font class=\\"jammer\\">& {) l3 X, B9 C8 T: U6 t# v</font><br /> s;<br /> <span style=\\"display:none\\">\\" X. h# s1 |) P; _2 X5 `</span> [clock()-t0]/1000;

结果：
24002271.69189384
, s# e1 c$ h7 A& [1.625
! M" H$ h% M9 s* m" `0 @/ N6 k
$ I; k9 T1 r( F8 mmatlab效率已不及Forcal。
: `+ P% k4 z8 J1 @* j2 M* C
) @% j" _- I. m8 B/ c' u提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
* t# E* _& A7 ^- G4 O0 |
# n7 g( Z7 t) L+ v( ]6 @1 ?Forcal代码：

2. !using["sys"];
   
3. / j% K, `8 ~( U: v3 J1 s) y
4. mvar:5 @1 t7 d# \( x! L
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组) {* X; P: a. p) v; k4 G
   
8. t0=clock(),
   
9. 9 ?\\" D9 v: |. X' d; s7 }+ U7 H
10. k=0, while{k<10000,
    
11. \\" A( [\\" _$ D& m0 E6 j
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值\\" n; F& {+ \) w7 z9 G! H0 |
    
14.   k++
    
15. 6 A: A$ j( o# `\\" u; O: V
16. },( x\\" R3 ]. O* X: h3 x\\" K: H\\" ]
    
18. [clock()-t0]/1000;5 q1 V4 I! m. }  b! o0 P- h
    
20. . S# t) o2 N- X% L) N8 H/ \
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组) ^2 S+ Y# |. T  ?' l
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值% R\\" I2 W\\" r6 o' i. U
    
26. t0=clock(),- J+ v1 f1 J0 c, j# V) ?
    
28. s=0,! N6 F  O& O$ X; M
    
30. k=0, while{k<10000,& J; v# \; W6 N+ ~: E
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. 1 `: R& s. V1 D- j
34.   k++\\" G2 V  X; w. l* H# _& c3 _/ `
    
36. },
    
37. 8 o5 u$ w, Z\\" H' h8 ^8 {6 q: `
38. s;
    
39. - e7 s$ D* [5 ]& N: M
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">/ j% K, `8 ~( U: v3 J1 s) y</span> mvar:<font class=\\"jammer\\">5 @1 t7 d# \( x! L</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">) {* X; P: a. p) v; k4 G</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">9 ?\\" D9 v: |. X' d; s7 }+ U7 H</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">\\" A( [\\" _$ D& m0 E6 j</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">\\" n; F& {+ \) w7 z9 G! H0 |</font><br />   k++<br /> <span style=\\"display:none\\">6 A: A$ j( o# `\\" u; O: V</span> },<font class=\\"jammer\\">( x\\" R3 ]. O* X: h3 x\\" K: H\\" ]</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">5 q1 V4 I! m. }  b! o0 P- h</font><br /> <font class=\\"jammer\\">. S# t) o2 N- X% L) N8 H/ \</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">) ^2 S+ Y# |. T  ?' l</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">% R\\" I2 W\\" r6 o' i. U</font><br /> t0=clock(),<font class=\\"jammer\\">- J+ v1 f1 J0 c, j# V) ?</font><br /> s=0,<font class=\\"jammer\\">! N6 F  O& O$ X; M</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">& J; v# \; W6 N+ ~: E</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">1 `: R& s. V1 D- j</span>   k++<font class=\\"jammer\\">\\" G2 V  X; w. l* H# _& c3 _/ `</font><br /> },<br /> <span style=\\"display:none\\">8 o5 u$ w, Z\\" H' h8 ^8 {6 q: `</span> s;<br /> <span style=\\"display:none\\">- e7 s$ D* [5 ]& N: M</span> [clock()-t0]/1000;

结果：
: Z# d. D; k" U7 ^% |. g- i. c9 q0.5620000000000001  //赋值用时间，秒
2 N  t2 x9 g# g, y1 d% I, }20000000.
0.579      //取值用时间，秒