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

forcal

1、简单循环
+ A( B+ e3 a* _0 N5 P# {
, p7 W' R* H* X! G7 nMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   & k$ E% Y2 U1 r- P0 U$ {
2. tic
   : f* b\" ^' S  n2 H# z- d. x
3. s=0;& n. H2 h/ o- X
   
4. for k=1:10000
   : _) a% j1 [& R( s. k- W7 c4 ]
5.   for i=1:1000
   ) w, A0 ~$ ^\" m. W& E- U\" n
6.     s=s+1;
   # _5 U, L& P$ T- R& _
7.   end) ~& t; r: n$ o# k7 [/ H# w
   
8. end
   5 \8 o8 H5 V) G+ `, }\" e
9. s; W/ G* C5 u6 x6 r
   
10. toc4 z' d5 e2 j3 y' Z+ H; R6 f
    
11. 
    . G\" ~3 f  V\" |+ B) u' G\" i\" ]8 }
12. s =* c( T- a$ u; P4 F
    
13. 5 S  a& }* J4 z% r$ N2 ?! X$ _4 h
    
14.     10000000. D7 q4 ~( R! x) [; g8 y
    
15. # I; E3 p) ]& h0 W
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. + h4 P( w1 X- @4 S7 T
4. mvar:
   
5. 1 O4 S7 Q- j\\" f! B6 W
6. t0=clock(),5 z1 u% d8 a& s! |5 }; e
   
8. s=0,
   
9. 3 z2 ~/ {& E: @
10. k=0, while{k<10000,% E% v8 k1 X\\" G5 J3 W
    
12.   i=0, while{i<1000,s=s+1, i++},
    
13. * w0 O& i9 s/ s6 \8 ~1 }
14.   k++' d( I+ P0 I7 t4 p; n: I
    
16. },$ @4 Y$ v; T* s9 c$ l* f2 b
    
18. s;6 c  N( r9 a+ `# Y; l/ L. ?4 \
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">+ h4 P( w1 X- @4 S7 T</span> mvar:<br /> <span style=\\"display:none\\">1 O4 S7 Q- j\\" f! B6 W</span> t0=clock(),<font class=\\"jammer\\">5 z1 u% d8 a& s! |5 }; e</font><br /> s=0,<br /> <span style=\\"display:none\\">3 z2 ~/ {& E: @</span> k=0, while{k<10000,<font class=\\"jammer\\">% E% v8 k1 X\\" G5 J3 W</font><br />   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">* w0 O& i9 s/ s6 \8 ~1 }</span>   k++<font class=\\"jammer\\">' d( I+ P0 I7 t4 p; n: I</font><br /> },<font class=\\"jammer\\">$ @4 Y$ v; T* s9 c$ l* f2 b</font><br /> s;<font class=\\"jammer\\">6 c  N( r9 a+ `# Y; l/ L. ?4 \</font><br /> [clock()-t0]/1000;

结果：
2 I& m' x' W  Q8 ?* g9 q- }10000000.
8 Z; c- A6 c0 E8 a0.704  秒
& H7 ?& h) u4 R. _
: ]! _# ~7 z% N) c二者运行效率相差7倍多。

* T9 S7 I; z. c" B' k/ X2、循环中稍微加点东西

Matlab代码：

1. clear all
   * d5 L4 r$ y) ^; v
2. tic+ o) R7 k. [  B$ R8 M* d
   
3. s=0;$ f. z\" [5 x  P( t
   
4. for k=1:1000  v' C- ~\" v! q. q8 v6 }
   
5.   for i=1:1000
   / e7 Q* b! _6 n, q$ F/ {! n
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);. I2 t# \$ @! S4 w# C/ R
   
7.   end
     X, i. G5 `# ^! l6 ^7 k
8. end
   $ A& K& ^( s& ?4 @' E: ^
9. s* u, s  P4 G/ G7 q
   
10. toc8 B& G, Q/ ?- B) n6 {
    
11. 
    ! I0 k' H$ ^8 o
12. s =& `2 a! \3 m& i. W/ z: z' y\" u
    
13. 
    6 Z2 k$ Z% |1 {3 m% y& d) ~* z4 r
14.   1.1009e+007
    1 [; t; P- ]9 ?8 @$ u
15. & ^& J. H7 B) a  B. q  q4 D
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. + \) e; _! Q$ H! k& {( X/ W
4. mvar:
   
5. 3 S* i! w+ n; U\\" x, T# ?# i
6. t0=clock(),
   
7. : i1 L. x$ R1 ^0 J  e
8. s=0,5 K; i! v0 g! R+ _9 u: D$ T: ^
   
10. k=1, while{k<=1000,9 d& U5 b; R# Y- x/ Z/ [
    
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++},* ?: _4 C/ L. Q0 X  _1 v
    
14.   k++\\" ^# y; g! H! F
    
16. },) M+ b! M) e\\" A. M# s: q, Z
    
18. s;
    
19. & v; o  ?, z: q9 r  D
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">+ \) e; _! Q$ H! k& {( X/ W</span> mvar:<br /> <span style=\\"display:none\\">3 S* i! w+ n; U\\" x, T# ?# i</span> t0=clock(),<br /> <span style=\\"display:none\\">: i1 L. x$ R1 ^0 J  e</span> s=0,<font class=\\"jammer\\">5 K; i! v0 g! R+ _9 u: D$ T: ^</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">9 d& U5 b; R# Y- x/ Z/ [</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\\">* ?: _4 C/ L. Q0 X  _1 v</font><br />   k++<font class=\\"jammer\\">\\" ^# y; g! H! F</font><br /> },<font class=\\"jammer\\">) M+ b! M) e\\" A. M# s: q, Z</font><br /> s;<br /> <span style=\\"display:none\\">& v; o  ?, z: q9 r  D</span> [clock()-t0]/1000;

结果：
11008953.66602346
1.031
% R& `1 O4 m( T3 t7 v
% w$ g  ]0 y3 H# e6 U1 `二者效率相差已比较小。
' ]. Y7 d! v+ f9 J1 Q, d
( f, Z3 R3 q9 l2 R3、再增加一个函数调用
$ \/ y! P7 B2 w
Matlab代码：

1. clear all
     p# r* H- `2 ?& @
2. f=@(x)x+1;2 D5 u  {' _) s$ r6 I
   
3. tic  N* f( W0 {! \2 A3 d) S
   
4. s=0;
   ! j3 [3 |. t. U6 J1 U
5. for k=1:1000: \! r: L7 u9 U0 ]4 S
   
6.   for i=1:1000
   ( D' N, i; O, G3 W
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 o9 b/ l6 A/ E2 J  z
8.   end3 t\" B! o& Z5 k! X* ?9 y+ m
   
9. end
   2 F( [( ?\" ^# B- B6 u$ u
10. s$ r4 w0 O# B4 Y& ?4 o% U
    
11. toc
    3 x; n! P: x9 K
12. & K  h8 ?  Y0 ~4 ?
    
13. s =
    # A  |0 s$ p2 Y
14. 7 p% G* P8 A0 U$ G' k% V8 c' {3 O
    
15.   5.1251e+008' y; p6 h' _& f% M3 v, s
    
16. . j1 r+ M0 Y& O9 Q8 f* d\" D0 [; `$ b, a
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. * R2 Y+ z/ J, m$ }1 e) a
4. f(x)=x+1;
   
5. : }\\" x0 R) F. p\\" P( o
6. mvar:+ T' d; v! {1 _3 \1 Y
   
8. t0=clock(),+ j  U9 h% x# t9 H+ @) q4 [9 F+ R
   
10. s=0,\\" a) I: Y: y% H/ Z/ W1 D) d- C
    
12. k=1, while{k<=1000,  e* D. e. A0 f2 Z. ~/ b; z  e1 E. L
    
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. 0 q% d\\" a* p  B
16.   k++! O\\" q; t6 q7 S! |5 _$ W
    
18. },# S8 ~0 @6 l# _7 Y- T% _4 z+ u
    
20. s;
    
21. 4 G- u: ^& X* W/ u- S
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">* R2 Y+ z/ J, m$ }1 e) a</span> f(x)=x+1;<br /> <span style=\\"display:none\\">: }\\" x0 R) F. p\\" P( o</span> mvar:<font class=\\"jammer\\">+ T' d; v! {1 _3 \1 Y</font><br /> t0=clock(),<font class=\\"jammer\\">+ j  U9 h% x# t9 H+ @) q4 [9 F+ R</font><br /> s=0,<font class=\\"jammer\\">\\" a) I: Y: y% H/ Z/ W1 D) d- C</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">  e* D. e. A0 f2 Z. ~/ b; z  e1 E. L</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\\">0 q% d\\" a* p  B</span>   k++<font class=\\"jammer\\">! O\\" q; t6 q7 S! |5 _$ W</font><br /> },<font class=\\"jammer\\"># S8 ~0 @6 l# _7 Y- T% _4 z+ u</font><br /> s;<br /> <span style=\\"display:none\\">4 G- u: ^& X* W/ u- S</span> [clock()-t0]/1000;

结果：
8 i$ x  g7 d0 Q( G+ r, v/ W  a512508953.6658993
# E5 y. F' [. ]* [% S1.078

0 N: c. u" b  M3 u4 y5 P可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

forcal

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

Forcal代码：

2. !using["sys"];1 h7 X  P5 B$ _5 z& h
   
4. mvar:- \& \) v% ]- v% Q+ g; X
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组. c8 v# O; R! e. h# s8 }
   
8. t0=clock(),
   
9. 4 O2 S  H, G5 K$ J
10. k=0, while{k<10000,
    
11. . c8 i) X) n1 C) ^2 x8 ~3 v# }
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值8 k4 U. ~$ a2 x3 L3 o
    
14.   k++8 d% E1 P/ A. L3 v
    
16. },
    
17. ( i$ r+ H* u\\" _1 E; l, s
18. [clock()-t0]/1000;
    
19. . M( K' Z& V. \8 i/ ]\\" L
20. + B3 y4 @4 j\\" s6 p) U1 l- }  g
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组' ~+ H6 c) H1 c1 Z- t; D4 K
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值5 z9 {1 x; A, f. ^7 S
    
26. t0=clock(),. Z4 [) c# P/ e$ z% H, h0 W\\" N
    
28. s=0,' R' [7 V0 R# I
    
30. k=0, while{k<10000,
    
31. \\" k0 c9 s6 d. E3 j
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值. v$ W\\" O5 j  d: o2 V
    
34.   k++4 S% K. [; T7 O# m5 ^6 E
    
36. },% w/ K- ^! H7 a+ _
    
38. s;
    
39. 7 V, l! i6 K) @) k
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">1 h7 X  P5 B$ _5 z& h</font><br /> mvar:<font class=\\"jammer\\">- \& \) v% ]- v% Q+ g; X</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">. c8 v# O; R! e. h# s8 }</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">4 O2 S  H, G5 K$ J</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">. c8 i) X) n1 C) ^2 x8 ~3 v# }</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">8 k4 U. ~$ a2 x3 L3 o</font><br />   k++<font class=\\"jammer\\">8 d% E1 P/ A. L3 v</font><br /> },<br /> <span style=\\"display:none\\">( i$ r+ H* u\\" _1 E; l, s</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">. M( K' Z& V. \8 i/ ]\\" L</span> <font class=\\"jammer\\">+ B3 y4 @4 j\\" s6 p) U1 l- }  g</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">' ~+ H6 c) H1 c1 Z- t; D4 K</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">5 z9 {1 x; A, f. ^7 S</font><br /> t0=clock(),<font class=\\"jammer\\">. Z4 [) c# P/ e$ z% H, h0 W\\" N</font><br /> s=0,<font class=\\"jammer\\">' R' [7 V0 R# I</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">\\" k0 c9 s6 d. E3 j</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">. v$ W\\" O5 j  d: o2 V</font><br />   k++<font class=\\"jammer\\">4 S% K. [; T7 O# m5 ^6 E</font><br /> },<font class=\\"jammer\\">% w/ K- ^! H7 a+ _</font><br /> s;<br /> <span style=\\"display:none\\">7 V, l! i6 K) @) k</span> [clock()-t0]/1000;

结果：
% D  I" J% r1 l$ a$ x0.5620000000000001  //赋值用时间，秒
20000000.
7 S/ M8 h; g+ K7 A9 i0.579      //取值用时间，秒

forcal

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

# g9 \6 Y% m5 C0 ZMatlab代码：

1. >> clear all
   7 _$ j& w0 Y* q$ h. C  d: m
2. a = linspace(1,1,1000);0 \. V0 B& w0 ]: I9 i7 g
   
3. tic
   . {\" o# a2 E& {5 ?* X- a: \
4. for k=1:10000/ d0 p( C: y& \8 b9 D. `
   
5.   for i=1:1000
   % a8 b# r3 e; T' x' S$ {+ q
6.     a(i)=2;$ J3 z3 y! ]2 l( _3 o# E% I
   
7.   end7 @  H% j/ A2 h3 q
   
8. end\" ~\" v7 G( }4 I( X7 I1 j
   
9. toc! n( {( Y) Q# i  ?/ }/ V
   
10. 
    ' M) f8 r9 Z' x. p: S
11. tic
      s$ K4 C- u( E, ?- q4 h
12. s=0;
    8 x, E# w$ P9 z) t) o6 N/ p4 d* Q/ e
13. for k=1:10000: u* l2 Z9 \7 u% F; K
    
14.   for i=1:1000  M1 o( ]4 [- P/ I- {! s8 ?
    
15.     s=s+a(i);2 F* G6 a\" T% v' H# r2 U4 N
    
16.   end
    \" f: J2 I# V/ A9 M! [2 X1 M
17. end
    & e# J! M# A( o$ j
18. s
    5 q4 Q3 `: R* Q, g0 O
19. toc* e7 }1 {9 ]# ?% C% a9 y  n
    
20. Elapsed time is 0.140658 seconds.
    4 z1 H5 Z; a. S& k2 Q3 f
21. 
    / T/ e8 @' S% N( {& G- g& ]8 `
22. s =
    : C' p% ~4 `8 Y* Q
23. 
    6 @' m/ Q0 J3 K\" w! S  g$ ?
24.     20000000
    % k% c0 m2 P7 |0 n' Y; k: R! I
25. 5 j& z! a. ]- n8 ?- V  \
    
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
& [% j. r$ z) E0 \. A
Forcal代码：

2. !using["math","sys"];
   
3. 8 i& c+ H2 m- M
4. mvar:2 T2 |# E% R* B* O: q
   
6. main0(::a)=
   
7. ' Y4 A7 e( [, o& }. `7 z, m3 U0 ]
8. a = linspace(1,1,1000),4 y' e/ i' d( t, O* p: b3 ]8 k6 D3 o
   
10. t0=clock(),) w( E' s: V, L7 A5 H
    
12. k=0, while{k<10000,4 v0 O* M5 r4 E2 i, G' y
    
14.   i=0, while{i<1000, a(i)=2, i++},7 d6 s- ], n% S5 s! T, h9 `, b
    
16.   k++3 H' e, Y2 Y4 M\\" g
    
18. },
    
19. % v- N; h5 A  {  p1 N3 K
20. [clock()-t0]/1000;
    
21. 8 _& z' M/ E$ b' |7 r
22. 
    
23. ) u- L* s- @+ J$ n$ s( l
24. main1(::a)=
    
25. + S# L6 i' [. T; z( K1 {& z& I- q
26. t0=clock(), s=0,! ?7 @' r1 ~' Z! ^2 n( ]
    
28. k=0, while{k<10000,/ b$ V8 E4 }% a0 E
    
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. 2 B  B% ~! d7 F) ~* D  p0 _$ b
32.   k++
    
33. ' n* u) e+ D0 |- Y6 C) ^
34. },% z' L4 e; T( L1 O# R/ z5 p6 V
    
36. s;
    
37. ! I. Z4 r( W* O8 s7 u! F
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">8 i& c+ H2 m- M</span> mvar:<font class=\\"jammer\\">2 T2 |# E% R* B* O: q</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">' Y4 A7 e( [, o& }. `7 z, m3 U0 ]</span> a = linspace(1,1,1000),<font class=\\"jammer\\">4 y' e/ i' d( t, O* p: b3 ]8 k6 D3 o</font><br /> t0=clock(),<font class=\\"jammer\\">) w( E' s: V, L7 A5 H</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">4 v0 O* M5 r4 E2 i, G' y</font><br />   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">7 d6 s- ], n% S5 s! T, h9 `, b</font><br />   k++<font class=\\"jammer\\">3 H' e, Y2 Y4 M\\" g</font><br /> },<br /> <span style=\\"display:none\\">% v- N; h5 A  {  p1 N3 K</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">8 _& z' M/ E$ b' |7 r</span> <br /> <span style=\\"display:none\\">) u- L* s- @+ J$ n$ s( l</span> main1(::a)=<br /> <span style=\\"display:none\\">+ S# L6 i' [. T; z( K1 {& z& I- q</span> t0=clock(), s=0,<font class=\\"jammer\\">! ?7 @' r1 ~' Z! ^2 n( ]</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">/ b$ V8 E4 }% a0 E</font><br />   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">2 B  B% ~! d7 F) ~* D  p0 _$ b</span>   k++<br /> <span style=\\"display:none\\">' n* u) e+ D0 |- Y6 C) ^</span> },<font class=\\"jammer\\">% z' L4 e; T( L1 O# R/ z5 p6 V</font><br /> s;<br /> <span style=\\"display:none\\">! I. Z4 r( W* O8 s7 u! F</span> [clock()-t0]/1000;

结果：
/ X( M+ h' f! {" C. R15.484  //存数据所用时间，秒
8 W, p9 O5 P" b* [20000000.
5.531   //取数据所用时间，秒

5 `1 }' n% Q4 Q－－－－－－－－－－
. a- r( G  C' P1 i; u. g
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. & Z/ {: a- t+ m- e\\" S/ @5 z( }3 A: _
4. mvar:/ U& B6 R; P& T* @
   
6. main0(::a)=, x/ Z7 z# C( P
   
8. a = linspace(1,1,1000),
   
9. 9 n- _. J0 ]2 c. z) w
10. t0=clock(),
    
11. 9 n% {. a; u+ H/ m& D
12. k=0, while{k<10000,; [9 m% l7 K' @) _\\" _
    
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. ! [; f( y. R. ^# a$ y1 |+ c
16.   k++' `3 b6 j! b, v! f1 k( [9 A
    
18. },! O% J7 w1 u1 o& {! f0 H! [; c
    
20. [clock()-t0]/1000;
    
21. # r! Z* O! T2 H5 N\\" N: R8 w& p
22. 
    
23. $ I: y4 b0 N\\" h\\" J+ h
24. main1(::a)=8 b, U) ~( e, Y5 E  v/ f
    
26. t0=clock(), s=0,
    
27. 3 r% U3 K9 T5 h
28. k=0, while{k<10000,! ^& m! j) p: p& `6 |% M# x, a' w
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. \\" ?$ C: j$ e; o: [9 p
32.   k++
    
33. * W* v\\" P( a5 G% `4 G2 o
34. },
    
35.   r6 T! Y7 U# y  c4 z
36. s;' n! D( K, ?+ b% z$ l' T
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">& Z/ {: a- t+ m- e\\" S/ @5 z( }3 A: _</span> mvar:<font class=\\"jammer\\">/ U& B6 R; P& T* @</font><br /> main0(::a)=<font class=\\"jammer\\">, x/ Z7 z# C( P</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">9 n- _. J0 ]2 c. z) w</span> t0=clock(),<br /> <span style=\\"display:none\\">9 n% {. a; u+ H/ m& D</span> k=0, while{k<10000,<font class=\\"jammer\\">; [9 m% l7 K' @) _\\" _</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">! [; f( y. R. ^# a$ y1 |+ c</span>   k++<font class=\\"jammer\\">' `3 b6 j! b, v! f1 k( [9 A</font><br /> },<font class=\\"jammer\\">! O% J7 w1 u1 o& {! f0 H! [; c</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\"># r! Z* O! T2 H5 N\\" N: R8 w& p</span> <br /> <span style=\\"display:none\\">$ I: y4 b0 N\\" h\\" J+ h</span> main1(::a)=<font class=\\"jammer\\">8 b, U) ~( e, Y5 E  v/ f</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">3 r% U3 K9 T5 h</span> k=0, while{k<10000,<font class=\\"jammer\\">! ^& m! j) p: p& `6 |% M# x, a' w</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">\\" ?$ C: j$ e; o: [9 p</span>   k++<br /> <span style=\\"display:none\\">* W* v\\" P( a5 G% `4 G2 o</span> },<br /> <span style=\\"display:none\\">  r6 T! Y7 U# y  c4 z</span> s;<font class=\\"jammer\\">' n! D( K, ?+ b% z$ l' T</font><br /> [clock()-t0]/1000;

结果：
10.469
20000000.
3.781

. Y6 }6 X9 i! c1 @2 a4 X6 m7 c4 G－－－－－－－

. S1 l5 L' z  z- r在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
# U" a* k) z, W  q$ i& ]0 p
* U% ]1 w/ F! w$ hForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
, C8 Q4 z8 E2 M& B. I( y
/ S/ O2 Q' a( g" e* g7 M! |5、在循环中加点东西

2 G3 Z  y6 v" K5 {$ V6 q7 Vmatlab代码：

1. >> clear all
   5 \3 B8 I& f3 S, m
2. f=@(x)x+1;8 G; y+ {! s4 P0 r! z\" O9 \
   
3. a = linspace(2,2,1000);' [. ^\" v. ^9 s; n! _+ x
   
4. tic
   & k5 H  s3 J\" x; [* J  `3 O
5. s=0;# _9 L3 O1 H% k/ {. c5 q
   
6. for k=1:1000
   ; h# T* W% D* N! ~- Y
7.   for i=1:1000
   . ]* k, `6 G5 d* D% m1 C
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\" C) D9 E/ B2 s6 {
9.   end
   9 v6 I% I% m9 T# @. W
10. end
    8 N5 |$ f4 {+ H: |$ u& q
11. s5 o* k8 h8 \! K( B; k' r& ]
    
12. toc6 z+ @3 ?( ^6 V8 I  ?9 u& A
    
13. 4 s1 Y: t9 @/ P: Z( @
    
14. s =
    , t' }4 M( ?; Q, Y8 Y2 y3 z
15. 
    3 T3 O& E, u0 R  ^7 `
16.   2.4002e+0077 p0 e9 }1 d4 @7 N8 l. [
    
17. 9 G6 S9 H0 A\" U; ^' D9 u
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];+ [0 j/ f* G2 T( Q6 i3 o% M) Y9 \
   
4. mvar:
   
5. 6 R% }6 s7 T/ k+ [/ [% R* t
6. f(x)=x+1;! g0 ]2 I' w, m1 k
   
8. main0(::a)=# @' `\\" t# C( T4 g* {
   
10. a = linspace(2,2,1000),- m( W! S( P) Q- |, w
    
12. t0=clock(), s=0,' L2 i0 z' b7 H
    
14. k=0, while{k<1000,
    
15. # S( a7 g5 o3 k. }6 E  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++},9 l8 ?) R: f- H8 T1 w7 T+ V( t$ T
    
18.   k++8 b5 ^- }1 \\\" V# k2 j4 s! N; ^
    
20. },
    
21. 0 S6 M4 X1 y$ `: [
22. s;4 e% F- a  z9 z4 S  T
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ [0 j/ f* G2 T( Q6 i3 o% M) Y9 \</font><br /> mvar:<br /> <span style=\\"display:none\\">6 R% }6 s7 T/ k+ [/ [% R* t</span> f(x)=x+1;<font class=\\"jammer\\">! g0 ]2 I' w, m1 k</font><br /> main0(::a)=<font class=\\"jammer\\"># @' `\\" t# C( T4 g* {</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">- m( W! S( P) Q- |, w</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">' L2 i0 z' b7 H</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\"># S( a7 g5 o3 k. }6 E  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\\">9 l8 ?) R: f- H8 T1 w7 T+ V( t$ T</font><br />   k++<font class=\\"jammer\\">8 b5 ^- }1 \\\" V# k2 j4 s! N; ^</font><br /> },<br /> <span style=\\"display:none\\">0 S6 M4 X1 y$ `: [</span> s;<font class=\\"jammer\\">4 e% F- a  z9 z4 S  T</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
5.172

! Q. [4 r$ s% y! _二者效率现在只有一倍差距。

/ p  P. g+ z: S% L－－－－－－－－－－
2 g8 l2 C9 O9 j0 e
$ T6 K" c) ^, _不过，显然matlab和Forcal的程序应该这样写：
, @7 @, i/ T" r+ z6 H3 `6 H
* {- Z+ @  j0 Z$ Z4 D5 T5 M/ D0 {/ lmatlab程序：

1. clear all
     T' a- R& H6 O4 L$ \
2. f=@(x)x+1;. h3 J2 f4 r2 x- v
   
3. a = linspace(2,2,1000);
   8 u+ E* A& L' \. {* \1 _
4. tic
   1 h* b\" `\" }. I0 O2 g
5. s=0;
   7 I& n0 F/ h\" b. ^# L2 {5 g7 s
6. for k=1:1000
   3 w; V4 [4 g6 ^: s( F; C
7.   for i=1:1000% X, T0 j  q  a0 F
   
8.     t=a(i);& X+ p4 w/ e1 ^6 R! D/ Z1 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);2 ?/ I# c0 ^* S
   
10.   end
    % H0 ?7 z/ ^- b
11. end  a8 g7 Z+ k& j6 g
    
12. s
    1 S  K* a\" ~# r
13. toc7 R\" W8 p% y7 A% L+ X
    
14. & I, d9 [) m: A2 v2 I
    
15. s =: t' i/ q\" Y3 i+ [5 Z7 B
    
16. 3 e\" [! R% F9 A' Y5 ^: s9 q
    
17.   2.4002e+007' p9 f2 q8 v/ L' |1 [( l
    
18. 6 Z& S- F0 }( U( T& X$ d3 I\" q
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];( T% @\\" t* k5 L5 N! ]4 _3 b8 s
   
4. mvar:
   
5. % W7 i2 D. p, D0 Q6 V7 e4 L( X
6. f(x)=x+1;
   
7.   W  H' x* j4 U! C, b
8. main0(::a)=
   
9. ; D- s+ b9 ^\\" J3 S) w  k, o
10. a = linspace(2,2,1000),
    
11. & d. y5 f& O; n9 b
12. t0=clock(), s=0,
    
13. : C  D% z4 I4 M$ I. L
14. k=0, while{k<1000,2 T8 l+ A5 i! ^# b
    
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. 0 z: J' L6 k- Q, |3 N$ c\\" B9 f* P: y
18.   k++
    
19. + r\\" |! E- Z9 r
20. },
    
21. 0 M, C) s4 m- n$ O  }5 h0 C7 p/ x
22. s;7 j' }' A$ ?; c; ]# C) @9 a
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">( T% @\\" t* k5 L5 N! ]4 _3 b8 s</font><br /> mvar:<br /> <span style=\\"display:none\\">% W7 i2 D. p, D0 Q6 V7 e4 L( X</span> f(x)=x+1;<br /> <span style=\\"display:none\\">  W  H' x* j4 U! C, b</span> main0(::a)=<br /> <span style=\\"display:none\\">; D- s+ b9 ^\\" J3 S) w  k, o</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">& d. y5 f& O; n9 b</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">: C  D% z4 I4 M$ I. L</span> k=0, while{k<1000,<font class=\\"jammer\\">2 T8 l+ A5 i! ^# b</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\\">0 z: J' L6 k- Q, |3 N$ c\\" B9 f* P: y</span>   k++<br /> <span style=\\"display:none\\">+ r\\" |! E- Z9 r</span> },<br /> <span style=\\"display:none\\">0 M, C) s4 m- n$ O  }5 h0 C7 p/ x</span> s;<font class=\\"jammer\\">7 j' }' A$ ?; c; ]# C) @9 a</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
2 K3 t0 z6 D+ T0 a4 B! h$ O7 ^, h1.625

matlab效率已不及Forcal。
" k7 L' l9 J3 Y* A! E+ l. z( n
/ I1 n; @( s5 C) M提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

海水