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

forcal

1、简单循环
9 ~! z5 r4 P2 o- D# u+ ~
3 `1 h) M0 G3 j% K, hMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   \" o, \) N8 @7 L) b* c3 m' @
2. tic9 q7 }7 U2 W, ]/ d9 |- h3 S
   
3. s=0;$ j- N' @& a/ a, a' b) n* ]+ J3 G* n
   
4. for k=1:100009 r# u5 c# F; ]+ e# D$ R9 N+ n/ A2 {
   
5.   for i=1:1000
   * l, u  F$ L6 `5 d8 A
6.     s=s+1;
   5 t4 b! J& b/ r) w, t) A
7.   end
   4 Z' m! c# e2 x& ^\" b9 }
8. end- M# ?% Y4 w0 Y! _. d3 U3 f5 \' m5 B
   
9. s
   * [. z8 ]\" ~: ?) x3 X  B% Y% n
10. toc
    % t: ~. X5 w; G, U
11. # R5 K6 m  s7 m  s, ?
    
12. s =
    , s, B7 U4 r0 O( k
13. 
    4 m/ Y6 S, P* C
14.     10000000
    , z: }% j* i) X* {5 f+ W; w
15. : \- v! B\" I; Z* U, \3 T1 b
    
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
* K. G3 [" ~* \" }, k" a3 x
8 L2 H$ w. X4 z  s/ hForcal没有JIT加速器，代码：

2. !using["sys"];) a0 E( o4 t7 R
   
4. mvar:( ]8 Y: e* N8 h0 G  B6 i
   
6. t0=clock(),3 `( k5 z* q- x: M9 H
   
8. s=0,7 G  f. Q/ D2 Q- H2 ]/ ~6 ^% m
   
10. k=0, while{k<10000,
    
11. / p! E* ?1 t' Q2 A# f! S1 J
12.   i=0, while{i<1000,s=s+1, i++},
    
13. : T1 h! E( h+ j\\" u- A3 C
14.   k++
    
15. - N$ u8 C+ k8 a& s\\" R  s1 U4 h+ [
16. },
    
17. + p$ T! j! N) z3 l% J+ x
18. s;
    
19. , ~. J- ?  z4 n3 R0 y8 V  v; z
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">) a0 E( o4 t7 R</font><br /> mvar:<font class=\\"jammer\\">( ]8 Y: e* N8 h0 G  B6 i</font><br /> t0=clock(),<font class=\\"jammer\\">3 `( k5 z* q- x: M9 H</font><br /> s=0,<font class=\\"jammer\\">7 G  f. Q/ D2 Q- H2 ]/ ~6 ^% m</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">/ p! E* ?1 t' Q2 A# f! S1 J</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">: T1 h! E( h+ j\\" u- A3 C</span>   k++<br /> <span style=\\"display:none\\">- N$ u8 C+ k8 a& s\\" R  s1 U4 h+ [</span> },<br /> <span style=\\"display:none\\">+ p$ T! j! N) z3 l% J+ x</span> s;<br /> <span style=\\"display:none\\">, ~. J- ?  z4 n3 R0 y8 V  v; z</span> [clock()-t0]/1000;

结果：
: L( V9 D6 {9 }! @0 C0 y6 l10000000.
! e2 M) G6 y  p# t4 B9 \0.704  秒
; q% c/ T+ k( x
二者运行效率相差7倍多。
! t- m6 D7 Q! Z2 v3 V! M$ ~3 W
- J6 i1 ?: R( \6 t7 f2、循环中稍微加点东西
+ N8 A5 G; U* H- s" h* }+ a
% ~( R. g2 }, q6 _6 C5 \Matlab代码：

1. clear all
   1 v# k7 B. E& ]8 b5 j7 z, [
2. tic4 `/ |2 w5 [: y1 P' `6 u
   
3. s=0;
   6 e5 K) z0 C+ ?. F8 K
4. for k=1:1000
   % Q1 `% @: _+ H/ I
5.   for i=1:1000. Z8 ~. p' b6 G6 i. ?, b& {8 J0 f% _
   
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);4 t\" p: N; C! i\" P7 U+ n1 B. G, d
   
7.   end
   % Z1 y: v' K5 _* b) p: @
8. end5 B/ ]9 z% t% E7 y
   
9. s
   ' x/ ^7 N$ ?- v
10. toc
    : f5 r1 {) k! g5 J! X7 I$ G
11. 9 E! e: U5 K3 ^) [+ w! i* P& f
    
12. s =) X8 ?+ D\" \. ]% c, R, M
    
13. 
    2 V/ @5 M/ b7 c3 b/ b( [$ d
14.   1.1009e+007, E: s+ e! c) K; Y3 p$ M$ t: l& o
    
15. 
    * \6 d4 s4 A+ R' Y- }6 U/ P/ I; H6 u
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];+ ?0 a1 `2 ^9 ^\\" q0 t# F2 R
   
4. mvar:\\" U; E! V; ]4 Q1 X\\" L7 [) g
   
6. t0=clock(),
   
7. ! r! F: F* P- i, @9 c; a( U
8. s=0,, x* R/ h/ ~7 _1 Y$ v- x  |
   
10. k=1, while{k<=1000,, ^6 w5 r% d: {5 M2 W
    
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++},! P. g5 I6 q( {+ I
    
14.   k++
    
15. 5 Z# d, d& r\\" G1 h0 a
16. },
    
17. ( C9 C\\" j  A) h: V
18. s;
    
19. 2 m& A4 F) Q1 Z
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">+ ?0 a1 `2 ^9 ^\\" q0 t# F2 R</font><br /> mvar:<font class=\\"jammer\\">\\" U; E! V; ]4 Q1 X\\" L7 [) g</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">! r! F: F* P- i, @9 c; a( U</span> s=0,<font class=\\"jammer\\">, x* R/ h/ ~7 _1 Y$ v- x  |</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">, ^6 w5 r% d: {5 M2 W</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++},<font class=\\"jammer\\">! P. g5 I6 q( {+ I</font><br />   k++<br /> <span style=\\"display:none\\">5 Z# d, d& r\\" G1 h0 a</span> },<br /> <span style=\\"display:none\\">( C9 C\\" j  A) h: V</span> s;<br /> <span style=\\"display:none\\">2 m& A4 F) Q1 Z</span> [clock()-t0]/1000;

结果：
11008953.66602346
* s6 N0 q- h5 Z+ `- n% k$ R, K1.031
6 t- f: `! L( S: d7 q
2 @8 Q6 a* j  F二者效率相差已比较小。

3、再增加一个函数调用

( |$ u$ p! d/ LMatlab代码：

1. clear all5 |, t# q: A( Y8 K\" y' Y5 V
   
2. f=@(x)x+1;6 w6 S7 a3 x; J\" E' G% E; M; g
   
3. tic
   * {- p( @- D* F/ r/ Z
4. s=0;1 B\" Z4 f, p) Z& g- m
   
5. for k=1:1000
   , b1 s: v! s$ ^3 k& `# w
6.   for i=1:1000' e. S- d; o& j& z3 G1 i7 T, A, u! Q
   
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);2 T, A\" n( o4 I\" _6 x6 W6 K
   
8.   end4 }3 C! b. q' p; m, h+ R8 I0 e1 w0 O
   
9. end' e7 k1 {7 X' \2 m
   
10. s: x5 G\" U5 {0 Y' p2 n/ y4 f
    
11. toc
    0 _6 _9 _3 ?/ v; @2 x9 n
12. 3 S5 q+ O$ N6 r4 o0 a( q* S
    
13. s =3 z% g. u* t/ l; ?
    
14. - }3 m' f\" q% r, l6 K2 g7 j) U
    
15.   5.1251e+008
    8 w' `* W5 h3 s5 @) }9 c* R( V* ^0 L
16. 
    3 M+ W; G/ R' z
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 3 v5 m$ H% k8 T/ z0 X
4. f(x)=x+1;( `4 C0 ]8 {7 O1 b5 S
   
6. mvar:
   
7. - x8 M& _, d! o0 c, h; _
8. t0=clock(),) W& ~1 w- y1 e' v, N& h2 z
   
10. s=0,
    
11. 3 I+ F, B; R% R# g$ D( t, @3 s\\" D
12. k=1, while{k<=1000,
    
13. $ T/ j, O: S% G& y! `8 S9 e+ ]
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++},5 `9 L, Y3 H$ \
    
16.   k++' M5 h5 y( E9 Z0 u- I) U1 E( E
    
18. },
    
19. ! l7 T/ n* [& [; g$ R6 n
20. s;
    
21. 2 h6 M+ l4 U0 Z5 [/ U
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">3 v5 m$ H% k8 T/ z0 X</span> f(x)=x+1;<font class=\\"jammer\\">( `4 C0 ]8 {7 O1 b5 S</font><br /> mvar:<br /> <span style=\\"display:none\\">- x8 M& _, d! o0 c, h; _</span> t0=clock(),<font class=\\"jammer\\">) W& ~1 w- y1 e' v, N& h2 z</font><br /> s=0,<br /> <span style=\\"display:none\\">3 I+ F, B; R% R# g$ D( t, @3 s\\" D</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">$ T/ j, O: S% G& y! `8 S9 e+ ]</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\\">5 `9 L, Y3 H$ \</font><br />   k++<font class=\\"jammer\\">' M5 h5 y( E9 Z0 u- I) U1 E( E</font><br /> },<br /> <span style=\\"display:none\\">! l7 T/ n* [& [; g$ R6 n</span> s;<br /> <span style=\\"display:none\\">2 h6 M+ l4 U0 Z5 [/ U</span> [clock()-t0]/1000;

结果：
2 [1 R, @) K$ `512508953.6658993
- Y# i+ {+ _) }0 j& C2 P9 ^; @0 {1.078

可以看出，Matlab效率已不及Forcal。
2 m9 D( X9 z8 F1 ~3 M% C故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
# S) {8 Z$ ^$ z* X' x+ v' A. v
) V9 x7 ?9 o2 ?Matlab代码：

1. >> clear all  z+ e+ q( ]' x7 W9 }
   
2. a = linspace(1,1,1000);- e4 b- l( T7 E\" r& g  @
   
3. tic$ ~4 V* T; l& w% f# Z$ E2 b: H
   
4. for k=1:10000\" b$ c/ d' D$ v* g& L
   
5.   for i=1:1000
   ; E' V\" v7 b+ L- D* }. T3 U) Y- h) z
6.     a(i)=2;
   $ s+ O) F+ h0 f5 A+ w
7.   end5 _6 t: y- _; c5 Y9 q- z: D/ f
   
8. end& t8 M- U\" c3 A\" \( [: g, S  ^
   
9. toc6 V2 h/ N7 f; U# E/ |
   
10. $ ^1 O7 l  b7 `
    
11. tic, Q. c4 q\" T3 Y/ K+ B
    
12. s=0;  s( r# V; Z. t: x
    
13. for k=1:10000* k  P: ^6 k! ?2 r! ?
    
14.   for i=1:1000- h- j/ u9 K; F) u, C4 {5 g( [6 ~
    
15.     s=s+a(i);\" d; `) M4 j# z( y% ^
    
16.   end
      r3 c# A5 e! j/ J
17. end
    3 T$ h: d- z\" I5 M% r
18. s
    % O: \: I8 t% W! a5 f6 W  J
19. toc
    \" }\" _3 v# r! S4 A9 E: I# g
20. Elapsed time is 0.140658 seconds.) R& h* j+ G6 {% n; t1 R
    
21. ' n2 J6 f  H: K
    
22. s =
    $ t) a3 C$ x  E* b( `; A+ Q( _
23. 
    ! L8 E; g! b' _( t
24.     20000000' L: v7 e7 l- i( P
    
25. 
    3 Z$ _4 a! k7 B
26. Elapsed time is 0.091794 seconds.

复制代码

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

  O7 ~/ S7 O0 `' n+ XForcal代码：

2. !using["math","sys"];
   
3. 7 [7 t8 r0 [2 y6 D
4. mvar:
   
5. ; X3 _# E! g) [: H4 u8 [
6. main0(::a)=\\" k# C, ]* _4 j
   
8. a = linspace(1,1,1000),- t9 [4 r' Y6 M8 \% r
   
10. t0=clock(),
    
11. 4 b# n2 S# E, A) D- x4 J
12. k=0, while{k<10000,) b& }/ q: E, U6 e$ i0 F
    
14.   i=0, while{i<1000, a(i)=2, i++},* b+ Y1 _. q9 H% w+ p$ @' |
    
16.   k++
    
17. , N3 s7 [, q) w6 s
18. },3 z. f. w4 q1 w9 P! c( L
    
20. [clock()-t0]/1000;
    
21. 8 w\\" j1 e& i: T) j- T' H- X% I
22. * h' R, |/ C4 d, h4 L  }9 x
    
24. main1(::a)=
    
25. 9 D2 r2 u5 o) ?
26. t0=clock(), s=0,/ z5 w) E  [: a$ ]* S! s7 ^
    
28. k=0, while{k<10000,3 j7 z# o$ g6 A5 |; t
    
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. * f4 G4 l- `( [
32.   k++
    
33. + X1 f- m( P1 C$ K( k+ [7 c; Q. v7 Z! Y9 C
34. },
    
35. $ r, H# Q0 T; D/ y
36. s;
    
37. + r0 v% d/ `3 j. x5 P; J. [  B
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">7 [7 t8 r0 [2 y6 D</span> mvar:<br /> <span style=\\"display:none\\">; X3 _# E! g) [: H4 u8 [</span> main0(::a)=<font class=\\"jammer\\">\\" k# C, ]* _4 j</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">- t9 [4 r' Y6 M8 \% r</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">4 b# n2 S# E, A) D- x4 J</span> k=0, while{k<10000,<font class=\\"jammer\\">) b& }/ q: E, U6 e$ i0 F</font><br />   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">* b+ Y1 _. q9 H% w+ p$ @' |</font><br />   k++<br /> <span style=\\"display:none\\">, N3 s7 [, q) w6 s</span> },<font class=\\"jammer\\">3 z. f. w4 q1 w9 P! c( L</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">8 w\\" j1 e& i: T) j- T' H- X% I</span> <font class=\\"jammer\\">* h' R, |/ C4 d, h4 L  }9 x</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">9 D2 r2 u5 o) ?</span> t0=clock(), s=0,<font class=\\"jammer\\">/ z5 w) E  [: a$ ]* S! s7 ^</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">3 j7 z# o$ g6 A5 |; t</font><br />   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">* f4 G4 l- `( [</span>   k++<br /> <span style=\\"display:none\\">+ X1 f- m( P1 C$ K( k+ [7 c; Q. v7 Z! Y9 C</span> },<br /> <span style=\\"display:none\\">$ r, H# Q0 T; D/ y</span> s;<br /> <span style=\\"display:none\\">+ r0 v% d/ `3 j. x5 P; J. [  B</span> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒

－－－－－－－－－－

" j+ R, G# N6 N8 S2 h4 E用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. ' ~% s2 c$ w! q& R
4. mvar:
   
5. $ o. B$ M! f) t\\" T0 x+ H2 Y3 Z
6. main0(::a)=
   
7. ! ?) N; ]7 x7 q' B! r. s3 V. y
8. a = linspace(1,1,1000),
   
9. ( t% Q, ^) W$ [2 ?4 `4 Q  X( V5 f& O
10. t0=clock(),+ T( V3 f. o+ ?# R# |, T+ ]3 Q8 j
    
12. k=0, while{k<10000,( C. u4 _2 A& I% G$ N
    
14.   i=0, while{i<1000, A(a,i)=2, i++},, d! }  h5 ^! |7 t\\" e# Y
    
16.   k++
    
17. 9 F/ |1 s! ]; S) v# N# I6 Y9 L) H
18. },
    
19. 5 N  A+ J  C  D+ J* _: ]5 L' t5 S
20. [clock()-t0]/1000;
    
21. \\" Z\\" Q0 w: ?2 a0 N\\" Y
22. ; M: v5 v4 h1 f
    
24. main1(::a)=' F& l) x! ?' }8 k) D
    
26. t0=clock(), s=0,
    
27. 9 y3 k0 K* d; l# O( H
28. k=0, while{k<10000,. }5 _' C6 R. ^; }2 x6 k8 h
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. ) i2 X+ {( L4 S0 m0 r) T! ?! A
32.   k++
    
33. ; n- N& m, q3 }) S; a/ k. a* E
34. },7 V4 k7 x7 N) z6 @
    
36. s;; ^+ @& t! {4 b7 y# l
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">' ~% s2 c$ w! q& R</span> mvar:<br /> <span style=\\"display:none\\">$ o. B$ M! f) t\\" T0 x+ H2 Y3 Z</span> main0(::a)=<br /> <span style=\\"display:none\\">! ?) N; ]7 x7 q' B! r. s3 V. y</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">( t% Q, ^) W$ [2 ?4 `4 Q  X( V5 f& O</span> t0=clock(),<font class=\\"jammer\\">+ T( V3 f. o+ ?# R# |, T+ ]3 Q8 j</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">( C. u4 _2 A& I% G$ N</font><br />   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">, d! }  h5 ^! |7 t\\" e# Y</font><br />   k++<br /> <span style=\\"display:none\\">9 F/ |1 s! ]; S) v# N# I6 Y9 L) H</span> },<br /> <span style=\\"display:none\\">5 N  A+ J  C  D+ J* _: ]5 L' t5 S</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">\\" Z\\" Q0 w: ?2 a0 N\\" Y</span> <font class=\\"jammer\\">; M: v5 v4 h1 f</font><br /> main1(::a)=<font class=\\"jammer\\">' F& l) x! ?' }8 k) D</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">9 y3 k0 K* d; l# O( H</span> k=0, while{k<10000,<font class=\\"jammer\\">. }5 _' C6 R. ^; }2 x6 k8 h</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">) i2 X+ {( L4 S0 m0 r) T! ?! A</span>   k++<br /> <span style=\\"display:none\\">; n- N& m, q3 }) S; a/ k. a* E</span> },<font class=\\"jammer\\">7 V4 k7 x7 N) z6 @</font><br /> s;<font class=\\"jammer\\">; ^+ @& t! {4 b7 y# l</font><br /> [clock()-t0]/1000;

结果：
) f, Y/ ^+ c  S- a% T6 w10.469
2 p) j/ u) l0 {1 }' p- S20000000.
$ a% F3 _$ j; i* i; A3.781

－－－－－－－

在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
. n5 G) ?: H( I9 K! C8 V" A在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
& I" s6 e( Q, o+ r5 p1 b# G! G) p9 Y" q
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

4 [' o5 @, z+ A; c+ S: jForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
, v( r5 ]! o5 \9 ?( b% e
  V, D1 k& `3 R+ A( ^5、在循环中加点东西

matlab代码：

1. >> clear all
   ; _$ P) T- T5 Q7 V* m
2. f=@(x)x+1;
   5 j, O! i. u, H7 o\" P
3. a = linspace(2,2,1000);8 F) @+ O  @/ q1 M: s: z
   
4. tic
   % y  M5 F( v% d# J' }& K5 K
5. s=0;1 Z, a4 o4 s4 c5 W
   
6. for k=1:1000
   + x* U4 t2 A9 V$ s5 c: n9 o' c
7.   for i=1:1000
   1 Z1 o  m, R$ W# `0 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));9 z; [; _7 e. x
   
9.   end
   ' [2 I\" H: y5 k$ x( X  u4 ]% |
10. end$ @5 D% |& b3 v# p( q) ?+ o% n
    
11. s
    0 e+ b% b7 w7 h1 @) a* H
12. toc( K( L8 z1 @2 R* O3 |3 d8 _
    
13. 
    $ ]) M4 q# K# y' q0 S: Q& H
14. s =
    9 R; w, H3 Q) u& q$ N5 X. Z/ v  S3 e2 u
15. 
    : }; l6 a! H8 `* d) ~& \\" {
16.   2.4002e+007- ]3 z4 q1 e2 l\" f! ?; C
    
17. 
    / b' F# i1 R: b' G$ S, k6 j, G2 t
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];, I( n! ^; d6 M: D  n6 M
   
4. mvar:
   
5. \\" _; l: M9 y3 w& \7 O  X
6. f(x)=x+1;
   
7. 5 ?9 a0 @& c8 J- k4 w& h! ~, G
8. main0(::a)=$ {+ ^5 _3 M) x- [2 c' H3 \
   
10. a = linspace(2,2,1000),' |' [$ ?) |9 R  M# o$ Q
    
12. t0=clock(), s=0,
    
13. $ c+ h\\" p; q1 r0 Y5 n3 e' Q9 K% e
14. k=0, while{k<1000,
    
15. ( V( V$ f+ \, d: a4 R
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++},; ^- o0 s, h7 H8 q
    
18.   k++! M0 c  k+ {) _0 n, c% T
    
20. },) y1 T/ O) J# c8 R! g
    
22. s;# ?; Z, D' E) p5 Q5 Z. O/ K( L$ E
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">, I( n! ^; d6 M: D  n6 M</font><br /> mvar:<br /> <span style=\\"display:none\\">\\" _; l: M9 y3 w& \7 O  X</span> f(x)=x+1;<br /> <span style=\\"display:none\\">5 ?9 a0 @& c8 J- k4 w& h! ~, G</span> main0(::a)=<font class=\\"jammer\\">$ {+ ^5 _3 M) x- [2 c' H3 \</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">' |' [$ ?) |9 R  M# o$ Q</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">$ c+ h\\" p; q1 r0 Y5 n3 e' Q9 K% e</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">( V( V$ f+ \, d: a4 R</span>   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\\">; ^- o0 s, h7 H8 q</font><br />   k++<font class=\\"jammer\\">! M0 c  k+ {) _0 n, c% T</font><br /> },<font class=\\"jammer\\">) y1 T/ O) J# c8 R! g</font><br /> s;<font class=\\"jammer\\"># ?; Z, D' E) p5 Q5 Z. O/ K( L$ E</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
5.172

二者效率现在只有一倍差距。

－－－－－－－－－－

不过，显然matlab和Forcal的程序应该这样写：

' ?1 e; H# x6 W5 J% c$ t; d* rmatlab程序：

1. clear all9 ]- @5 x0 c- V! y' ~$ F: b
   
2. f=@(x)x+1;* d5 [* j9 Y0 D, O$ I; G
   
3. a = linspace(2,2,1000);
   ; _1 H; R* j; H\" q# \4 D; u- @
4. tic7 A( @5 y# e& S- R' ?6 K+ g
   
5. s=0;
   - f; b' s6 ~1 i  T( Z# T
6. for k=1:1000
   ! A; ?: ~' f0 k
7.   for i=1:1000
   6 L  A/ K- }# ^! f& e* t
8.     t=a(i);
   : z0 C, ?# L5 s/ Z' B- C  @) p
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);) G& O3 M4 m$ K  E: a3 j6 `' t7 x
   
10.   end
    0 p8 M0 I7 N, Z0 A- P8 b) N
11. end6 I9 L0 `2 e6 q& u) `: v
    
12. s
    % n$ t+ w0 V& l. e: V
13. toc5 W% u% @& Z' [! t2 w
    
14. 
    7 @' X  m6 }5 o
15. s =: I% e$ Z* m. H0 b
    
16. 3 a, ^: [% N% ~: T% G
    
17.   2.4002e+0071 H) n3 ^7 Y5 M% p% [
    
18. 
    . O! i$ |: k2 ^. z% i* {
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. % s  H. w5 d5 \( ]3 F\\" ?0 y
4. mvar:\\" \\\" r' M8 b4 e; v; R2 ^
   
6. f(x)=x+1;9 ]. X* s# u+ ]; f6 [
   
8. main0(::a)=- S% D; ]0 w: ^\\" l$ |3 V. y
   
10. a = linspace(2,2,1000),0 i8 V! a* a1 X7 b( M- d+ d
    
12. t0=clock(), s=0,
    
13.   w& j' }3 W; P) i& }; l1 X2 |' _
14. k=0, while{k<1000,6 z\\" r* o% z* ?! e9 a
    
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. 9 W  n' U) H3 N+ |* M6 c9 ?- p9 ~
18.   k++
    
19. & P9 {5 v! n) d+ |/ F2 n0 e, @
20. },# i8 C( [! z( ~\\" q2 e
    
22. s;
    
23. 9 f\\" R8 m2 b: [
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">% s  H. w5 d5 \( ]3 F\\" ?0 y</span> mvar:<font class=\\"jammer\\">\\" \\\" r' M8 b4 e; v; R2 ^</font><br /> f(x)=x+1;<font class=\\"jammer\\">9 ]. X* s# u+ ]; f6 [</font><br /> main0(::a)=<font class=\\"jammer\\">- S% D; ]0 w: ^\\" l$ |3 V. y</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">0 i8 V! a* a1 X7 b( M- d+ d</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">  w& j' }3 W; P) i& }; l1 X2 |' _</span> k=0, while{k<1000,<font class=\\"jammer\\">6 z\\" r* o% z* ?! e9 a</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\\">9 W  n' U) H3 N+ |* M6 c9 ?- p9 ~</span>   k++<br /> <span style=\\"display:none\\">& P9 {5 v! n) d+ |/ F2 n0 e, @</span> },<font class=\\"jammer\\"># i8 C( [! z( ~\\" q2 e</font><br /> s;<br /> <span style=\\"display:none\\">9 f\\" R8 m2 b: [</span> [clock()-t0]/1000;

结果：
24002271.69189384
: o+ T; l" a1 n1.625

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

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
2 B. D' f+ j% U2 b3 M) q- C
: [' }' U) A/ x; _Forcal代码：

2. !using["sys"];
   
3. 6 r: {4 i$ v% s0 \, }
4. mvar:8 o! m' J1 p( T' n9 D' c/ Y3 F
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组( ?, T* R5 J! U: u: y  V6 U
   
8. t0=clock(),4 r* X  R. ^5 }0 ]0 C1 Z
   
10. k=0, while{k<10000,
    
11. ; a2 ?& ^+ u( ]1 i
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. ' G8 N( J+ D4 r- X+ b
14.   k++
    
15. ( O8 v6 K4 K+ }) K
16. },  Y+ h+ S& C/ a7 x% E6 h) o8 R\\" f
    
18. [clock()-t0]/1000;& E' R& J\\" X9 b\\" Y4 z( }5 i* @
    
20. : z% R9 X0 y, a4 u& f8 W
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组$ I; m. `! s$ A; i- ]4 d. V; L! w- m
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值, J4 O# D' D- R; O
    
26. t0=clock(),
    
27. 2 W0 ?4 ?, E# T; R- [3 F& q* L8 z
28. s=0,) J( e6 N& I0 l1 e5 P$ f9 Q: c
    
30. k=0, while{k<10000,# ~/ V/ x0 O; [& [
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值8 h5 t2 s( Z* z: f) j! x
    
34.   k++
    
35. # ?- a% `1 y* q, d% q* w  g
36. },
    
37. . E' g3 B5 z$ l9 Q* I7 ]
38. s;, q; r- D* V: o' J! O. z
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">6 r: {4 i$ v% s0 \, }</span> mvar:<font class=\\"jammer\\">8 o! m' J1 p( T' n9 D' c/ Y3 F</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">( ?, T* R5 J! U: u: y  V6 U</font><br /> t0=clock(),<font class=\\"jammer\\">4 r* X  R. ^5 }0 ]0 C1 Z</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">; a2 ?& ^+ u( ]1 i</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">' G8 N( J+ D4 r- X+ b</span>   k++<br /> <span style=\\"display:none\\">( O8 v6 K4 K+ }) K</span> },<font class=\\"jammer\\">  Y+ h+ S& C/ a7 x% E6 h) o8 R\\" f</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">& E' R& J\\" X9 b\\" Y4 z( }5 i* @</font><br /> <font class=\\"jammer\\">: z% R9 X0 y, a4 u& f8 W</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">$ I; m. `! s$ A; i- ]4 d. V; L! w- m</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">, J4 O# D' D- R; O</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">2 W0 ?4 ?, E# T; R- [3 F& q* L8 z</span> s=0,<font class=\\"jammer\\">) J( e6 N& I0 l1 e5 P$ f9 Q: c</font><br /> k=0, while{k<10000,<font class=\\"jammer\\"># ~/ V/ x0 O; [& [</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">8 h5 t2 s( Z* z: f) j! x</font><br />   k++<br /> <span style=\\"display:none\\"># ?- a% `1 y* q, d% q* w  g</span> },<br /> <span style=\\"display:none\\">. E' g3 B5 z$ l9 Q* I7 ]</span> s;<font class=\\"jammer\\">, q; r- D* V: o' J! O. z</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
20000000.
0.579      //取值用时间，秒