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

forcal

1、简单循环

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

1. clear all\" S7 R6 |& |  a, z  R
   
2. tic
   ! y\" e/ e$ ]# q; o1 j8 ?
3. s=0;+ r1 L6 ]6 I: K
   
4. for k=1:10000  N+ ]7 [8 F+ Z, p& E& C
   
5.   for i=1:1000  ~; `' P- g3 }; A
   
6.     s=s+1;
   0 n\" i) O0 P. m. C. s& t& x
7.   end
   0 U. j/ U( S+ B; w& g
8. end
   / N3 D- f& t& e! @. u% B2 ~
9. s
   ; \; E\" M+ N% J# u& W
10. toc' Y0 S4 E  {/ ]. d6 W& t+ }4 I3 y
    
11. 
    . i7 ^% F- {6 s* y  d
12. s =
    5 W# L( G6 c\" F2 R) g
13. 
    1 u& t' u2 F6 j, A; q
14.     100000004 r: Z7 Y* P% i, B  b
    
15. 
    + |7 Q! k$ I* O% o0 M! s8 I; Y
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
+ w5 N  k' J1 J2 g: W- N& s
Forcal没有JIT加速器，代码：

2. !using["sys"];( T8 ~# B\\" k9 E& T\\" h6 O: c+ H
   
4. mvar:
   
5. 3 ]& c) h( O( ?1 i$ h& O  i
6. t0=clock(),
   
7. 5 i+ p& e, w' [. Q5 o# j  y! s! g1 X
8. s=0,/ |! Y! v+ R/ N2 S) C. G0 s8 O
   
10. k=0, while{k<10000,
    
11. + B9 _8 j  ^  h& W4 X$ ?2 M! o& a
12.   i=0, while{i<1000,s=s+1, i++},
    
13. / r\\" w4 W' K' K- z5 S3 ~  w2 T0 q
14.   k++9 {0 y\\" ~- d* s/ q6 [
    
16. },$ x% P; j7 a$ l
    
18. s;
    
19. + G: u2 z4 S  Y3 x! G
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">( T8 ~# B\\" k9 E& T\\" h6 O: c+ H</font><br /> mvar:<br /> <span style=\\"display:none\\">3 ]& c) h( O( ?1 i$ h& O  i</span> t0=clock(),<br /> <span style=\\"display:none\\">5 i+ p& e, w' [. Q5 o# j  y! s! g1 X</span> s=0,<font class=\\"jammer\\">/ |! Y! v+ R/ N2 S) C. G0 s8 O</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">+ B9 _8 j  ^  h& W4 X$ ?2 M! o& a</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">/ r\\" w4 W' K' K- z5 S3 ~  w2 T0 q</span>   k++<font class=\\"jammer\\">9 {0 y\\" ~- d* s/ q6 [</font><br /> },<font class=\\"jammer\\">$ x% P; j7 a$ l</font><br /> s;<br /> <span style=\\"display:none\\">+ G: u2 z4 S  Y3 x! G</span> [clock()-t0]/1000;

结果：
" M) f" P& K+ j! Z' K8 r) _1 Z10000000.
5 u! e* S9 n! ?0.704  秒
, Y$ s& U2 ?8 m7 V3 z5 v# \
$ i7 a1 |4 Q) O% a二者运行效率相差7倍多。

3 V; p: l) M$ H" `: e+ a" V8 \2、循环中稍微加点东西

8 u+ r; }" b: r+ \: b  x, HMatlab代码：

1. clear all
   & E6 Y. U7 h\" f3 t, H
2. tic7 a8 A$ F: I- |8 X8 N/ J$ _$ `$ c
   
3. s=0;
   . G) l# X) @: z0 J* N
4. for k=1:1000) V0 O! u5 ^& Z# C/ o
   
5.   for i=1:1000
   ' U) O- I/ \( k) [: S* T1 a
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);
   / U4 v) Z* k! _4 x/ Y\" L
7.   end2 I9 N/ p0 J/ l. C# H5 @
   
8. end- y0 z1 B; R; N+ T/ R- M7 |\" |
   
9. s1 d7 e2 y( G6 T) I. w0 V\" i
   
10. toc# h# b; {; c6 {( \6 @$ w1 o- _, r: d
    
11. 7 q  m5 r* R+ B2 t
    
12. s =0 k* S# ?+ ~9 i) j
    
13. 
    8 L  z2 _0 L' U5 }$ c4 k% }
14.   1.1009e+007
    7 c7 n7 h5 ?( r
15. ; U% \) w/ j- O5 @; M
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 7 }  Y3 T( }; ^# H2 T# p# I
4. mvar:4 c+ e- i7 a, Q- b! _6 Y
   
6. t0=clock(),2 M  U. @. f) H  C
   
8. s=0,
   
9. : b8 a. D, t\\" |- p1 \' U) v
10. k=1, while{k<=1000,
    
11. % e% j. N/ q* O
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. 6 L9 |+ r5 o4 x' Q3 J/ Q\\" c  E\\" M
14.   k++) \) J7 q2 O5 S. {8 ~
    
16. },
    
17. ) j. D( r, f# w; \- N* b4 T- Z
18. s;0 g1 R$ V5 u4 K0 z1 F! v
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">7 }  Y3 T( }; ^# H2 T# p# I</span> mvar:<font class=\\"jammer\\">4 c+ e- i7 a, Q- b! _6 Y</font><br /> t0=clock(),<font class=\\"jammer\\">2 M  U. @. f) H  C</font><br /> s=0,<br /> <span style=\\"display:none\\">: b8 a. D, t\\" |- p1 \' U) v</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">% e% j. N/ q* O</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\\">6 L9 |+ r5 o4 x' Q3 J/ Q\\" c  E\\" M</span>   k++<font class=\\"jammer\\">) \) J7 q2 O5 S. {8 ~</font><br /> },<br /> <span style=\\"display:none\\">) j. D( r, f# w; \- N* b4 T- Z</span> s;<font class=\\"jammer\\">0 g1 R$ V5 u4 K0 z1 F! v</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
1.031
7 Z. Q2 g2 \* @7 \# w
二者效率相差已比较小。

3、再增加一个函数调用

: J7 Z& ?9 X6 B  F$ PMatlab代码：

1. clear all
   0 |; v9 J( n4 H2 l
2. f=@(x)x+1;
   + G, B% t: N; h' l8 C) h
3. tic  X1 r3 j, D) u/ b
   
4. s=0;1 J' P\" b; q& C: r: r
   
5. for k=1:1000) X$ ~( a2 A8 p. R
   
6.   for i=1:1000* F3 `\" I5 i* v7 M2 O3 l5 A7 X
   
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 d- M! r/ L2 \  v4 {: w& f
   
8.   end5 P7 z' Y( {; O. O6 F6 W# z
   
9. end
   / s% _0 J7 g7 X: O# J
10. s
    ) r& y6 F' ]- g* H\" j9 M# f. R
11. toc4 s& l1 k4 b0 o3 X% {& ]% d% ^
    
12. # c% r% j9 e! w& ?0 Z0 r1 A$ y
    
13. s =/ w4 t) f# t$ E
    
14. & G\" r8 W. c! ~9 `
    
15.   5.1251e+0085 k$ K$ k$ s! x, r
    
16. 
    5 E) O$ L' f% D, n
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];  I% Y+ E6 e4 C: N1 v; E4 A$ p& h
   
4. f(x)=x+1;2 m6 N9 P3 G- v7 w
   
6. mvar:! W+ T9 v! y6 B' R* f\\" b& W3 g) K
   
8. t0=clock(),; j$ f5 w* E1 Z; s
   
10. s=0,
    
11. 4 G( g8 t3 W. J8 M2 R, x& J
12. k=1, while{k<=1000,
    
13. ' k, J( k\\" ]' U3 m6 h: z# @
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. - k+ N( v8 V) {* w* `! ]
16.   k++( K5 C\\" Q. Q5 T% O; |4 H
    
18. },( c! e( H+ X; g1 G: x
    
20. s;3 n& c$ W$ y6 y* p  F+ J' n
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">  I% Y+ E6 e4 C: N1 v; E4 A$ p& h</font><br /> f(x)=x+1;<font class=\\"jammer\\">2 m6 N9 P3 G- v7 w</font><br /> mvar:<font class=\\"jammer\\">! W+ T9 v! y6 B' R* f\\" b& W3 g) K</font><br /> t0=clock(),<font class=\\"jammer\\">; j$ f5 w* E1 Z; s</font><br /> s=0,<br /> <span style=\\"display:none\\">4 G( g8 t3 W. J8 M2 R, x& J</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">' k, J( k\\" ]' U3 m6 h: z# @</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\\">- k+ N( v8 V) {* w* `! ]</span>   k++<font class=\\"jammer\\">( K5 C\\" Q. Q5 T% O; |4 H</font><br /> },<font class=\\"jammer\\">( c! e( H+ X; g1 G: x</font><br /> s;<font class=\\"jammer\\">3 n& c$ W$ y6 y* p  F+ J' n</font><br /> [clock()-t0]/1000;

结果：
512508953.6658993
1.078
3 G" {# w  n2 X1 {
/ a$ l" d. s7 P$ S4 ~& N; v$ [可以看出，Matlab效率已不及Forcal。
8 O: H) `2 q& ]% a故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

Matlab代码：

1. >> clear all6 i( ^8 G: W. i4 g5 f) o& k
   
2. a = linspace(1,1,1000);
   4 R- `7 o- V& j# s
3. tic
   ; R6 k6 h9 ?5 j, V4 m& d2 N
4. for k=1:10000) J- y! b1 p* _6 f; j
   
5.   for i=1:1000
   9 @0 Y4 p! E$ R0 \6 U% i2 e' W/ C
6.     a(i)=2;\" J! P' q8 v: _  K
   
7.   end, O; y) P: U* v$ l# j2 X
   
8. end7 I7 K% v3 d6 z1 c
   
9. toc7 q( [\" j$ u% O+ w+ F
   
10. 6 M7 s* R1 H6 Y7 ?4 V) V( x
    
11. tic
    6 _2 J& e9 L! s5 e  O: x+ z- z- u% G6 G
12. s=0;
    # Y* l6 {# G' N: G: O
13. for k=1:10000. M2 ?: J9 p! G3 z3 O
    
14.   for i=1:1000
    / X& ~7 k  u+ w5 h) l  G1 M  h
15.     s=s+a(i);
      X# i: W: ~, S* m+ M
16.   end$ D1 @- w# b- q\" f! Q
    
17. end% c; |4 n* J- s: F& x
    
18. s* k8 m2 E, ~4 B' V8 e
    
19. toc
    : U. D2 X& R- N+ V/ J4 A
20. Elapsed time is 0.140658 seconds.
    $ `6 d( n* }3 q+ P% Z
21. 
    7 u* r* z% R  P( @2 p
22. s =
    / N0 q8 f+ Q% q* N5 T3 T4 Q3 }
23. + m/ h2 w5 p5 ]
    
24.     20000000: z  Z4 ?3 a1 @) ^
    
25. ! a' o\" K' R8 w6 ]4 Z6 f9 q4 S, s6 X
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];
   
3. 0 W, s' \5 y5 O\\" S( N% G: O3 @
4. mvar:, }\\" o8 a( m& F& W0 t) L
   
6. main0(::a)=
   
7. + C8 d9 _. Z+ }
8. a = linspace(1,1,1000),
   
9. ( x, O6 [7 N, C* I4 G3 _
10. t0=clock(),
    
11. # U- g1 [8 G\\" {& m+ j2 q
12. k=0, while{k<10000,
    
13. : J4 `5 d* @& b9 @4 \; C
14.   i=0, while{i<1000, a(i)=2, i++},! W  Q- a* Q4 S6 a  w
    
16.   k++
    
17. 8 {9 R$ u$ y- N5 Y9 J  F
18. },
    
19. & ]! x; p* ^' G  u) i
20. [clock()-t0]/1000;4 K$ S( f. h, b! c+ |
    
22. 
    
23. 5 J  |( V# K& L0 h
24. main1(::a)=
    
25. 9 c, n3 n$ A  J: F: ]$ P  b9 q0 Y# \' S
26. t0=clock(), s=0,. q/ Z2 x- @! E, D
    
28. k=0, while{k<10000,* K/ ^. ]! q$ _: ~( J' F
    
30.   i=0, while{i<1000, s=s+a(i), i++},\\" G, `* a, P/ H! x& S4 t6 s& y2 X
    
32.   k++: B- S4 k* D2 X+ C% z\\" B# [
    
34. },
    
35. 9 i0 C% L0 E3 O/ I
36. s;* Y: g/ Z- b' c- m4 C
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">0 W, s' \5 y5 O\\" S( N% G: O3 @</span> mvar:<font class=\\"jammer\\">, }\\" o8 a( m& F& W0 t) L</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">+ C8 d9 _. Z+ }</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">( x, O6 [7 N, C* I4 G3 _</span> t0=clock(),<br /> <span style=\\"display:none\\"># U- g1 [8 G\\" {& m+ j2 q</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">: J4 `5 d* @& b9 @4 \; C</span>   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">! W  Q- a* Q4 S6 a  w</font><br />   k++<br /> <span style=\\"display:none\\">8 {9 R$ u$ y- N5 Y9 J  F</span> },<br /> <span style=\\"display:none\\">& ]! x; p* ^' G  u) i</span> [clock()-t0]/1000;<font class=\\"jammer\\">4 K$ S( f. h, b! c+ |</font><br /> <br /> <span style=\\"display:none\\">5 J  |( V# K& L0 h</span> main1(::a)=<br /> <span style=\\"display:none\\">9 c, n3 n$ A  J: F: ]$ P  b9 q0 Y# \' S</span> t0=clock(), s=0,<font class=\\"jammer\\">. q/ Z2 x- @! E, D</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">* K/ ^. ]! q$ _: ~( J' F</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">\\" G, `* a, P/ H! x& S4 t6 s& y2 X</font><br />   k++<font class=\\"jammer\\">: B- S4 k* D2 X+ C% z\\" B# [</font><br /> },<br /> <span style=\\"display:none\\">9 i0 C% L0 E3 O/ I</span> s;<font class=\\"jammer\\">* Y: g/ Z- b' c- m4 C</font><br /> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒

－－－－－－－－－－
& p1 B& V" i/ L) s. \% i7 j
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. 1 h9 _; s( g\\" @8 M) a
4. mvar:1 [$ X2 v+ L% S3 C* A
   
6. main0(::a)=
   
7. 8 l* c* i4 w& @, o+ t( Y' t) C- |
8. a = linspace(1,1,1000),  \! K9 l: q/ d
   
10. t0=clock(),3 [0 O: B' a$ _9 ]9 z$ K
    
12. k=0, while{k<10000,$ c* _1 K$ |2 f
    
14.   i=0, while{i<1000, A(a,i)=2, i++},- u; ?( w% M6 P6 G$ k% Z$ _- a: w4 L+ W
    
16.   k++! B) N1 }. }4 {, J. A+ s
    
18. },
    
19. ! c/ x7 D( L' t4 j
20. [clock()-t0]/1000;! L5 M& t$ b! z/ ?4 H- ^
    
22. 
    
23. 4 Z0 H& w6 p; a
24. main1(::a)=; q  Z) W( u; F) Y
    
26. t0=clock(), s=0,5 \+ q7 z% J, r, O9 G! X: Q
    
28. k=0, while{k<10000,/ [* N2 v- Q$ C  U! ^1 T4 t+ \
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},; O1 F\\" L, ]- l! \0 P
    
32.   k++' p, W7 {( d# u* ^
    
34. },. H! e, X! p% J& J0 B3 V  `) G
    
36. s;
    
37. 2 N( r4 F2 c* w# n, v
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">1 h9 _; s( g\\" @8 M) a</span> mvar:<font class=\\"jammer\\">1 [$ X2 v+ L% S3 C* A</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">8 l* c* i4 w& @, o+ t( Y' t) C- |</span> a = linspace(1,1,1000),<font class=\\"jammer\\">  \! K9 l: q/ d</font><br /> t0=clock(),<font class=\\"jammer\\">3 [0 O: B' a$ _9 ]9 z$ K</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">$ c* _1 K$ |2 f</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">- u; ?( w% M6 P6 G$ k% Z$ _- a: w4 L+ W</font><br />   k++<font class=\\"jammer\\">! B) N1 }. }4 {, J. A+ s</font><br /> },<br /> <span style=\\"display:none\\">! c/ x7 D( L' t4 j</span> [clock()-t0]/1000;<font class=\\"jammer\\">! L5 M& t$ b! z/ ?4 H- ^</font><br /> <br /> <span style=\\"display:none\\">4 Z0 H& w6 p; a</span> main1(::a)=<font class=\\"jammer\\">; q  Z) W( u; F) Y</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">5 \+ q7 z% J, r, O9 G! X: Q</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">/ [* N2 v- Q$ C  U! ^1 T4 t+ \</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">; O1 F\\" L, ]- l! \0 P</font><br />   k++<font class=\\"jammer\\">' p, W7 {( d# u* ^</font><br /> },<font class=\\"jammer\\">. H! e, X! p% J& J0 B3 V  `) G</font><br /> s;<br /> <span style=\\"display:none\\">2 N( r4 F2 c* w# n, v</span> [clock()-t0]/1000;

结果：
10.469
; B& t! V0 ~" [9 E20000000.
3.781
1 L5 n8 b7 k/ e5 m
－－－－－－－

" I" a6 C5 T# Y# p! a在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
# y" d; @) Z4 _& x% G( L在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

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

2 S; d- Z5 @0 W, i& G" ]; C5、在循环中加点东西

matlab代码：

1. >> clear all
   4 C) z/ s9 S$ y0 Z9 v9 h3 I
2. f=@(x)x+1;
   ; L$ S3 C: z0 K' p4 ~
3. a = linspace(2,2,1000);; h: |. D3 ^( n/ g6 W5 f6 X
   
4. tic/ W1 u/ X) I; p* q: d, ]
   
5. s=0;
   ! i& @' N- I\" O8 ]5 f- Y
6. for k=1:1000
   3 f3 j& C! j$ ~- }$ {1 {\" U+ g* |
7.   for i=1:1000
   * a: n0 h$ n! t
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));$ R7 V: P8 ^# N& u% H$ W  |9 r5 H
   
9.   end* j- e* h3 o+ S* O. ^3 T. B4 C
   
10. end# z: W% d, T; U
    
11. s
    4 U4 W1 E( |; b1 [8 ^- y4 s
12. toc
    - q- O0 m+ Y0 v; q' p
13. 
    6 C  d  r6 ]4 {
14. s =
    ( y+ E' I9 S7 M5 j! S! L4 ]$ x7 t
15. * P6 y5 @2 n4 s% \' J' D. ]
    
16.   2.4002e+007' ~4 I. {8 d7 U$ z1 |8 A9 O+ `
    
17. 
    8 I/ u9 E3 x5 g4 f. [) f/ h: Y
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];4 j0 V: X8 Y. d- x1 }) @- F
   
4. mvar:7 u/ R9 e) E+ `; R3 {! ]! y3 s
   
6. f(x)=x+1;
   
7. 8 u; Y/ Z- Q2 \$ W# s1 V
8. main0(::a)=
   
9. # x2 B\\" A5 }$ L
10. a = linspace(2,2,1000),; q0 \$ E  v7 q3 \: ?1 {* F7 R
    
12. t0=clock(), s=0,
    
13.   V, i\\" q  [  B% b, g0 o
14. k=0, while{k<1000,
    
15. ' F  t7 k$ }2 O
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. - s8 g\\" ~. |\\" g5 c
18.   k++2 C\\" _( M# U  ^, v* M- Z
    
20. },
    
21. 1 r, F9 l! ~, x' _2 ?! t
22. s;1 m+ K# [; u: x$ w% N
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">4 j0 V: X8 Y. d- x1 }) @- F</font><br /> mvar:<font class=\\"jammer\\">7 u/ R9 e) E+ `; R3 {! ]! y3 s</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">8 u; Y/ Z- Q2 \$ W# s1 V</span> main0(::a)=<br /> <span style=\\"display:none\\"># x2 B\\" A5 }$ L</span> a = linspace(2,2,1000),<font class=\\"jammer\\">; q0 \$ E  v7 q3 \: ?1 {* F7 R</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">  V, i\\" q  [  B% b, g0 o</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">' F  t7 k$ }2 O</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\\">- s8 g\\" ~. |\\" g5 c</span>   k++<font class=\\"jammer\\">2 C\\" _( M# U  ^, v* M- Z</font><br /> },<br /> <span style=\\"display:none\\">1 r, F9 l! ~, x' _2 ?! t</span> s;<font class=\\"jammer\\">1 m+ K# [; u: x$ w% N</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
5.172

; J2 i+ }3 r) [6 t# q2 F$ e5 h二者效率现在只有一倍差距。
) A& H5 A; m, _1 A2 y8 K, k
* ^+ j5 d: E. [4 [# c  ^. b9 U－－－－－－－－－－
' ^$ V% G; r+ [0 x" d2 W' Y& Z
不过，显然matlab和Forcal的程序应该这样写：

5 u3 X% B0 k1 d  p5 m$ w, Q* |) ematlab程序：

1. clear all
   / g% t% ?6 [8 x; K
2. f=@(x)x+1;2 W, l9 P( M3 z1 G! f9 Z0 Y2 P
   
3. a = linspace(2,2,1000);
   % G9 N\" U. V$ m- N* ]
4. tic* U# d3 k# Z! H+ h1 u1 E$ _# W1 z
   
5. s=0;) M2 d3 w: C& z1 f* x8 v
   
6. for k=1:10006 a& ?2 A* e3 O! Q2 \: L2 \  j
   
7.   for i=1:1000$ C  h) t  H0 e! v$ b7 u) M
   
8.     t=a(i);$ x) k5 V9 L1 J. {6 k1 I
   
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);( E\" G3 ~& Y3 y: P, K4 D1 N
   
10.   end
    : g( o: S2 |- x2 R
11. end\" }+ o: }2 F  W# h- u
    
12. s% ~4 d% A) u5 T. q5 B
    
13. toc* ], v$ ?9 `+ ]6 _' |8 O+ h& d
    
14. 1 {7 T4 E, P; r+ p$ {: g' V
    
15. s =8 \3 ~7 D6 Z6 F
    
16. 
    : U\" F8 P4 p8 Q; x% g+ U' L
17.   2.4002e+007
    : Q5 V6 Q( ?! j! @
18. . c5 }4 @' J  C, @7 R
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];/ J. c1 b0 e2 S0 ^1 n/ C
   
4. mvar:
   
5. 8 a+ i* w9 v1 S+ U
6. f(x)=x+1;! p$ Y# Y- O( ^) n8 z
   
8. main0(::a)=1 \( `* O, v9 }$ @2 w
   
10. a = linspace(2,2,1000),/ p# k: A1 n1 k: z  b
    
12. t0=clock(), s=0,% }3 x2 r\\" I( H, P
    
14. k=0, while{k<1000,
    
15. $ m! S5 k$ k+ i* x7 `\\" k
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++},& W8 h5 C' y; I( k
    
18.   k++
    
19. 4 y& U* D+ ^# B. m\\" Q+ |
20. },
    
21. ; T0 M! k& s, g$ W8 r
22. s;
    
23. 7 o' a- O: l5 U\\" z9 H1 U
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">/ J. c1 b0 e2 S0 ^1 n/ C</font><br /> mvar:<br /> <span style=\\"display:none\\">8 a+ i* w9 v1 S+ U</span> f(x)=x+1;<font class=\\"jammer\\">! p$ Y# Y- O( ^) n8 z</font><br /> main0(::a)=<font class=\\"jammer\\">1 \( `* O, v9 }$ @2 w</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">/ p# k: A1 n1 k: z  b</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">% }3 x2 r\\" I( H, P</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">$ m! S5 k$ k+ i* x7 `\\" k</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++},<font class=\\"jammer\\">& W8 h5 C' y; I( k</font><br />   k++<br /> <span style=\\"display:none\\">4 y& U* D+ ^# B. m\\" Q+ |</span> },<br /> <span style=\\"display:none\\">; T0 M! k& s, g$ W8 r</span> s;<br /> <span style=\\"display:none\\">7 o' a- O: l5 U\\" z9 H1 U</span> [clock()-t0]/1000;

结果：
24002271.69189384
" x5 {5 t( _9 u4 n1.625

4 l2 ^* m$ b& Z  b4 A% X5 Rmatlab效率已不及Forcal。

% R8 J4 K9 v5 e( }, z" q' S3 S提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
  ^2 x. m, ^6 S+ k7 s" F
. D" I. u' R' O* z4 C, ^Forcal代码：

2. !using["sys"];( W8 q% ^& R& J. e. C# C6 J
   
4. mvar:
   
5. 9 R\\" u8 z$ m2 q+ T  ]
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. \\" ^: e$ v6 l8 F\\" B. A. A9 O. O
8. t0=clock(),
   
9. 8 ?- b7 ~& H+ d& E1 @9 J# R1 A' p# ~9 w
10. k=0, while{k<10000,
    
11.   R9 ?$ i/ a0 t4 ^' \5 ]2 [/ s8 g5 y
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值: E  X( w2 x+ @4 p2 x0 N\\" {& C; v# E% ~
    
14.   k++
    
15. ' P) ^: n9 E2 |( A, E* g# Q2 U' B
16. },4 Q& s% h4 C( l3 s
    
18. [clock()-t0]/1000;
    
19. / y8 ]( b3 Z; R0 i. L4 b
20. 
    
21. ) W+ I2 @9 o6 V# a9 ]- s1 a' d
22. a="\&\[1000]"/4,    //用字符串模拟静态数组& t! M: M2 [! D  R1 c9 w/ Q
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值6 V7 T, @& @1 U
    
26. t0=clock(),
    
27. 3 D- B$ c; G. b& N3 B\\" M
28. s=0,7 R- G, Q% Q8 @. m\\" A
    
30. k=0, while{k<10000,, G: e9 E! L4 h
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. % o& v; s7 U$ z& l- d\\" L
34.   k++' W2 e% v4 k* N  m4 f
    
36. },
    
37.   e4 l- q  ?2 _( v3 e5 O+ p
38. s;
    
39. , h1 L\\" `\\" a) M
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">( W8 q% ^& R& J. e. C# C6 J</font><br /> mvar:<br /> <span style=\\"display:none\\">9 R\\" u8 z$ m2 q+ T  ]</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">\\" ^: e$ v6 l8 F\\" B. A. A9 O. O</span> t0=clock(),<br /> <span style=\\"display:none\\">8 ?- b7 ~& H+ d& E1 @9 J# R1 A' p# ~9 w</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">  R9 ?$ i/ a0 t4 ^' \5 ]2 [/ s8 g5 y</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">: E  X( w2 x+ @4 p2 x0 N\\" {& C; v# E% ~</font><br />   k++<br /> <span style=\\"display:none\\">' P) ^: n9 E2 |( A, E* g# Q2 U' B</span> },<font class=\\"jammer\\">4 Q& s% h4 C( l3 s</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">/ y8 ]( b3 Z; R0 i. L4 b</span> <br /> <span style=\\"display:none\\">) W+ I2 @9 o6 V# a9 ]- s1 a' d</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">& t! M: M2 [! D  R1 c9 w/ Q</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">6 V7 T, @& @1 U</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">3 D- B$ c; G. b& N3 B\\" M</span> s=0,<font class=\\"jammer\\">7 R- G, Q% Q8 @. m\\" A</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">, G: e9 E! L4 h</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">% o& v; s7 U$ z& l- d\\" L</span>   k++<font class=\\"jammer\\">' W2 e% v4 k* N  m4 f</font><br /> },<br /> <span style=\\"display:none\\">  e4 l- q  ?2 _( v3 e5 O+ p</span> s;<br /> <span style=\\"display:none\\">, h1 L\\" `\\" a) M</span> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
20000000.
9 i3 ]/ x4 u& S% x; K2 a0 M0.579      //取值用时间，秒