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

forcal

1、简单循环

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

1. clear all3 B0 L1 _$ Y& I% U5 A2 p* |3 |! L
   
2. tic3 B0 u, H- P! c. a0 C! E5 L0 W
   
3. s=0;5 ]- X0 ]* b( p* I0 b3 Z
   
4. for k=1:100003 O6 r; ~& Q* n7 F9 x
   
5.   for i=1:1000/ Z0 g! g. Y5 R0 T
   
6.     s=s+1;
   7 Q\" M. C+ C/ N6 L# j& \$ z
7.   end/ V+ i. N1 s/ K
   
8. end& s$ s; R$ J) Z; {( {% k
   
9. s( Y5 D, L/ A# l
   
10. toc+ i\" R1 r  l5 |4 B# Z* H$ k
    
11. 
    ! L0 N1 N: j0 n5 z- i, W: P0 F8 `
12. s =1 g% |& t2 y- n5 t4 Z( J9 I
    
13. 
    \" x+ Y2 T& n$ V\" g
14.     10000000
    , ~1 h( T6 C9 f3 j* _* f
15. ! p$ J6 h) x: ]5 d. q; [) ?! r
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

: d7 v% H8 c9 v; A: V& LForcal没有JIT加速器，代码：

2. !using["sys"];
   
3. ' ]% T8 |! ^\\" g/ Y) ?
4. mvar:
   
5. 9 q1 U1 a( r9 ]9 k$ f9 _
6. t0=clock(),
   
7. 0 g3 w. B9 n  @3 b
8. s=0,$ x- l3 ?- }# r2 M8 w
   
10. k=0, while{k<10000,
    
11. ) \6 x. r1 E. E# ]+ j. U3 ~9 E- C
12.   i=0, while{i<1000,s=s+1, i++},
    
13. 6 y$ v& E) x* L' j! D  o
14.   k++' w. n, D6 u5 e: z5 p% Q/ q
    
16. },
    
17. - F& \4 ~0 Q# R1 {3 Z  t' O
18. s;
    
19. $ O+ ^8 G+ r5 R# e8 l7 P
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">' ]% T8 |! ^\\" g/ Y) ?</span> mvar:<br /> <span style=\\"display:none\\">9 q1 U1 a( r9 ]9 k$ f9 _</span> t0=clock(),<br /> <span style=\\"display:none\\">0 g3 w. B9 n  @3 b</span> s=0,<font class=\\"jammer\\">$ x- l3 ?- }# r2 M8 w</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">) \6 x. r1 E. E# ]+ j. U3 ~9 E- C</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">6 y$ v& E) x* L' j! D  o</span>   k++<font class=\\"jammer\\">' w. n, D6 u5 e: z5 p% Q/ q</font><br /> },<br /> <span style=\\"display:none\\">- F& \4 ~0 Q# R1 {3 Z  t' O</span> s;<br /> <span style=\\"display:none\\">$ O+ ^8 G+ r5 R# e8 l7 P</span> [clock()-t0]/1000;

结果：
10000000.
0.704  秒

二者运行效率相差7倍多。

6 v1 v  @6 G, W; T, B$ K4 q( V2、循环中稍微加点东西
0 w5 ]5 C4 i2 K
9 i( e) G/ l9 d) `% Q) `Matlab代码：

1. clear all
   8 ~0 A( e5 E4 u8 E+ e
2. tic9 d2 }' u& ?( H- F( |
   
3. s=0;
   5 ]\" D: A& k: `& p# H# _2 X
4. for k=1:1000
   ( y: j. Y\" |' P1 u\" Y' E
5.   for i=1:10002 n1 I; A) U* H, h- X: M
   
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);7 z* X' M\" C+ V5 A
   
7.   end
   - B. T; @& @) T. s
8. end; w9 ]( Q' K7 ~( A# v2 X
   
9. s7 t9 C: v) v2 J
   
10. toc7 N) ?4 p9 }9 J$ c+ }
    
11. 6 ^% d  R- z5 ?9 t. x, T
    
12. s =$ [# w9 t$ y6 I! _
    
13. 
    - q; i+ J& L1 c' ^  E
14.   1.1009e+0075 f3 [- z6 O\" B) S1 |0 |- T; c
    
15. ) |& P& G& A  j/ O$ ]' Y
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];* n2 j8 T) J8 B, C
   
4. mvar:
   
5. # W+ u6 H$ z1 [4 @9 i1 H9 r& f
6. t0=clock(),5 I! @5 z4 ?3 e9 k4 R: P& y: [7 M
   
8. s=0,
   
9. * a7 Z! @+ s, g* }: x4 G\\" _
10. k=1, while{k<=1000,
    
11. . j! ^  H3 V\\" ~2 p$ 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++},# X2 ]+ m6 I) N
    
14.   k++
    
15. $ }& L& L1 _5 ~* ?4 s. T) p
16. },
    
17. ) {: {' F; o+ Y* n
18. s;
    
19. ( K! ~4 R& o$ [1 M
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">* n2 j8 T) J8 B, C</font><br /> mvar:<br /> <span style=\\"display:none\\"># W+ u6 H$ z1 [4 @9 i1 H9 r& f</span> t0=clock(),<font class=\\"jammer\\">5 I! @5 z4 ?3 e9 k4 R: P& y: [7 M</font><br /> s=0,<br /> <span style=\\"display:none\\">* a7 Z! @+ s, g* }: x4 G\\" _</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">. j! ^  H3 V\\" ~2 p$ 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++},<font class=\\"jammer\\"># X2 ]+ m6 I) N</font><br />   k++<br /> <span style=\\"display:none\\">$ }& L& L1 _5 ~* ?4 s. T) p</span> },<br /> <span style=\\"display:none\\">) {: {' F; o+ Y* n</span> s;<br /> <span style=\\"display:none\\">( K! ~4 R& o$ [1 M</span> [clock()-t0]/1000;

结果：
11008953.66602346
2 [# _4 R0 J0 E; D1.031

二者效率相差已比较小。

3、再增加一个函数调用

2 {2 m) [, K2 X  h3 Q5 C7 EMatlab代码：

1. clear all0 j$ [5 ^1 i& F
   
2. f=@(x)x+1;
   , d2 N2 {9 b# G
3. tic3 u: D% m: c8 ^8 G9 i
   
4. s=0;
   - m: Q/ `* B6 l, x7 ]6 [
5. for k=1:1000
   / e2 ]) h0 X. y6 {' P' y
6.   for i=1:10000 T2 _1 V$ \& @1 D  v
   
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);
   : S9 H& ]2 g1 j
8.   end
   \" t. ~. [+ u* `5 r
9. end7 x\" F; Q9 Y& A' M# B3 [8 U  q1 J
   
10. s  {  [1 A- _# x# }! d# g
    
11. toc
    / r- m; |: l9 M4 t3 A& [
12. 4 `\" T\" ~( i. Z4 h' q0 w
    
13. s =5 F! s# N- {  o3 O2 p9 g
    
14. 2 g2 N8 a$ f$ M/ ^7 l) `0 V( D) b  f% g
    
15.   5.1251e+008. [: e. V( u, q' r
    
16. , ^4 [+ M! A- a7 j/ @
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];# }  n5 @3 y% k+ F5 U5 x
   
4. f(x)=x+1;
   
5. . D7 J1 U; _6 E. k( _2 ~
6. mvar:- y2 [2 k8 Q. P$ s\\" D  U
   
8. t0=clock(),
   
9. + m) w4 Q& t3 U
10. s=0,* _4 W4 n8 w\\" u2 z; ]# p
    
12. k=1, while{k<=1000,
    
13. . E* s! z0 n- L; c3 `
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++},. o5 q: y$ H; q4 e' k/ {7 U
    
16.   k++5 p& p4 ]8 ?\\" J\\" s+ i- t. j
    
18. },  `) g1 X4 S& T. K7 b
    
20. s;! C! T9 Q\\" v7 T+ T# A
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\"># }  n5 @3 y% k+ F5 U5 x</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">. D7 J1 U; _6 E. k( _2 ~</span> mvar:<font class=\\"jammer\\">- y2 [2 k8 Q. P$ s\\" D  U</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">+ m) w4 Q& t3 U</span> s=0,<font class=\\"jammer\\">* _4 W4 n8 w\\" u2 z; ]# p</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">. E* s! z0 n- L; c3 `</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\\">. o5 q: y$ H; q4 e' k/ {7 U</font><br />   k++<font class=\\"jammer\\">5 p& p4 ]8 ?\\" J\\" s+ i- t. j</font><br /> },<font class=\\"jammer\\">  `) g1 X4 S& T. K7 b</font><br /> s;<font class=\\"jammer\\">! C! T9 Q\\" v7 T+ T# A</font><br /> [clock()-t0]/1000;

结果：
512508953.6658993
1.078

2 U+ m* N. {5 g' `' e+ A可以看出，Matlab效率已不及Forcal。
1 N8 L  @( p. g: @- n9 X故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

Matlab代码：

1. >> clear all
   ' A  v5 Q, n\" Z2 P* J
2. a = linspace(1,1,1000);
   % W' X' B1 t  h! Q9 r& B) k; _; k1 F3 _
3. tic6 p: r2 b& O4 l( j: |
   
4. for k=1:10000% Y6 c4 [8 X# n/ y0 b
   
5.   for i=1:10002 v: Y% G  z6 W) Z4 y0 ?
   
6.     a(i)=2;) w, L; I1 F/ x
   
7.   end
   6 i; y% _( r: z+ V! P+ z0 D
8. end# z8 K\" B( l# w8 h  {8 R' Z
   
9. toc
   & r+ q# t5 T: a/ ?: e. H: X8 b
10. % ?3 U) k) W) {( u$ p* W  L: L% d
    
11. tic* X7 v. W; G- _6 B) a
    
12. s=0;
    4 B\" k) Y: H5 f0 }' u
13. for k=1:100001 h0 d) C1 ?. L. q/ s( O' A
    
14.   for i=1:1000
    # k9 g2 @\" W/ `+ K5 _
15.     s=s+a(i);/ X6 g( O: s! u  h! c9 f0 V2 u
    
16.   end+ T; r4 C$ P7 p0 ^. W
    
17. end4 n. w) E' G- ^1 }( ?9 L! m
    
18. s6 T5 k7 ~6 s# f1 j( g  z
    
19. toc
    : U! C9 |0 o% k; N  M5 z
20. Elapsed time is 0.140658 seconds.
    - W2 C8 G4 P& X. `/ t- @
21.   i! j( I$ |2 y8 o3 L1 N
    
22. s =
    3 M: y* X6 l* {2 f3 ^
23. 
    ) d. X1 S\" w! U\" X+ r
24.     200000000 o) C+ N# X3 D$ I
    
25. 9 O% V# h& X; w3 i: H
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

/ ~0 f( N& R2 N+ z& cForcal代码：

2. !using["math","sys"];: E5 N9 n; e1 j& E' w4 ^
   
4. mvar:
   
5. , g( y; M/ ]! q# P7 w& E2 K
6. main0(::a)=
   
7. - a5 ~& E$ J1 N. i. i( v# f
8. a = linspace(1,1,1000),
   
9. + u' U' X2 k\\" P! \, d1 ?
10. t0=clock(),
    
11. 4 k8 }( h; J1 s7 D# i
12. k=0, while{k<10000,+ k+ I* Q+ ?/ H# v
    
14.   i=0, while{i<1000, a(i)=2, i++},5 M- ~& l$ i3 D1 q\\" _
    
16.   k++
    
17. 8 x* U. {5 ^' X9 z- s. N
18. },
    
19. / \' s+ |3 s+ F+ g/ B0 y% i% i
20. [clock()-t0]/1000;
    
21. - ]; N0 H4 [2 H7 X7 {; L% L
22. 
    
23. , G  u; i  z4 z/ ~
24. main1(::a)=
    
25. 6 w$ z$ n) a+ K1 Z
26. t0=clock(), s=0,# G, l4 W; U- ]( N
    
28. k=0, while{k<10000,
    
29. 6 Y! @- L# C4 b$ r0 b8 a
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. . l9 t! h# G+ w- x) C$ h# x; e9 J
32.   k++
    
33. ) a( P( g9 F# a/ X; o2 A2 w. h6 [: V
34. },
    
35. 5 y& `8 j8 Q' n7 E/ T
36. s;$ c4 b' Y3 l; h1 }  A% N+ N
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">: E5 N9 n; e1 j& E' w4 ^</font><br /> mvar:<br /> <span style=\\"display:none\\">, g( y; M/ ]! q# P7 w& E2 K</span> main0(::a)=<br /> <span style=\\"display:none\\">- a5 ~& E$ J1 N. i. i( v# f</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">+ u' U' X2 k\\" P! \, d1 ?</span> t0=clock(),<br /> <span style=\\"display:none\\">4 k8 }( h; J1 s7 D# i</span> k=0, while{k<10000,<font class=\\"jammer\\">+ k+ I* Q+ ?/ H# v</font><br />   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">5 M- ~& l$ i3 D1 q\\" _</font><br />   k++<br /> <span style=\\"display:none\\">8 x* U. {5 ^' X9 z- s. N</span> },<br /> <span style=\\"display:none\\">/ \' s+ |3 s+ F+ g/ B0 y% i% i</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">- ]; N0 H4 [2 H7 X7 {; L% L</span> <br /> <span style=\\"display:none\\">, G  u; i  z4 z/ ~</span> main1(::a)=<br /> <span style=\\"display:none\\">6 w$ z$ n) a+ K1 Z</span> t0=clock(), s=0,<font class=\\"jammer\\"># G, l4 W; U- ]( N</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">6 Y! @- L# C4 b$ r0 b8 a</span>   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">. l9 t! h# G+ w- x) C$ h# x; e9 J</span>   k++<br /> <span style=\\"display:none\\">) a( P( g9 F# a/ X; o2 A2 w. h6 [: V</span> },<br /> <span style=\\"display:none\\">5 y& `8 j8 Q' n7 E/ T</span> s;<font class=\\"jammer\\">$ c4 b' Y3 l; h1 }  A% N+ N</font><br /> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
20000000.
2 W3 r( G4 r( r% L$ r6 ~+ @+ ~5.531   //取数据所用时间，秒

" S" x" S2 L8 q# F, W! y－－－－－－－－－－
/ h/ d( r2 p1 r( g$ ?$ v; }8 n
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];# [( S, _; t& v- ~8 W
   
4. mvar:7 C6 @- i+ M! p\\" A- b/ y0 F
   
6. main0(::a)=
   
7. 4 }3 h\\" ~% u( z, Q* u6 a
8. a = linspace(1,1,1000),
   
9. : E( O! C2 h\\" k& H5 Y, }
10. t0=clock(),
    
11. - C  u( O8 o6 b
12. k=0, while{k<10000,3 y2 \* u7 h* p2 a. b9 K\\" T  L8 T
    
14.   i=0, while{i<1000, A(a,i)=2, i++},  M% f7 u: v& W0 [6 O, g) L- n
    
16.   k++$ q3 ]1 `3 \( y+ A
    
18. },
    
19. # N; @$ H9 J8 [
20. [clock()-t0]/1000;* v3 ]0 P2 H+ l: b/ @: c
    
22. 2 s& L+ w2 o! \# `) f( F
    
24. main1(::a)=7 F, N2 r! i2 o' l6 q7 p
    
26. t0=clock(), s=0,
    
27. 4 X! j: a- \& J' \
28. k=0, while{k<10000,
    
29. ( `. z: C3 {) c
30.   i=0, while{i<1000, s=s+A(a,i), i++},! g9 }- J% X) P9 Z5 D
    
32.   k++
    
33. 8 w( I( u: g3 m
34. },
    
35. * s# ?\\" ]5 g2 ^3 {
36. s;
    
37. * V7 v' c: g0 u1 ~% |
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\"># [( S, _; t& v- ~8 W</font><br /> mvar:<font class=\\"jammer\\">7 C6 @- i+ M! p\\" A- b/ y0 F</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">4 }3 h\\" ~% u( z, Q* u6 a</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">: E( O! C2 h\\" k& H5 Y, }</span> t0=clock(),<br /> <span style=\\"display:none\\">- C  u( O8 o6 b</span> k=0, while{k<10000,<font class=\\"jammer\\">3 y2 \* u7 h* p2 a. b9 K\\" T  L8 T</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">  M% f7 u: v& W0 [6 O, g) L- n</font><br />   k++<font class=\\"jammer\\">$ q3 ]1 `3 \( y+ A</font><br /> },<br /> <span style=\\"display:none\\"># N; @$ H9 J8 [</span> [clock()-t0]/1000;<font class=\\"jammer\\">* v3 ]0 P2 H+ l: b/ @: c</font><br /> <font class=\\"jammer\\">2 s& L+ w2 o! \# `) f( F</font><br /> main1(::a)=<font class=\\"jammer\\">7 F, N2 r! i2 o' l6 q7 p</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">4 X! j: a- \& J' \</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">( `. z: C3 {) c</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">! g9 }- J% X) P9 Z5 D</font><br />   k++<br /> <span style=\\"display:none\\">8 w( I( u: g3 m</span> },<br /> <span style=\\"display:none\\">* s# ?\\" ]5 g2 ^3 {</span> s;<br /> <span style=\\"display:none\\">* V7 v' c: g0 u1 ~% |</span> [clock()-t0]/1000;

结果：
10.469
20000000.
3.781

－－－－－－－

在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
* C8 e/ m8 _* N& m在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

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

Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
' c2 z% \3 a6 t+ ~
5、在循环中加点东西

matlab代码：

1. >> clear all3 n8 x# B* [* y7 t
   
2. f=@(x)x+1;
   : ~$ J# V0 u8 M; @0 N0 Q1 M
3. a = linspace(2,2,1000);; ^) ~- h  E2 ~\" Z- R6 C\" S( _% Q
   
4. tic. j2 b\" C7 E# ]' ]
   
5. s=0;3 U6 o- b. i4 W# H
   
6. for k=1:1000: d$ O0 Q6 N\" |\" Q, Q0 g0 V& p+ q
   
7.   for i=1:1000
   , m& H* u1 T# ?
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));
   / q, Y4 z4 z+ f. |5 Q7 ^
9.   end; Z: g: ?2 l  z* ]\" S
   
10. end
    - H7 e7 V) u. ^  |; T% A9 K
11. s
      l* z3 s1 ]\" P
12. toc  q# Z% i4 E  N- n) Z% H
    
13. # W: M% h1 O. m( [8 f
    
14. s =. c  g2 M: m) \' j
    
15. 
    6 B+ F7 b; J* f
16.   2.4002e+007
    ) c0 M- ~! k+ R0 k3 N& B9 \
17. . ^\" z* n: n% Y% ]9 R6 G# n# O
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];3 |2 j\\" W- D9 X7 j, i+ B
   
4. mvar:) c6 n# @7 M/ Y) v( A! C
   
6. f(x)=x+1;
   
7. + |  \\\" ?. G& T7 ?5 r  E2 e
8. main0(::a)=
   
9. - {6 M% W, u- z- T\\" y
10. a = linspace(2,2,1000),: }, @! Q. ]' a5 E. T( w2 i# t( O6 m
    
12. t0=clock(), s=0,
    
13. 9 M* l* z2 e0 ~% o& ?: m7 g
14. k=0, while{k<1000,& \  @7 {6 m& |4 B6 o) U4 T
    
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. ! F+ M- B5 O0 {! F* k8 G: J
18.   k++
    
19. , s+ m1 @6 t: b* \% m. P7 u  J
20. },
    
21. ) d\\" P6 H2 a+ o  c! C0 R
22. s;
    
23. % a8 n# c  k) B6 ]0 y  o5 @
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">3 |2 j\\" W- D9 X7 j, i+ B</font><br /> mvar:<font class=\\"jammer\\">) c6 n# @7 M/ Y) v( A! C</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">+ |  \\\" ?. G& T7 ?5 r  E2 e</span> main0(::a)=<br /> <span style=\\"display:none\\">- {6 M% W, u- z- T\\" y</span> a = linspace(2,2,1000),<font class=\\"jammer\\">: }, @! Q. ]' a5 E. T( w2 i# t( O6 m</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">9 M* l* z2 e0 ~% o& ?: m7 g</span> k=0, while{k<1000,<font class=\\"jammer\\">& \  @7 {6 m& |4 B6 o) U4 T</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\\">! F+ M- B5 O0 {! F* k8 G: J</span>   k++<br /> <span style=\\"display:none\\">, s+ m1 @6 t: b* \% m. P7 u  J</span> },<br /> <span style=\\"display:none\\">) d\\" P6 H2 a+ o  c! C0 R</span> s;<br /> <span style=\\"display:none\\">% a8 n# c  k) B6 ]0 y  o5 @</span> [clock()-t0]/1000;

结果：
24002271.69189384
5.172

$ f! l5 ]# U: }: y二者效率现在只有一倍差距。
, L( R; n+ k4 R9 o+ \" N
3 }# z; m7 H0 \( u& l－－－－－－－－－－
- T( Y' G' f8 o8 f6 [+ ], I
. ^* r5 n2 Q$ l& P9 x不过，显然matlab和Forcal的程序应该这样写：

& m$ V% f2 M; Xmatlab程序：

1. clear all
   2 v' i3 I' W7 W$ f& ~6 A3 P. u
2. f=@(x)x+1;! I, O. E4 N8 Q9 N, n6 Q# [
   
3. a = linspace(2,2,1000);
   3 Z; y. m# n/ Q\" W. g$ |8 J# C) ~4 y
4. tic
   0 i0 `\" H; s; S5 x& \- Y
5. s=0;
   3 Y5 O4 n7 B  U  i
6. for k=1:1000
   & E9 W1 {9 @# i# ]\" S5 p
7.   for i=1:1000
   ) \( R9 E7 F6 ~( g/ o: c9 y/ @
8.     t=a(i);; P\" }* s3 d\" f) S\" N/ N# \
   
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);
   3 j8 f\" k/ U5 b. b' m) s: ]/ x. x
10.   end
    8 c( \! \9 q) N$ j
11. end8 i  w+ X\" t* G8 y9 \& M! b, `2 S' M8 T
    
12. s
    ( k1 n; B* ~6 v9 `' h$ q0 g
13. toc  A4 X. l) w, ]
    
14. 
    \" B7 g9 }7 o2 m( U
15. s =* k. `6 L6 b, X. p8 {
    
16. 
    - _5 G( K8 G& j8 J
17.   2.4002e+007
    ( B' D. ?8 Y7 u' T* b3 S/ v8 }$ B
18. * N\" n# d9 X4 Y4 ]\" H
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. - v$ v; N1 p9 M% O6 y( H
4. mvar:4 Z\\" }+ F( p4 X$ ^' R
   
6. f(x)=x+1;
   
7. * ]9 `8 c+ H7 T; N
8. main0(::a)=, q8 C4 L+ j- Z\\" [' U, V6 r! M
   
10. a = linspace(2,2,1000),6 Q; G  {( o7 b1 P. H: ]0 N! O
    
12. t0=clock(), s=0,+ k* M\\" r$ i3 e8 X! a; V$ i/ P
    
14. k=0, while{k<1000,
    
15. 2 a% x/ H/ `) Y: {. G9 K
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++},8 n( Z  G' A& a0 _
    
18.   k++: c. }2 q  V% {  E/ G
    
20. },! m( J( V6 l0 ~' Y# K, E  C# |
    
22. s;
    
23. 2 L7 R, o: r: G( e0 S4 m4 Z
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">- v$ v; N1 p9 M% O6 y( H</span> mvar:<font class=\\"jammer\\">4 Z\\" }+ F( p4 X$ ^' R</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">* ]9 `8 c+ H7 T; N</span> main0(::a)=<font class=\\"jammer\\">, q8 C4 L+ j- Z\\" [' U, V6 r! M</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">6 Q; G  {( o7 b1 P. H: ]0 N! O</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">+ k* M\\" r$ i3 e8 X! a; V$ i/ P</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">2 a% x/ H/ `) Y: {. G9 K</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\\">8 n( Z  G' A& a0 _</font><br />   k++<font class=\\"jammer\\">: c. }2 q  V% {  E/ G</font><br /> },<font class=\\"jammer\\">! m( J( V6 l0 ~' Y# K, E  C# |</font><br /> s;<br /> <span style=\\"display:none\\">2 L7 R, o: r: G( e0 S4 m4 Z</span> [clock()-t0]/1000;

结果：
8 ~$ s, b! P9 W. I1 M* N24002271.69189384
1.625

matlab效率已不及Forcal。

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

forcal

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

Forcal代码：

2. !using["sys"];
   
3. ) [9 {0 O/ m( E' s7 F& h; A
4. mvar:/ Z' K  Z5 X1 y
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. 6 U: p4 Y+ n0 F; Z
8. t0=clock(),
   
9. ! ?1 c1 I8 k: [
10. k=0, while{k<10000,
    
11. . m& b. k, M% d8 |# t
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. 4 T2 M: O: `2 r* o, B3 m( ~% @* G  \
14.   k++
    
15. 1 f) U; r6 u* T
16. },. i* _4 O1 ~. t) T9 x! I+ ?. i
    
18. [clock()-t0]/1000;
    
19. 9 b; X% e& U  n! e+ G2 o' W) f
20. 
    
21. * J. E- B9 s. K4 e3 Y
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. 7 g3 C+ A* q4 {% C: O- I
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. 3 R9 @8 g1 M! d; ^
26. t0=clock(),
    
27. 1 H% c$ A0 ?1 d1 |
28. s=0,
    
29. 0 \6 b! U# U( w9 |\\" D1 I
30. k=0, while{k<10000,
    
31. ) K  ^& P0 B. L
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. & s& c5 @) v- d
34.   k++) n2 m% I2 a* _0 [- X8 M5 H
    
36. },7 Y! _) o! x6 _6 w* m/ o
    
38. s;6 U! R; U$ L3 t
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">) [9 {0 O/ m( E' s7 F& h; A</span> mvar:<font class=\\"jammer\\">/ Z' K  Z5 X1 y</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">6 U: p4 Y+ n0 F; Z</span> t0=clock(),<br /> <span style=\\"display:none\\">! ?1 c1 I8 k: [</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">. m& b. k, M% d8 |# t</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">4 T2 M: O: `2 r* o, B3 m( ~% @* G  \</span>   k++<br /> <span style=\\"display:none\\">1 f) U; r6 u* T</span> },<font class=\\"jammer\\">. i* _4 O1 ~. t) T9 x! I+ ?. i</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">9 b; X% e& U  n! e+ G2 o' W) f</span> <br /> <span style=\\"display:none\\">* J. E- B9 s. K4 e3 Y</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">7 g3 C+ A* q4 {% C: O- I</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">3 R9 @8 g1 M! d; ^</span> t0=clock(),<br /> <span style=\\"display:none\\">1 H% c$ A0 ?1 d1 |</span> s=0,<br /> <span style=\\"display:none\\">0 \6 b! U# U( w9 |\\" D1 I</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">) K  ^& P0 B. L</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">& s& c5 @) v- d</span>   k++<font class=\\"jammer\\">) n2 m% I2 a* _0 [- X8 M5 H</font><br /> },<font class=\\"jammer\\">7 Y! _) o! x6 _6 w* m/ o</font><br /> s;<font class=\\"jammer\\">6 U! R; U$ L3 t</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
20000000.
9 _( y- i0 D7 D0.579      //取值用时间，秒