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

forcal

1、简单循环
3 V3 Q: s5 R) B' @: ?& ~3 ]( G
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   4 u. w+ }0 l* n+ [; V
2. tic
   0 w6 y  V) w7 H' D, B
3. s=0;
     u: o8 k) l\" z
4. for k=1:10000* z2 o% D0 q0 ~9 [9 H
   
5.   for i=1:1000
     H7 ?, Y. R( @+ B% @
6.     s=s+1;; d) _- H/ h* h. P; f. V
   
7.   end
   ; u' b6 U  R( ?$ ~\" a  z8 d6 ]
8. end
   - U9 S1 ?( a' v/ G& E
9. s  Q0 D1 D& R1 [
   
10. toc
    ! Z7 r8 _/ H9 J( c% ~
11. ( {' _: T( R0 L0 K2 K) c% `4 d. P9 N
    
12. s =- @) W2 G9 l$ c6 k  }1 h
    
13. 1 z2 \; h\" m$ ]; L' j' N, x* H
    
14.     100000002 @7 w' r+ Q7 M- b' J2 P
    
15. 
    8 z: W5 M  T( I- }
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
- I/ K6 `- N4 D
Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. / G) ?8 r$ c1 @1 `6 h( W+ q# H
4. mvar:
   
5. 2 R! B9 |( ^$ g1 \
6. t0=clock(),
   
7. ' L  q  a6 m\\" o) O
8. s=0,5 ^4 j5 }$ j* C7 f& A3 A5 f  F
   
10. k=0, while{k<10000,' M# W5 W3 z7 q; i8 y9 h3 L
    
12.   i=0, while{i<1000,s=s+1, i++},$ I! Y, G1 T* c* B7 L7 D
    
14.   k++
    
15. / a8 J7 z6 g/ o' W\\" O1 e
16. },
    
17. 6 y4 }* D% `+ k\\" J# f
18. s;3 F9 H( a, G4 H9 g+ m/ f' q
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">/ G) ?8 r$ c1 @1 `6 h( W+ q# H</span> mvar:<br /> <span style=\\"display:none\\">2 R! B9 |( ^$ g1 \</span> t0=clock(),<br /> <span style=\\"display:none\\">' L  q  a6 m\\" o) O</span> s=0,<font class=\\"jammer\\">5 ^4 j5 }$ j* C7 f& A3 A5 f  F</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">' M# W5 W3 z7 q; i8 y9 h3 L</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">$ I! Y, G1 T* c* B7 L7 D</font><br />   k++<br /> <span style=\\"display:none\\">/ a8 J7 z6 g/ o' W\\" O1 e</span> },<br /> <span style=\\"display:none\\">6 y4 }* D% `+ k\\" J# f</span> s;<font class=\\"jammer\\">3 F9 H( a, G4 H9 g+ m/ f' q</font><br /> [clock()-t0]/1000;

结果：
10000000.
0.704  秒

) C: L9 _( Q9 z$ M3 d3 O# V二者运行效率相差7倍多。
; P4 V; c+ D- Q* o
" q; G9 m" x, r0 D2、循环中稍微加点东西
- c# P/ s# R4 `  E+ c, w
6 U; x# ~+ o+ O. |; tMatlab代码：

1. clear all
   ! w4 ?7 P9 D8 v: `
2. tic
   ; u, h- ^: B& I  M. a5 O
3. s=0;, v8 l% k4 A& i9 q3 B( R
   
4. for k=1:1000
   * S9 _9 m5 S  h& A* ~0 W
5.   for i=1:1000  g3 C: Q\" V! X$ S' ]$ H! e. ]5 c: T
   
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);) k% c, ~# |3 E1 t: M1 T( A
   
7.   end
   0 r- I$ x! q' f- `# m. a8 a) _, G& N
8. end
   # _5 F# u9 }/ g7 r0 u; E! A8 ?
9. s
   ' }9 j  [0 e7 \, w6 ]! u
10. toc. H% o0 S) l) o) r
    
11. 
    4 L1 t8 A+ E$ k* J$ O6 e$ R8 t
12. s =
    7 t' e+ G6 x  G
13. - S% p9 Q& p. J# x% v# ?
    
14.   1.1009e+007
    ' ^$ W/ I8 @; A* B: }: J' X7 [' C/ t
15. 
    1 |8 }- P: b* l7 w; ?
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 3 g2 P( l, j1 R  H
4. mvar:
   
5. - m. R1 v+ [. {% H  D# A7 I
6. t0=clock(),; ?3 Z) F4 c% |4 T; _. X* i
   
8. s=0,' a' b9 Z- a& a- z$ ?) N
   
10. k=1, while{k<=1000,\\" v+ w3 \) ]' ]
    
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. 9 w( u3 X& u\\" q
14.   k++
    
15. % V\\" l( i) _\\" E- q6 v: C
16. },
    
17. # f% Z% ]* t& ^, |: B
18. s;
    
19. 7 _* Q) a9 y/ l0 f' M# W
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">3 g2 P( l, j1 R  H</span> mvar:<br /> <span style=\\"display:none\\">- m. R1 v+ [. {% H  D# A7 I</span> t0=clock(),<font class=\\"jammer\\">; ?3 Z) F4 c% |4 T; _. X* i</font><br /> s=0,<font class=\\"jammer\\">' a' b9 Z- a& a- z$ ?) N</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">\\" v+ w3 \) ]' ]</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\\">9 w( u3 X& u\\" q</span>   k++<br /> <span style=\\"display:none\\">% V\\" l( i) _\\" E- q6 v: C</span> },<br /> <span style=\\"display:none\\"># f% Z% ]* t& ^, |: B</span> s;<br /> <span style=\\"display:none\\">7 _* Q) a9 y/ l0 f' M# W</span> [clock()-t0]/1000;

结果：
/ N. e" u7 o( E2 M11008953.66602346
1.031

/ a* r& s, V( n8 a- Y- U二者效率相差已比较小。
9 t) c% ^( m- G
0 Z9 S' ~+ V6 V! }( @3、再增加一个函数调用

Matlab代码：

1. clear all\" ~/ a/ V\" V: F
   
2. f=@(x)x+1;
   / Y* o5 D$ I2 S2 m+ L; L9 Z
3. tic1 u. u, W, O) {. Z; g- u
   
4. s=0;8 Z9 |, |3 `: n5 \* e) A
   
5. for k=1:1000
   6 ?1 q& V2 m5 \. ?* o  y, N  Z0 I
6.   for i=1:10008 D, J! w  a5 H/ L
   
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);: R# O3 Q0 N) C; B2 ]4 Z
   
8.   end\" n( n) O# T2 x3 }6 Z) o# ~
   
9. end
   0 j+ G9 K0 |/ Q* {0 [4 j! Z
10. s' I8 {) o. T+ x) q5 N
    
11. toc! B. B: V$ R  S* |0 `6 o* K
    
12. 
    ' H4 u: Q7 N. \% o  m- `
13. s =# d! O3 f- ]- U: X9 S2 T
    
14. 3 q+ l4 a( r3 H2 e) E
    
15.   5.1251e+008
    * g4 d0 t7 {0 [# Z, Z
16. \" ~+ P) i8 _- A7 q
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. : ~\\" p9 D, x$ \
4. f(x)=x+1;
   
5. 9 a8 g2 b0 n# N+ }8 M& f
6. mvar:3 M4 q) Q( r% i$ T# i/ W+ S
   
8. t0=clock(),
   
9. & m$ D2 s4 Y3 z8 i: J  a* U
10. s=0,7 x) X1 u6 V/ `( Y# j
    
12. k=1, while{k<=1000,
    
13. ( f7 k7 i\\" r  j
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. / c9 E, ^% V, N. u
16.   k++
    
17. 9 M1 s: b' ^+ E
18. },
    
19. ) A/ }4 K2 F/ H  b5 N
20. s;
    
21.   |. r8 L0 \4 T7 D
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">: ~\\" p9 D, x$ \</span> f(x)=x+1;<br /> <span style=\\"display:none\\">9 a8 g2 b0 n# N+ }8 M& f</span> mvar:<font class=\\"jammer\\">3 M4 q) Q( r% i$ T# i/ W+ S</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">& m$ D2 s4 Y3 z8 i: J  a* U</span> s=0,<font class=\\"jammer\\">7 x) X1 u6 V/ `( Y# j</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">( f7 k7 i\\" r  j</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\\">/ c9 E, ^% V, N. u</span>   k++<br /> <span style=\\"display:none\\">9 M1 s: b' ^+ E</span> },<br /> <span style=\\"display:none\\">) A/ }4 K2 F/ H  b5 N</span> s;<br /> <span style=\\"display:none\\">  |. r8 L0 \4 T7 D</span> [clock()-t0]/1000;

结果：
* g" H# \7 B* q# [/ W7 H512508953.6658993
8 ]$ C- p4 o; g' P+ x% J1.078

可以看出，Matlab效率已不及Forcal。
% h# J3 _1 ?# b' N故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
9 H. k2 e  Q) `5 Z+ n; \# b0 X3 R7 s
8 F" N& g# k9 L. i  c% D9 [Matlab代码：

1. >> clear all; b; q) Q$ z+ F, O6 Y. V/ b
   
2. a = linspace(1,1,1000);
   4 I' e& z2 h, J) O
3. tic
   * ?3 W- w4 w( Y0 p& W
4. for k=1:10000% u% w4 L6 N! y
   
5.   for i=1:1000$ c\" x1 T6 L2 {1 D9 r/ e
   
6.     a(i)=2;
   & D2 K& j4 y8 F
7.   end
   # p+ z% u6 M/ D' M6 n7 O2 w+ K
8. end5 V7 J0 }* _$ ]1 L  \\" K
   
9. toc
   * {9 @; N$ T6 a1 O
10. 
    # ^6 z7 ~# S. @8 N8 J' f) D
11. tic
    ; e/ }\" R/ q7 M+ P
12. s=0;
    2 E- x3 ]) A' j$ D2 q* s2 l% b7 p
13. for k=1:10000
    ! s* N+ {% T* d' d9 ?
14.   for i=1:1000
    ) e3 j4 n3 M3 D% Y: D
15.     s=s+a(i);\" @! D/ ^4 M+ o4 J# n4 ?4 u6 q
    
16.   end
    * w3 {; O$ H9 U/ M$ N' [' Z
17. end\" l! H5 M; Y! Y: e3 C
    
18. s6 U4 s\" [% D3 O. ~
    
19. toc+ I: W3 z4 G\" b! C$ a9 _. s
    
20. Elapsed time is 0.140658 seconds.
    3 h3 {- b* T# Z+ g
21. 
    2 n; Q4 Q, u- ?
22. s =
    & z8 g' F0 I$ e, G- l6 s\" J  |
23. 5 c5 z  r9 X! _7 B
    
24.     20000000/ `4 X/ W  U- h1 K. q* W
    
25. + h* d! Z8 g$ ~& {; @4 d
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];, L4 I4 s8 n# G- H: B: D) @; O& `
   
4. mvar:7 n+ y\\" d1 y+ T+ `
   
6. main0(::a)=
   
7. & Y- U+ ?! |8 b\\" h  R% N! |
8. a = linspace(1,1,1000),# ]; D. I. T. i+ P9 k8 a
   
10. t0=clock(),
    
11. 9 \& x) {' [$ v( g9 h/ l- y# F
12. k=0, while{k<10000,8 \3 a6 S: G# E% [8 r: x
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. , s2 N' }2 D% t1 ?: f
16.   k++/ j/ s8 Y2 T# O# e/ j
    
18. },; p4 P  N- B! F5 F7 k
    
20. [clock()-t0]/1000;, P; I; d  \5 d
    
22. , l1 y: w: d. I1 P; T
    
24. main1(::a)=) l+ q. q. f5 {# c* T
    
26. t0=clock(), s=0,
    
27. 9 Y+ Y% [* v. R. F
28. k=0, while{k<10000,
    
29. 3 b  v. }\\" H! d
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. ! T' {; L. ]8 l, [' X
32.   k++
    
33. + w7 p+ r' b) w6 K/ o
34. },
    
35. 6 p# K, H+ C  v3 d\\" V+ B  `
36. s;
    
37. % f  i\\" H: `( |8 b% U0 @7 S& J
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">, L4 I4 s8 n# G- H: B: D) @; O& `</font><br /> mvar:<font class=\\"jammer\\">7 n+ y\\" d1 y+ T+ `</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">& Y- U+ ?! |8 b\\" h  R% N! |</span> a = linspace(1,1,1000),<font class=\\"jammer\\"># ]; D. I. T. i+ P9 k8 a</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">9 \& x) {' [$ v( g9 h/ l- y# F</span> k=0, while{k<10000,<font class=\\"jammer\\">8 \3 a6 S: G# E% [8 r: x</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">, s2 N' }2 D% t1 ?: f</span>   k++<font class=\\"jammer\\">/ j/ s8 Y2 T# O# e/ j</font><br /> },<font class=\\"jammer\\">; p4 P  N- B! F5 F7 k</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">, P; I; d  \5 d</font><br /> <font class=\\"jammer\\">, l1 y: w: d. I1 P; T</font><br /> main1(::a)=<font class=\\"jammer\\">) l+ q. q. f5 {# c* T</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">9 Y+ Y% [* v. R. F</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">3 b  v. }\\" H! d</span>   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">! T' {; L. ]8 l, [' X</span>   k++<br /> <span style=\\"display:none\\">+ w7 p+ r' b) w6 K/ o</span> },<br /> <span style=\\"display:none\\">6 p# K, H+ C  v3 d\\" V+ B  `</span> s;<br /> <span style=\\"display:none\\">% f  i\\" H: `( |8 b% U0 @7 S& J</span> [clock()-t0]/1000;

结果：
5 r, O- |! u" b) g0 @5 \15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒
( l  p' M0 O: k$ a
( ?9 l9 R3 M/ B5 `+ l& W－－－－－－－－－－

6 [6 c5 M+ T# d; u用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];8 P8 n3 T0 l0 `/ T- j+ ?  j! o, w
   
4. mvar:5 U7 h! Z+ A3 L' U# z
   
6. main0(::a)=# e% \$ ^' o9 c) W\\" {  T; @
   
8. a = linspace(1,1,1000),
   
9. 7 Y- e4 j8 J/ t0 U5 d2 M7 f: u
10. t0=clock(),- M; g' @; b7 F) U# ]
    
12. k=0, while{k<10000,\\" C; C+ Q! T6 n' W\\" R
    
14.   i=0, while{i<1000, A(a,i)=2, i++},$ N: q  Y* F. M+ e8 U4 ^7 G
    
16.   k++1 G  r$ q9 ~9 y$ _: z5 N
    
18. },, t0 ]% _) m- t2 [- o) A
    
20. [clock()-t0]/1000;2 H% K# L- S5 f( q& h
    
22. 4 d% R! P3 T' D5 Q4 c( _
    
24. main1(::a)=
    
25. , x\\" ^\\" H3 f* B\\" j% O
26. t0=clock(), s=0,
    
27. 3 c! C1 A! X\\" _3 E8 L
28. k=0, while{k<10000,1 w- g2 x+ \* y0 o8 C
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. 0 n1 V- k# Z0 \: i
32.   k++
    
33. + e3 r5 P  G: i! x' Z
34. },
    
35. 8 O. V7 z; }# s+ {9 |( b8 E
36. s;
    
37. ; [- T1 E1 `$ n6 T6 g0 |+ j
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">8 P8 n3 T0 l0 `/ T- j+ ?  j! o, w</font><br /> mvar:<font class=\\"jammer\\">5 U7 h! Z+ A3 L' U# z</font><br /> main0(::a)=<font class=\\"jammer\\"># e% \$ ^' o9 c) W\\" {  T; @</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">7 Y- e4 j8 J/ t0 U5 d2 M7 f: u</span> t0=clock(),<font class=\\"jammer\\">- M; g' @; b7 F) U# ]</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">\\" C; C+ Q! T6 n' W\\" R</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">$ N: q  Y* F. M+ e8 U4 ^7 G</font><br />   k++<font class=\\"jammer\\">1 G  r$ q9 ~9 y$ _: z5 N</font><br /> },<font class=\\"jammer\\">, t0 ]% _) m- t2 [- o) A</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">2 H% K# L- S5 f( q& h</font><br /> <font class=\\"jammer\\">4 d% R! P3 T' D5 Q4 c( _</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">, x\\" ^\\" H3 f* B\\" j% O</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">3 c! C1 A! X\\" _3 E8 L</span> k=0, while{k<10000,<font class=\\"jammer\\">1 w- g2 x+ \* y0 o8 C</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">0 n1 V- k# Z0 \: i</span>   k++<br /> <span style=\\"display:none\\">+ e3 r5 P  G: i! x' Z</span> },<br /> <span style=\\"display:none\\">8 O. V7 z; }# s+ {9 |( b8 E</span> s;<br /> <span style=\\"display:none\\">; [- T1 E1 `$ n6 T6 g0 |+ j</span> [clock()-t0]/1000;

结果：
! N! h6 O5 V. j" Y  S4 G5 m10.469
20000000.
3.781
3 j) ?) |/ H# C6 v5 M. N$ C1 m8 k
* V3 K' }5 n8 s3 I7 c3 {" I－－－－－－－

& Y! Z) I7 E" g. [# n' D在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
+ v3 G6 h$ j* f4 C7 c! l: g; f, Z9 r在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

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

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

5、在循环中加点东西
0 p. n9 D8 w' b! T  u
matlab代码：

1. >> clear all
   ' r. N; J\" v% O0 _
2. f=@(x)x+1;
   \" U9 G4 U. S* W
3. a = linspace(2,2,1000);+ y% P# H\" h& O
   
4. tic
   * @9 e4 q\" ^7 T
5. s=0;
   / h, x! T9 w7 j: G) R, a; i
6. for k=1:10004 c* z  J$ Z+ Y& j: }( I0 D. m* D
   
7.   for i=1:1000# \2 u6 Y; e' 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));- c' O2 \0 A$ \# t1 e( W& p1 g
   
9.   end
   , |  v4 }, F\" s% j2 ]& ~. A
10. end
    7 n; D* [2 m, _# j- K+ U( ^$ z
11. s5 }! n% M2 R: |3 Z9 i' Q
    
12. toc
    2 ?. [' v1 u7 k* L
13. 
    ( m% v3 P0 ^  g) s; a3 o+ T
14. s =4 K  i; e0 [  z0 x; x
    
15. $ I3 J* N1 H, l% [7 y- b% K6 |; T
    
16.   2.4002e+007% z- F) _( c+ ^! R/ u
    
17. 
    1 J\" q- U0 l& s1 }. @0 T6 g7 |
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. 9 @, l% Z$ f# a' x/ a
4. mvar:% A3 d' O7 I% [$ z7 z7 `
   
6. f(x)=x+1;
   
7. , F' G$ p. R# j
8. main0(::a)=+ `4 o  Y/ J5 u\\" Z8 V& u
   
10. a = linspace(2,2,1000),$ L! t  p  I  ~, O5 q$ A+ y& c) j7 D
    
12. t0=clock(), s=0,
    
13. ' W2 ]) M4 Y8 W2 u. n
14. k=0, while{k<1000,
    
15. ' Q0 R8 O$ h2 u  G7 e5 w6 Z0 I$ H
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++},& j5 L& _7 e9 a! B- c- n* g
    
18.   k++
    
19. ( b/ v5 D, h4 }- l
20. },
    
21. ! i- f* j4 U- B- A' j- x
22. s;8 `8 `# I- }4 C% ~) @4 S  n* \$ g3 A
    
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">9 @, l% Z$ f# a' x/ a</span> mvar:<font class=\\"jammer\\">% A3 d' O7 I% [$ z7 z7 `</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">, F' G$ p. R# j</span> main0(::a)=<font class=\\"jammer\\">+ `4 o  Y/ J5 u\\" Z8 V& u</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">$ L! t  p  I  ~, O5 q$ A+ y& c) j7 D</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">' W2 ]) M4 Y8 W2 u. n</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">' Q0 R8 O$ h2 u  G7 e5 w6 Z0 I$ H</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\\">& j5 L& _7 e9 a! B- c- n* g</font><br />   k++<br /> <span style=\\"display:none\\">( b/ v5 D, h4 }- l</span> },<br /> <span style=\\"display:none\\">! i- f* j4 U- B- A' j- x</span> s;<font class=\\"jammer\\">8 `8 `# I- }4 C% ~) @4 S  n* \$ g3 A</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
( m+ @8 a2 N* V$ p) j+ d5.172
4 r; d5 d( u4 x( y6 d+ q: }
二者效率现在只有一倍差距。

, C) m1 k5 G' U6 ?, w/ P－－－－－－－－－－

不过，显然matlab和Forcal的程序应该这样写：

$ A8 ]9 |5 L# P  U; `4 h) amatlab程序：

1. clear all* q: q8 U0 B\" N
   
2. f=@(x)x+1;
   5 {9 H8 t0 V0 e* i+ J
3. a = linspace(2,2,1000);& r# y0 G. h3 I; S; X5 |+ G
   
4. tic0 J# j3 ^\" M+ I' q2 {& u4 \
   
5. s=0;
   \" `& k8 ]/ V$ s
6. for k=1:1000, ^. X\" z9 ?+ Q) X2 O$ B8 r
   
7.   for i=1:10001 E6 a4 J8 l2 o/ O\" o5 r
   
8.     t=a(i);! y: |9 W5 g: x8 m5 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);
   2 }$ R- b' ~$ ~+ w( n. F
10.   end
    $ \2 [' V7 C' }  s2 y5 [
11. end
    ' o7 G2 H3 r1 H/ j5 J2 V% @
12. s
    \" ]1 c/ O6 ?- g\" l. w\" m& ^$ ?
13. toc' H) x) D+ c\" C
    
14. 1 R; @7 S6 S\" `! A3 m3 o9 K0 R/ H
    
15. s =
    ' H2 d/ C) D* l8 x
16. ) y+ o, V+ K* k* P% f! d\" y! o% k
    
17.   2.4002e+007
    : h\" [1 q\" x; {* ^5 K
18. . D! i8 g: g1 L7 X, \' c
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];. }7 i$ B$ ?! f  U. |1 W- ?
   
4. mvar:
   
5. 0 ^  _  m1 S1 I- p# J8 T5 c
6. f(x)=x+1;
   
7. , U; n( l& Z5 n' b6 k  s: J! I6 D
8. main0(::a)=
   
9. 9 D/ I6 p. s) e6 L! W
10. a = linspace(2,2,1000),2 Z6 P/ k* T( l- \) t0 b
    
12. t0=clock(), s=0,
    
13. & s5 a- d! U, p8 p$ f) A! B0 D
14. k=0, while{k<1000,7 _) {( z' z& l# q5 |2 L
    
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++},* l' j$ Y$ Y3 }& t) i
    
18.   k++
    
19. ) h: G7 q& j1 W2 N. |+ E* D
20. },+ Z7 m' |0 A\\" M5 M\\" G( w. o, I1 G: U
    
22. s;
    
23. / ], J  Y! u) L, }( d\\" U3 U
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">. }7 i$ B$ ?! f  U. |1 W- ?</font><br /> mvar:<br /> <span style=\\"display:none\\">0 ^  _  m1 S1 I- p# J8 T5 c</span> f(x)=x+1;<br /> <span style=\\"display:none\\">, U; n( l& Z5 n' b6 k  s: J! I6 D</span> main0(::a)=<br /> <span style=\\"display:none\\">9 D/ I6 p. s) e6 L! W</span> a = linspace(2,2,1000),<font class=\\"jammer\\">2 Z6 P/ k* T( l- \) t0 b</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">& s5 a- d! U, p8 p$ f) A! B0 D</span> k=0, while{k<1000,<font class=\\"jammer\\">7 _) {( z' z& l# q5 |2 L</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++},<font class=\\"jammer\\">* l' j$ Y$ Y3 }& t) i</font><br />   k++<br /> <span style=\\"display:none\\">) h: G7 q& j1 W2 N. |+ E* D</span> },<font class=\\"jammer\\">+ Z7 m' |0 A\\" M5 M\\" G( w. o, I1 G: U</font><br /> s;<br /> <span style=\\"display:none\\">/ ], J  Y! u) L, }( d\\" U3 U</span> [clock()-t0]/1000;

结果：
24002271.69189384
9 L- C; k1 h8 {7 ?- h1.625
8 ^7 i6 b, X4 D4 Z+ `% r$ `
7 J. s& p" ^: T" nmatlab效率已不及Forcal。

( n9 w$ N0 }  o8 @: e; k提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
$ O+ R- A. @5 K* u9 O# n( _
9 l5 j* n! l  G5 _- cForcal代码：

2. !using["sys"];
   
3. 0 Z2 g2 @\\" {& Y# l/ D& R  E& V
4. mvar:7 R7 y# [  a7 V3 H. Q\\" c. d
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. 8 V& ^5 ^0 U# [' G# z
8. t0=clock(),! u5 j& L9 g7 f6 u3 B5 P. A
   
10. k=0, while{k<10000,3 V5 |. \4 o$ H; ~; X+ O! ?
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. 4 B* f6 k* v, ]5 [) \: O
14.   k++7 S7 L) \6 a' T) X
    
16. },
    
17.   Y\\" L7 c' r! x# s# W
18. [clock()-t0]/1000;: W3 Z, I+ V! T! N0 {, }/ Y7 p
    
20. 
    
21. 4 D( a3 S7 R8 q' ~$ |& t6 K
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. , E% [, z  S( p& d- s2 E/ D6 l
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. % e  S9 m5 a+ _5 U\\" |; _
26. t0=clock(),0 n7 }, z( f: V6 B4 M& Y
    
28. s=0,
    
29. 2 `  R6 N$ p* q
30. k=0, while{k<10000,
    
31. + J  o5 K; C% T7 V5 q& t4 p7 l
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. ' m' y/ w. t9 t. ^4 N
34.   k++5 [) f. S2 D, J# I' n! q, @! f% k
    
36. },: K\\" g( z  E7 ?
    
38. s;
    
39. 0 K# n\\" u! B0 B. B
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">0 Z2 g2 @\\" {& Y# l/ D& R  E& V</span> mvar:<font class=\\"jammer\\">7 R7 y# [  a7 V3 H. Q\\" c. d</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">8 V& ^5 ^0 U# [' G# z</span> t0=clock(),<font class=\\"jammer\\">! u5 j& L9 g7 f6 u3 B5 P. A</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">3 V5 |. \4 o$ H; ~; X+ O! ?</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">4 B* f6 k* v, ]5 [) \: O</span>   k++<font class=\\"jammer\\">7 S7 L) \6 a' T) X</font><br /> },<br /> <span style=\\"display:none\\">  Y\\" L7 c' r! x# s# W</span> [clock()-t0]/1000;<font class=\\"jammer\\">: W3 Z, I+ V! T! N0 {, }/ Y7 p</font><br /> <br /> <span style=\\"display:none\\">4 D( a3 S7 R8 q' ~$ |& t6 K</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">, E% [, z  S( p& d- s2 E/ D6 l</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">% e  S9 m5 a+ _5 U\\" |; _</span> t0=clock(),<font class=\\"jammer\\">0 n7 }, z( f: V6 B4 M& Y</font><br /> s=0,<br /> <span style=\\"display:none\\">2 `  R6 N$ p* q</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">+ J  o5 K; C% T7 V5 q& t4 p7 l</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">' m' y/ w. t9 t. ^4 N</span>   k++<font class=\\"jammer\\">5 [) f. S2 D, J# I' n! q, @! f% k</font><br /> },<font class=\\"jammer\\">: K\\" g( z  E7 ?</font><br /> s;<br /> <span style=\\"display:none\\">0 K# n\\" u! B0 B. B</span> [clock()-t0]/1000;

结果：
5 _8 Y# d, a5 ^/ I% B0.5620000000000001  //赋值用时间，秒
# X" E* y* f8 o3 {. b6 a- E20000000.
: E+ w+ [$ X6 }' _0.579      //取值用时间，秒