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

forcal

1、简单循环

  v& r2 F* P: R& h, ?/ OMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   ; U1 t6 r6 r1 i2 H1 O- ~
2. tic: N7 M/ q3 s+ y* R( ]6 k
   
3. s=0;0 Y  d; X& `1 F9 l, M# ]
   
4. for k=1:10000
   1 l5 n, }7 C) s- ~
5.   for i=1:1000% F; o2 q  C+ o5 ^# s
   
6.     s=s+1;
   5 ~& p$ {, c3 h& z* _6 l; Z. a8 X
7.   end
   ( z) t8 ^& V; q+ B% \- t
8. end$ V- r3 i% `, D: ~
   
9. s6 c. s/ m4 O3 X' m0 x7 U- M3 ^4 X
   
10. toc
    % M: s2 T4 W0 b* y) I6 b
11. 
    6 i  ]\" h$ g' N/ |
12. s =
    ! |8 K' E: e$ Q6 w4 c1 {
13. 
    2 e; ]! Q/ h% o$ G+ G2 k
14.     10000000
    : ^- m$ S& J! Y! W
15. 
    3 |' p; N8 R! C- G
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
& T6 T" H$ m, X  g  {
& J1 q/ B4 ?9 y5 B2 j* c# M, j0 NForcal没有JIT加速器，代码：

2. !using["sys"];
   
3. % \, m7 v, P! H. `0 J8 O4 c
4. mvar:3 V% ~$ y9 W. t9 U, |6 @
   
6. t0=clock(),
   
7. & M\\" V7 u8 }, {5 ~2 A/ c+ m
8. s=0,$ a0 A3 o% j$ }  |) ?( [
   
10. k=0, while{k<10000,\\" w6 d$ X1 Q+ U* x% Z2 d# |  e
    
12.   i=0, while{i<1000,s=s+1, i++},
    
13. ) j0 q- F: I! h8 E$ W
14.   k++
    
15. ( a& q. M7 D+ [6 o% O$ R7 `$ w- U
16. },# L9 D3 |- J( ~  g& G  q5 K4 I
    
18. s;\\" p% c3 N4 y) g1 p: d4 ?6 G+ Z0 x& o* _
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">% \, m7 v, P! H. `0 J8 O4 c</span> mvar:<font class=\\"jammer\\">3 V% ~$ y9 W. t9 U, |6 @</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">& M\\" V7 u8 }, {5 ~2 A/ c+ m</span> s=0,<font class=\\"jammer\\">$ a0 A3 o% j$ }  |) ?( [</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">\\" w6 d$ X1 Q+ U* x% Z2 d# |  e</font><br />   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">) j0 q- F: I! h8 E$ W</span>   k++<br /> <span style=\\"display:none\\">( a& q. M7 D+ [6 o% O$ R7 `$ w- U</span> },<font class=\\"jammer\\"># L9 D3 |- J( ~  g& G  q5 K4 I</font><br /> s;<font class=\\"jammer\\">\\" p% c3 N4 y) g1 p: d4 ?6 G+ Z0 x& o* _</font><br /> [clock()-t0]/1000;

结果：
6 n' K' y& ^, B/ U8 q, q! H10000000.
# B) e0 ]* y6 k, G5 g0.704  秒

$ S# d. J) r8 x/ n" [二者运行效率相差7倍多。

2、循环中稍微加点东西

Matlab代码：

1. clear all
   . t+ W* }2 }8 z; t$ @
2. tic3 y\" i3 |1 H9 G
   
3. s=0;
   3 q+ p( V: K* B2 P2 I+ W7 D3 m
4. for k=1:10001 ?; ?) y. W! ]8 M9 d+ ]
   
5.   for i=1:1000
   ) h& Z& g: c& V* `
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! V: j+ c3 x$ n! H
   
7.   end
   ' M& |5 _# x4 S- P% ]6 {
8. end; R7 u2 W- D: ], a$ w5 f( P& p
   
9. s2 \# O9 y% F/ u% p9 ?
   
10. toc
    2 U3 @/ b8 g0 g5 z\" g/ O( a
11. 2 I  K$ ~+ E4 J\" ~3 W7 C. _
    
12. s =
    \" N  M3 R  `) g$ M' g
13. $ i2 i6 D8 F& X  B& s, ?
    
14.   1.1009e+007- @9 L! O5 P0 Y0 o! h. m6 ?( S
    
15. 2 d\" t\" c2 k8 K1 l7 u: ~
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];6 S3 C% B0 v% ]$ A' L3 p# J\\" A
   
4. mvar:\\" e7 A\\" N' ~; D# t( x
   
6. t0=clock(),\\" Q- H% T& z/ B* s  t2 }
   
8. s=0,\\" Z- ^- R8 B: [3 f
   
10. k=1, while{k<=1000,
    
11. ! ?  x: c6 t0 T+ Y* G- ~
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. # D: c1 G' t! d1 i& W3 e
14.   k++
    
15. ( n$ p% P* R- V0 n+ R5 S
16. },
    
17. 5 {3 F) D* V! V6 R& h\\" w\\" T: k
18. s;. O6 k, Y' W0 x; [6 u# ]\\" c) |
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">6 S3 C% B0 v% ]$ A' L3 p# J\\" A</font><br /> mvar:<font class=\\"jammer\\">\\" e7 A\\" N' ~; D# t( x</font><br /> t0=clock(),<font class=\\"jammer\\">\\" Q- H% T& z/ B* s  t2 }</font><br /> s=0,<font class=\\"jammer\\">\\" Z- ^- R8 B: [3 f</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">! ?  x: c6 t0 T+ Y* G- ~</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\\"># D: c1 G' t! d1 i& W3 e</span>   k++<br /> <span style=\\"display:none\\">( n$ p% P* R- V0 n+ R5 S</span> },<br /> <span style=\\"display:none\\">5 {3 F) D* V! V6 R& h\\" w\\" T: k</span> s;<font class=\\"jammer\\">. O6 k, Y' W0 x; [6 u# ]\\" c) |</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
1.031
" {2 D7 v) l4 w3 X+ Y
二者效率相差已比较小。
3 r! G( d4 l6 |7 C6 }- A& b
3、再增加一个函数调用

0 Q" O$ K/ g7 W) wMatlab代码：

1. clear all$ t4 ]$ D* X: ?$ N; ^4 l- K* {1 U
   
2. f=@(x)x+1;
   7 f/ s6 [# d' N* r. I$ T
3. tic2 m4 q& @/ O8 f9 m* I
   
4. s=0;
   7 C$ Y$ ~8 a9 T  q3 A
5. for k=1:1000
   \" ]8 w. b2 @! x3 H( m
6.   for i=1:1000/ x0 N$ F* a6 _1 J* n5 w% C# b- 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);
   1 I\" ^0 B; W: |' h$ E
8.   end' O* K; B- P6 [( P
   
9. end7 @/ X9 \* b% s
   
10. s
    8 J' V1 I1 H8 O6 C5 f1 I1 ^
11. toc: U2 M; F$ B6 ?% d, @
    
12. 8 l0 W& T+ r8 C\" ]$ O( b
    
13. s =
    \" m* r- R; X' [& M! \- T0 z
14. / G  f+ h! z! w\" s7 l+ i
    
15.   5.1251e+008, U3 Y; m; E9 v
    
16. 
    \" W2 h+ c+ P/ m5 w7 Y' }
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];& `' P2 P; L( S2 Y6 p
   
4. f(x)=x+1;
   
5. 9 I+ V& o2 ?, ^% |( t/ C
6. mvar:$ i$ X8 M2 o/ E& s* q0 \4 J1 E1 G
   
8. t0=clock(),7 b\\" o, P- K. x. A1 Z) k7 \
   
10. s=0,4 g\\" h2 k! E: F4 Y' D
    
12. k=1, while{k<=1000,4 o/ ^5 q% D; 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++},* h9 L& `2 l! C$ w* e/ X! i& O9 {& s
    
16.   k++& \& t0 x- \6 Y
    
18. },) k* i9 w1 w$ o3 A
    
20. s;& k8 Q; A0 i4 w* W* p6 B! c+ D
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">& `' P2 P; L( S2 Y6 p</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">9 I+ V& o2 ?, ^% |( t/ C</span> mvar:<font class=\\"jammer\\">$ i$ X8 M2 o/ E& s* q0 \4 J1 E1 G</font><br /> t0=clock(),<font class=\\"jammer\\">7 b\\" o, P- K. x. A1 Z) k7 \</font><br /> s=0,<font class=\\"jammer\\">4 g\\" h2 k! E: F4 Y' D</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">4 o/ ^5 q% D; F</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++},<font class=\\"jammer\\">* h9 L& `2 l! C$ w* e/ X! i& O9 {& s</font><br />   k++<font class=\\"jammer\\">& \& t0 x- \6 Y</font><br /> },<font class=\\"jammer\\">) k* i9 w1 w$ o3 A</font><br /> s;<font class=\\"jammer\\">& k8 Q; A0 i4 w* W* p6 B! c+ D</font><br /> [clock()-t0]/1000;

结果：
512508953.6658993
1.078
+ t3 e, R0 x% S: Q- D) |% U
可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
/ A" F* ~7 v4 d0 K
  E: F% t8 R% g* D! `Matlab代码：

1. >> clear all
   9 D# R) }* S. r8 B/ @
2. a = linspace(1,1,1000);
   0 ^# Z+ F) Y; V. P9 D% C6 S
3. tic
   : [% M) b4 `9 `& J, g; P
4. for k=1:100009 t5 K+ U% i$ `' V1 ~% w8 }
   
5.   for i=1:1000% q7 K2 x5 U' q4 }
   
6.     a(i)=2;- N; A; {# y& b4 H2 p/ y( }9 J
   
7.   end: _9 u& |% N& K3 F1 l\" C
   
8. end
   % T4 m) o\" }$ z
9. toc5 k! V& J9 r4 K& q
   
10. 
    . z: X( H% ~7 M2 G
11. tic) X1 K) S* v+ x9 h
    
12. s=0;* W! \* Y+ E! C6 E. b( m% D! q: ]
    
13. for k=1:10000
    8 f* P; F5 F7 q0 A
14.   for i=1:1000! z\" B# Z* e( `7 J: X
    
15.     s=s+a(i);; U5 o* _8 s. V2 V
    
16.   end& v0 j5 W7 ]' q) g4 D# w+ P- j) O
    
17. end& |0 g# w( G2 a& e
    
18. s! ~& A: D; p6 o' \( n- R% c: b4 z
    
19. toc/ G/ ~1 `( y1 X( e2 ~0 I
    
20. Elapsed time is 0.140658 seconds.
    2 f. ]2 Z6 Z% h7 N& j$ K
21. + x$ x0 n. `5 V
    
22. s =3 B4 x3 ~: K6 @
    
23. 
    ; r) j$ Z$ Z( o7 o7 W) k
24.     200000006 \' e* ^( v6 z7 \4 |
    
25. 
    \" p* o6 c\" V1 b/ P3 ~: ?. ?1 @
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];
   
3. ; b  p% e\\" P, d! G7 }4 i* d
4. mvar:
   
5. 6 g8 y* V6 ^/ c8 [* y
6. main0(::a)=
   
7. * Z9 V8 C5 r3 }' y$ c3 O2 Y
8. a = linspace(1,1,1000),
   
9. . c7 N$ {8 n6 D\\" P. P7 U
10. t0=clock(),% r4 F* L+ n# g( t8 P
    
12. k=0, while{k<10000,
    
13. : P  X\\" n9 ~8 n. Q( _\\" d0 z
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. ! k' i9 ^& {+ R& R
16.   k+++ b3 \$ i! [% E- S7 x* s' Q
    
18. },
    
19. $ Y  K8 ~, J/ Q  @$ g& W1 \
20. [clock()-t0]/1000;
    
21. 8 T' x9 V! [+ `3 f& w% r
22. 
    
23. . }& Q5 J- H& B, u* {1 j; m
24. main1(::a)=
    
25. 1 t: U\\" c4 Q/ x$ n2 A
26. t0=clock(), s=0,
    
27. - a# m% M$ @; Z0 _5 g+ X+ C6 J
28. k=0, while{k<10000,
    
29. ( r: u6 [' t8 E) G! p6 r7 W# E0 e# j
30.   i=0, while{i<1000, s=s+a(i), i++},, }* r5 F# e+ G- ~
    
32.   k++. V# y; X9 s/ F& k9 x
    
34. },1 ~9 W$ }) D  h0 [% p' \; F5 n
    
36. s;
    
37. * P9 D, e+ P( }
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">; b  p% e\\" P, d! G7 }4 i* d</span> mvar:<br /> <span style=\\"display:none\\">6 g8 y* V6 ^/ c8 [* y</span> main0(::a)=<br /> <span style=\\"display:none\\">* Z9 V8 C5 r3 }' y$ c3 O2 Y</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">. c7 N$ {8 n6 D\\" P. P7 U</span> t0=clock(),<font class=\\"jammer\\">% r4 F* L+ n# g( t8 P</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">: P  X\\" n9 ~8 n. Q( _\\" d0 z</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">! k' i9 ^& {+ R& R</span>   k++<font class=\\"jammer\\">+ b3 \$ i! [% E- S7 x* s' Q</font><br /> },<br /> <span style=\\"display:none\\">$ Y  K8 ~, J/ Q  @$ g& W1 \</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">8 T' x9 V! [+ `3 f& w% r</span> <br /> <span style=\\"display:none\\">. }& Q5 J- H& B, u* {1 j; m</span> main1(::a)=<br /> <span style=\\"display:none\\">1 t: U\\" c4 Q/ x$ n2 A</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">- a# m% M$ @; Z0 _5 g+ X+ C6 J</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">( r: u6 [' t8 E) G! p6 r7 W# E0 e# j</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">, }* r5 F# e+ G- ~</font><br />   k++<font class=\\"jammer\\">. V# y; X9 s/ F& k9 x</font><br /> },<font class=\\"jammer\\">1 ~9 W$ }) D  h0 [% p' \; F5 n</font><br /> s;<br /> <span style=\\"display:none\\">* P9 D, e+ P( }</span> [clock()-t0]/1000;

结果：
% g' \7 \$ A. u0 Q: G- L6 o  ?15.484  //存数据所用时间，秒
20000000.
( N& _- N/ n8 [5 p9 u8 |8 a; @3 r5.531   //取数据所用时间，秒

－－－－－－－－－－
2 T3 {  d* b. u/ _1 F
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. ! N7 x+ @; h2 H- ^5 Y\\" B+ p- n
4. mvar:
   
5. \\" `  t2 h4 L' E/ d. X/ y1 X+ Q
6. main0(::a)=
   
7. 8 {, V( x! ?$ x+ o$ l4 c) R* @) {
8. a = linspace(1,1,1000),
   
9. ' s: {: w0 a  P
10. t0=clock(),3 W5 N; e/ ?9 ^! ?+ g# i\\" p% D
    
12. k=0, while{k<10000,
    
13. \\" v( {: w3 s5 C7 x, Q& h
14.   i=0, while{i<1000, A(a,i)=2, i++},& q4 D) S# \! @! b/ C( Z& R\\" H
    
16.   k++
    
17. 5 _% \; i- |& g4 j2 }
18. },
    
19. ) n( b3 Y* H\\" Q, t' B
20. [clock()-t0]/1000;0 R. J6 q$ r, l
    
22. 
    
23. 8 z, a8 l7 V, Z$ L8 a
24. main1(::a)=  q6 x\\" f+ |& y, s$ Q: p% m; [
    
26. t0=clock(), s=0,  p1 M6 I  j% }* j/ \+ B\\" l
    
28. k=0, while{k<10000,
    
29. ( G2 n* O: X. }. f5 |: L8 V
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. - I% |& Q. S0 X
32.   k++3 i4 B& x) G5 t0 P
    
34. },
    
35. ; L: n9 W7 e7 m5 B8 l; ]( S
36. s;
    
37. 9 n6 T2 N) c9 j+ }\\" B6 o
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">! N7 x+ @; h2 H- ^5 Y\\" B+ p- n</span> mvar:<br /> <span style=\\"display:none\\">\\" `  t2 h4 L' E/ d. X/ y1 X+ Q</span> main0(::a)=<br /> <span style=\\"display:none\\">8 {, V( x! ?$ x+ o$ l4 c) R* @) {</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">' s: {: w0 a  P</span> t0=clock(),<font class=\\"jammer\\">3 W5 N; e/ ?9 ^! ?+ g# i\\" p% D</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">\\" v( {: w3 s5 C7 x, Q& h</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">& q4 D) S# \! @! b/ C( Z& R\\" H</font><br />   k++<br /> <span style=\\"display:none\\">5 _% \; i- |& g4 j2 }</span> },<br /> <span style=\\"display:none\\">) n( b3 Y* H\\" Q, t' B</span> [clock()-t0]/1000;<font class=\\"jammer\\">0 R. J6 q$ r, l</font><br /> <br /> <span style=\\"display:none\\">8 z, a8 l7 V, Z$ L8 a</span> main1(::a)=<font class=\\"jammer\\">  q6 x\\" f+ |& y, s$ Q: p% m; [</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">  p1 M6 I  j% }* j/ \+ B\\" l</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">( G2 n* O: X. }. f5 |: L8 V</span>   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">- I% |& Q. S0 X</span>   k++<font class=\\"jammer\\">3 i4 B& x) G5 t0 P</font><br /> },<br /> <span style=\\"display:none\\">; L: n9 W7 e7 m5 B8 l; ]( S</span> s;<br /> <span style=\\"display:none\\">9 n6 T2 N) c9 j+ }\\" B6 o</span> [clock()-t0]/1000;

结果：
10.469
20000000.
3.781

－－－－－－－

* G: H( B( R' o; w9 l! L' R在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
% [8 y8 E; i% F3 q
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
- l8 J0 e, K5 y7 G
3 G! N9 T0 l2 ]6 @3 [* |, ]Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

8 A  t) a2 U! L# W8 X, R5、在循环中加点东西

matlab代码：

1. >> clear all
   $ {; B; Q4 t* s0 f! @: C
2. f=@(x)x+1;  |4 W( D7 ]. ?\" u7 I+ h- K& \$ G3 s
   
3. a = linspace(2,2,1000);
   9 I) `\" `  b  D+ S5 _% Z9 Q
4. tic
   % j' G. |5 ~. T5 A
5. s=0;6 O+ P, E7 {* G+ Y) c
   
6. for k=1:1000  p( v) l4 m& ~* m
   
7.   for i=1:1000
   6 [. @, b8 b) h, L
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));
   ! f1 ^2 E3 F# q/ T) t) |
9.   end- s1 q* Z5 q9 ~! i0 p% ^- o
   
10. end& i3 ]9 x0 G+ {# `9 H  o  X
    
11. s* n4 Q0 ^# v' O# o5 s8 i
    
12. toc; J9 w% \3 C8 s' s  F. p8 W& k
    
13. : ^- g- V% b) U2 ]4 P# a9 f
    
14. s =
    ' P1 a  |0 _6 h# h7 C
15. 5 ^% j6 t6 G0 m( T7 L
    
16.   2.4002e+0078 J- b3 O) {8 H
    
17. 9 _4 z: o2 L. B6 ?
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];4 }\\" i\\" K\\" R$ G1 e2 P
   
4. mvar:
   
5. ! T\\" K0 |; V+ h: z
6. f(x)=x+1;
   
7.   T/ T$ ~* i; z
8. main0(::a)=- Y& q3 d, r! y0 Z
   
10. a = linspace(2,2,1000),3 l* v+ |6 n3 k% I% l0 _& ~
    
12. t0=clock(), s=0,
    
13. 9 s3 Z1 F' t: B; f- g
14. k=0, while{k<1000,  h4 E0 P$ C4 [% @/ o
    
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++},: ~+ R+ M! O; T; K; n
    
18.   k++
    
19. & ^6 a  j# ^& v' ^/ l
20. },3 _2 @. F9 r3 ]
    
22. s;, \* a. K8 g* A
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">4 }\\" i\\" K\\" R$ G1 e2 P</font><br /> mvar:<br /> <span style=\\"display:none\\">! T\\" K0 |; V+ h: z</span> f(x)=x+1;<br /> <span style=\\"display:none\\">  T/ T$ ~* i; z</span> main0(::a)=<font class=\\"jammer\\">- Y& q3 d, r! y0 Z</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">3 l* v+ |6 n3 k% I% l0 _& ~</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">9 s3 Z1 F' t: B; f- g</span> k=0, while{k<1000,<font class=\\"jammer\\">  h4 E0 P$ C4 [% @/ o</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\\">: ~+ R+ M! O; T; K; n</font><br />   k++<br /> <span style=\\"display:none\\">& ^6 a  j# ^& v' ^/ l</span> },<font class=\\"jammer\\">3 _2 @. F9 r3 ]</font><br /> s;<font class=\\"jammer\\">, \* a. K8 g* A</font><br /> [clock()-t0]/1000;

结果：
- v' v  l, a( i# H: `( }$ B24002271.69189384
5.172

3 `: i5 t# |" t) j2 @! z  N二者效率现在只有一倍差距。

－－－－－－－－－－
4 C+ u: v7 K+ x3 ]
不过，显然matlab和Forcal的程序应该这样写：
8 {' a( @9 q( E, U
matlab程序：

1. clear all& U/ R# W& S( S2 A7 v7 w: F
   
2. f=@(x)x+1;
   5 W) m' f2 T$ N' s9 {
3. a = linspace(2,2,1000);
   7 Y$ d  V0 U! I- I
4. tic
   : ]# ~2 n. t! P+ X
5. s=0;8 t2 ]$ Z5 b. j# \  b
   
6. for k=1:1000
   + O! A, ^5 Y& H\" L4 T( s
7.   for i=1:1000
   % m2 L/ r2 e& }
8.     t=a(i);8 B\" ^0 n\" o+ G( t
   
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);
   % y& {  S* N5 ?: @
10.   end
    \" h+ v* v3 y1 G( {; _
11. end6 W# l6 g- {0 M! \$ g: q
    
12. s
    2 M$ c5 n9 ?) q7 ]; ~
13. toc
    7 S6 I: }5 I6 w( s: z
14. 
    ( E; L! l5 T6 `, z
15. s =' r, Q/ ?  f; b0 G- b- V% H
    
16. 
    1 O# r) ]\" G& {6 |: X2 b: h6 v
17.   2.4002e+007! D# T- O' Q  z8 @  ^* C' l( {
    
18. 
    8 q7 }1 u! r- A+ f; W; l
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];: P! X# V# p\\" a4 {0 B9 _
   
4. mvar:5 n7 ?. b( @/ H5 p, ?3 A
   
6. f(x)=x+1;
   
7. ( f; I1 [+ R3 i  O: e5 J3 X
8. main0(::a)=! N# S, n- ?0 e) i# x2 t) c
   
10. a = linspace(2,2,1000),
    
11. & U\\" m) h. M; h
12. t0=clock(), s=0,
    
13. ! P# U% _* I- y1 V( B3 G
14. k=0, while{k<1000,
    
15. \\" K$ U\\" h5 H1 W, b  W/ X% v
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. . y% i1 T5 q$ b' V1 z+ f* s2 Z2 X* Y
18.   k++
    
19. & p) ?, ]. f3 ^  I; j- ~
20. },
    
21. ( m  }6 M9 A2 a% d
22. s;4 p, I& G5 s4 n3 H4 R0 {. Y* z
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">: P! X# V# p\\" a4 {0 B9 _</font><br /> mvar:<font class=\\"jammer\\">5 n7 ?. b( @/ H5 p, ?3 A</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">( f; I1 [+ R3 i  O: e5 J3 X</span> main0(::a)=<font class=\\"jammer\\">! N# S, n- ?0 e) i# x2 t) c</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">& U\\" m) h. M; h</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">! P# U% _* I- y1 V( B3 G</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">\\" K$ U\\" h5 H1 W, b  W/ X% v</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++},<br /> <span style=\\"display:none\\">. y% i1 T5 q$ b' V1 z+ f* s2 Z2 X* Y</span>   k++<br /> <span style=\\"display:none\\">& p) ?, ]. f3 ^  I; j- ~</span> },<br /> <span style=\\"display:none\\">( m  }6 M9 A2 a% d</span> s;<font class=\\"jammer\\">4 p, I& G5 s4 n3 H4 R0 {. Y* z</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
1.625
8 ~/ ]6 n. Z: V- l3 a/ \
matlab效率已不及Forcal。
8 R% H- s. Q% i0 c7 x& D8 }" W! m4 k/ l
提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
" H" }; S# F" H+ K2 p6 ?% k
Forcal代码：

2. !using["sys"];
   
3. . ?6 G# c; X% O; e) u) R$ a
4. mvar:% ?3 e  A+ m# ~
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. 1 Z$ p  R/ w' i5 L
8. t0=clock(),8 K5 D4 ^. R* T  g- M; o, Q  y3 G
   
10. k=0, while{k<10000,1 {  Y+ I0 X* c9 x; P
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. $ F5 B  g: Q. N: V
14.   k++
    
15. - q3 K8 [+ |# f
16. },6 Q6 d9 i4 S0 o\\" ?
    
18. [clock()-t0]/1000;
    
19. 9 H. R0 m\\" U) f9 }* S7 g
20. 7 |\\" D: z/ v4 V6 s: `+ g
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组2 b  l5 \- K0 e0 u5 D
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值7 C1 L% Z\\" g# p1 m! R% Y3 b
    
26. t0=clock(),
    
27.   _; h# o0 R& I
28. s=0,
    
29. 7 g4 X\\" {0 K+ |+ U! E: M) z
30. k=0, while{k<10000,; N0 d) b: j& d9 g& T8 d  ~
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值! T# q% S\\" b+ v) V) m1 a
    
34.   k++
    
35. \\" \+ B. K\\" I! P# O2 K3 v
36. },3 f' z4 M/ A4 X0 w, K
    
38. s;- f! C( u& o3 h
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">. ?6 G# c; X% O; e) u) R$ a</span> mvar:<font class=\\"jammer\\">% ?3 e  A+ m# ~</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">1 Z$ p  R/ w' i5 L</span> t0=clock(),<font class=\\"jammer\\">8 K5 D4 ^. R* T  g- M; o, Q  y3 G</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">1 {  Y+ I0 X* c9 x; P</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">$ F5 B  g: Q. N: V</span>   k++<br /> <span style=\\"display:none\\">- q3 K8 [+ |# f</span> },<font class=\\"jammer\\">6 Q6 d9 i4 S0 o\\" ?</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">9 H. R0 m\\" U) f9 }* S7 g</span> <font class=\\"jammer\\">7 |\\" D: z/ v4 V6 s: `+ g</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">2 b  l5 \- K0 e0 u5 D</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">7 C1 L% Z\\" g# p1 m! R% Y3 b</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">  _; h# o0 R& I</span> s=0,<br /> <span style=\\"display:none\\">7 g4 X\\" {0 K+ |+ U! E: M) z</span> k=0, while{k<10000,<font class=\\"jammer\\">; N0 d) b: j& d9 g& T8 d  ~</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">! T# q% S\\" b+ v) V) m1 a</font><br />   k++<br /> <span style=\\"display:none\\">\\" \+ B. K\\" I! P# O2 K3 v</span> },<font class=\\"jammer\\">3 f' z4 M/ A4 X0 w, K</font><br /> s;<font class=\\"jammer\\">- f! C( u& o3 h</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
20000000.
- B) K1 ?: H( }0 O$ l5 b. ^0.579      //取值用时间，秒