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

forcal

1、简单循环

4 [' ]& p! w! g6 ?5 UMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all  H, B, c! L% d& \4 Q9 [  j
   
2. tic: G  K2 l) M* Z0 `# ~0 I+ e
   
3. s=0;
     ^& K# s8 w% X( R
4. for k=1:10000
   ) \  w4 o+ I9 O& R* R9 ^
5.   for i=1:1000
   + ]7 p\" }9 z4 U1 Y\" M
6.     s=s+1;
   $ o6 Q* R9 ^  ?6 O6 O( e
7.   end& c1 o) T: L0 T+ c9 {/ M
   
8. end
   . m$ C+ O6 _# h
9. s& _& Q9 T* w5 s' C1 R
   
10. toc
    % y' ~! T4 V# b
11. 
    / M5 p: g6 e* ?& g% Y7 B' x
12. s =
    9 y  I- E3 e  q, `
13. 
    - P( ~1 d3 `$ R7 ?/ \0 b
14.     10000000% f# {& J1 |2 s* Y- U( u
    
15. 
    . g8 a7 P0 ^\" b9 ^4 \& s
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
: I4 v' _, T7 e- B) _5 g
Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. - V% b2 z) i. ?  P9 a9 A, L3 |
4. mvar:
   
5. ) w' E1 R  Z  W$ B- d
6. t0=clock(),
   
7. ! @) |: y, z$ s  T2 W. v% E2 d
8. s=0,5 [. ?4 A8 F/ W) t
   
10. k=0, while{k<10000,
    
11. ! h7 L+ m6 Q7 e% k* u* L* R4 h
12.   i=0, while{i<1000,s=s+1, i++},) O0 [3 G( G- @# o& Y) u+ c
    
14.   k++
    
15. 2 g3 d6 }; z  u0 e# z5 \: X
16. },
    
17. ' D9 O  H# B& ^8 r, F6 S' g* k6 C  |
18. s;
    
19. ; m\\" N: Z1 ^; m; s$ b% ~5 a2 o
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">- V% b2 z) i. ?  P9 a9 A, L3 |</span> mvar:<br /> <span style=\\"display:none\\">) w' E1 R  Z  W$ B- d</span> t0=clock(),<br /> <span style=\\"display:none\\">! @) |: y, z$ s  T2 W. v% E2 d</span> s=0,<font class=\\"jammer\\">5 [. ?4 A8 F/ W) t</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">! h7 L+ m6 Q7 e% k* u* L* R4 h</span>   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">) O0 [3 G( G- @# o& Y) u+ c</font><br />   k++<br /> <span style=\\"display:none\\">2 g3 d6 }; z  u0 e# z5 \: X</span> },<br /> <span style=\\"display:none\\">' D9 O  H# B& ^8 r, F6 S' g* k6 C  |</span> s;<br /> <span style=\\"display:none\\">; m\\" N: Z1 ^; m; s$ b% ~5 a2 o</span> [clock()-t0]/1000;

结果：
10000000.
0.704  秒
" m' V* E5 r$ x- z8 D
二者运行效率相差7倍多。
' A6 z9 P) @+ k: m- Y
" ]' ?, i; |7 n, m2、循环中稍微加点东西

Matlab代码：

1. clear all. M! {! ?! {/ L, a: m6 [
   
2. tic5 x! t1 s, }$ _0 c, r
   
3. s=0;; I& t. W$ T0 S+ y# U! p
   
4. for k=1:1000
   4 n! ]) n\" u! [\" p  m
5.   for i=1:1000
   . t* u; I/ O- _2 R
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);
   6 K* @8 z4 f* l! y
7.   end
   ; k- |8 t% D  Z8 j
8. end2 X5 A+ P4 J8 r# v) Q3 }
   
9. s
   * K  N4 C% o% P& P, N2 x
10. toc
    ! O( k( `& ]! j9 t9 I
11. 
    1 H7 f( L8 H/ o, z$ f0 V. j6 x
12. s =. R& B2 M% l6 ^( B, ~: Q
    
13. 
    - v7 {\" o+ y  W  t
14.   1.1009e+007
      L$ {( ?  c3 ]# P+ v0 a8 l
15. 
    ! R( U7 P8 {4 y3 F. e7 k
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 8 h0 ^0 ?# E, U8 N
4. mvar:; L+ m4 L# j$ o1 I' B
   
6. t0=clock(),, Y# |- l7 ^. s; M* M$ T
   
8. s=0,. n\\" ^: f; j) P. [
   
10. k=1, while{k<=1000,
    
11.   f* N/ q; h' f9 B2 |' s
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. 7 T5 ?3 h6 `+ I! `& ~- `$ P: s
14.   k++/ V, K, M6 h( m1 t; O
    
16. },
    
17. 3 v$ s  ^- [' x8 r
18. s;
    
19. + ^\\" Z  u# y\\" Y- H9 i2 a# o
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">8 h0 ^0 ?# E, U8 N</span> mvar:<font class=\\"jammer\\">; L+ m4 L# j$ o1 I' B</font><br /> t0=clock(),<font class=\\"jammer\\">, Y# |- l7 ^. s; M* M$ T</font><br /> s=0,<font class=\\"jammer\\">. n\\" ^: f; j) P. [</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">  f* N/ q; h' f9 B2 |' s</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\\">7 T5 ?3 h6 `+ I! `& ~- `$ P: s</span>   k++<font class=\\"jammer\\">/ V, K, M6 h( m1 t; O</font><br /> },<br /> <span style=\\"display:none\\">3 v$ s  ^- [' x8 r</span> s;<br /> <span style=\\"display:none\\">+ ^\\" Z  u# y\\" Y- H9 i2 a# o</span> [clock()-t0]/1000;

结果：
& Y' w7 ]8 L7 D1 u11008953.66602346
1.031
) c5 k8 _* `9 G! T
: w* p- z3 y$ ~7 G% z& `2 d3 r二者效率相差已比较小。

3、再增加一个函数调用

4 Q/ ~) v, q; x) ?7 RMatlab代码：

1. clear all; k/ j5 s+ c0 V/ Y
   
2. f=@(x)x+1;% C, H! f5 `& y7 h+ c
   
3. tic
   8 L  Y1 p( b6 f8 q0 A
4. s=0;
   0 e& e% |- G) ~/ r+ L\" U
5. for k=1:1000- x6 N, q* B- r0 l% [
   
6.   for i=1:1000
   ! ^! G* c7 J( {0 ]. f9 j
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);
   ! |* z3 f4 x& K
8.   end  b) C+ ?3 D5 U# J% W
   
9. end
     L$ O. |4 _3 [& F) I$ E3 b5 V
10. s0 o( j0 i; @. y, I' M( r6 y
    
11. toc
    6 f: b4 Q% u) B: k6 X# y
12. 
    6 f4 o  K; G; X
13. s =
    # u/ N+ `8 S! p# b7 d  f
14. 9 C0 S0 b/ D% f
    
15.   5.1251e+008\" Y* ]8 ~6 R5 u+ ]
    
16. : ~. e( j2 b% s( ?; N$ Y
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 2 }0 n. P3 |- @& ?
4. f(x)=x+1;& w\\" M7 j+ e2 w- I
   
6. mvar:9 U0 _; F7 i: q
   
8. t0=clock(),
   
9. + a( c. R3 ~9 }+ M! Z
10. s=0,9 T0 v5 }( d! x/ i/ {8 N
    
12. k=1, while{k<=1000,
    
13. % N% g5 s# P% \0 E- A; y
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. ! W; Q# C4 W( z) \\\" Y7 i
16.   k++! }8 @7 G! s; X) q/ t0 @- t\\" x8 _
    
18. },- _1 `+ L3 Z4 A\\" y
    
20. s;
    
21. ! f- D4 z2 ]: x5 {2 b
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">2 }0 n. P3 |- @& ?</span> f(x)=x+1;<font class=\\"jammer\\">& w\\" M7 j+ e2 w- I</font><br /> mvar:<font class=\\"jammer\\">9 U0 _; F7 i: q</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">+ a( c. R3 ~9 }+ M! Z</span> s=0,<font class=\\"jammer\\">9 T0 v5 }( d! x/ i/ {8 N</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">% N% g5 s# P% \0 E- A; y</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\\">! W; Q# C4 W( z) \\\" Y7 i</span>   k++<font class=\\"jammer\\">! }8 @7 G! s; X) q/ t0 @- t\\" x8 _</font><br /> },<font class=\\"jammer\\">- _1 `+ L3 Z4 A\\" y</font><br /> s;<br /> <span style=\\"display:none\\">! f- D4 z2 ]: x5 {2 b</span> [clock()-t0]/1000;

结果：
512508953.6658993
0 d: A& w* T; i2 Y' J6 ~1.078

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

海水

forcal

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

) P' J7 u0 R+ t3 P+ V) p, u; NMatlab代码：

1. >> clear all/ ?+ o& _\" p; V
   
2. a = linspace(1,1,1000);
   - S/ {3 A) {# L/ g/ F& N) b6 _
3. tic  C1 X$ B# h/ \
   
4. for k=1:10000
   9 ^* k8 f  ~. U( d2 q
5.   for i=1:1000( X1 k. Z9 b* d6 k7 L  Z8 S: V
   
6.     a(i)=2;9 \+ ^) l* P% F# r0 W; Z6 @
   
7.   end
   ! h# b, @& w4 o; \
8. end
     @& g& I- v# E7 J$ f- x
9. toc( w7 ^\" A: `1 e- Q
   
10. 
    & B' v  ^: W% t3 y$ e0 ?- c7 t1 L
11. tic
    . v2 H6 D' ]; L' `. l! P3 r( o) p$ E
12. s=0;2 c4 ?8 B) b! ]
    
13. for k=1:100006 |6 }. i0 p3 b$ T
    
14.   for i=1:10007 ~2 A1 z3 {6 a: @9 k. G
    
15.     s=s+a(i);- W: u( E, q. |8 B
    
16.   end6 n* `0 b: Z* }. x! X( j
    
17. end
      Q$ J% ^: h0 r2 }5 {
18. s! }  g# ?5 o, d. X  z4 v
    
19. toc
    6 C\" [4 }$ \& }3 [# y: C! [
20. Elapsed time is 0.140658 seconds.5 O8 Y* i3 [5 x* J7 a; P$ o. X* m
    
21. 
    6 z  I( c+ b) L$ ~  E' y
22. s =  V( \9 y: a6 m7 a: c
    
23. . _' \7 S1 e6 l9 T% Z
    
24.     20000000
    4 H3 t2 n4 o% }- i+ K! ?
25. , U7 n  i: S9 y6 X& L+ h
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];6 j* s4 B4 d+ ^- W9 f
   
4. mvar:8 |9 ?& P\\" E  P4 e
   
6. main0(::a)=
   
7. % T! w0 u9 R  \' ~* [
8. a = linspace(1,1,1000),. ^  I$ X+ A/ Q  M
   
10. t0=clock(),' U/ }4 S. h4 J; v
    
12. k=0, while{k<10000,
    
13. ! j2 H) e7 f1 _
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. 9 J/ Q( k' l4 N8 c6 y6 p* y0 T
16.   k++
    
17. / S8 r8 @2 M) A+ Q+ e* G
18. },
    
19. 9 t* p* W/ q+ X% ?% T) }  y' x
20. [clock()-t0]/1000;
    
21. * U* }2 k4 t5 k) f7 V/ j
22. 
    
23. 9 ?( n- u* Z% I7 ]( u  \; N
24. main1(::a)=9 c- n5 d- W; j. [
    
26. t0=clock(), s=0,4 C2 s& \5 L' e# w3 w
    
28. k=0, while{k<10000,/ E. A- _% M2 j
    
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. $ x1 n/ F8 ^\\" O2 H* [7 U# c
32.   k++
    
33. , [* _: F& ~5 I6 Q! |
34. },$ S' Y# X, P8 E: n% t- ?\\" f
    
36. s;( K) E% X$ q# [  @
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">6 j* s4 B4 d+ ^- W9 f</font><br /> mvar:<font class=\\"jammer\\">8 |9 ?& P\\" E  P4 e</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">% T! w0 u9 R  \' ~* [</span> a = linspace(1,1,1000),<font class=\\"jammer\\">. ^  I$ X+ A/ Q  M</font><br /> t0=clock(),<font class=\\"jammer\\">' U/ }4 S. h4 J; v</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">! j2 H) e7 f1 _</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">9 J/ Q( k' l4 N8 c6 y6 p* y0 T</span>   k++<br /> <span style=\\"display:none\\">/ S8 r8 @2 M) A+ Q+ e* G</span> },<br /> <span style=\\"display:none\\">9 t* p* W/ q+ X% ?% T) }  y' x</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">* U* }2 k4 t5 k) f7 V/ j</span> <br /> <span style=\\"display:none\\">9 ?( n- u* Z% I7 ]( u  \; N</span> main1(::a)=<font class=\\"jammer\\">9 c- n5 d- W; j. [</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">4 C2 s& \5 L' e# w3 w</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">/ E. A- _% M2 j</font><br />   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">$ x1 n/ F8 ^\\" O2 H* [7 U# c</span>   k++<br /> <span style=\\"display:none\\">, [* _: F& ~5 I6 Q! |</span> },<font class=\\"jammer\\">$ S' Y# X, P8 E: n% t- ?\\" f</font><br /> s;<font class=\\"jammer\\">( K) E% X$ q# [  @</font><br /> [clock()-t0]/1000;

结果：
) P4 a5 G( J+ x  `6 [- ?15.484  //存数据所用时间，秒
+ [' X) i; Y( I20000000.
5.531   //取数据所用时间，秒
  i6 z$ j" O" s6 O- c& ]
7 k- I+ U: r; E% {/ [3 b－－－－－－－－－－
* }7 D: |2 a6 A7 |. l, p
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. 2 w3 f6 A: j5 k  J) w
4. mvar:
   
5. % v; L* z0 A( K, Q8 A% B
6. main0(::a)=$ w  Q- Q- O- ?  A+ W
   
8. a = linspace(1,1,1000),
   
9. 2 d9 m0 p7 ]) S. u$ f# p: ?5 `, @
10. t0=clock(),
    
11. 2 w; \% V4 f9 p0 q2 G) O1 a* @
12. k=0, while{k<10000,* S8 s* @$ p7 v% m
    
14.   i=0, while{i<1000, A(a,i)=2, i++},$ m* j0 W8 q4 v3 G% `
    
16.   k++1 k* W. v  g% ^# P; J, `\\" v
    
18. },; [; c- n\\" ], o. [- D* H
    
20. [clock()-t0]/1000;% u: Y, I/ l7 Y; o3 A/ P
    
22. : h& F: I8 f( Y( u
    
24. main1(::a)=
    
25. . g* v0 Q: M' Z0 O. V) ]' a, Y
26. t0=clock(), s=0,6 ^) X7 r( w\\" r
    
28. k=0, while{k<10000,
    
29. , _5 P5 J/ W: O' \! [
30.   i=0, while{i<1000, s=s+A(a,i), i++},3 S6 a+ F* O( E, v, S
    
32.   k++3 q8 B8 Z0 F; U; \/ b: p; O
    
34. },- ?) L/ f1 K; L2 Z$ p6 h7 U& j' w9 k( y
    
36. s;: @8 J$ z% ~# N* d3 Q  ]4 y- U
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">2 w3 f6 A: j5 k  J) w</span> mvar:<br /> <span style=\\"display:none\\">% v; L* z0 A( K, Q8 A% B</span> main0(::a)=<font class=\\"jammer\\">$ w  Q- Q- O- ?  A+ W</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">2 d9 m0 p7 ]) S. u$ f# p: ?5 `, @</span> t0=clock(),<br /> <span style=\\"display:none\\">2 w; \% V4 f9 p0 q2 G) O1 a* @</span> k=0, while{k<10000,<font class=\\"jammer\\">* S8 s* @$ p7 v% m</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">$ m* j0 W8 q4 v3 G% `</font><br />   k++<font class=\\"jammer\\">1 k* W. v  g% ^# P; J, `\\" v</font><br /> },<font class=\\"jammer\\">; [; c- n\\" ], o. [- D* H</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">% u: Y, I/ l7 Y; o3 A/ P</font><br /> <font class=\\"jammer\\">: h& F: I8 f( Y( u</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">. g* v0 Q: M' Z0 O. V) ]' a, Y</span> t0=clock(), s=0,<font class=\\"jammer\\">6 ^) X7 r( w\\" r</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">, _5 P5 J/ W: O' \! [</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">3 S6 a+ F* O( E, v, S</font><br />   k++<font class=\\"jammer\\">3 q8 B8 Z0 F; U; \/ b: p; O</font><br /> },<font class=\\"jammer\\">- ?) L/ f1 K; L2 Z$ p6 h7 U& j' w9 k( y</font><br /> s;<font class=\\"jammer\\">: @8 J$ z% ~# N* d3 Q  ]4 y- U</font><br /> [clock()-t0]/1000;

结果：
10.469
" B7 X: n0 z+ ?3 p" {, d20000000.
3.781
6 G5 @$ ]3 H' j0 s1 o
－－－－－－－
/ N. X* ]/ L/ c! V% `% |* ~
在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
) y( s" L' \# I1 C+ n# i( o7 S# M
5 _0 H- n) W( j' ^matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
. m; c4 t: {- V& D$ }
Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

  G; Q$ u/ v- v) {; k* j1 \1 D9 m# J5、在循环中加点东西
1 m% p: |" {/ A2 H+ \, A3 q
6 f" z# b$ b) X6 y0 u3 Z' zmatlab代码：

1. >> clear all# n% S$ _! E- v' f\" F# q& j& o0 b
   
2. f=@(x)x+1;
     `! f4 }6 p& k9 e' K
3. a = linspace(2,2,1000);/ B! K# I6 z2 N- v
   
4. tic( t, [6 k* H* D$ c8 L
   
5. s=0;
   2 x( R* {) h+ w+ t7 ]( d7 Y1 n$ k3 [
6. for k=1:1000( W! ~7 _; H$ r
   
7.   for i=1:1000
   ; ]9 \6 I6 K# s& l- F7 S6 K& A
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));2 ?& t: D% J- V% j8 T4 K% v
   
9.   end6 W5 J4 s% H* X
   
10. end
    * ]0 q9 ^+ }# c; A+ {+ }* k\" }\" ~  f
11. s/ i5 v3 g2 o2 \9 u7 t, Q5 \
    
12. toc
    ( K7 J$ Z5 z6 w1 z: P
13. 
    ( |: R* T! F% h
14. s =
    3 k( ]7 K% ^3 J  @7 y8 B
15. 
    8 O, N, U6 U2 U8 ^0 U% A
16.   2.4002e+007
    6 M) S% q7 o6 a! L! O1 E6 j6 E
17. % ^! a\" M7 F7 o* l' G6 F
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];+ j( v3 w/ o) Y$ t! S# d5 n: [6 q: n
   
4. mvar:. Q) W( q5 @- E
   
6. f(x)=x+1;
   
7. % E% k' x: v& Y& U\\" s
8. main0(::a)=# u4 k6 Q  i1 _3 M: n! i9 ^
   
10. a = linspace(2,2,1000),6 q\\" I: w- G# o5 ^6 |
    
12. t0=clock(), s=0,' W/ T: Z7 e) W2 j( V# Q; G
    
14. k=0, while{k<1000,
    
15. . \9 }8 T* Q6 b& \) A* 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++},2 s% t$ l) U$ G; a3 _4 P' m
    
18.   k+++ ^5 ~1 I  u6 T2 y% t# a8 j8 m4 A* P
    
20. },# u6 P, k# H, c- s
    
22. s;
    
23. % h6 q/ C, Q. I) [9 S' \# X9 g7 y. l
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ j( v3 w/ o) Y$ t! S# d5 n: [6 q: n</font><br /> mvar:<font class=\\"jammer\\">. Q) W( q5 @- E</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">% E% k' x: v& Y& U\\" s</span> main0(::a)=<font class=\\"jammer\\"># u4 k6 Q  i1 _3 M: n! i9 ^</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">6 q\\" I: w- G# o5 ^6 |</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">' W/ T: Z7 e) W2 j( V# Q; G</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">. \9 }8 T* Q6 b& \) A* g</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\\">2 s% t$ l) U$ G; a3 _4 P' m</font><br />   k++<font class=\\"jammer\\">+ ^5 ~1 I  u6 T2 y% t# a8 j8 m4 A* P</font><br /> },<font class=\\"jammer\\"># u6 P, k# H, c- s</font><br /> s;<br /> <span style=\\"display:none\\">% h6 q/ C, Q. I) [9 S' \# X9 g7 y. l</span> [clock()-t0]/1000;

结果：
# y# n6 E5 {0 c24002271.69189384
5.172
6 [1 \3 W! i1 S- f: K
( z$ C# S* U5 J' r8 N( v二者效率现在只有一倍差距。

－－－－－－－－－－
: v/ U, r% V4 H5 |- Y1 v
5 C6 t! X5 Y# D" ~不过，显然matlab和Forcal的程序应该这样写：
" |) l2 _' W" m; ~
matlab程序：

1. clear all
   / p\" e3 b) C& i' _( V
2. f=@(x)x+1;; \3 H2 v$ w+ {. b% ~1 L1 u0 Y
   
3. a = linspace(2,2,1000);
   4 i\" q. J. D; @7 p* w  U
4. tic
   5 m. C\" h\" {2 Q9 {4 T0 D$ g
5. s=0;0 U3 z+ O5 K7 m  ~, d) O4 S/ s3 q
   
6. for k=1:1000, Z& W! a. l' g/ L6 H
   
7.   for i=1:1000
   3 I* N( J5 ^' S! s( d
8.     t=a(i);& \3 O6 c1 u5 _* X
   
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& a# P9 j4 _+ p
10.   end
    4 m  a' N8 e5 \
11. end$ @\" |/ g$ g% E. U
    
12. s
    ) J( B+ Y& h' d6 C; m; `6 J1 j
13. toc
    $ H8 d. l6 x) w* C  B
14. 
    3 |- s: u; \0 [8 d6 j. Y
15. s =
    : m2 c\" {$ v3 v5 V9 T
16. & r5 d' o6 ^3 r2 [
    
17.   2.4002e+007
    0 X5 ], v, V$ g# w
18. 
    * x\" l' K% {7 U: N
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. 7 G& t; N# N0 [: ~  F
4. mvar:- L) ?; F* D  M; C9 d& V  R
   
6. f(x)=x+1;
   
7. 5 Q- z1 D9 p5 i8 P
8. main0(::a)=
   
9. # R  C' A$ y2 ~- a\\" ^+ e
10. a = linspace(2,2,1000),( |# ?, D& m8 A. ~7 R9 D1 o
    
12. t0=clock(), s=0,$ W# O  n: K' `) m2 \) w( z
    
14. k=0, while{k<1000,/ j9 s; o& ?6 L* y) z
    
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. . P/ T( V# K1 U( O
18.   k++
    
19. . O5 E5 K3 w/ w
20. },
    
21. 2 |\\" J' C2 W; N0 a6 i
22. s;  T! {9 ^& q- o4 j( c8 M
    
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">7 G& t; N# N0 [: ~  F</span> mvar:<font class=\\"jammer\\">- L) ?; F* D  M; C9 d& V  R</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">5 Q- z1 D9 p5 i8 P</span> main0(::a)=<br /> <span style=\\"display:none\\"># R  C' A$ y2 ~- a\\" ^+ e</span> a = linspace(2,2,1000),<font class=\\"jammer\\">( |# ?, D& m8 A. ~7 R9 D1 o</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">$ W# O  n: K' `) m2 \) w( z</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">/ j9 s; o& ?6 L* y) z</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\\">. P/ T( V# K1 U( O</span>   k++<br /> <span style=\\"display:none\\">. O5 E5 K3 w/ w</span> },<br /> <span style=\\"display:none\\">2 |\\" J' C2 W; N0 a6 i</span> s;<font class=\\"jammer\\">  T! {9 ^& q- o4 j( c8 M</font><br /> [clock()-t0]/1000;

结果：
0 Q8 P2 Y: M7 z$ B24002271.69189384
1.625
" }1 ^/ q/ T) V" Y- C$ s* f2 l9 F
/ O# L, ^7 M! l. I. ]3 j5 L" nmatlab效率已不及Forcal。
% I% c8 F6 Y% F) w, I6 r! e1 }
提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

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

Forcal代码：

2. !using["sys"];' Q5 G7 ^: m\\" B) n$ ^; A
   
4. mvar:3 v( F& {# t% g  t4 O
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. 0 }6 q: P  v- ?8 c+ y$ q
8. t0=clock(),
   
9. 8 k& u# X# X+ f7 B% d* C
10. k=0, while{k<10000,( J6 R9 T7 q- K' s+ Q# d
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值* |6 K/ x3 i  _. }: D; @
    
14.   k++$ B/ c1 E/ u* y' a3 }: l; E
    
16. },
    
17. \\" W* w' X5 v3 a8 h- L: X+ i: m
18. [clock()-t0]/1000;
    
19. . O3 l$ V9 f6 Q/ C0 ^9 b. ]
20. 
    
21. 5 c# M/ q$ G& m
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. # `/ I5 e5 H/ d) v( ?' x
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. - g6 l! ~9 h$ J: U* G6 {2 C
26. t0=clock(),; X7 Y! X% m% l  N
    
28. s=0,& u; r, j/ G/ d9 Q. m
    
30. k=0, while{k<10000,( P6 F( c6 \* n2 y5 [
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值8 {( {9 w# q1 U
    
34.   k++6 u2 x$ K/ g/ l\\" x
    
36. },
    
37.   e; C9 f' m- U$ A
38. s;. u' K3 }2 X# ~\\" T% B
    
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">' Q5 G7 ^: m\\" B) n$ ^; A</font><br /> mvar:<font class=\\"jammer\\">3 v( F& {# t% g  t4 O</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">0 }6 q: P  v- ?8 c+ y$ q</span> t0=clock(),<br /> <span style=\\"display:none\\">8 k& u# X# X+ f7 B% d* C</span> k=0, while{k<10000,<font class=\\"jammer\\">( J6 R9 T7 q- K' s+ Q# d</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">* |6 K/ x3 i  _. }: D; @</font><br />   k++<font class=\\"jammer\\">$ B/ c1 E/ u* y' a3 }: l; E</font><br /> },<br /> <span style=\\"display:none\\">\\" W* w' X5 v3 a8 h- L: X+ i: m</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">. O3 l$ V9 f6 Q/ C0 ^9 b. ]</span> <br /> <span style=\\"display:none\\">5 c# M/ q$ G& m</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\"># `/ I5 e5 H/ d) v( ?' x</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">- g6 l! ~9 h$ J: U* G6 {2 C</span> t0=clock(),<font class=\\"jammer\\">; X7 Y! X% m% l  N</font><br /> s=0,<font class=\\"jammer\\">& u; r, j/ G/ d9 Q. m</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">( P6 F( c6 \* n2 y5 [</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">8 {( {9 w# q1 U</font><br />   k++<font class=\\"jammer\\">6 u2 x$ K/ g/ l\\" x</font><br /> },<br /> <span style=\\"display:none\\">  e; C9 f' m- U$ A</span> s;<font class=\\"jammer\\">. u' K3 }2 X# ~\\" T% B</font><br /> [clock()-t0]/1000;

结果：
4 o6 R1 ?2 U' J& I# x0.5620000000000001  //赋值用时间，秒
20000000.
/ \, F$ G3 ^) H" i. G0.579      //取值用时间，秒