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

forcal

1、简单循环

Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   ( I3 ~3 H$ D! w3 r* g2 V0 T% H( Y
2. tic6 Y; f7 p. U& h5 O( w- n
   
3. s=0;
   ) g, o' w1 O/ S4 g% _; P
4. for k=1:10000+ k. M* t- ?) Y3 v; T6 P' _3 Y
   
5.   for i=1:10002 @9 f8 }  y- [! |( |
   
6.     s=s+1;
   ' H+ n  o, r1 d# C8 v
7.   end
   1 I1 E! V) E\" y7 g- d( ]
8. end
   , \( j1 v& J6 a) a# ^: U; S
9. s
   . Y, x% \- F% b* _+ }, a! @
10. toc8 N6 \3 R, c. N  ^3 t2 z
    
11. 
    * `4 v& X3 w; @/ [$ `
12. s =
    & I& B+ P/ f( f& q7 o. V1 c9 A
13. / C+ M6 t' w+ K& \4 h$ \
    
14.     10000000( `/ h8 c8 Z+ w8 T9 j
    
15. : I, x2 l1 i/ ]: H! e2 C\" G' w+ ?
    
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
! N' K- i; U& f2 Q' Y/ f) M
Forcal没有JIT加速器，代码：

2. !using["sys"];- Q1 H. l  n( ~9 L% Y
   
4. mvar:' }/ L0 l) y( m0 l, _# ?  k4 X
   
6. t0=clock(),# W% I- O\\" l4 D$ e
   
8. s=0,; K; ^7 d\\" Q5 o( ~
   
10. k=0, while{k<10000,
    
11. & K1 g2 a/ F: O+ S7 B) X# j4 O
12.   i=0, while{i<1000,s=s+1, i++},
    
13. ( c9 V0 I! n8 O& k- U; |( H! B
14.   k++
    
15. $ e5 G) m+ m\\" s8 z& n
16. },
    
17. # B- }5 |0 i8 ?4 h8 o2 {: w; A( ^: B( w
18. s;! ]3 [/ L/ B, J  w
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">- Q1 H. l  n( ~9 L% Y</font><br /> mvar:<font class=\\"jammer\\">' }/ L0 l) y( m0 l, _# ?  k4 X</font><br /> t0=clock(),<font class=\\"jammer\\"># W% I- O\\" l4 D$ e</font><br /> s=0,<font class=\\"jammer\\">; K; ^7 d\\" Q5 o( ~</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">& K1 g2 a/ F: O+ S7 B) X# j4 O</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">( c9 V0 I! n8 O& k- U; |( H! B</span>   k++<br /> <span style=\\"display:none\\">$ e5 G) m+ m\\" s8 z& n</span> },<br /> <span style=\\"display:none\\"># B- }5 |0 i8 ?4 h8 o2 {: w; A( ^: B( w</span> s;<font class=\\"jammer\\">! ]3 [/ L/ B, J  w</font><br /> [clock()-t0]/1000;

结果：
" G: X% ]9 c) R10000000.
0.704  秒
. ]& t2 d4 B! o5 `. }" U0 S% x. \; y
% g9 L1 a% [) K7 M, h: c+ ^二者运行效率相差7倍多。
2 Z- u; a/ I' }, H5 |
2、循环中稍微加点东西
) y! C8 K& O. Z6 T7 x+ a
3 p8 B. r! I. [  p& W, J# S/ sMatlab代码：

1. clear all
   ' s. w# x7 G! ]$ Q. Y6 F
2. tic, U: P5 U1 e0 M, o1 R- r) a7 c
   
3. s=0;$ c) x+ ]0 c4 N: v1 h6 y1 O
   
4. for k=1:1000! J' I9 D0 Z: u1 _' g) U
   
5.   for i=1:1000! h* R+ T) P4 `4 u
   
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);
   : a9 w2 o* g\" D# K! E& `- q
7.   end5 S* t& Q' [4 N' f
   
8. end\" ]# }( B. _- l* p2 E
   
9. s
   4 F( x& \  M% Y* s
10. toc4 _& }8 a2 D3 g9 G
    
11. % g% ^) E8 ]6 G9 l/ g! W4 m
    
12. s =0 d\" Y3 c) I' k) C
    
13. 5 t; q5 n; V. F: L+ |, r! M* i0 q
    
14.   1.1009e+007
    3 m2 l+ s: X) s7 s# \5 b' ^8 C
15. 6 U: u0 |3 }\" x7 X- D: F- o0 e* Q
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];$ x# O6 F0 J+ y) z  p% Z( D1 ?4 V
   
4. mvar:* i, q' G6 o! P9 n% {: r; N
   
6. t0=clock(),- D1 y8 M- t! Y1 O$ c
   
8. s=0,9 H$ e& R$ {  b
   
10. k=1, while{k<=1000,
    
11. 6 E# C+ E* k4 Y- J\\" j
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++},# M! J0 D2 V; k/ v( A$ l' a' _
    
14.   k++
    
15. 4 V2 N2 J+ O5 _1 \. [: k
16. },
    
17. : a, o, S! M, H
18. s;
    
19. # Z/ c0 m- M+ Z5 o& _& X
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">$ x# O6 F0 J+ y) z  p% Z( D1 ?4 V</font><br /> mvar:<font class=\\"jammer\\">* i, q' G6 o! P9 n% {: r; N</font><br /> t0=clock(),<font class=\\"jammer\\">- D1 y8 M- t! Y1 O$ c</font><br /> s=0,<font class=\\"jammer\\">9 H$ e& R$ {  b</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">6 E# C+ E* k4 Y- J\\" j</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\\"># M! J0 D2 V; k/ v( A$ l' a' _</font><br />   k++<br /> <span style=\\"display:none\\">4 V2 N2 J+ O5 _1 \. [: k</span> },<br /> <span style=\\"display:none\\">: a, o, S! M, H</span> s;<br /> <span style=\\"display:none\\"># Z/ c0 m- M+ Z5 o& _& X</span> [clock()-t0]/1000;

结果：
11008953.66602346
$ ?4 G! E3 ?2 L+ j; s. A1.031

二者效率相差已比较小。

4 v! c/ _* j% s3、再增加一个函数调用

Matlab代码：

1. clear all\" K& w/ b9 i, p0 {; V, {
   
2. f=@(x)x+1;! l! T* p2 z3 p. g4 r. z5 ^
   
3. tic2 m' g! c* y\" `* C2 @, h
   
4. s=0;
   * s% b7 q* V7 {
5. for k=1:1000# E  Z8 o! P1 E. s
   
6.   for i=1:1000% T+ {& ^# H8 U! K8 n& E, n9 a8 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);
   7 h1 g* c/ [* v6 `
8.   end
   9 c( _- V4 J5 B/ {! Q$ s8 I' J\" ?
9. end
   , V- K7 P1 U- H, |$ x; a$ ^* [
10. s- X5 h; j* N/ o# x( M# j1 r
    
11. toc  k2 S8 l\" _$ R2 ^' r8 K* Q
    
12. 
    * I8 ?1 [- N. o& A/ f
13. s =- e8 S/ x3 X4 w9 K
    
14. ! g! ^$ v% Z: e4 h7 K8 t
    
15.   5.1251e+0085 H4 x: k. S& ~! I! f9 b\" I8 L0 e3 `
    
16. . f% q$ t0 L: D; q5 L
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];' ^' _6 S, [) J
   
4. f(x)=x+1;
   
5. ! d& Q% t3 z  S9 e
6. mvar:1 O; l2 V\\" x$ x0 Y  Y# m4 m# _
   
8. t0=clock(),
   
9. 2 u* a\\" T. K% Z6 f
10. s=0,. }, ~5 g3 {- _* t
    
12. k=1, while{k<=1000,9 A' ^: ?( X+ R0 G+ ~4 ~% M9 c* B
    
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. 7 U9 l5 u\\" t) ^/ s; N
16.   k++* ^# j6 J5 f. D* \3 {
    
18. },
    
19. $ }5 d3 w, L: c  e5 I4 A
20. s;
    
21.   T% m+ F1 p5 G/ v3 h0 b2 Z, o3 {
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">' ^' _6 S, [) J</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">! d& Q% t3 z  S9 e</span> mvar:<font class=\\"jammer\\">1 O; l2 V\\" x$ x0 Y  Y# m4 m# _</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">2 u* a\\" T. K% Z6 f</span> s=0,<font class=\\"jammer\\">. }, ~5 g3 {- _* t</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">9 A' ^: ?( X+ R0 G+ ~4 ~% M9 c* B</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++},<br /> <span style=\\"display:none\\">7 U9 l5 u\\" t) ^/ s; N</span>   k++<font class=\\"jammer\\">* ^# j6 J5 f. D* \3 {</font><br /> },<br /> <span style=\\"display:none\\">$ }5 d3 w, L: c  e5 I4 A</span> s;<br /> <span style=\\"display:none\\">  T% m+ F1 p5 G/ v3 h0 b2 Z, o3 {</span> [clock()-t0]/1000;

结果：
512508953.6658993
4 @& N0 }( u* P2 W3 y; }1.078

& \% Y. G9 P+ k. h9 z/ l; p可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

9 T4 ]; g  ^! O1 P2 u6 N8 m8 k; W2 xMatlab代码：

1. >> clear all
   % ?7 [9 X5 z2 A) y+ K! n( X
2. a = linspace(1,1,1000);2 j! V- V: l. i/ t' T# e# L
   
3. tic- s; z- L# ]/ T7 H( \: X5 B
   
4. for k=1:10000, J) g\" o4 M9 I( Z* ^' E
   
5.   for i=1:1000
   ! M5 l\" D: z/ P* p
6.     a(i)=2;* r  `, B\" y, A! _9 o
   
7.   end
   1 B, y0 r  Y' k* C6 {6 E
8. end
     a6 Q* K! }; v  r
9. toc
   / f5 ^: Q- o9 L2 m
10. % K$ x, @0 I2 `) X
    
11. tic4 W# M3 r2 l- ]# F
    
12. s=0;
    1 w% y0 E$ R. l7 \( e
13. for k=1:10000) o+ q4 p3 x& D1 Y3 W$ l5 E; R
    
14.   for i=1:1000
    : I# u; C- L' `  r6 |0 O( d
15.     s=s+a(i);$ O3 F\" u8 O1 H& F$ c
    
16.   end\" ]1 L% Y7 _$ `6 G
    
17. end% w- x\" W1 c' }  u: y
    
18. s) b$ z  @4 W) P
    
19. toc
    ' u/ d% N. h* x3 t
20. Elapsed time is 0.140658 seconds.* Y' ~+ @2 m\" Z  ^7 b8 y$ ?4 ]
    
21. 
    0 k% }4 u# ~- \8 [
22. s =
    $ J+ `7 b, i. c( p
23. 
    0 a5 u8 i: N! @
24.     20000000
    ( @& v7 {\" Z, P7 \% Y2 {# c6 x+ M
25. 
    7 p6 V( z% d1 N5 p; n7 N& Y
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];
   
3. ' @' A. F7 s! ~& ^3 v9 f- ~
4. mvar:* m) x2 B3 n- o* N$ a
   
6. main0(::a)=+ n' C0 |$ [1 N5 k
   
8. a = linspace(1,1,1000),1 u6 |$ h6 d8 t) J
   
10. t0=clock(),
    
11. % Z: y) [+ {: Q4 P; v* d9 F
12. k=0, while{k<10000,: j\\" e' k! ^# G7 C
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. + Z: o4 X- z# ~9 e4 b5 K) y
16.   k++
    
17. 6 O' B- R5 P. P! L  s! G
18. },* Z: m* Q' J\\" u7 `- y/ ]4 v; u5 C
    
20. [clock()-t0]/1000;
    
21. & a6 A! ^& z+ c3 f8 O
22. \\" I$ h6 {3 [5 U! b% B
    
24. main1(::a)=
    
25. ) H' P9 S1 C\\" ^1 ], L, t
26. t0=clock(), s=0,5 ]/ y0 x' a, V+ w/ i. s0 Y' Q
    
28. k=0, while{k<10000,
    
29. & r2 [: V\\" ]# @\\" Z
30.   i=0, while{i<1000, s=s+a(i), i++},- E& x0 l! c% ?& o& l
    
32.   k++3 D% w& m& Y6 P# Q9 ^* A( ~
    
34. },$ Q& g/ N2 [: I+ G' Y\\" [
    
36. s;
    
37. & Y\\" o0 f4 F9 n  }6 E
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">' @' A. F7 s! ~& ^3 v9 f- ~</span> mvar:<font class=\\"jammer\\">* m) x2 B3 n- o* N$ a</font><br /> main0(::a)=<font class=\\"jammer\\">+ n' C0 |$ [1 N5 k</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">1 u6 |$ h6 d8 t) J</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">% Z: y) [+ {: Q4 P; v* d9 F</span> k=0, while{k<10000,<font class=\\"jammer\\">: j\\" e' k! ^# G7 C</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">+ Z: o4 X- z# ~9 e4 b5 K) y</span>   k++<br /> <span style=\\"display:none\\">6 O' B- R5 P. P! L  s! G</span> },<font class=\\"jammer\\">* Z: m* Q' J\\" u7 `- y/ ]4 v; u5 C</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">& a6 A! ^& z+ c3 f8 O</span> <font class=\\"jammer\\">\\" I$ h6 {3 [5 U! b% B</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">) H' P9 S1 C\\" ^1 ], L, t</span> t0=clock(), s=0,<font class=\\"jammer\\">5 ]/ y0 x' a, V+ w/ i. s0 Y' Q</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">& r2 [: V\\" ]# @\\" Z</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">- E& x0 l! c% ?& o& l</font><br />   k++<font class=\\"jammer\\">3 D% w& m& Y6 P# Q9 ^* A( ~</font><br /> },<font class=\\"jammer\\">$ Q& g/ N2 [: I+ G' Y\\" [</font><br /> s;<br /> <span style=\\"display:none\\">& Y\\" o0 f4 F9 n  }6 E</span> [clock()-t0]/1000;

结果：
3 C/ D0 g6 K8 C) j2 E5 g- X4 [' a3 p15.484  //存数据所用时间，秒
5 q* h$ c* h$ N20000000.
5.531   //取数据所用时间，秒
1 I) S' V2 ^8 u) u9 l" d( Q
－－－－－－－－－－

- Q) s% r* c  \9 G用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];5 A% o& S6 c\\" G. k
   
4. mvar:
   
5. \\" x: n' S7 B6 Z
6. main0(::a)=: I6 L& T, t7 [
   
8. a = linspace(1,1,1000),
   
9. / H( @& v; \, E9 s& w& L: t
10. t0=clock(),  [( K2 I, ~* o6 O8 d
    
12. k=0, while{k<10000,
    
13. * ~/ W( n2 W, i! S; H+ K
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. 8 R\\" }6 j9 U5 [; j& v- U! d& q
16.   k++
    
17. % V, _  S' x9 ]+ Y
18. },
    
19. 6 d0 M# k8 {0 G( Q: d7 b2 |
20. [clock()-t0]/1000;5 @5 V5 [. O/ o) E' v& A/ \
    
22. 
    
23. 6 L2 |% b9 R7 P4 y6 `
24. main1(::a)=1 X! H0 B- W1 e! C
    
26. t0=clock(), s=0,
    
27. ! h\\" G/ i/ a  ?4 Y1 ^
28. k=0, while{k<10000,/ Q! w' A8 I4 C: e
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. . p# K0 n5 ^' M9 C0 E
32.   k++
    
33. * @# U( q4 I\\" P
34. },\\" S0 L  C0 W1 {! |/ q0 k0 S
    
36. s;
    
37. + q  t9 G6 }) }\\" x0 k: V7 L\\" }
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">5 A% o& S6 c\\" G. k</font><br /> mvar:<br /> <span style=\\"display:none\\">\\" x: n' S7 B6 Z</span> main0(::a)=<font class=\\"jammer\\">: I6 L& T, t7 [</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">/ H( @& v; \, E9 s& w& L: t</span> t0=clock(),<font class=\\"jammer\\">  [( K2 I, ~* o6 O8 d</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">* ~/ W( n2 W, i! S; H+ K</span>   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">8 R\\" }6 j9 U5 [; j& v- U! d& q</span>   k++<br /> <span style=\\"display:none\\">% V, _  S' x9 ]+ Y</span> },<br /> <span style=\\"display:none\\">6 d0 M# k8 {0 G( Q: d7 b2 |</span> [clock()-t0]/1000;<font class=\\"jammer\\">5 @5 V5 [. O/ o) E' v& A/ \</font><br /> <br /> <span style=\\"display:none\\">6 L2 |% b9 R7 P4 y6 `</span> main1(::a)=<font class=\\"jammer\\">1 X! H0 B- W1 e! C</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">! h\\" G/ i/ a  ?4 Y1 ^</span> k=0, while{k<10000,<font class=\\"jammer\\">/ Q! w' A8 I4 C: e</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">. p# K0 n5 ^' M9 C0 E</span>   k++<br /> <span style=\\"display:none\\">* @# U( q4 I\\" P</span> },<font class=\\"jammer\\">\\" S0 L  C0 W1 {! |/ q0 k0 S</font><br /> s;<br /> <span style=\\"display:none\\">+ q  t9 G6 }) }\\" x0 k: V7 L\\" }</span> [clock()-t0]/1000;

结果：
10.469
20000000.
0 r' `& |& k) j1 Y7 H; |0 a8 h3.781

( }- I4 H2 x3 E) J# Y$ ?1 i' {2 }! S－－－－－－－

3 V3 t2 E% h( C- I4 M, V在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
8 j) r+ j' z9 y" B8 z- M
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
) L- _6 E$ d9 X/ @
8 n2 t. z# F2 q4 \7 S" s" q8 |Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
3 w5 {  @, {, m$ ^) L/ i+ i
5、在循环中加点东西

3 K, }( c: y6 q6 ^3 Mmatlab代码：

1. >> clear all* ?& k& C1 z4 J) y0 P
   
2. f=@(x)x+1;# R: w& w+ C% a% t; `+ i
   
3. a = linspace(2,2,1000);
   , N6 b! m# I* [, o
4. tic
   + h- L. l, n+ x  R- p& g
5. s=0;\" \$ z4 a4 B) \1 x
   
6. for k=1:10000 u5 R- z7 g' _
   
7.   for i=1:1000
   0 i$ r* K8 X) s5 ^5 P% {  J  m
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));
   & ?9 y: X9 Q& o5 d. U  ^
9.   end! ?; Y% y  c- I  x/ ?
   
10. end
    , s7 m8 s! U9 t+ {3 Y, i
11. s8 d) j9 i2 A! ^
    
12. toc
    % i# g. R% C2 w
13. % p4 ~\" ]4 _, I7 D* I, u
    
14. s =
    : G% t3 [$ X  T7 I\" U# V& B
15. 
    0 _4 R# @, c- z. n- J! y. Z: \( a
16.   2.4002e+007
    / u4 k2 v6 z  B( h0 @
17. 
    , ^8 s$ D: N$ o# h8 p' J$ ?
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];7 P9 \% J  W/ s$ I
   
4. mvar:# x9 C) q7 i( z! a1 P+ ^
   
6. f(x)=x+1;' H& f! s& R6 K
   
8. main0(::a)=. U/ L2 ]3 `! S\\" p: L9 E* F7 f% H
   
10. a = linspace(2,2,1000),
    
11. * {% W\\" F6 x3 c# x. \2 X
12. t0=clock(), s=0,
    
13. 9 Q1 ]( S6 `8 t0 ~  l9 t* M, }) x
14. k=0, while{k<1000,1 T\\" e6 V\\" X8 l4 F8 }5 m
    
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++},6 m6 Y\\" v' |3 \$ a\\" }5 q' G7 \
    
18.   k++
    
19. ' x+ k9 W6 ?1 H# T& d) n9 h8 t
20. },
    
21. 7 h$ M$ f' j- ~0 s
22. s;) @7 c+ k\\" d, S2 c6 U
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">7 P9 \% J  W/ s$ I</font><br /> mvar:<font class=\\"jammer\\"># x9 C) q7 i( z! a1 P+ ^</font><br /> f(x)=x+1;<font class=\\"jammer\\">' H& f! s& R6 K</font><br /> main0(::a)=<font class=\\"jammer\\">. U/ L2 ]3 `! S\\" p: L9 E* F7 f% H</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">* {% W\\" F6 x3 c# x. \2 X</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">9 Q1 ]( S6 `8 t0 ~  l9 t* M, }) x</span> k=0, while{k<1000,<font class=\\"jammer\\">1 T\\" e6 V\\" X8 l4 F8 }5 m</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\\">6 m6 Y\\" v' |3 \$ a\\" }5 q' G7 \</font><br />   k++<br /> <span style=\\"display:none\\">' x+ k9 W6 ?1 H# T& d) n9 h8 t</span> },<br /> <span style=\\"display:none\\">7 h$ M$ f' j- ~0 s</span> s;<font class=\\"jammer\\">) @7 c+ k\\" d, S2 c6 U</font><br /> [clock()-t0]/1000;

结果：
7 l2 X7 t6 [6 C  s8 X24002271.69189384
5.172

) V8 _9 }# Y6 p+ ^# H0 P, q" O二者效率现在只有一倍差距。
  R/ ]5 t' j* Q9 @/ L! W) M% q
－－－－－－－－－－
) q* j1 j& F7 J7 S0 m5 O4 ]
不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

1. clear all
   4 h7 e- m# q8 Y0 c  K) E
2. f=@(x)x+1;; h2 B! m3 r% a
   
3. a = linspace(2,2,1000);
   1 H4 t' C* B1 B$ C8 e. @
4. tic
   / q$ @) S' ^4 t8 [
5. s=0;
   % z$ B. F' P- C% j
6. for k=1:10007 f+ g& N1 W: T+ m, T; g
   
7.   for i=1:1000% h: Q$ n3 G. K2 I
   
8.     t=a(i);! Y, M& M\" y) [
   
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);
   , y1 s2 U/ `% T; t! K8 c2 R
10.   end& `6 A4 [9 w' U0 f7 N
    
11. end
    . {8 X. B\" Y7 L  W& D
12. s# c& D+ n% i9 W+ h2 M4 X, f9 c
    
13. toc# |  o% U$ k9 H; L$ v
    
14. 
    7 v9 @5 C0 I+ q. L  E8 ?! A( ?
15. s =0 {+ t; A& [6 }: ]
    
16. , I- \5 z, c/ V! J* _
    
17.   2.4002e+007# k. m6 ^( U2 |\" K& N% {
    
18. 
    $ v\" v, G, l. r
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. ) J' `- F\\" @4 r( V- I7 m
4. mvar:
   
5. 2 m6 @8 P3 d4 l
6. f(x)=x+1;
   
7. 1 `( V6 A/ L: P
8. main0(::a)=
   
9.   b9 u; v  w9 J2 I; E' }' l
10. a = linspace(2,2,1000),, d2 N+ P% s9 ]: c4 {2 @- B
    
12. t0=clock(), s=0,- a8 Z, C- O+ W0 `
    
14. k=0, while{k<1000,
    
15. / s\\" G) b/ j) Z' \, @' C! h
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& h; K2 J; \
18.   k++
    
19. & k2 c7 I( z& Y( ~\\" t, n. x
20. },1 Q5 a. s; S) X- P
    
22. s;. O6 Z1 L5 y! ~\\" K
    
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">) J' `- F\\" @4 r( V- I7 m</span> mvar:<br /> <span style=\\"display:none\\">2 m6 @8 P3 d4 l</span> f(x)=x+1;<br /> <span style=\\"display:none\\">1 `( V6 A/ L: P</span> main0(::a)=<br /> <span style=\\"display:none\\">  b9 u; v  w9 J2 I; E' }' l</span> a = linspace(2,2,1000),<font class=\\"jammer\\">, d2 N+ P% s9 ]: c4 {2 @- B</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">- a8 Z, C- O+ W0 `</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">/ s\\" G) b/ j) Z' \, @' C! h</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\\">+ P& h; K2 J; \</span>   k++<br /> <span style=\\"display:none\\">& k2 c7 I( z& Y( ~\\" t, n. x</span> },<font class=\\"jammer\\">1 Q5 a. s; S) X- P</font><br /> s;<font class=\\"jammer\\">. O6 Z1 L5 y! ~\\" K</font><br /> [clock()-t0]/1000;

结果：
% G- }1 K# o2 ]9 X0 V) L: f24002271.69189384
1.625
: p3 X/ D) i) ^+ W/ x8 z  S
6 ], u9 g" D/ S$ X) Gmatlab效率已不及Forcal。
* j& Z) _$ y  c0 V/ l3 x% b
提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
* J. x5 L( r/ I8 I# u
Forcal代码：

2. !using["sys"];4 \: r; G5 m- y( w4 Y3 o1 L& m
   
4. mvar:, d( U3 O9 ^$ }
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组/ S5 r: S, g2 e) W9 _: J1 M& }% O9 J
   
8. t0=clock(),  V  J0 J, b, J! v2 l
   
10. k=0, while{k<10000,
    
11. $ j0 g+ S$ ^* G+ s6 w4 ]0 o+ w9 k
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. 5 S\\" I+ K4 d) P7 |$ ?$ E) |. s9 E/ h
14.   k++
    
15. 7 p9 ~  U) r  Z( Q6 M
16. },
    
17. 0 ^( n% g$ F$ ?! \5 H\\" s: k* s
18. [clock()-t0]/1000;
    
19. 5 u% J8 u) ?* L# }4 }, _6 E0 H. ^) k
20. 
    
21. 5 }; y6 J$ ?6 F; @
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. 6 K\\" A- d. l: H4 g3 E: [
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. . [2 x% Q. O$ ]2 S# o! [
26. t0=clock(),
    
27. 0 k/ W& i\\" ?+ P
28. s=0,5 A: ~; |3 i2 h1 ^- d\\" |% v+ U0 o+ O
    
30. k=0, while{k<10000,0 g6 o) r\\" W# U5 n6 Y+ d6 c
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值! e8 i# b/ e/ T2 R
    
34.   k++
    
35. 1 h  s: W( \8 s/ f2 o! W. X
36. },
    
37. # D$ V* j1 j/ d3 r6 M1 h
38. s;
    
39. - S! x' a2 \' l- z\\" Q; P\\" h2 y
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">4 \: r; G5 m- y( w4 Y3 o1 L& m</font><br /> mvar:<font class=\\"jammer\\">, d( U3 O9 ^$ }</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">/ S5 r: S, g2 e) W9 _: J1 M& }% O9 J</font><br /> t0=clock(),<font class=\\"jammer\\">  V  J0 J, b, J! v2 l</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">$ j0 g+ S$ ^* G+ s6 w4 ]0 o+ w9 k</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">5 S\\" I+ K4 d) P7 |$ ?$ E) |. s9 E/ h</span>   k++<br /> <span style=\\"display:none\\">7 p9 ~  U) r  Z( Q6 M</span> },<br /> <span style=\\"display:none\\">0 ^( n% g$ F$ ?! \5 H\\" s: k* s</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">5 u% J8 u) ?* L# }4 }, _6 E0 H. ^) k</span> <br /> <span style=\\"display:none\\">5 }; y6 J$ ?6 F; @</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">6 K\\" A- d. l: H4 g3 E: [</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">. [2 x% Q. O$ ]2 S# o! [</span> t0=clock(),<br /> <span style=\\"display:none\\">0 k/ W& i\\" ?+ P</span> s=0,<font class=\\"jammer\\">5 A: ~; |3 i2 h1 ^- d\\" |% v+ U0 o+ O</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">0 g6 o) r\\" W# U5 n6 Y+ d6 c</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">! e8 i# b/ e/ T2 R</font><br />   k++<br /> <span style=\\"display:none\\">1 h  s: W( \8 s/ f2 o! W. X</span> },<br /> <span style=\\"display:none\\"># D$ V* j1 j/ d3 r6 M1 h</span> s;<br /> <span style=\\"display:none\\">- S! x' a2 \' l- z\\" Q; P\\" h2 y</span> [clock()-t0]/1000;

结果：
% Q' c, f- E6 a: X1 u0.5620000000000001  //赋值用时间，秒
/ Y% `* a9 P8 N8 [% \20000000.
& H& Y+ d3 n* n( K( K0.579      //取值用时间，秒