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

forcal

1、简单循环
* ?8 Y2 s  i4 @% d3 D
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   . [+ F7 U2 h% B\" z2 b- l; s
2. tic
   : \& E& p. Z) L! {; Y
3. s=0;0 C\" e& d\" a& {* h$ p
   
4. for k=1:100000 C; {5 X/ j. Z( {! Y  H
   
5.   for i=1:1000
   5 ?0 r& ?3 N( c+ Q0 u0 u- _8 B
6.     s=s+1;/ t, r; y3 b3 ]
   
7.   end
   / n4 x- R' E/ U# _! Q) v
8. end
   9 u* T) R1 o; ~* M! g1 S, \\" B
9. s
   + d9 C! O! U  K# F5 n\" t* o
10. toc+ t# k2 c7 \' o9 b4 c* I9 B
    
11. 
    ( O7 T1 b2 V$ b  _
12. s =
    & m0 ]: l; u7 j: R\" u* M! j
13. 
    2 r9 Z- ?5 L. u, X% s3 Z
14.     100000002 {& a* V- S\" c9 K4 p
    
15. ; h. \2 K* ]! J6 O* r
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

* M: k0 F) H( o: [; }Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. $ H+ i; Q; r0 G; m
4. mvar:3 K* w& L: V, L$ I' x4 E' w
   
6. t0=clock(),# u. q% S9 F! Y/ |2 M' ^
   
8. s=0,
   
9. / L3 O  ?, ^  I( \9 G
10. k=0, while{k<10000,& Q- J\\" {' P  x' I
    
12.   i=0, while{i<1000,s=s+1, i++},\\" m# f) q! n2 @# A
    
14.   k++
    
15. ' `* N- V5 Y) }3 o0 N
16. },- m, V) _% X& O5 {\\" w
    
18. s;/ m8 J2 E5 D* f) _5 V
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">$ H+ i; Q; r0 G; m</span> mvar:<font class=\\"jammer\\">3 K* w& L: V, L$ I' x4 E' w</font><br /> t0=clock(),<font class=\\"jammer\\"># u. q% S9 F! Y/ |2 M' ^</font><br /> s=0,<br /> <span style=\\"display:none\\">/ L3 O  ?, ^  I( \9 G</span> k=0, while{k<10000,<font class=\\"jammer\\">& Q- J\\" {' P  x' I</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">\\" m# f) q! n2 @# A</font><br />   k++<br /> <span style=\\"display:none\\">' `* N- V5 Y) }3 o0 N</span> },<font class=\\"jammer\\">- m, V) _% X& O5 {\\" w</font><br /> s;<font class=\\"jammer\\">/ m8 J2 E5 D* f) _5 V</font><br /> [clock()-t0]/1000;

结果：
10000000.
) q( u, e+ e5 y3 n0.704  秒

二者运行效率相差7倍多。
: T" [) S- q) ]
2、循环中稍微加点东西
4 F& H- B3 y0 J8 ]7 L# C/ v1 a2 i( C' @
Matlab代码：

1. clear all
   6 @5 q9 R1 x2 f* m1 M- {# ^
2. tic( F7 b' \% \/ V4 B1 P
   
3. s=0;4 g) R0 K$ D* Z6 e0 p
   
4. for k=1:1000, ~- H4 C\" y! e1 k4 K% p  `/ o1 b) b8 [
   
5.   for i=1:10009 }, j6 }4 H* h& 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);' I1 L8 n* O6 Q% y
   
7.   end* T+ X3 t  z: c* m7 Y
   
8. end
   7 Q6 b# x$ `$ E& \& u0 p
9. s& D2 \# s2 a4 O8 Q
   
10. toc* _4 }, [2 R) h$ U- ]\" b
    
11. 
    ( T& W1 o7 a. _% U
12. s =
    \" b' _* f( s7 g1 M- V\" |
13. 
    - R7 l3 k& N  x( O
14.   1.1009e+007# q( h( v! I% Y; w
    
15. ' U6 t. p% |# o% C/ S0 ~* Q8 C
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];; K$ S; f! ~  V* S( `  G
   
4. mvar:
   
5. + W/ Q; A% y7 m+ k6 h' d1 J  e
6. t0=clock(),
   
7. 6 a9 Z1 }. G& ~/ |6 f# Q; ~& Z
8. s=0,* {\\" a, Z6 ^- G. A
   
10. k=1, while{k<=1000,
    
11. 6 a3 a- b% h* }5 Y* k
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. ! {( g# o, |5 J
14.   k++
    
15. ' p8 W7 k2 F' \  F* C9 x% m: S- {& m
16. },: F  G  y, C\\" j1 D9 D6 s! f, r
    
18. s;5 C/ i' \4 k, s, c7 Y
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">; K$ S; f! ~  V* S( `  G</font><br /> mvar:<br /> <span style=\\"display:none\\">+ W/ Q; A% y7 m+ k6 h' d1 J  e</span> t0=clock(),<br /> <span style=\\"display:none\\">6 a9 Z1 }. G& ~/ |6 f# Q; ~& Z</span> s=0,<font class=\\"jammer\\">* {\\" a, Z6 ^- G. A</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">6 a3 a- b% h* }5 Y* k</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\\">! {( g# o, |5 J</span>   k++<br /> <span style=\\"display:none\\">' p8 W7 k2 F' \  F* C9 x% m: S- {& m</span> },<font class=\\"jammer\\">: F  G  y, C\\" j1 D9 D6 s! f, r</font><br /> s;<font class=\\"jammer\\">5 C/ i' \4 k, s, c7 Y</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
  g" g& L# t1 I" \( V* W1.031

二者效率相差已比较小。

3、再增加一个函数调用
) k9 X1 i1 I8 \3 g5 J3 R. I
Matlab代码：

1. clear all
   ; U) |( E! a% j- b5 R* ]9 }& k( @
2. f=@(x)x+1;
   ; F7 W4 e0 O/ U  ?
3. tic+ q! s* {) g' T* y- p
   
4. s=0;2 n: ]1 B' f: V: C3 g( m* Q) h& N6 \! o
   
5. for k=1:1000
   ! _: ?# b' N) [1 y5 b0 F
6.   for i=1:10007 o6 D. |% J! C- V, y
   
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);
   ) H4 n0 N) e0 g) T
8.   end
   $ H& N$ N5 Q7 z7 q$ a
9. end
   4 [1 }( K5 G' \; D/ ]4 x# t
10. s1 Y' u, p/ a0 _/ `& D4 N\" \
    
11. toc- ]% N+ C0 ~  j; S8 Y\" r; G  v; C
    
12. - `2 N: l0 d7 L' R7 t
    
13. s =
    : f) I8 m. H' m4 a7 T
14. 1 z1 \7 F1 q$ F% R
    
15.   5.1251e+008
    . {7 N* u% d) b: r
16. \" P  W; t\" x3 k' P8 |6 Z
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 4 ~% [+ F  W1 w- C
4. f(x)=x+1;; X& a\\" L, W0 A, N5 x0 D% f
   
6. mvar:
   
7. % n8 Y& y) E* U( u\\" G
8. t0=clock(),
   
9. & u0 w! Y' K8 G. h' n  {% S
10. s=0,
    
11. 8 C% d8 D  ^! Q& i2 T8 z! ?: [* i
12. k=1, while{k<=1000,
    
13. . z- f8 S  Y! b1 V\\" H1 [+ g
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++},  |6 L- ]9 i8 W5 w* i4 L
    
16.   k++
    
17. # v( f) Z5 e& K$ N
18. },+ o, Z6 S4 X) r$ W9 d9 {6 M
    
20. s;' R8 V* d: d& G/ q2 J9 J
    
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">4 ~% [+ F  W1 w- C</span> f(x)=x+1;<font class=\\"jammer\\">; X& a\\" L, W0 A, N5 x0 D% f</font><br /> mvar:<br /> <span style=\\"display:none\\">% n8 Y& y) E* U( u\\" G</span> t0=clock(),<br /> <span style=\\"display:none\\">& u0 w! Y' K8 G. h' n  {% S</span> s=0,<br /> <span style=\\"display:none\\">8 C% d8 D  ^! Q& i2 T8 z! ?: [* i</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">. z- f8 S  Y! b1 V\\" H1 [+ g</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\\">  |6 L- ]9 i8 W5 w* i4 L</font><br />   k++<br /> <span style=\\"display:none\\"># v( f) Z5 e& K$ N</span> },<font class=\\"jammer\\">+ o, Z6 S4 X) r$ W9 d9 {6 M</font><br /> s;<font class=\\"jammer\\">' R8 V* d: d& G/ q2 J9 J</font><br /> [clock()-t0]/1000;

结果：
512508953.6658993
3 o* L0 q3 h! x# u# u: O: c1.078

可以看出，Matlab效率已不及Forcal。
: R& Y8 p  J6 D$ z故JIT加速器虽好，但效果似乎并不明显。

forcal

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

Forcal代码：

2. !using["sys"];
   
3. / j, w; Z: u, \4 R5 W
4. mvar:, ?- S! G) w! |' z
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组: W/ v1 Q& D6 g- ^
   
8. t0=clock(),
   
9. & i1 ?. \4 q5 {/ y  v% b\\" Y5 S
10. k=0, while{k<10000,4 W# S+ b+ E+ @8 m
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. ' {5 ]1 @0 c* ?
14.   k++4 A& m& r\\" j: ~$ H4 g
    
16. },% @' o( M. C% T3 o
    
18. [clock()-t0]/1000;  w; c5 V8 K' Z4 N- O* w* d: B
    
20. 
    
21. 9 w! g- u# F& h8 F' |3 [! U' s\\" a
22. a="\&\[1000]"/4,    //用字符串模拟静态数组1 ^: R* I; e, T1 F
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. ! |* m+ s5 N3 x( a
26. t0=clock(),$ O( R2 M0 U$ k1 M' e+ ?. H
    
28. s=0,
    
29. ( d$ s7 `0 r9 x8 _' W
30. k=0, while{k<10000,' h9 l: T3 Q. \$ X% u  V
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值% J8 M' o% N3 u
    
34.   k++
    
35. / O: L8 h0 |3 }) w( e
36. },
    
37. \\" ~* x8 [; H3 _) ~3 K) t5 K
38. s;2 n' |: t5 l7 B; P
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">/ j, w; Z: u, \4 R5 W</span> mvar:<font class=\\"jammer\\">, ?- S! G) w! |' z</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">: W/ v1 Q& D6 g- ^</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">& i1 ?. \4 q5 {/ y  v% b\\" Y5 S</span> k=0, while{k<10000,<font class=\\"jammer\\">4 W# S+ b+ E+ @8 m</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">' {5 ]1 @0 c* ?</span>   k++<font class=\\"jammer\\">4 A& m& r\\" j: ~$ H4 g</font><br /> },<font class=\\"jammer\\">% @' o( M. C% T3 o</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">  w; c5 V8 K' Z4 N- O* w* d: B</font><br /> <br /> <span style=\\"display:none\\">9 w! g- u# F& h8 F' |3 [! U' s\\" a</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">1 ^: R* I; e, T1 F</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">! |* m+ s5 N3 x( a</span> t0=clock(),<font class=\\"jammer\\">$ O( R2 M0 U$ k1 M' e+ ?. H</font><br /> s=0,<br /> <span style=\\"display:none\\">( d$ s7 `0 r9 x8 _' W</span> k=0, while{k<10000,<font class=\\"jammer\\">' h9 l: T3 Q. \$ X% u  V</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">% J8 M' o% N3 u</font><br />   k++<br /> <span style=\\"display:none\\">/ O: L8 h0 |3 }) w( e</span> },<br /> <span style=\\"display:none\\">\\" ~* x8 [; H3 _) ~3 K) t5 K</span> s;<font class=\\"jammer\\">2 n' |: t5 l7 B; P</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
. w+ g' Y% u0 C" g4 o20000000.
( @( @+ H4 K3 ^5 C0.579      //取值用时间，秒

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
4 \# R* k) ~8 [( B% F8 y0 M5 {
9 Y' O! L) y" V8 VMatlab代码：

1. >> clear all) n% Y- _! T7 b+ e6 |, K! |
   
2. a = linspace(1,1,1000);
   % P( T; v; c$ o' O
3. tic0 d- {9 \- u\" d
   
4. for k=1:10000
   ! [: f8 `9 K$ W. b# u0 j; U- ?
5.   for i=1:1000: ^. K! M# _; O1 Q8 c
   
6.     a(i)=2;
   , F  n$ A' L- q\" u$ ]0 \& O( z4 x
7.   end% M1 l# S9 N) Y+ s
   
8. end2 J, ~5 b2 Z* d
   
9. toc
   7 M/ n, V  `$ ?
10. + s$ P8 ?8 r\" P! G7 x7 r! l
    
11. tic
    / a$ K1 w. u( e3 X9 K4 q- i
12. s=0;
    6 w8 k; }' Q; @! {+ F\" p
13. for k=1:10000: u5 U9 O5 ~6 }( |
    
14.   for i=1:1000& `1 U+ {+ u# R/ `. G0 Q
    
15.     s=s+a(i);
    6 m0 m: X/ _  ?# z% g
16.   end
    1 {+ `4 ]& d: P\" I1 V( E
17. end
    6 U0 y7 ]4 D8 i\" v3 g' v, D
18. s$ ~3 H5 M4 X' O% e# D
    
19. toc; _# W' |5 L. ]) c5 U* B
    
20. Elapsed time is 0.140658 seconds.
    : O& K3 S6 Y: X2 {4 [8 X\" N
21. 
    ; ^6 G9 N# `. s% f2 P
22. s =\" D, s. K6 ]! \' O- g
    
23. 
    . n& f' H- c) U5 Z
24.     200000001 p5 W6 Z8 r0 p/ p( m
    
25.   D5 @8 @9 }5 A
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

& g/ y; C# p  R" V( e% UForcal代码：

2. !using["math","sys"];
   
3. # B7 B2 }+ ?5 u. I- E
4. mvar:5 Y\\" e3 @\\" J2 S8 u, Q\\" R
   
6. main0(::a)=
   
7. $ |  m5 \7 _8 |, g5 o' a# d$ }
8. a = linspace(1,1,1000),8 n) Y$ N, B/ Q2 h2 R+ M1 n
   
10. t0=clock(),
    
11.   x1 h9 R9 |3 y: F+ w( S4 l+ W
12. k=0, while{k<10000,; W( @5 t, G  m+ r! l# |* [: N2 f
    
14.   i=0, while{i<1000, a(i)=2, i++},2 H8 W' @+ n7 I7 |. A  p, C9 P$ L
    
16.   k++
    
17. ; J- S9 w- _0 r6 _& g2 P- u
18. },\\" e/ u. R* C9 E4 P\\" y6 O$ F( h\\" n$ @
    
20. [clock()-t0]/1000;3 T/ o. _( U- f) s0 \: j& h
    
22. 
    
23. % W/ d5 l. p8 A& h( V- l) B
24. main1(::a)=
    
25. - v, q/ f$ `* L1 l$ @
26. t0=clock(), s=0,
    
27. 4 @  H) b0 j- M2 r, S! v
28. k=0, while{k<10000,
    
29. ( F2 i) _7 x2 I0 R, ?* x  x
30.   i=0, while{i<1000, s=s+a(i), i++},; c0 h8 w2 `\\" t
    
32.   k++8 T& r7 F* c. Z8 m/ u: e
    
34. },
    
35. 7 B' D$ w5 @8 b
36. s;# W8 H  y2 B8 G2 N$ D4 w8 }
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\"># B7 B2 }+ ?5 u. I- E</span> mvar:<font class=\\"jammer\\">5 Y\\" e3 @\\" J2 S8 u, Q\\" R</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">$ |  m5 \7 _8 |, g5 o' a# d$ }</span> a = linspace(1,1,1000),<font class=\\"jammer\\">8 n) Y$ N, B/ Q2 h2 R+ M1 n</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">  x1 h9 R9 |3 y: F+ w( S4 l+ W</span> k=0, while{k<10000,<font class=\\"jammer\\">; W( @5 t, G  m+ r! l# |* [: N2 f</font><br />   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">2 H8 W' @+ n7 I7 |. A  p, C9 P$ L</font><br />   k++<br /> <span style=\\"display:none\\">; J- S9 w- _0 r6 _& g2 P- u</span> },<font class=\\"jammer\\">\\" e/ u. R* C9 E4 P\\" y6 O$ F( h\\" n$ @</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">3 T/ o. _( U- f) s0 \: j& h</font><br /> <br /> <span style=\\"display:none\\">% W/ d5 l. p8 A& h( V- l) B</span> main1(::a)=<br /> <span style=\\"display:none\\">- v, q/ f$ `* L1 l$ @</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">4 @  H) b0 j- M2 r, S! v</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">( F2 i) _7 x2 I0 R, ?* x  x</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">; c0 h8 w2 `\\" t</font><br />   k++<font class=\\"jammer\\">8 T& r7 F* c. Z8 m/ u: e</font><br /> },<br /> <span style=\\"display:none\\">7 B' D$ w5 @8 b</span> s;<font class=\\"jammer\\"># W8 H  y2 B8 G2 N$ D4 w8 }</font><br /> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
" D: ]3 k5 N' s5 g20000000.
5.531   //取数据所用时间，秒
" T0 [% h0 |8 F/ V, z5 l
9 `; o1 t5 V2 B* g: h－－－－－－－－－－
2 N5 T% e$ L9 N2 f5 B/ h: ?4 }; q, |
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. 0 w  i$ x4 h; Q  N& o# S. L\\" M
4. mvar:& Y- @$ i  H! p
   
6. main0(::a)=( P( C) j& P: {5 S% j
   
8. a = linspace(1,1,1000),: Q, Y1 V  ]: Z. P; S/ E
   
10. t0=clock(),
    
11. ( X. _, Q  x\\" ~2 h5 |+ `
12. k=0, while{k<10000,
    
13. & n( h! |* w8 P2 K1 q# p: ]1 L
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. , D/ V2 ?! |6 c, g\\" h3 J
16.   k+++ P  R! }) B1 e+ q; Y
    
18. },; j6 G. O7 F' J# f0 f/ L0 L
    
20. [clock()-t0]/1000;. }& {* H) a* R# f$ n7 B4 g! B8 r
    
22. 9 M% |# o6 m) y/ |4 p& v; o
    
24. main1(::a)=
    
25. 5 y* v1 P2 ]8 S4 l* [$ t( e
26. t0=clock(), s=0,% {! F+ j! }4 U+ N0 W
    
28. k=0, while{k<10000,5 c: v$ ^5 E( t
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},! M+ `/ V4 d* B, q& v8 s
    
32.   k++- u- n; }3 B\\" c- F
    
34. },! [; G* a+ L% E
    
36. s;+ K9 i4 Q& D# e* r
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">0 w  i$ x4 h; Q  N& o# S. L\\" M</span> mvar:<font class=\\"jammer\\">& Y- @$ i  H! p</font><br /> main0(::a)=<font class=\\"jammer\\">( P( C) j& P: {5 S% j</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">: Q, Y1 V  ]: Z. P; S/ E</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">( X. _, Q  x\\" ~2 h5 |+ `</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">& n( h! |* w8 P2 K1 q# p: ]1 L</span>   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">, D/ V2 ?! |6 c, g\\" h3 J</span>   k++<font class=\\"jammer\\">+ P  R! }) B1 e+ q; Y</font><br /> },<font class=\\"jammer\\">; j6 G. O7 F' J# f0 f/ L0 L</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">. }& {* H) a* R# f$ n7 B4 g! B8 r</font><br /> <font class=\\"jammer\\">9 M% |# o6 m) y/ |4 p& v; o</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">5 y* v1 P2 ]8 S4 l* [$ t( e</span> t0=clock(), s=0,<font class=\\"jammer\\">% {! F+ j! }4 U+ N0 W</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">5 c: v$ ^5 E( t</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">! M+ `/ V4 d* B, q& v8 s</font><br />   k++<font class=\\"jammer\\">- u- n; }3 B\\" c- F</font><br /> },<font class=\\"jammer\\">! [; G* a+ L% E</font><br /> s;<font class=\\"jammer\\">+ K9 i4 Q& D# e* r</font><br /> [clock()-t0]/1000;

结果：
10.469
20000000.
3.781

6 |" @- g8 d+ T$ k' q－－－－－－－

; D9 J$ S- L! _- V在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
9 o  w/ f. w5 P. V在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
& g( @3 X) ^4 b7 ]( b+ P1 ^  x
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
8 J! e8 K( J6 `# a6 p* p
Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

+ f- I* o) l0 T; W6 Y, v5、在循环中加点东西

matlab代码：

1. >> clear all) R! ^, I  v: P' t
   
2. f=@(x)x+1;9 k0 _, w9 B/ v
   
3. a = linspace(2,2,1000);
   ( d9 V9 H+ _8 ~* z' Z; E! u
4. tic1 r& t( s5 b2 o6 h
   
5. s=0;
   ( [/ e\" m3 a/ Z3 L0 a, ?
6. for k=1:10005 X& }: T: {' S; d8 e7 {
   
7.   for i=1:1000
   * i- J/ ]) A( B
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));; M% C+ I3 w! r\" f( U; g  \1 z$ y8 ]6 P5 y
   
9.   end) B, g# R( a% K0 a# u# [
   
10. end
    % A9 S# A* w- S) O
11. s
    + B: ]  A5 @% ]
12. toc
      `0 @& W' a9 N3 N/ W
13. # [9 \% |* D: w  x) \3 O
    
14. s =
    ; O# c& L3 S! m4 N/ v
15. ; v2 B8 |! @( w  f3 J
    
16.   2.4002e+007$ ]% K. T5 b: z9 u' W, w! t2 t2 N
    
17. 
    6 [' N& r& L7 {3 l! X2 d7 ~( y
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];, G0 B8 m4 P' F& m0 o
   
4. mvar:
   
5. % @. u  ~6 H8 A2 I
6. f(x)=x+1;9 q9 z$ T+ |# M- E5 @. b9 r! I
   
8. main0(::a)=8 n, z! ?# D* R, }) A
   
10. a = linspace(2,2,1000),
    
11. & {0 X\\" }; Y  _
12. t0=clock(), s=0,
    
13. ( h9 O/ _/ n- D0 j; ~2 p) u: ~
14. k=0, while{k<1000,9 l& B) W) C( Q2 F: @7 F
    
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++},; |& D: p8 w- t& o7 ^
    
18.   k+++ U8 F6 b5 I# J9 E0 ~5 e
    
20. },1 y% m3 X+ N  I. y( }7 B8 y6 a- K
    
22. s;( v, U+ [& K3 l# c: o! g
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">, G0 B8 m4 P' F& m0 o</font><br /> mvar:<br /> <span style=\\"display:none\\">% @. u  ~6 H8 A2 I</span> f(x)=x+1;<font class=\\"jammer\\">9 q9 z$ T+ |# M- E5 @. b9 r! I</font><br /> main0(::a)=<font class=\\"jammer\\">8 n, z! ?# D* R, }) A</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">& {0 X\\" }; Y  _</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">( h9 O/ _/ n- D0 j; ~2 p) u: ~</span> k=0, while{k<1000,<font class=\\"jammer\\">9 l& B) W) C( Q2 F: @7 F</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\\">; |& D: p8 w- t& o7 ^</font><br />   k++<font class=\\"jammer\\">+ U8 F6 b5 I# J9 E0 ~5 e</font><br /> },<font class=\\"jammer\\">1 y% m3 X+ N  I. y( }7 B8 y6 a- K</font><br /> s;<font class=\\"jammer\\">( v, U+ [& K3 l# c: o! g</font><br /> [clock()-t0]/1000;

结果：
8 ^) v) d* V5 E$ k/ v! d24002271.69189384
5.172

5 Y  t# l6 Y) p* w/ F9 i; \二者效率现在只有一倍差距。
# `9 \  b  G5 T7 m4 O! q
－－－－－－－－－－
- S" U' c  b% u0 B
不过，显然matlab和Forcal的程序应该这样写：

4 Z7 ?% }. l4 i1 \/ I1 ^matlab程序：

1. clear all
   2 X0 d# i. I1 H( R5 X
2. f=@(x)x+1;  ?8 e9 {3 ]& S8 z9 o
   
3. a = linspace(2,2,1000);
   0 t% S\" H- v& w% T& V
4. tic
   # B; |* U- F1 a# u3 X7 n\" r
5. s=0;( d; ~  N; B% C0 O2 n
   
6. for k=1:10005 u# j( w% X9 J. I; t# Q& y
   
7.   for i=1:1000
   : b' m! O3 l! b% Q& Y6 C3 d
8.     t=a(i);
     q. F+ G$ K7 p6 X. G+ E
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);6 x3 K- V$ `5 z+ X$ l
   
10.   end
    ' t. X, {3 G/ o; o$ g1 z
11. end0 p+ ~) n8 d, L6 g3 H
    
12. s
    * |  g\" w3 ]2 ]$ g$ e' @. L
13. toc, I# _4 r; Q% h- _0 K
    
14. ) T. z  a. l! l0 s3 x6 ^
    
15. s =0 a0 Q% v/ C\" c7 p' C/ n9 o4 v9 O4 n
    
16. 2 u2 Q, b+ @+ i4 g, N
    
17.   2.4002e+0075 s3 F4 L5 }* w- H6 I
    
18. 5 P# z! p  L: Y! f& s- M
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. $ G1 H7 |( h3 J* f
4. mvar:
   
5. ( C- H% g/ M/ P! ], Q1 F. ^
6. f(x)=x+1;$ b: _$ u* u: i! E* r
   
8. main0(::a)=+ X- T9 ~9 b( A! d% |
   
10. a = linspace(2,2,1000),
    
11. 5 Y! I; Z2 S6 ^
12. t0=clock(), s=0,
    
13. 9 @: Z) U9 x6 q: D: w
14. k=0, while{k<1000,
    
15. 0 [4 x4 z8 I6 H: F% r
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++},5 q5 u  S. ?& M/ q, w
    
18.   k++! ^- X6 V6 m+ x+ d6 c, G9 x$ ^
    
20. },
    
21. * x7 |2 `! t  i, @# n7 l0 m' |7 \
22. s;
    
23. ; `7 M# J& ~* q; q7 Z7 a4 a$ J
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">$ G1 H7 |( h3 J* f</span> mvar:<br /> <span style=\\"display:none\\">( C- H% g/ M/ P! ], Q1 F. ^</span> f(x)=x+1;<font class=\\"jammer\\">$ b: _$ u* u: i! E* r</font><br /> main0(::a)=<font class=\\"jammer\\">+ X- T9 ~9 b( A! d% |</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">5 Y! I; Z2 S6 ^</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">9 @: Z) U9 x6 q: D: w</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">0 [4 x4 z8 I6 H: F% r</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\\">5 q5 u  S. ?& M/ q, w</font><br />   k++<font class=\\"jammer\\">! ^- X6 V6 m+ x+ d6 c, G9 x$ ^</font><br /> },<br /> <span style=\\"display:none\\">* x7 |2 `! t  i, @# n7 l0 m' |7 \</span> s;<br /> <span style=\\"display:none\\">; `7 M# J& ~* q; q7 Z7 a4 a$ J</span> [clock()-t0]/1000;

结果：
+ s+ w6 @! n( A9 n- h2 }24002271.69189384
- e  e* ?7 y* Q  _  o1.625

matlab效率已不及Forcal。

1 N8 W- r( I+ h  x' N1 X6 V% G- m' M提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

海水