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

forcal

1、简单循环

) f) l/ M( W* P. L& |- o/ u6 uMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all1 n) J* c% t9 F4 F
   
2. tic
   $ k, z% A+ X5 M# j
3. s=0;3 X+ h: Z: Q3 y  |2 q8 J: t\" I
   
4. for k=1:100002 V1 ^& v! K. E/ y/ x! C! ?
   
5.   for i=1:1000
   ( o/ s5 b! ^9 J, `7 Z0 I\" \
6.     s=s+1;; b+ \$ K' z7 s+ i2 W' }6 M6 ?6 @% G
   
7.   end
   0 g/ S) c+ T5 j, m9 U* e7 c
8. end# X& S0 r% @) b9 _7 q/ g
   
9. s
   : ?& c1 ~\" E\" G
10. toc
    - k/ h: ~- T* f9 s, w4 h3 O
11. 
    * ]' I\" c$ S) g
12. s =
    , ]. W6 l; C; H) O
13. 
    # ^. ?/ \\" {: K' ?2 x
14.     10000000. @+ }2 D% E: E, K# H
    
15. 7 Y+ A- S2 n. Z) A. l- b0 P0 `
    
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
) _7 H$ c; B8 D7 e
Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. 1 G( R1 E' l% X/ A: H, T' @* j* ?
4. mvar:1 E0 c8 \- [7 k) s
   
6. t0=clock(),5 S/ O- V$ t% F3 }& ~
   
8. s=0,/ ?. A. \8 f! r& Q3 ]% v) o
   
10. k=0, while{k<10000,* _$ R' K2 b3 _7 S, W: r) O% Z& j
    
12.   i=0, while{i<1000,s=s+1, i++},6 P. V3 C* }; Z( P
    
14.   k++
    
15. 8 e# m5 ?6 t; g
16. },! Z- X$ L5 [7 j$ p7 T* {
    
18. s;
    
19.   ~# K6 a/ J# `4 `4 L\\" j7 s6 x  W
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">1 G( R1 E' l% X/ A: H, T' @* j* ?</span> mvar:<font class=\\"jammer\\">1 E0 c8 \- [7 k) s</font><br /> t0=clock(),<font class=\\"jammer\\">5 S/ O- V$ t% F3 }& ~</font><br /> s=0,<font class=\\"jammer\\">/ ?. A. \8 f! r& Q3 ]% v) o</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">* _$ R' K2 b3 _7 S, W: r) O% Z& j</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">6 P. V3 C* }; Z( P</font><br />   k++<br /> <span style=\\"display:none\\">8 e# m5 ?6 t; g</span> },<font class=\\"jammer\\">! Z- X$ L5 [7 j$ p7 T* {</font><br /> s;<br /> <span style=\\"display:none\\">  ~# K6 a/ J# `4 `4 L\\" j7 s6 x  W</span> [clock()-t0]/1000;

结果：
10000000.
0.704  秒
" O8 H+ E  B$ y- b
% ?) R. n2 e) j2 Y二者运行效率相差7倍多。

2、循环中稍微加点东西

* g. X& c* T/ z, |* cMatlab代码：

1. clear all, C; E% `. C8 S/ L, y& s
   
2. tic
   ; r. B- x6 n1 v% L+ M6 H( a8 b
3. s=0;' C! q- L. Y% D% |) _, _
   
4. for k=1:1000
   ; ^: @* @+ B# _  U- e9 ^
5.   for i=1:10002 U* Z$ v& t! n# F. C: F. w. g
   
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);
   * w4 a\" H0 w( Q9 J7 Q
7.   end
   . t1 Z9 q! t3 h( i  @
8. end
   ' [; K, f% l$ W9 H  c- `7 L( |
9. s6 |- H% b2 z5 r4 y\" i- m# r! g' ^; C( D
   
10. toc
    6 F$ x  T\" I1 I
11. 
    ' q2 H, I$ _% W) A6 t) ~
12. s =
      C) F' W5 ]& a( {( x/ Q  E2 F
13. 
    ; G4 d4 X! j) ~2 F
14.   1.1009e+007
    1 K! U9 h\" b% j
15. 
    9 u* v6 o5 W, f2 i0 ]
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];$ c. J- O1 b, e4 J
   
4. mvar:3 Z5 ]  i2 W+ v, s
   
6. t0=clock(),
   
7. ( h: n' q  h, u5 o9 Z\\" l* d
8. s=0,
   
9. % N9 ?$ k  _0 u2 j0 A  j
10. k=1, while{k<=1000,1 r( I, t0 _1 ~0 d+ O4 A\\" |2 w; @
    
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++},/ [$ l8 V! A7 A8 K% Q* E' u
    
14.   k++$ I  P5 l) D, v; G) j) c3 j
    
16. },
    
17. * q0 d' y4 n% m7 I) q* _+ ^
18. s;
    
19. 8 ~- z2 R& i) ?5 ]
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">$ c. J- O1 b, e4 J</font><br /> mvar:<font class=\\"jammer\\">3 Z5 ]  i2 W+ v, s</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">( h: n' q  h, u5 o9 Z\\" l* d</span> s=0,<br /> <span style=\\"display:none\\">% N9 ?$ k  _0 u2 j0 A  j</span> k=1, while{k<=1000,<font class=\\"jammer\\">1 r( I, t0 _1 ~0 d+ O4 A\\" |2 w; @</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++},<font class=\\"jammer\\">/ [$ l8 V! A7 A8 K% Q* E' u</font><br />   k++<font class=\\"jammer\\">$ I  P5 l) D, v; G) j) c3 j</font><br /> },<br /> <span style=\\"display:none\\">* q0 d' y4 n% m7 I) q* _+ ^</span> s;<br /> <span style=\\"display:none\\">8 ~- z2 R& i) ?5 ]</span> [clock()-t0]/1000;

结果：
9 h6 t) J! {" J' K; s8 g11008953.66602346
1.031

二者效率相差已比较小。
; I; f7 W+ S4 p8 ]+ V9 u
/ B* N& w: A, r; d4 D3、再增加一个函数调用
" o  O& O3 q; M( i+ d7 _
7 g6 k' I/ G% B5 ~! N  mMatlab代码：

1. clear all
   * Y. m9 c* X) S5 _5 A
2. f=@(x)x+1;+ \9 S+ P& G\" f! y* S\" ~3 J. T
   
3. tic- j* T: v- i) C
   
4. s=0;* f6 @9 h0 l7 v! D! e8 K
   
5. for k=1:1000
   ( U\" J1 z/ M, R: K* v# m6 y. d
6.   for i=1:1000
   7 Z& p5 \9 t9 H\" c$ F9 U- R; v
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);
   : m2 ^. F\" n. j4 |3 N4 j
8.   end
   3 @4 {& e! S# |4 U7 E% i1 h* {+ X8 i
9. end! t  \7 |& A  G. X' I0 b! y. R$ H. A
   
10. s
    % h% t6 l. S  R! G
11. toc
    4 U; w) f& B; _. V& M1 b) u3 u
12. 
    ! y* z+ d7 j1 o/ u# u; T6 b4 `+ H
13. s =
    . w0 U. [7 }4 _5 G: k; O
14. 
    + p+ _, ?) Z$ F\" g3 o+ q
15.   5.1251e+008: L# g3 I+ \/ B( r/ q2 D
    
16. 
    5 H1 w6 |; t\" l& M
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];  J9 i  _1 s$ z\\" f. j; C4 x$ a0 i
   
4. f(x)=x+1;
   
5. : b# f\\" _, [0 G' N& {' ]5 T. Q
6. mvar:* P' \\\" ~5 b9 W% U
   
8. t0=clock(),
   
9. % H& U! \5 ]2 r. s+ ]# U% a\\" n
10. s=0,9 g' P3 ~8 A1 b\\" r5 j0 X( q
    
12. k=1, while{k<=1000,
    
13. + }2 a- i8 Y; x0 e6 x
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++},; V! x4 B\\" Y1 a! n8 n
    
16.   k++  C) e8 {4 ?8 @% V! K! s; q
    
18. },
    
19. \\" \, q' W( l4 E& D4 X
20. s;7 j4 u$ i4 {( L3 t8 z2 ~6 w7 Z* Y! g
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">  J9 i  _1 s$ z\\" f. j; C4 x$ a0 i</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">: b# f\\" _, [0 G' N& {' ]5 T. Q</span> mvar:<font class=\\"jammer\\">* P' \\\" ~5 b9 W% U</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">% H& U! \5 ]2 r. s+ ]# U% a\\" n</span> s=0,<font class=\\"jammer\\">9 g' P3 ~8 A1 b\\" r5 j0 X( q</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">+ }2 a- i8 Y; x0 e6 x</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\\">; V! x4 B\\" Y1 a! n8 n</font><br />   k++<font class=\\"jammer\\">  C) e8 {4 ?8 @% V! K! s; q</font><br /> },<br /> <span style=\\"display:none\\">\\" \, q' W( l4 E& D4 X</span> s;<font class=\\"jammer\\">7 j4 u$ i4 {( L3 t8 z2 ~6 w7 Z* Y! g</font><br /> [clock()-t0]/1000;

结果：
512508953.6658993
1.078

可以看出，Matlab效率已不及Forcal。
9 x' \! J6 I) A7 }$ X' s: p3 a故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
) Y% y6 G5 N) b. z. a/ P0 r' C
Matlab代码：

1. >> clear all+ T: v1 ?9 P, j0 K\" F
   
2. a = linspace(1,1,1000);
     Q1 I1 p' b5 z
3. tic- ~  ?/ b/ x! r. K$ z7 Y
   
4. for k=1:10000# ^9 y. S6 x5 i4 Z6 ^5 A
   
5.   for i=1:10009 u  p1 m! [. v
   
6.     a(i)=2;6 N( X\" J- V; \: Z, Q
   
7.   end
   & `5 I8 e& Z. S6 g( G8 m3 ~
8. end+ S; k# m+ u- H5 v- a, Q
   
9. toc
   ( y* E! X4 C3 c, g5 }
10. % O* Q2 n7 X) X2 e/ C5 w
    
11. tic, e. G4 E: x0 K. L( _& C
    
12. s=0;4 G# T8 P' s6 p/ I
    
13. for k=1:10000, G* m- K; h. u! C  U\" F  q
    
14.   for i=1:1000* @0 ?# l9 A' n) z/ B8 R  P* U
    
15.     s=s+a(i);
    - n3 Z: P8 w- L( d
16.   end3 _- m3 s% G5 J& n+ X/ d) }! E
    
17. end
    0 J3 e1 o7 U* Q' Y& \
18. s/ ~: p- |, [# N7 h
    
19. toc& b\" K, z! U  k# X4 u\" W
    
20. Elapsed time is 0.140658 seconds.
    ) W$ A2 D) U% `
21. \" d, ]6 a& P6 j5 C, C( f7 F3 k
    
22. s =
    : x9 J( Z# s9 _! ?7 B
23. 
    4 q2 P  |0 z( Y7 \& _
24.     20000000; V. N* {\" z3 R/ w2 k, ^
    
25. 7 v! \$ L# h0 K8 e. [6 f7 w
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];
   
3. 6 ?3 w8 P/ L! {( h& F
4. mvar:$ b3 V* Z$ P; Q) C8 [2 r8 r
   
6. main0(::a)=
   
7. ' U/ ?* l% M' ]/ T( b
8. a = linspace(1,1,1000),
   
9. , K8 [' z\\" G$ C3 o6 I3 `, `# ^
10. t0=clock(),
    
11. 3 ?& I1 ]: e. s- S1 N0 D
12. k=0, while{k<10000,
    
13. & V1 K; Y% S, C5 Z. l3 Z
14.   i=0, while{i<1000, a(i)=2, i++},+ T, a$ S, j1 `6 T( |/ e) D
    
16.   k++3 x* l- Y( a9 B- K4 s3 Y
    
18. },9 Z5 Q# V# I9 a! p# t
    
20. [clock()-t0]/1000;8 T0 V. {3 L! \6 [4 b
    
22. 
    
23. 3 B2 L  P* j- E/ Z
24. main1(::a)=
    
25. ' C+ ?# L. O6 v. w
26. t0=clock(), s=0,
    
27. 0 `6 z  O8 }( o5 K# v0 x, ^) P
28. k=0, while{k<10000,
    
29. # }\\" o) K' P1 S  K% B' Q# n
30.   i=0, while{i<1000, s=s+a(i), i++},- x5 B% N5 o% n# d8 F5 F+ J\\" N1 b
    
32.   k++. Q5 s* Q) `$ F* s7 K: j
    
34. },2 R# _5 p4 b( J
    
36. s;4 ~$ c; E: R2 H\\" {8 E
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">6 ?3 w8 P/ L! {( h& F</span> mvar:<font class=\\"jammer\\">$ b3 V* Z$ P; Q) C8 [2 r8 r</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">' U/ ?* l% M' ]/ T( b</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">, K8 [' z\\" G$ C3 o6 I3 `, `# ^</span> t0=clock(),<br /> <span style=\\"display:none\\">3 ?& I1 ]: e. s- S1 N0 D</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">& V1 K; Y% S, C5 Z. l3 Z</span>   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">+ T, a$ S, j1 `6 T( |/ e) D</font><br />   k++<font class=\\"jammer\\">3 x* l- Y( a9 B- K4 s3 Y</font><br /> },<font class=\\"jammer\\">9 Z5 Q# V# I9 a! p# t</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">8 T0 V. {3 L! \6 [4 b</font><br /> <br /> <span style=\\"display:none\\">3 B2 L  P* j- E/ Z</span> main1(::a)=<br /> <span style=\\"display:none\\">' C+ ?# L. O6 v. w</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">0 `6 z  O8 }( o5 K# v0 x, ^) P</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\"># }\\" o) K' P1 S  K% B' Q# n</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">- x5 B% N5 o% n# d8 F5 F+ J\\" N1 b</font><br />   k++<font class=\\"jammer\\">. Q5 s* Q) `$ F* s7 K: j</font><br /> },<font class=\\"jammer\\">2 R# _5 p4 b( J</font><br /> s;<font class=\\"jammer\\">4 ~$ c; E: R2 H\\" {8 E</font><br /> [clock()-t0]/1000;

结果：
" C* M7 x# E8 r7 F15.484  //存数据所用时间，秒
( p" t4 h4 z, F! W20000000.
- i4 j. z: u0 }/ E5 R! m* ^5.531   //取数据所用时间，秒
# a' [' k4 i- s. h  [& {6 D
( Y3 j1 y2 V. P- n－－－－－－－－－－
/ t7 \  w! e1 o# g9 c
( @. t6 s; Z7 R  H3 ?用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. $ ~7 U# o6 G/ ?
4. mvar:
   
5. ; n8 ]  c8 g4 a, @( K9 |2 p
6. main0(::a)=* n6 ?6 u, {$ n! V$ k, H( {
   
8. a = linspace(1,1,1000),
   
9. : u( k7 m4 C! C( y( Q2 n3 |1 ?
10. t0=clock(),
    
11. % B+ p5 {- j) d5 X, D# g$ q2 [
12. k=0, while{k<10000,% P+ e2 t; g6 q8 L! z
    
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. / q' m( Q1 A+ p9 D2 O% J\\" w, m& a2 ?8 v
16.   k++: u0 @% [  w! M- K, K
    
18. },
    
19. / O: r8 ~4 ?  r) Z$ C' {
20. [clock()-t0]/1000;+ g% w$ X\\" I  B
    
22. + f8 f0 O# T; N4 v4 F& y% w3 K
    
24. main1(::a)=
    
25. 6 X* E\\" Y$ F& A; C7 p
26. t0=clock(), s=0,\\" r; o% P$ b\\" m3 a
    
28. k=0, while{k<10000,! y% Q/ l! q4 g4 h' ?
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},; m) J* r; b  A/ e' s
    
32.   k++
    
33. 3 v9 K0 C5 q. s6 }. a! _
34. },' q) M: z2 c2 W$ n: f6 C- Y
    
36. s;& V7 k8 D  r4 I+ S\\" x9 J
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">$ ~7 U# o6 G/ ?</span> mvar:<br /> <span style=\\"display:none\\">; n8 ]  c8 g4 a, @( K9 |2 p</span> main0(::a)=<font class=\\"jammer\\">* n6 ?6 u, {$ n! V$ k, H( {</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">: u( k7 m4 C! C( y( Q2 n3 |1 ?</span> t0=clock(),<br /> <span style=\\"display:none\\">% B+ p5 {- j) d5 X, D# g$ q2 [</span> k=0, while{k<10000,<font class=\\"jammer\\">% P+ e2 t; g6 q8 L! z</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">/ q' m( Q1 A+ p9 D2 O% J\\" w, m& a2 ?8 v</span>   k++<font class=\\"jammer\\">: u0 @% [  w! M- K, K</font><br /> },<br /> <span style=\\"display:none\\">/ O: r8 ~4 ?  r) Z$ C' {</span> [clock()-t0]/1000;<font class=\\"jammer\\">+ g% w$ X\\" I  B</font><br /> <font class=\\"jammer\\">+ f8 f0 O# T; N4 v4 F& y% w3 K</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">6 X* E\\" Y$ F& A; C7 p</span> t0=clock(), s=0,<font class=\\"jammer\\">\\" r; o% P$ b\\" m3 a</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">! y% Q/ l! q4 g4 h' ?</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">; m) J* r; b  A/ e' s</font><br />   k++<br /> <span style=\\"display:none\\">3 v9 K0 C5 q. s6 }. a! _</span> },<font class=\\"jammer\\">' q) M: z2 c2 W$ n: f6 C- Y</font><br /> s;<font class=\\"jammer\\">& V7 k8 D  r4 I+ S\\" x9 J</font><br /> [clock()-t0]/1000;

结果：
1 g5 {/ ?: B# i% {/ ]7 s5 X& [10.469
20000000.
3.781

－－－－－－－
: |& U) n: K3 A1 X% L5 j& E
. I  m5 i  Q* C& P% P- v在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
/ c3 V- V/ Y+ h) z. h
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
. _9 x3 A. C- n8 C) ]+ S; }
( B6 m: P% E! Z, I+ oForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

5 E8 K, m# i6 G0 o6 `7 k3 \5、在循环中加点东西
0 z5 e) {6 ^0 Z7 e1 u: q& T" w4 w
matlab代码：

1. >> clear all
   ; C$ m. x  ~$ K: @# G
2. f=@(x)x+1;: k) {1 s' U3 d2 x1 X0 o4 \
   
3. a = linspace(2,2,1000);9 W( G! Y+ L8 ?' E! E
   
4. tic1 p9 ^4 ^7 y# F1 ?
   
5. s=0;: L0 D/ ^8 f0 G$ X2 f  m9 |3 R
   
6. for k=1:10006 B3 E7 Y0 S7 u1 M' _; d
   
7.   for i=1:1000
   2 _$ L3 e; g5 |1 Z: q3 c% z  D  ^
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));
   5 z( T\" d\" ^1 |# T/ o& G1 k
9.   end  k0 _* M1 H1 n
   
10. end! L+ D, m% w  N1 G# b% a
    
11. s# |1 B( M\" T' r
    
12. toc
    6 q. l$ P: I2 F\" K; g1 i( N
13.   h* n$ d* v  u2 W( b8 v+ a
    
14. s =
    . b# A1 }) C7 h
15. ) ^2 V' \0 @* F0 ]
    
16.   2.4002e+007
    7 i# E6 c# ?: l2 o1 Z( P
17. 
    6 J: l) r+ ^# \. f
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. 8 `. S6 F  @5 o9 q
4. mvar:9 r% O- N. B, }9 S
   
6. f(x)=x+1;5 P6 H; e8 i- E$ t, i
   
8. main0(::a)=
   
9. 2 d1 c# \9 Z) ^: c9 j
10. a = linspace(2,2,1000),8 ]5 X8 k& [. s  X
    
12. t0=clock(), s=0,
    
13. * i8 ]* Q! B  m1 W7 j. |1 C9 Z
14. k=0, while{k<1000,
    
15. 5 w* u6 S1 h  X  f$ X
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++},
    
17. $ l6 ?- W% }! e  I
18.   k++, i1 g' b7 C& z5 {& T+ j
    
20. },3 Z' X\\" E( Y) `- u
    
22. s;
    
23. 4 Y( [% S7 s* y7 {2 _- A# A
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">8 `. S6 F  @5 o9 q</span> mvar:<font class=\\"jammer\\">9 r% O- N. B, }9 S</font><br /> f(x)=x+1;<font class=\\"jammer\\">5 P6 H; e8 i- E$ t, i</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">2 d1 c# \9 Z) ^: c9 j</span> a = linspace(2,2,1000),<font class=\\"jammer\\">8 ]5 X8 k& [. s  X</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">* i8 ]* Q! B  m1 W7 j. |1 C9 Z</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">5 w* u6 S1 h  X  f$ X</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++},<br /> <span style=\\"display:none\\">$ l6 ?- W% }! e  I</span>   k++<font class=\\"jammer\\">, i1 g' b7 C& z5 {& T+ j</font><br /> },<font class=\\"jammer\\">3 Z' X\\" E( Y) `- u</font><br /> s;<br /> <span style=\\"display:none\\">4 Y( [% S7 s* y7 {2 _- A# A</span> [clock()-t0]/1000;

结果：
24002271.69189384
/ k( c0 a3 A1 X# X3 ]/ S" I5.172

) g6 l$ C! _% H- {* e1 T% h" c: J: f二者效率现在只有一倍差距。

, Y( ^$ Z4 o7 f# u8 |－－－－－－－－－－
" x1 [2 b4 @' s$ z7 I4 Z
不过，显然matlab和Forcal的程序应该这样写：
; {" @! k8 U4 l& r* A9 O" g6 Z0 N
+ [: B: }3 F" B6 h) qmatlab程序：

1. clear all5 \, a4 N  `/ v3 ~3 i1 l
   
2. f=@(x)x+1;
   5 s* d' [4 B) c: B; ]6 f: J
3. a = linspace(2,2,1000);9 p( |# n: [; K. J- \, ~\" T
   
4. tic& |7 E8 T. L4 G7 a5 j
   
5. s=0;
   ) E3 U3 z. l  I- f
6. for k=1:1000: w7 \7 j1 Z6 q# u
   
7.   for i=1:1000
   7 m/ A\" ?* C2 I: V  d1 G* N& k7 [4 P
8.     t=a(i);2 l. {% ]- z\" Z& k; O1 t/ G
   
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);
   , ~) l2 ]: J6 O) b8 x8 T3 k
10.   end
      t8 n6 ~' K% h' f' _$ Y
11. end3 O8 x! ~  p) z+ j9 `- W
    
12. s
    ; c5 s4 F6 Y. c; U
13. toc
    % ~5 y& k) h4 C  c7 |
14. 
    1 l, o3 J7 Z. m* E' B
15. s =
    \" h: n( |7 Z* A* x- g8 c
16. 9 @, O\" i- C4 `& W  H2 [# s9 }$ P
    
17.   2.4002e+0075 ^$ J$ i1 i( k* c8 X: M
    
18. 
    + a- |5 B7 ?% k( a
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];+ S8 c* E8 h0 T
   
4. mvar:
   
5. \\" V% D5 S5 ~+ Y* X* y
6. f(x)=x+1;$ o' |4 e3 @) `
   
8. main0(::a)=: C5 k$ Y9 Y; _0 E* }
   
10. a = linspace(2,2,1000),' M( ]+ C. Y: g, }
    
12. t0=clock(), s=0,
    
13. ! x5 i0 d. G! J- l1 c
14. k=0, while{k<1000,2 j* g& m8 U5 ^
    
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. , v* o; t8 L  r! Z3 s/ c5 e5 y
18.   k++
    
19. 9 Z* C4 ^  g) x7 p/ k5 ?5 O$ K8 \' V  T
20. },
    
21. ; U6 i  x0 B, u. @+ h, N
22. s;
    
23. : m2 y% q: f/ _: p, F( f9 R, `8 o4 z
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ S8 c* E8 h0 T</font><br /> mvar:<br /> <span style=\\"display:none\\">\\" V% D5 S5 ~+ Y* X* y</span> f(x)=x+1;<font class=\\"jammer\\">$ o' |4 e3 @) `</font><br /> main0(::a)=<font class=\\"jammer\\">: C5 k$ Y9 Y; _0 E* }</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">' M( ]+ C. Y: g, }</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">! x5 i0 d. G! J- l1 c</span> k=0, while{k<1000,<font class=\\"jammer\\">2 j* g& m8 U5 ^</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\\">, v* o; t8 L  r! Z3 s/ c5 e5 y</span>   k++<br /> <span style=\\"display:none\\">9 Z* C4 ^  g) x7 p/ k5 ?5 O$ K8 \' V  T</span> },<br /> <span style=\\"display:none\\">; U6 i  x0 B, u. @+ h, N</span> s;<br /> <span style=\\"display:none\\">: m2 y% q: f/ _: p, F( f9 R, `8 o4 z</span> [clock()-t0]/1000;

结果：
5 {8 x. r8 p. S' x: t. K) g24002271.69189384
1.625

matlab效率已不及Forcal。

) Y- N3 p: x& l) }提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
( h0 ]' G; d! |: |+ }; ?9 {
% X/ C  ~8 P; }1 NForcal代码：

2. !using["sys"];
   
3. & r& `) O$ ^: J4 N! O
4. mvar:1 }. i% h9 J* f. Y8 H( {/ U
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组  y4 \+ T  W1 v, x# a/ A3 o4 J' u  H
   
8. t0=clock(),6 ~0 f4 A$ X: n3 f# c
   
10. k=0, while{k<10000,/ L: @; h5 g/ k4 ]; [5 ^/ U
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. 9 ?- a# b/ A( q, I1 g8 {
14.   k++8 m6 i4 F; R) V
    
16. },, @% V/ z8 t; [- b2 Y- S
    
18. [clock()-t0]/1000;
    
19. / P; ]$ e, b% ^\\" D5 J
20. 
    
21. # l9 n7 [8 w, S: `
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. ( h7 f/ T/ h+ M0 _6 J
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. % |\\" K, d5 Y# S, L8 P\\" t' Q
26. t0=clock(),. r& t1 T5 o6 ?2 E( Z
    
28. s=0,
    
29. 8 X$ @. w. _6 ]
30. k=0, while{k<10000,7 l3 K0 ]+ f7 f3 w% i
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. ! ?6 b! J: b5 u8 {9 \
34.   k++4 q0 d5 d1 D. t  Z. l5 f
    
36. },6 E! t  v& ?& g\\" v& [* ^
    
38. s;9 M* W1 W. p% R8 s
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">& r& `) O$ ^: J4 N! O</span> mvar:<font class=\\"jammer\\">1 }. i% h9 J* f. Y8 H( {/ U</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">  y4 \+ T  W1 v, x# a/ A3 o4 J' u  H</font><br /> t0=clock(),<font class=\\"jammer\\">6 ~0 f4 A$ X: n3 f# c</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">/ L: @; h5 g/ k4 ]; [5 ^/ U</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">9 ?- a# b/ A( q, I1 g8 {</span>   k++<font class=\\"jammer\\">8 m6 i4 F; R) V</font><br /> },<font class=\\"jammer\\">, @% V/ z8 t; [- b2 Y- S</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">/ P; ]$ e, b% ^\\" D5 J</span> <br /> <span style=\\"display:none\\"># l9 n7 [8 w, S: `</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">( h7 f/ T/ h+ M0 _6 J</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">% |\\" K, d5 Y# S, L8 P\\" t' Q</span> t0=clock(),<font class=\\"jammer\\">. r& t1 T5 o6 ?2 E( Z</font><br /> s=0,<br /> <span style=\\"display:none\\">8 X$ @. w. _6 ]</span> k=0, while{k<10000,<font class=\\"jammer\\">7 l3 K0 ]+ f7 f3 w% i</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">! ?6 b! J: b5 u8 {9 \</span>   k++<font class=\\"jammer\\">4 q0 d5 d1 D. t  Z. l5 f</font><br /> },<font class=\\"jammer\\">6 E! t  v& ?& g\\" v& [* ^</font><br /> s;<font class=\\"jammer\\">9 M* W1 W. p% R8 s</font><br /> [clock()-t0]/1000;

结果：
* c5 b. r% {+ Y- Y0.5620000000000001  //赋值用时间，秒
( O' }$ g8 E+ p) T, u4 k$ x. L$ Y20000000.
0.579      //取值用时间，秒