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

forcal

1、简单循环

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

1. clear all
   - A$ J: k  Q7 J3 w/ j1 N
2. tic
   $ y. c7 l\" _% Z. K+ s/ n( K, W: X
3. s=0;
   ) R\" L4 K7 ]% e9 `9 k6 h9 m2 Z
4. for k=1:10000/ m% `: Y\" z, n; a6 }
   
5.   for i=1:1000
   5 h+ k8 [1 \6 M7 ~, O2 ?8 c+ N
6.     s=s+1;
   ! c* \2 z1 z+ e) `) W6 L/ t
7.   end\" [9 p* b% @) g4 e1 P
   
8. end$ }\" h4 o& u/ C. `1 m
   
9. s\" F, b7 P8 ~5 t/ r% e7 ?. x
   
10. toc( R8 W8 I3 Q6 T0 r
    
11. 
    1 U. @7 ~! O$ I9 p4 \- I
12. s =/ O) [+ t8 C\" Q5 T\" q
    
13. 
    0 Y! v+ ?: F\" e$ l: X1 ]
14.     10000000
    / j7 V& B! N( }; V/ \
15. . y# M$ ~5 f  Z1 `: T/ |6 k
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. 5 m2 D. P0 Y/ R; x( B  E# H
4. mvar:
   
5. 1 H: q' {' L( z0 ^0 q\\" }* v( C
6. t0=clock(),3 H  r- X\\" P, L' I( e* \1 i
   
8. s=0,
   
9.   }' S- g\\" c1 k* N. K$ R6 ]
10. k=0, while{k<10000,- F$ R& A/ Z7 s, G$ n8 m
    
12.   i=0, while{i<1000,s=s+1, i++},! ^- T9 ?% k$ S1 W, d
    
14.   k++) ]0 B\\" }$ Q5 [+ k
    
16. },
    
17. * h2 d8 U# v. r+ X8 s+ R: B
18. s;
    
19. ; n% P/ J$ |+ J; }; s
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">5 m2 D. P0 Y/ R; x( B  E# H</span> mvar:<br /> <span style=\\"display:none\\">1 H: q' {' L( z0 ^0 q\\" }* v( C</span> t0=clock(),<font class=\\"jammer\\">3 H  r- X\\" P, L' I( e* \1 i</font><br /> s=0,<br /> <span style=\\"display:none\\">  }' S- g\\" c1 k* N. K$ R6 ]</span> k=0, while{k<10000,<font class=\\"jammer\\">- F$ R& A/ Z7 s, G$ n8 m</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">! ^- T9 ?% k$ S1 W, d</font><br />   k++<font class=\\"jammer\\">) ]0 B\\" }$ Q5 [+ k</font><br /> },<br /> <span style=\\"display:none\\">* h2 d8 U# v. r+ X8 s+ R: B</span> s;<br /> <span style=\\"display:none\\">; n% P/ J$ |+ J; }; s</span> [clock()-t0]/1000;

结果：
10000000.
. d3 n) r% }$ v# |: B0.704  秒
& k- C: \" A8 x0 j* o9 I
二者运行效率相差7倍多。

9 K; N6 z2 i0 b* k: K# \. R2、循环中稍微加点东西
; z, z  n; Q8 ]* n. n& I1 N
Matlab代码：

1. clear all
   $ U1 |% F- O, K& B& S* ?. R
2. tic, `+ p% K! i1 n5 H: r
   
3. s=0;
   # C, ?2 o9 A\" |, Z  Q
4. for k=1:1000, u3 A+ n$ \0 v' s
   
5.   for i=1:1000
   ) l! |6 u0 B; m0 E$ n' E( S
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);) z# m0 `  \0 i- d4 v4 m
   
7.   end
   7 C2 n' A0 e8 }1 w& l( A2 @
8. end
   \" u8 a# }+ j% A8 w/ M# o
9. s  N4 `0 Z5 U$ W% t\" W9 E+ W
   
10. toc0 x+ A* ~\" U% P' f$ E
    
11. 
    8 W% F/ q# |2 Y' v2 L; [/ q
12. s =
    0 ^8 ], U; T8 L0 J( I
13. ) t, h8 F6 R2 U; z( Q
    
14.   1.1009e+007\" r! ~! w  x; U2 ]/ B
    
15. \" Q9 L( D+ k6 B' X
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];/ F4 L8 j# ]' i  K8 a) M
   
4. mvar:\\" V- p3 Q; r  E/ l! `  {* E% a
   
6. t0=clock(),
   
7. ! D/ J5 }; g6 J. D
8. s=0,/ [. y7 @# Z; M( x0 ~8 `
   
10. k=1, while{k<=1000,8 R; v- R  \0 h) }( r% H/ 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++},
    
13. \\" g. Q1 J3 X) d' K, M\\" T
14.   k++
    
15. 9 [$ v$ o. _3 @6 X4 N/ c# j
16. },
    
17. 4 m, {2 P- j, Y6 @6 q  {
18. s;5 C* M' v. d2 S( M
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">/ F4 L8 j# ]' i  K8 a) M</font><br /> mvar:<font class=\\"jammer\\">\\" V- p3 Q; r  E/ l! `  {* E% a</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">! D/ J5 }; g6 J. D</span> s=0,<font class=\\"jammer\\">/ [. y7 @# Z; M( x0 ~8 `</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">8 R; v- R  \0 h) }( r% H/ 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++},<br /> <span style=\\"display:none\\">\\" g. Q1 J3 X) d' K, M\\" T</span>   k++<br /> <span style=\\"display:none\\">9 [$ v$ o. _3 @6 X4 N/ c# j</span> },<br /> <span style=\\"display:none\\">4 m, {2 P- j, Y6 @6 q  {</span> s;<font class=\\"jammer\\">5 C* M' v. d2 S( M</font><br /> [clock()-t0]/1000;

结果：
4 A- ^  ~# }: S+ ~( Z) Q7 Q11008953.66602346
1.031
3 [5 P: \% B8 M* V1 B. d
二者效率相差已比较小。

( A, @6 ~6 `7 W  [, o* F3、再增加一个函数调用
$ O0 B$ g  G. ~( }( }
4 J" M+ V  R# T/ AMatlab代码：

1. clear all% g, y1 ^* H' d  F( D7 E8 ]
   
2. f=@(x)x+1;: r8 x, h6 `/ F- c: j0 N\" t  Y( U
   
3. tic
   9 [( c+ N$ @1 W. J. K\" O
4. s=0;
   ) @: v3 ^9 J5 v4 y
5. for k=1:1000
   9 x  E/ M9 ?% K- \6 {. Q
6.   for i=1:1000
   ! z9 M+ ]; l9 ?8 |4 H; [- `; B
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);$ }5 M$ W  H( }+ f3 L
   
8.   end0 R* |9 D, K+ t% K9 c0 i  w% B
   
9. end/ m  \2 A. t4 H( _& }5 R
   
10. s9 S4 [! @4 J* e7 T
    
11. toc+ ^  Z9 ~- h7 Z, `$ }+ i. G
    
12. ' b% `) ]\" @+ j
    
13. s =
    + F+ u& W) K) ~0 G0 m5 U9 x
14. 1 a# `4 \3 V5 {
    
15.   5.1251e+008$ {( e& a7 X9 @, M8 o
    
16. 4 G* M' R; x4 \3 x
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];/ B6 a2 C3 A\\" S2 c0 _, M4 ?! j- k
   
4. f(x)=x+1;& R5 O2 G# A, @% @$ A
   
6. mvar:6 ^+ _4 H% H2 @, u9 r
   
8. t0=clock(),
   
9. 5 x3 g! w2 b, `1 W
10. s=0,& u7 X4 o2 U  M. A
    
12. k=1, while{k<=1000,
    
13. ; T  q7 y) p# r+ i  ~9 F
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++},, B. S0 y; l\\" P0 {8 h' C' j7 n
    
16.   k++/ G& Q8 I+ S' ~' }2 R+ {
    
18. },
    
19. 3 \8 T: n1 s\\" l$ P- X5 N& T
20. s;
    
21. 3 P' V! e* }) t0 w$ x% }: x8 e
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">/ B6 a2 C3 A\\" S2 c0 _, M4 ?! j- k</font><br /> f(x)=x+1;<font class=\\"jammer\\">& R5 O2 G# A, @% @$ A</font><br /> mvar:<font class=\\"jammer\\">6 ^+ _4 H% H2 @, u9 r</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">5 x3 g! w2 b, `1 W</span> s=0,<font class=\\"jammer\\">& u7 X4 o2 U  M. A</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">; T  q7 y) p# r+ i  ~9 F</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\\">, B. S0 y; l\\" P0 {8 h' C' j7 n</font><br />   k++<font class=\\"jammer\\">/ G& Q8 I+ S' ~' }2 R+ {</font><br /> },<br /> <span style=\\"display:none\\">3 \8 T: n1 s\\" l$ P- X5 N& T</span> s;<br /> <span style=\\"display:none\\">3 P' V! e* }) t0 w$ x% }: x8 e</span> [clock()-t0]/1000;

结果：
512508953.6658993
1.078
; K' ]# D' O! H' z  P
/ F& H0 |/ n% W$ D可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

forcal

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

Forcal代码：

2. !using["sys"];: U( f* A9 t8 ]0 w( f  X
   
4. mvar:
   
5. 2 P' N- b9 Q! [5 K% {3 `
6. a="\&\[1000]"/4,    //用字符串模拟静态数组6 }! {8 H  ?/ k- ^, e
   
8. t0=clock(),
   
9. 1 W. _) _) ]3 Y' A, H5 n
10. k=0, while{k<10000,0 s+ D* k/ c\\" V* u3 j$ @2 f
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. , O' @5 m4 F# f2 d1 z& Q* i% d
14.   k++7 j\\" v/ v* t+ `6 P+ S
    
16. },
    
17. # Z# n2 q# ~; X3 i* x
18. [clock()-t0]/1000;) O& M4 N, q% Z* b, |3 c7 j7 i$ d% u
    
20. , L3 l4 ^4 D2 i: J( a9 L
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. 8 ]& t7 c% w/ Z, X. I3 J7 Q
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值; g0 c, b2 |: Z; x* B' M
    
26. t0=clock(),
    
27. 7 Q3 p; U- f$ w2 c
28. s=0,; J8 H) k: Z1 {\\" N8 m\\" h, O; ^
    
30. k=0, while{k<10000,
    
31. & \\\" h\\" [; H5 b8 V& a$ G* S
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值. p( N9 }% r6 t2 W0 b5 ~! X
    
34.   k+++ f' a. s6 m# U0 \
    
36. },/ d. z7 C\\" P1 W' K
    
38. s;
    
39. \\" Z2 y  {9 g: O1 C! F! y0 W
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">: U( f* A9 t8 ]0 w( f  X</font><br /> mvar:<br /> <span style=\\"display:none\\">2 P' N- b9 Q! [5 K% {3 `</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">6 }! {8 H  ?/ k- ^, e</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">1 W. _) _) ]3 Y' A, H5 n</span> k=0, while{k<10000,<font class=\\"jammer\\">0 s+ D* k/ c\\" V* u3 j$ @2 f</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">, O' @5 m4 F# f2 d1 z& Q* i% d</span>   k++<font class=\\"jammer\\">7 j\\" v/ v* t+ `6 P+ S</font><br /> },<br /> <span style=\\"display:none\\"># Z# n2 q# ~; X3 i* x</span> [clock()-t0]/1000;<font class=\\"jammer\\">) O& M4 N, q% Z* b, |3 c7 j7 i$ d% u</font><br /> <font class=\\"jammer\\">, L3 l4 ^4 D2 i: J( a9 L</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">8 ]& t7 c% w/ Z, X. I3 J7 Q</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">; g0 c, b2 |: Z; x* B' M</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">7 Q3 p; U- f$ w2 c</span> s=0,<font class=\\"jammer\\">; J8 H) k: Z1 {\\" N8 m\\" h, O; ^</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">& \\\" h\\" [; H5 b8 V& a$ G* S</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">. p( N9 }% r6 t2 W0 b5 ~! X</font><br />   k++<font class=\\"jammer\\">+ f' a. s6 m# U0 \</font><br /> },<font class=\\"jammer\\">/ d. z7 C\\" P1 W' K</font><br /> s;<br /> <span style=\\"display:none\\">\\" Z2 y  {9 g: O1 C! F! y0 W</span> [clock()-t0]/1000;

结果：
; N- {) a1 x7 Q: e0 G+ q6 B5 o" B0.5620000000000001  //赋值用时间，秒
. n: ^$ C+ `6 m1 u2 Q4 y20000000.
0.579      //取值用时间，秒

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
" F' v5 H6 u; y( K* F$ M6 `- p
Matlab代码：

1. >> clear all2 O6 w( d' u3 v( x: C
   
2. a = linspace(1,1,1000);* M* k9 U9 t- J
   
3. tic
   % D1 p0 M) k& i# \2 F; {
4. for k=1:10000
   3 {) C; m) u) P/ y5 A& V1 u' }7 s) F
5.   for i=1:10000 u1 J# s# e! L1 _
   
6.     a(i)=2;% ]& f  f6 {' U1 y5 |7 ~. ]3 p
   
7.   end/ ]4 l6 O: z9 ~$ B  n
   
8. end7 g! \$ Q& P5 J9 C  a
   
9. toc
   ( M. k$ q% T6 H\" _
10. 
    ( A- Z) t7 Y2 \+ v9 Y' P* b
11. tic7 f! F2 R4 b# i  j/ c5 C. Q
    
12. s=0;* Z7 v( X6 }) F. P\" q, x
    
13. for k=1:10000
    8 s1 F, y7 V6 ]
14.   for i=1:1000' Q$ D3 \4 L. Z# M% w' a
    
15.     s=s+a(i);
    + I8 o) V, p8 J( [6 D
16.   end
    7 A5 Y% K\" j& q2 ^6 X
17. end
    \" b* T0 f6 [* O! {9 B8 N
18. s
    0 C  U2 i) `0 h- r( Y% f1 _! Z& _
19. toc
    - ^; ^6 v* M& s, Z3 m
20. Elapsed time is 0.140658 seconds.; B2 r- x: o. O$ @% P9 M
    
21. * I5 P. W0 \+ U1 g
    
22. s =
    7 }\" P' T, M% T- _\" c
23. 2 ^. ~0 v$ }, k! t4 B! _
    
24.     20000000
    * K  m8 N; A\" }3 R
25. , _; t7 J& i7 Q# Q+ R/ b; L
    
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
7 G0 C% Q; N0 s1 \. v. q7 J
# F+ d- H$ Y1 n8 }0 \9 n* hForcal代码：

2. !using["math","sys"];6 L3 a# ]- d  r3 y% k1 n
   
4. mvar:
   
5. $ @! J# s* \8 r' O( ^8 E8 I
6. main0(::a)=
   
7. , m7 ?$ F3 _9 x, o) ]1 a' a3 u
8. a = linspace(1,1,1000),
   
9. ' `. `, W0 L2 M* G( ]
10. t0=clock(),6 R& z( }- E+ p0 l* R
    
12. k=0, while{k<10000,: h. A; t7 x: i
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. + _) d# F6 _5 H5 D) H  B7 |  m) F, ^
16.   k++7 z6 C2 Y; ~: t  }3 _2 M
    
18. },- X  P# H3 j! N
    
20. [clock()-t0]/1000;
    
21. % k2 E* o7 T. P! k: s9 }: Y  z
22. 
    
23. 2 v: o& j\\" w! I1 r
24. main1(::a)=
    
25. . z+ f  H+ R. u8 e. u
26. t0=clock(), s=0,* k  J, C) l; h/ \/ |3 O
    
28. k=0, while{k<10000,0 v1 r+ V1 W9 I. L\\" c
    
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. ! o  D5 b3 s, k2 w: O# x: @
32.   k++9 x. h( }9 B9 F. s/ S$ f$ z; I, W
    
34. },
    
35. % R9 p3 O5 ^4 c6 {) L* y; w3 N- U9 ?% e
36. s;
    
37. ' D2 J5 w; K* ~2 N+ f  s
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">6 L3 a# ]- d  r3 y% k1 n</font><br /> mvar:<br /> <span style=\\"display:none\\">$ @! J# s* \8 r' O( ^8 E8 I</span> main0(::a)=<br /> <span style=\\"display:none\\">, m7 ?$ F3 _9 x, o) ]1 a' a3 u</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">' `. `, W0 L2 M* G( ]</span> t0=clock(),<font class=\\"jammer\\">6 R& z( }- E+ p0 l* R</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">: h. A; t7 x: i</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">+ _) d# F6 _5 H5 D) H  B7 |  m) F, ^</span>   k++<font class=\\"jammer\\">7 z6 C2 Y; ~: t  }3 _2 M</font><br /> },<font class=\\"jammer\\">- X  P# H3 j! N</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">% k2 E* o7 T. P! k: s9 }: Y  z</span> <br /> <span style=\\"display:none\\">2 v: o& j\\" w! I1 r</span> main1(::a)=<br /> <span style=\\"display:none\\">. z+ f  H+ R. u8 e. u</span> t0=clock(), s=0,<font class=\\"jammer\\">* k  J, C) l; h/ \/ |3 O</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">0 v1 r+ V1 W9 I. L\\" c</font><br />   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">! o  D5 b3 s, k2 w: O# x: @</span>   k++<font class=\\"jammer\\">9 x. h( }9 B9 F. s/ S$ f$ z; I, W</font><br /> },<br /> <span style=\\"display:none\\">% R9 p3 O5 ^4 c6 {) L* y; w3 N- U9 ?% e</span> s;<br /> <span style=\\"display:none\\">' D2 J5 w; K* ~2 N+ f  s</span> [clock()-t0]/1000;

结果：
+ g0 L, }( r7 V2 \: A1 D. ~0 F15.484  //存数据所用时间，秒
; V) m3 s4 ?/ R( f( i% B20000000.
5.531   //取数据所用时间，秒
5 y' U4 G' L2 R. u" `9 w& ^2 u
－－－－－－－－－－

用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];6 y5 Q  ?9 E0 m\\" l
   
4. mvar:
   
5. # [% N- Q2 G/ `  m. h$ @
6. main0(::a)=
   
7. ; U( x7 ]5 u9 ?6 |  i, }# ^
8. a = linspace(1,1,1000),
   
9. + m( n5 G% v1 m& M: O8 M/ |
10. t0=clock(),( G9 s; D( ^  f2 [: V
    
12. k=0, while{k<10000,
    
13. 7 [! u1 ], X1 n% I8 B2 X
14.   i=0, while{i<1000, A(a,i)=2, i++},. ]* J* y0 m( d: P
    
16.   k++
    
17. / [. }. q9 K7 Y) r3 q
18. },
    
19. 2 t7 B3 ?. M& F& c+ C
20. [clock()-t0]/1000;; x\\" y9 @: l5 w! P+ F
    
22. : z/ V5 l\\" u% T4 H' ~( A' h* f; D
    
24. main1(::a)=
    
25. 6 T: F& ~% i+ }4 }, H' b6 g\\" y
26. t0=clock(), s=0,
    
27. . T! O1 b/ j& f; }: z/ W\\" Q+ ]# F
28. k=0, while{k<10000,
    
29. 1 d) j( p+ V* ?* D
30.   i=0, while{i<1000, s=s+A(a,i), i++},' x$ r. V3 t$ u1 ?. K8 b: l
    
32.   k++
    
33. 5 J+ o6 l\\" V, \: Y) n
34. },
    
35. ' ?* n8 g6 q, {) k; N4 F6 f
36. s;
    
37. & j6 k5 j! O( N  X1 P
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">6 y5 Q  ?9 E0 m\\" l</font><br /> mvar:<br /> <span style=\\"display:none\\"># [% N- Q2 G/ `  m. h$ @</span> main0(::a)=<br /> <span style=\\"display:none\\">; U( x7 ]5 u9 ?6 |  i, }# ^</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">+ m( n5 G% v1 m& M: O8 M/ |</span> t0=clock(),<font class=\\"jammer\\">( G9 s; D( ^  f2 [: V</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">7 [! u1 ], X1 n% I8 B2 X</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">. ]* J* y0 m( d: P</font><br />   k++<br /> <span style=\\"display:none\\">/ [. }. q9 K7 Y) r3 q</span> },<br /> <span style=\\"display:none\\">2 t7 B3 ?. M& F& c+ C</span> [clock()-t0]/1000;<font class=\\"jammer\\">; x\\" y9 @: l5 w! P+ F</font><br /> <font class=\\"jammer\\">: z/ V5 l\\" u% T4 H' ~( A' h* f; D</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">6 T: F& ~% i+ }4 }, H' b6 g\\" y</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">. T! O1 b/ j& f; }: z/ W\\" Q+ ]# F</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">1 d) j( p+ V* ?* D</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">' x$ r. V3 t$ u1 ?. K8 b: l</font><br />   k++<br /> <span style=\\"display:none\\">5 J+ o6 l\\" V, \: Y) n</span> },<br /> <span style=\\"display:none\\">' ?* n8 g6 q, {) k; N4 F6 f</span> s;<br /> <span style=\\"display:none\\">& j6 k5 j! O( N  X1 P</span> [clock()-t0]/1000;

结果：
9 A* L2 U6 j  q10.469
20000000.
1 W6 s, S- U- w3.781
( V1 U0 N4 `8 I  F
3 \* p$ L8 ~( [4 \& i－－－－－－－

在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
, @" u2 V. D7 I在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

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

Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
4 B' s( g9 t4 `4 B0 l* B
5、在循环中加点东西

. K% z  m% |" e" kmatlab代码：

1. >> clear all  O: B- q' k) C  ~. |* {. W4 Y
   
2. f=@(x)x+1;9 F( I# d4 H( u( U# |+ y. R
   
3. a = linspace(2,2,1000);
   - o: U# a: i* M* e; b$ A
4. tic
   $ n& i+ C, n; {) t
5. s=0;
   , P9 _- B6 Z& T
6. for k=1:1000
   ( J6 u% P* V: F8 k5 s
7.   for i=1:1000# L' U. P7 M  d( h\" j% P* o0 Y: g
   
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));0 [1 {0 p6 G* C5 w9 c\" j
   
9.   end( Y( _# m2 e- Q+ ?& G
   
10. end, c! r) n3 U) B! h
    
11. s
    4 I* ~8 Y7 |1 `9 e0 F! W
12. toc4 F0 Y2 g% h1 J* G
    
13. 
    & W; y2 Q: }: D  k1 e
14. s =$ a3 B$ ?. q' ^+ c$ T' c: r
    
15. 
    + T: y# q# C3 l& C9 x! g
16.   2.4002e+007
    , ?4 v1 _* a: P( G\" ], r! Y( I
17. 
    : w/ V) `/ \8 K6 g4 s$ l3 u8 h
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];8 n; F) v! A( j6 P2 `
   
4. mvar:! D2 v9 l, W1 }! p, l
   
6. f(x)=x+1;
   
7. ) W$ }, A* ?0 C- {
8. main0(::a)=* T7 z' I\\" T! R  ^+ Y4 L  q
   
10. a = linspace(2,2,1000),/ d! e3 c! w: ]/ z, _\\" I
    
12. t0=clock(), s=0,7 Y( [/ P5 E& f0 w
    
14. k=0, while{k<1000,  L$ B, g7 ]6 k& g5 C3 l
    
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++},* C2 W, ]\\" K- ^& Q: r2 r  _7 j
    
18.   k++$ V9 S8 _( m4 W. \# [
    
20. },- s; m+ N2 C) V: G' Z+ t
    
22. s;% ]; h# ]* j\\" m
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">8 n; F) v! A( j6 P2 `</font><br /> mvar:<font class=\\"jammer\\">! D2 v9 l, W1 }! p, l</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">) W$ }, A* ?0 C- {</span> main0(::a)=<font class=\\"jammer\\">* T7 z' I\\" T! R  ^+ Y4 L  q</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">/ d! e3 c! w: ]/ z, _\\" I</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">7 Y( [/ P5 E& f0 w</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">  L$ B, g7 ]6 k& g5 C3 l</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\\">* C2 W, ]\\" K- ^& Q: r2 r  _7 j</font><br />   k++<font class=\\"jammer\\">$ V9 S8 _( m4 W. \# [</font><br /> },<font class=\\"jammer\\">- s; m+ N2 C) V: G' Z+ t</font><br /> s;<font class=\\"jammer\\">% ]; h# ]* j\\" m</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
5 q4 G3 q/ K! P( j5.172
  S5 Q+ M* [) h1 L# Z; J
二者效率现在只有一倍差距。

－－－－－－－－－－

2 h& W7 `1 t+ H1 K不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

1. clear all9 G\" h! [) X1 G' i% A% {
   
2. f=@(x)x+1;
   / F# w9 k: {% D! K( i3 q
3. a = linspace(2,2,1000);6 o- T% c6 D1 `* s$ h
   
4. tic. [: C5 F& l, S0 j6 h9 t
   
5. s=0;$ T9 }( x' U4 M0 x  S
   
6. for k=1:1000
     v4 U5 w- d- i, ~8 r; R
7.   for i=1:10001 [- x: v& m* b; p
   
8.     t=a(i);/ O; S; U/ k) u
   
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);
   # x, L3 n% Y# I
10.   end
    \" n2 Q( G3 @$ k( w6 |. J( z2 i
11. end
    / Y/ z+ F8 M5 W, {+ r7 M# z' ]
12. s' ^# k0 y+ ^3 Q
    
13. toc
    : _, O3 _1 \: }: m$ Z, ^\" _; i
14. 
    & a5 @5 m  ]- x/ x! ]* P
15. s =; |) S6 L' [& ?! V
    
16. ' T) ]$ H* |7 C9 L$ {
    
17.   2.4002e+007
    4 w6 D\" }! m3 g- ?% I
18. 
    \" N. [8 H  b) H7 Q
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. 9 B# u6 q- @\\" w5 m; z' \
4. mvar:/ M; _7 n5 R: ~- @# u! s; c8 Z
   
6. f(x)=x+1;
   
7. \\" \, [9 e* {+ S4 I
8. main0(::a)=8 r# o3 Q- z9 K5 h4 m& I9 o4 Z
   
10. a = linspace(2,2,1000),
    
11. ' Q) e! c+ `6 p/ \! e3 ]
12. t0=clock(), s=0,
    
13. ) h: b# e7 h' i5 H+ r
14. k=0, while{k<1000,9 g2 N- C; v; z1 Y
    
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++},* N* C5 v1 ?2 X& N: m
    
18.   k++% {2 L: m  N& ~5 F4 S- @% B% k
    
20. },
    
21. + R; f. E& w; r6 x
22. s;
    
23. / d/ C: q' U+ N! _' p4 Z) l
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">9 B# u6 q- @\\" w5 m; z' \</span> mvar:<font class=\\"jammer\\">/ M; _7 n5 R: ~- @# u! s; c8 Z</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">\\" \, [9 e* {+ S4 I</span> main0(::a)=<font class=\\"jammer\\">8 r# o3 Q- z9 K5 h4 m& I9 o4 Z</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">' Q) e! c+ `6 p/ \! e3 ]</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">) h: b# e7 h' i5 H+ r</span> k=0, while{k<1000,<font class=\\"jammer\\">9 g2 N- C; v; z1 Y</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\\">* N* C5 v1 ?2 X& N: m</font><br />   k++<font class=\\"jammer\\">% {2 L: m  N& ~5 F4 S- @% B% k</font><br /> },<br /> <span style=\\"display:none\\">+ R; f. E& w; r6 x</span> s;<br /> <span style=\\"display:none\\">/ d/ C: q' U+ N! _' p4 Z) l</span> [clock()-t0]/1000;

结果：
1 I; i, }. `% _24002271.69189384
8 _8 o8 ?+ Z" J: u! {1.625

' i) y! p7 x) ]1 `matlab效率已不及Forcal。

# d( R  w: y9 v$ e7 y提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

海水