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

forcal

1、简单循环
8 ]% W( d1 |+ f) W# q" Q
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all. a! V( ?0 O* }
   
2. tic! m2 _$ y, `! A  [! @) ^
   
3. s=0;
   5 d5 W# T$ \* j2 ]/ F3 v
4. for k=1:10000
   7 S0 v: A! q/ u- R) {
5.   for i=1:10002 l* [. Y. j0 r% A, U0 i# Y6 X/ W
   
6.     s=s+1;0 D4 U' v8 M$ g# i& O
   
7.   end5 u/ k2 m7 q7 G1 X5 X/ ?! k$ E
   
8. end, u1 t- ?2 p; y( m9 S8 U7 v
   
9. s
   % _& |9 D; C( ~
10. toc& O2 x* w1 P) ]/ q5 A* D4 u
    
11. % Z: e& `# z+ J# q# z& w
    
12. s =3 G. T: v- C' {7 U5 k
    
13. 4 d3 ?6 Z3 a. D; B\" V0 f
    
14.     10000000
    * P7 n% s' D2 w- f+ H, b
15. 
    ( f* C- p, P* `3 n* \1 G
16. Elapsed time is 0.092808 seconds.

复制代码

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

0 _# T( H4 R* Z. P7 U1 lForcal没有JIT加速器，代码：

2. !using["sys"];
   
3. \\" L2 g9 e2 K+ a# V/ k5 l1 J( @  h
4. mvar:
   
5. 6 k4 e) q  t2 B6 O2 r. Z, M! T
6. t0=clock(),
   
7. % k) h0 K! p1 e* J
8. s=0,
   
9. # Q8 w0 y4 m+ v/ }6 I
10. k=0, while{k<10000,# _' K' L, W& `$ f; ?
    
12.   i=0, while{i<1000,s=s+1, i++},7 P& D9 d* B7 c) W( F6 W0 y
    
14.   k++) n- b( u' h\\" m+ b, O+ A; E
    
16. },
    
17. 5 q/ Y! t* G1 J! C/ r\\" b8 d
18. s;1 ~. H6 k, [: O  g! z+ i' t! @
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">\\" L2 g9 e2 K+ a# V/ k5 l1 J( @  h</span> mvar:<br /> <span style=\\"display:none\\">6 k4 e) q  t2 B6 O2 r. Z, M! T</span> t0=clock(),<br /> <span style=\\"display:none\\">% k) h0 K! p1 e* J</span> s=0,<br /> <span style=\\"display:none\\"># Q8 w0 y4 m+ v/ }6 I</span> k=0, while{k<10000,<font class=\\"jammer\\"># _' K' L, W& `$ f; ?</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">7 P& D9 d* B7 c) W( F6 W0 y</font><br />   k++<font class=\\"jammer\\">) n- b( u' h\\" m+ b, O+ A; E</font><br /> },<br /> <span style=\\"display:none\\">5 q/ Y! t* G1 J! C/ r\\" b8 d</span> s;<font class=\\"jammer\\">1 ~. H6 k, [: O  g! z+ i' t! @</font><br /> [clock()-t0]/1000;

结果：
10000000.
8 E2 C3 }) \$ o/ D( `0.704  秒

二者运行效率相差7倍多。
) W2 \5 ^: e+ D6 r+ b8 R! G8 x
. R1 i" S* N0 j+ o' R2、循环中稍微加点东西
' Y. H: k8 x  H' G
2 J1 e4 b# p/ H5 U) U/ lMatlab代码：

1. clear all( D6 K\" r1 |$ y
   
2. tic- x: V8 a. j6 l5 l
   
3. s=0;
   $ v  X- |) h# p- n0 g: i1 W
4. for k=1:1000
   9 ]  {' t4 m, Z/ {
5.   for i=1:10009 T# T( b7 `\" o2 O
   
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);- m& L- A. x# L: s
   
7.   end
   7 Q# x* b# U1 T\" @/ r
8. end
   3 z* v$ r1 V7 {5 o
9. s
   1 u; n  F. L! Y* W& O
10. toc
    - }5 R  Q/ D. D7 S
11. 
    5 I# J! h  ^. C4 V; G# e
12. s =& Z+ s- `; C1 s/ f( f5 S; P
    
13. 
    - `/ m; F, i9 K' ?
14.   1.1009e+007
    ' ?; T$ Y$ |. t+ t) ^
15. 
    & [% U( H  |7 m- A# _$ S5 C% Z
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. \\" {4 {' L\\" [# t; k7 u( `# f
4. mvar:/ M. G& ?\\" D3 s# u2 [: F% s
   
6. t0=clock(),
   
7. ' Z' y. H/ O. c! q3 d& k
8. s=0,0 X& L  w  B\\" b9 `2 `: a# n
   
10. k=1, while{k<=1000,
    
11. 8 J2 e6 b3 f2 o. r
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 h6 S2 |. G: K
14.   k++
    
15. 8 {$ J3 ~0 e* v3 s# {3 ^) p
16. },1 q$ c1 P- q% L5 C, w\\" L$ G
    
18. s;
    
19. % [: b# j$ Q7 m+ q6 ?7 ~2 O# _/ }9 B
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">\\" {4 {' L\\" [# t; k7 u( `# f</span> mvar:<font class=\\"jammer\\">/ M. G& ?\\" D3 s# u2 [: F% s</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">' Z' y. H/ O. c! q3 d& k</span> s=0,<font class=\\"jammer\\">0 X& L  w  B\\" b9 `2 `: a# n</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">8 J2 e6 b3 f2 o. r</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 h6 S2 |. G: K</span>   k++<br /> <span style=\\"display:none\\">8 {$ J3 ~0 e* v3 s# {3 ^) p</span> },<font class=\\"jammer\\">1 q$ c1 P- q% L5 C, w\\" L$ G</font><br /> s;<br /> <span style=\\"display:none\\">% [: b# j$ Q7 m+ q6 ?7 ~2 O# _/ }9 B</span> [clock()-t0]/1000;

结果：
11008953.66602346
1.031

" S% g/ x  W8 ?3 h: y% ^: t二者效率相差已比较小。
" X( `: l. t0 x7 x
, |0 M6 E3 [& Q" T3、再增加一个函数调用
  K+ ?6 j$ Z. Y. I
Matlab代码：

1. clear all% w( M  H& f% C* O! p
   
2. f=@(x)x+1;
   9 q: H\" K( q! C2 L
3. tic: V9 V- m# c! z8 B6 S
   
4. s=0;) D8 R# g: V' W. y; W\" R6 c
   
5. for k=1:1000- H& k2 r0 d2 k! ~
   
6.   for i=1:1000
   # I) y1 u* F. d) i0 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);+ h+ n. j* U+ ?
   
8.   end* ]6 f( M9 T( S9 R; R: h
   
9. end
   : F, `) T9 o7 ?6 N# O
10. s, T: B* r3 \% d! A4 w6 H
    
11. toc
      b5 v$ F0 x! i& G1 P8 X& P! A: L  m
12. 1 Z# d4 N; i5 X2 O7 M
    
13. s =/ W6 K1 |; L- \  Q/ G, B- W
    
14. 1 T* Y7 R7 |- ]$ C2 l! o
    
15.   5.1251e+008
    : x2 _! f8 H3 Z\" w
16. 4 N( X; w: `, ~* [' I, |+ e
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];* ]- a2 }/ F8 R
   
4. f(x)=x+1;0 ^$ N, z! Q7 R5 Q, `3 C  P
   
6. mvar:8 G7 E9 m+ m. Q; W\\" d
   
8. t0=clock(),* Y\\" }- Q\\" g% Y
   
10. s=0,
    
11. / J2 @9 L* q% s- X
12. k=1, while{k<=1000,
    
13. 2 M* e9 z1 P5 Z! F3 y: o' I
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++},& o* K$ h2 H6 F3 i3 J. A8 x1 j
    
16.   k++
    
17. 3 }1 [/ C1 P# X5 v5 B) a
18. },  D% U' P8 U. m+ \* B; [, C
    
20. s;
    
21. ! u0 i4 k; J! j4 t
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">* ]- a2 }/ F8 R</font><br /> f(x)=x+1;<font class=\\"jammer\\">0 ^$ N, z! Q7 R5 Q, `3 C  P</font><br /> mvar:<font class=\\"jammer\\">8 G7 E9 m+ m. Q; W\\" d</font><br /> t0=clock(),<font class=\\"jammer\\">* Y\\" }- Q\\" g% Y</font><br /> s=0,<br /> <span style=\\"display:none\\">/ J2 @9 L* q% s- X</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">2 M* e9 z1 P5 Z! F3 y: o' I</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\\">& o* K$ h2 H6 F3 i3 J. A8 x1 j</font><br />   k++<br /> <span style=\\"display:none\\">3 }1 [/ C1 P# X5 v5 B) a</span> },<font class=\\"jammer\\">  D% U' P8 U. m+ \* B; [, C</font><br /> s;<br /> <span style=\\"display:none\\">! u0 i4 k; J! j4 t</span> [clock()-t0]/1000;

结果：
512508953.6658993
1.078
/ e. s+ r, V! k3 t  \  Z  P; h; w
5 W8 \6 T7 K) H5 }. q可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
7 J  ~% c0 _3 V  g: n, q# z
Matlab代码：

1. >> clear all1 g8 h' N, |; `% v+ Y
   
2. a = linspace(1,1,1000);- @& F. x, \3 R8 Q- B# \0 |
   
3. tic
   0 A6 n9 r. T% r! w. U' u/ ^
4. for k=1:10000  z% X, [7 ~0 h! [- A
   
5.   for i=1:1000! u  Y: `  T( z. j, D! d/ F: M$ X
   
6.     a(i)=2;
   - P- c! [9 P- D\" o# m
7.   end
   - {8 l2 u1 q  ~+ A/ f  ^0 e
8. end, C* x\" S6 [7 b
   
9. toc# F6 I& V, C: B7 ]  `
   
10. 
    : ~1 W1 L: h* @: ]2 h+ {2 |; d( \  x
11. tic
    0 {' b\" K\" p, b+ B8 h\" r
12. s=0;
    % d; m  U% ?# b+ i  K
13. for k=1:100001 p3 j( G& J( R; c6 B- D
    
14.   for i=1:1000) {1 S' O6 v! O/ V8 j+ Z
    
15.     s=s+a(i);( s$ |2 p. Y3 o3 v, c4 L# R
    
16.   end; h( Q- ^1 g1 P
    
17. end
    \" N1 C& Y3 H5 o8 G1 {1 D
18. s
    5 T9 c! W+ ^: g8 V- f
19. toc
    6 ]2 J. E\" z1 w$ H6 K' M
20. Elapsed time is 0.140658 seconds.
    ! Z) ?, _' \% u, i) U& X' }! T4 N
21. 9 G/ v- y6 l0 x$ u5 W\" |
    
22. s =! i1 |7 Q\" V- N
    
23. , ]6 ~9 J& \! f. ~4 B- K; {0 I
    
24.     20000000
    7 k! }/ v* ^: T, P: ?4 e
25. 5 Y  [+ \4 U$ F+ I
    
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
# l5 b: f" I0 j% T
Forcal代码：

2. !using["math","sys"];& q* s$ a  {7 B6 Q. o- i2 f
   
4. mvar:/ G; d+ A! p& Z4 u5 k
   
6. main0(::a)=/ }  k# T% w1 ?
   
8. a = linspace(1,1,1000),9 U& Z2 k( Z/ N/ V
   
10. t0=clock(),
    
11. , z; R- p# }1 T# s: {\\" s\\" }
12. k=0, while{k<10000,$ i5 O! ~3 i\\" V3 i8 z
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. 2 V9 d! ]$ W: L\\" J- ?, J
16.   k++- A7 M5 V+ ]8 T; a
    
18. },
    
19. ( c( ]4 J# p( o' S# z+ a1 h4 M0 f
20. [clock()-t0]/1000;
    
21. ! e  V4 w/ g6 q6 @  e$ p9 h7 w  a
22. 
    
23. 7 y% [. S1 \# k4 W$ z5 t% _
24. main1(::a)=, O0 C$ P' v% u6 I0 `) ^$ W1 V) q/ p
    
26. t0=clock(), s=0,
    
27. $ y+ o5 \  z) k\\" R* k' ]. Z
28. k=0, while{k<10000,
    
29. # Z, G) j/ f* `5 z* T, m
30.   i=0, while{i<1000, s=s+a(i), i++},& T7 @: V; X& D+ J
    
32.   k++
    
33. ! T  B; h) `: a5 s# [( D% j
34. },
    
35. 6 Z4 O+ h3 N  y\\" ]
36. s;  s4 j0 M+ c9 q3 L$ d
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">& q* s$ a  {7 B6 Q. o- i2 f</font><br /> mvar:<font class=\\"jammer\\">/ G; d+ A! p& Z4 u5 k</font><br /> main0(::a)=<font class=\\"jammer\\">/ }  k# T% w1 ?</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">9 U& Z2 k( Z/ N/ V</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">, z; R- p# }1 T# s: {\\" s\\" }</span> k=0, while{k<10000,<font class=\\"jammer\\">$ i5 O! ~3 i\\" V3 i8 z</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">2 V9 d! ]$ W: L\\" J- ?, J</span>   k++<font class=\\"jammer\\">- A7 M5 V+ ]8 T; a</font><br /> },<br /> <span style=\\"display:none\\">( c( ]4 J# p( o' S# z+ a1 h4 M0 f</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">! e  V4 w/ g6 q6 @  e$ p9 h7 w  a</span> <br /> <span style=\\"display:none\\">7 y% [. S1 \# k4 W$ z5 t% _</span> main1(::a)=<font class=\\"jammer\\">, O0 C$ P' v% u6 I0 `) ^$ W1 V) q/ p</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">$ y+ o5 \  z) k\\" R* k' ]. Z</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\"># Z, G) j/ f* `5 z* T, m</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">& T7 @: V; X& D+ J</font><br />   k++<br /> <span style=\\"display:none\\">! T  B; h) `: a5 s# [( D% j</span> },<br /> <span style=\\"display:none\\">6 Z4 O+ h3 N  y\\" ]</span> s;<font class=\\"jammer\\">  s4 j0 M+ c9 q3 L$ d</font><br /> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
; U9 C0 m& @$ A20000000.
5.531   //取数据所用时间，秒
4 ^& f0 l! p" Q$ Q. v
% [9 F/ }" R# }8 s; p" e+ K% J9 T6 Y－－－－－－－－－－

& w/ J- S/ ]- w; G% W; H/ i用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. % r4 \- O: Q& R) D0 e! @
4. mvar:
   
5. , \7 N4 L# I3 M% D0 p
6. main0(::a)=# |2 Z7 i1 M# X
   
8. a = linspace(1,1,1000),& I: w+ s- S0 u* o) o+ s
   
10. t0=clock(),! D1 X/ B8 V  H( i4 b) L( D
    
12. k=0, while{k<10000,' Q- `( L- X3 [2 Q) z/ M
    
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. ; c\\" S1 A5 C! a/ a, [' [! {7 O
16.   k++
    
17. 8 \4 T0 T5 E. \3 Z
18. },
    
19. ; W8 R8 K) `9 ]3 W1 ?( I5 s, D3 q
20. [clock()-t0]/1000;1 [# h! l, T6 ?- r/ m$ |9 D2 {
    
22. 
    
23. / |' e6 g) A, E0 M( d! |
24. main1(::a)=
    
25.   |6 i4 d# g* d1 e\\" K4 v
26. t0=clock(), s=0,8 o5 o/ `) `# }2 X/ o7 ~, t0 o
    
28. k=0, while{k<10000,
    
29. ' N. K% l\\" y; }+ I
30.   i=0, while{i<1000, s=s+A(a,i), i++},: `, u: e- ^1 U
    
32.   k++
    
33. 0 k0 v. k! V\\" ^
34. },6 X+ W/ D( u+ B6 U2 X2 ^3 ~
    
36. s;
    
37. 0 O9 H\\" Y; U4 t1 o- {  _% t8 I
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">% r4 \- O: Q& R) D0 e! @</span> mvar:<br /> <span style=\\"display:none\\">, \7 N4 L# I3 M% D0 p</span> main0(::a)=<font class=\\"jammer\\"># |2 Z7 i1 M# X</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">& I: w+ s- S0 u* o) o+ s</font><br /> t0=clock(),<font class=\\"jammer\\">! D1 X/ B8 V  H( i4 b) L( D</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">' Q- `( L- X3 [2 Q) z/ M</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">; c\\" S1 A5 C! a/ a, [' [! {7 O</span>   k++<br /> <span style=\\"display:none\\">8 \4 T0 T5 E. \3 Z</span> },<br /> <span style=\\"display:none\\">; W8 R8 K) `9 ]3 W1 ?( I5 s, D3 q</span> [clock()-t0]/1000;<font class=\\"jammer\\">1 [# h! l, T6 ?- r/ m$ |9 D2 {</font><br /> <br /> <span style=\\"display:none\\">/ |' e6 g) A, E0 M( d! |</span> main1(::a)=<br /> <span style=\\"display:none\\">  |6 i4 d# g* d1 e\\" K4 v</span> t0=clock(), s=0,<font class=\\"jammer\\">8 o5 o/ `) `# }2 X/ o7 ~, t0 o</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">' N. K% l\\" y; }+ I</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">: `, u: e- ^1 U</font><br />   k++<br /> <span style=\\"display:none\\">0 k0 v. k! V\\" ^</span> },<font class=\\"jammer\\">6 X+ W/ D( u+ B6 U2 X2 ^3 ~</font><br /> s;<br /> <span style=\\"display:none\\">0 O9 H\\" Y; U4 t1 o- {  _% t8 I</span> [clock()-t0]/1000;

结果：
10.469
20000000.
5 Y4 ^: m4 y+ u. ~6 i# d7 K4 }. ?4 B3.781

" O8 h0 I9 M5 V8 l/ ?－－－－－－－
4 m/ b5 V! C6 W) Q/ `
5 \, ], \! }; K: Y在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
) k! Q/ e8 N6 Q/ l7 i; c
, p- O. R$ }4 F$ d$ i# b& p; z" @matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

# C- t) Y  L- x' [  c) |& ]Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
  y2 J: i3 i" r0 a
2 o' [3 F9 ~. s5、在循环中加点东西
# @. C' [" o( L& F7 {
7 W8 O3 N4 S" [matlab代码：

1. >> clear all
   + D: @/ R4 S& C# P# x
2. f=@(x)x+1;( C0 O' M+ U. K) C/ ]' ]' A& o9 j
   
3. a = linspace(2,2,1000);
   5 x9 [+ C5 v, A, i7 T# N
4. tic/ B1 s% ~4 ~% `5 l: P& M
   
5. s=0;% P, t4 V: J' ^+ c# ]: R5 @
   
6. for k=1:1000
   0 Z, r8 `3 q, J) L7 s
7.   for i=1:1000/ k* V. _4 K* A* ]7 E
   
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/ O9 R\" G/ T$ O  {
   
9.   end
   \" d+ ^9 i* _- }) v# \4 d4 G
10. end
    3 f9 ~+ C' T! ~' N% G8 }7 u
11. s\" O! c% q7 }' w- S3 K
    
12. toc
    / Q6 O4 N  [$ D1 d% i) S
13. 1 O1 ^& Y. ~$ G5 C; t% c
    
14. s =
    # \& i( l/ N' P  z) P. A
15. / F  G' j: G6 G+ o\" H$ J# N
    
16.   2.4002e+007
    6 `- |  B4 t6 S
17. 
    ( N2 Z) X7 E\" b- K2 R4 T5 k
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];' e1 ^; K) T- @
   
4. mvar:0 \  |/ I/ S8 @$ G( }; A% t1 d
   
6. f(x)=x+1;5 F1 |. y7 d9 G: j% S
   
8. main0(::a)=) j$ E' U; H) B2 e* R. @
   
10. a = linspace(2,2,1000),
    
11. 8 d# W\\" }6 N$ B% y. y& V. L\\" B8 D
12. t0=clock(), s=0,
    
13. / ]! w+ q# o6 l0 R! c7 W
14. k=0, while{k<1000,4 g3 h. t# \  f# K' }/ D) N3 C  n
    
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. 2 l/ X\\" k/ P( a  {  a$ w
18.   k++! p3 c0 h9 C1 T6 r5 W& }5 _
    
20. },
    
21. $ s+ [2 j0 Y\\" V4 D
22. s;
    
23. 1 I  N8 B5 ]4 F- {/ H
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">' e1 ^; K) T- @</font><br /> mvar:<font class=\\"jammer\\">0 \  |/ I/ S8 @$ G( }; A% t1 d</font><br /> f(x)=x+1;<font class=\\"jammer\\">5 F1 |. y7 d9 G: j% S</font><br /> main0(::a)=<font class=\\"jammer\\">) j$ E' U; H) B2 e* R. @</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">8 d# W\\" }6 N$ B% y. y& V. L\\" B8 D</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">/ ]! w+ q# o6 l0 R! c7 W</span> k=0, while{k<1000,<font class=\\"jammer\\">4 g3 h. t# \  f# K' }/ D) N3 C  n</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++},<br /> <span style=\\"display:none\\">2 l/ X\\" k/ P( a  {  a$ w</span>   k++<font class=\\"jammer\\">! p3 c0 h9 C1 T6 r5 W& }5 _</font><br /> },<br /> <span style=\\"display:none\\">$ s+ [2 j0 Y\\" V4 D</span> s;<br /> <span style=\\"display:none\\">1 I  N8 B5 ]4 F- {/ H</span> [clock()-t0]/1000;

结果：
9 E0 Z0 E7 U4 @* h: v1 A' x24002271.69189384
5.172

二者效率现在只有一倍差距。
- o, P; V7 T- Y) Q
$ d; ]$ u$ C1 v" R－－－－－－－－－－
, \; l4 e3 Z: f6 v, H
不过，显然matlab和Forcal的程序应该这样写：
3 `- N; N5 d5 ?7 r: Z5 p5 X
' k. N8 y  B# ~7 M" smatlab程序：

1. clear all; G0 p  U$ X# G; M/ J
   
2. f=@(x)x+1;/ a' v4 @; c\" T) u
   
3. a = linspace(2,2,1000);
   ! U: g* n7 ]\" }7 @
4. tic; B/ f+ t( C5 y9 @8 f9 P; a) c. f
   
5. s=0;
   ( w0 X, t0 W! P# f
6. for k=1:1000  f; f5 m0 K; }
   
7.   for i=1:10007 m8 ]3 D/ c& _- q& S% s! j: h
   
8.     t=a(i);
   $ S3 S7 A& Q/ B
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);9 ^: Q0 s) ^; }0 X/ J7 \$ u$ f9 `
   
10.   end
    # x  d3 S5 v\" a
11. end
    , D  ?: h, l+ ?2 M7 i
12. s$ A- ~1 s3 I6 Z- N  {$ z/ r' y
    
13. toc1 b0 Y- k& s: U' N\" t
    
14. ) {; X6 R9 O- ?+ |
    
15. s =; J( N0 n) v) S9 U. u
    
16. . B. u9 j! O9 ?8 i
    
17.   2.4002e+007
    , c6 i0 P. j( m
18. & g5 ?/ d2 s  ]' w5 [  q
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. 9 L3 ^, w) V& W4 M! [( r+ a7 B
4. mvar:$ F4 o- ?1 K; }* K9 d0 O7 A9 @' }
   
6. f(x)=x+1;9 y' |4 P4 b9 f\\" W  ~
   
8. main0(::a)=/ o1 d/ V2 s. |
   
10. a = linspace(2,2,1000),% i9 I( \0 t+ _3 Z
    
12. t0=clock(), s=0,
    
13. 0 s/ G( i+ t' L5 Y
14. k=0, while{k<1000,; g- ^\\" s( }5 O% @, i9 p
    
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. - M1 p, C1 V2 U. \( N3 x0 f% I8 ^
18.   k++
    
19.   D6 R$ R$ I  |2 O
20. },
    
21. / O. [\\" x, c- {4 G& x
22. s;\\" G& L8 v5 a6 X0 `, g
    
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">9 L3 ^, w) V& W4 M! [( r+ a7 B</span> mvar:<font class=\\"jammer\\">$ F4 o- ?1 K; }* K9 d0 O7 A9 @' }</font><br /> f(x)=x+1;<font class=\\"jammer\\">9 y' |4 P4 b9 f\\" W  ~</font><br /> main0(::a)=<font class=\\"jammer\\">/ o1 d/ V2 s. |</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">% i9 I( \0 t+ _3 Z</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">0 s/ G( i+ t' L5 Y</span> k=0, while{k<1000,<font class=\\"jammer\\">; g- ^\\" s( }5 O% @, i9 p</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\\">- M1 p, C1 V2 U. \( N3 x0 f% I8 ^</span>   k++<br /> <span style=\\"display:none\\">  D6 R$ R$ I  |2 O</span> },<br /> <span style=\\"display:none\\">/ O. [\\" x, c- {4 G& x</span> s;<font class=\\"jammer\\">\\" G& L8 v5 a6 X0 `, g</font><br /> [clock()-t0]/1000;

结果：
$ Z7 a" l  @( G3 y: W5 T24002271.69189384
% t1 H: e& W& f9 u# Q4 _1.625

matlab效率已不及Forcal。

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

forcal

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

* _5 R, q& o2 n" ?Forcal代码：

2. !using["sys"];
   
3.   N# m. J* }1 F' ~6 x2 \' [# ^
4. mvar:
   
5. 5 R2 ], {  H7 q7 Y; [/ H
6. a="\&\[1000]"/4,    //用字符串模拟静态数组9 y' j& ]& Q# w0 i6 `
   
8. t0=clock(),
   
9. 0 L0 ]% J$ c& g& W
10. k=0, while{k<10000,, p5 J2 U; L5 n
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. ; i9 c# w: g9 N, k$ w4 J+ ^
14.   k++
    
15. - n- R+ r9 K( z7 z6 E5 A9 g
16. },
    
17. 9 h+ ]9 T3 [% C\\" Z3 H
18. [clock()-t0]/1000;6 O: K# ~& S! ~& D' V5 z# v
    
20. + z  x6 q' x& L8 K, y) g0 K
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. 6 W1 ?+ O! j& h, a3 p
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值5 k1 h0 Q, R5 C\\" |. f
    
26. t0=clock(),
    
27. $ V5 _% A2 H1 E' B1 Q& e6 X
28. s=0,
    
29. 8 q6 ?/ A+ M+ R0 ?\\" r
30. k=0, while{k<10000,
    
31. 5 J4 g0 R/ r% m5 J- H: k2 ~) e
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值7 e1 h* ?/ Y+ a  k; i
    
34.   k++/ h9 s7 l( R\\" \
    
36. },
    
37. - R$ }. k  P3 W( Q
38. s;7 k' m7 v5 h5 r+ T( `7 C$ j
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">  N# m. J* }1 F' ~6 x2 \' [# ^</span> mvar:<br /> <span style=\\"display:none\\">5 R2 ], {  H7 q7 Y; [/ H</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">9 y' j& ]& Q# w0 i6 `</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">0 L0 ]% J$ c& g& W</span> k=0, while{k<10000,<font class=\\"jammer\\">, p5 J2 U; L5 n</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">; i9 c# w: g9 N, k$ w4 J+ ^</span>   k++<br /> <span style=\\"display:none\\">- n- R+ r9 K( z7 z6 E5 A9 g</span> },<br /> <span style=\\"display:none\\">9 h+ ]9 T3 [% C\\" Z3 H</span> [clock()-t0]/1000;<font class=\\"jammer\\">6 O: K# ~& S! ~& D' V5 z# v</font><br /> <font class=\\"jammer\\">+ z  x6 q' x& L8 K, y) g0 K</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">6 W1 ?+ O! j& h, a3 p</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">5 k1 h0 Q, R5 C\\" |. f</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">$ V5 _% A2 H1 E' B1 Q& e6 X</span> s=0,<br /> <span style=\\"display:none\\">8 q6 ?/ A+ M+ R0 ?\\" r</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">5 J4 g0 R/ r% m5 J- H: k2 ~) e</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">7 e1 h* ?/ Y+ a  k; i</font><br />   k++<font class=\\"jammer\\">/ h9 s7 l( R\\" \</font><br /> },<br /> <span style=\\"display:none\\">- R$ }. k  P3 W( Q</span> s;<font class=\\"jammer\\">7 k' m7 v5 h5 r+ T( `7 C$ j</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
20000000.
7 Q& m! A7 K0 f2 L- k0.579      //取值用时间，秒