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

forcal

1、简单循环

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

1. clear all8 ?; L- q1 G* I) v0 d7 f0 U# T
   
2. tic
   $ W# b4 p* t1 A# T
3. s=0;
   # M9 I; z8 s, \5 J\" L5 h) j\" T
4. for k=1:10000+ o0 U  x. p7 V# b& `$ y+ F
   
5.   for i=1:1000
   \" q6 e7 A7 H: X+ m4 ~, R1 e. X6 I2 Y
6.     s=s+1;
   ' g$ p) Z2 g4 h' n) }
7.   end- b, a) g! V4 J9 k1 |
   
8. end2 B& Y$ \) E! w6 l: b5 X3 G
   
9. s
   : }! e3 T2 T: h- a\" i( P1 S4 O
10. toc
    + f1 k, q. i. i7 J) b: _% p
11. 
    0 ?# F) r( U7 o3 m$ v+ r\" s
12. s =/ s; ~5 g5 z, p, j2 L) V
    
13. 
    0 e/ Q\" g  P( B0 @% B0 j\" l! j
14.     10000000- v$ o! r: W. w  W
    
15. 0 w; J. C2 V! v* E
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. 7 D6 `' W- {. N8 E
4. mvar:
   
5. $ t1 h9 h* p5 X) p* W4 Y; a, A# h
6. t0=clock(),
   
7. # J3 Y/ ?1 s  d) `) x
8. s=0,( ^5 F, ?6 U! S- V. B! v$ P
   
10. k=0, while{k<10000,
    
11. 6 d9 |/ [, }1 I1 N) Q- A0 |3 o% l( U
12.   i=0, while{i<1000,s=s+1, i++},
    
13. 3 \, @  J8 }: v! D4 f! [
14.   k++
    
15. 8 J  n3 Z8 O9 D! K
16. },\\" K/ p9 [( e7 U, H
    
18. s;  R& z  M8 D7 M8 D$ c/ @, }* ]: C
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">7 D6 `' W- {. N8 E</span> mvar:<br /> <span style=\\"display:none\\">$ t1 h9 h* p5 X) p* W4 Y; a, A# h</span> t0=clock(),<br /> <span style=\\"display:none\\"># J3 Y/ ?1 s  d) `) x</span> s=0,<font class=\\"jammer\\">( ^5 F, ?6 U! S- V. B! v$ P</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">6 d9 |/ [, }1 I1 N) Q- A0 |3 o% l( U</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">3 \, @  J8 }: v! D4 f! [</span>   k++<br /> <span style=\\"display:none\\">8 J  n3 Z8 O9 D! K</span> },<font class=\\"jammer\\">\\" K/ p9 [( e7 U, H</font><br /> s;<font class=\\"jammer\\">  R& z  M8 D7 M8 D$ c/ @, }* ]: C</font><br /> [clock()-t0]/1000;

结果：
5 J# h" L0 M, c" o" M$ u! @10000000.
" q" }% V8 t( T7 q! N5 u; O0.704  秒

& L2 F6 J% z0 `二者运行效率相差7倍多。

2、循环中稍微加点东西
. d1 s* S4 Y1 o  k, p
Matlab代码：

1. clear all
   4 [! U* x$ F  v% \
2. tic
   ; n6 k  u* w6 ~$ z2 v( v+ R9 z
3. s=0;
   ( o7 R\" X& \) H
4. for k=1:1000
   $ G1 K9 q# o, Y
5.   for i=1:10008 j\" u' R. p! H; A6 k$ a' u2 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);7 {3 e2 K: c. ]9 w! v- {8 x6 q% Z1 `
   
7.   end
   , t9 f- F' S2 Z2 ^+ ^
8. end1 O8 J+ U1 a4 J) B0 L' ^) K
   
9. s
   # Q' _4 U6 z4 T  w0 n/ h
10. toc6 _% k0 P6 O* q7 a$ I& A# @4 ~
    
11. 9 ^/ W! \9 q+ V\" b: \# J( z, s
    
12. s =7 v9 k. t/ s% Q
    
13. 
    ! O\" M  O9 @' l( z0 M) i% F8 a) O
14.   1.1009e+007- x. b  F. Y: v$ N  B' Y
    
15. 
    4 o7 B6 N& Q& e- J: M\" I& T( X9 @
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. + z\\" T2 g# J* Z# f4 ^
4. mvar:
   
5. ! r# ]0 H/ K9 N4 y
6. t0=clock(),
   
7. $ J; S5 l. H1 I* B3 L* u  ^7 s; V
8. s=0,; v. z7 A) |, S) @# A
   
10. k=1, while{k<=1000,
    
11. . I8 a' U. `9 u8 o\\" k- ^
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. & |$ l7 I/ `; s3 r! {$ A& Q
14.   k++
    
15. * U( @% X2 F/ F& |, ~7 {
16. },
    
17. 0 B3 L3 u\\" j' Y9 i- m4 q: s
18. s;2 A) |9 y0 R1 _$ A# N! D# f
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">+ z\\" T2 g# J* Z# f4 ^</span> mvar:<br /> <span style=\\"display:none\\">! r# ]0 H/ K9 N4 y</span> t0=clock(),<br /> <span style=\\"display:none\\">$ J; S5 l. H1 I* B3 L* u  ^7 s; V</span> s=0,<font class=\\"jammer\\">; v. z7 A) |, S) @# A</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">. I8 a' U. `9 u8 o\\" k- ^</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\\">& |$ l7 I/ `; s3 r! {$ A& Q</span>   k++<br /> <span style=\\"display:none\\">* U( @% X2 F/ F& |, ~7 {</span> },<br /> <span style=\\"display:none\\">0 B3 L3 u\\" j' Y9 i- m4 q: s</span> s;<font class=\\"jammer\\">2 A) |9 y0 R1 _$ A# N! D# f</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
1.031

二者效率相差已比较小。
. j3 j1 z$ d. e; c+ G* `" w# l
1 G2 Z+ T3 o$ M; H! _! k) X3、再增加一个函数调用
# K# ~+ z& @" H, ^* c
  o1 C$ w/ l& T& IMatlab代码：

1. clear all
   7 L$ r& R% \6 K# T* Q- z
2. f=@(x)x+1;( Y# R$ e) x( c+ y
   
3. tic\" Z1 o  v( X) f
   
4. s=0;0 ~: _  s. {: H( u
   
5. for k=1:10008 i0 b\" i# V  H, f
   
6.   for i=1:1000( R# b' k# s( P( A9 j' O5 W
   
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);
   0 l' p- ?+ |: H
8.   end
   ! j: K% f6 D( A! N% V7 D
9. end. o! Y' z3 A4 z6 l. b9 O
   
10. s+ n1 Q& @, s4 g6 v7 L
    
11. toc  V5 N+ h1 ^! s\" P4 h0 G2 ?
    
12. 
    - `& i2 G( n9 q4 z$ U5 @/ w# c
13. s =8 ^$ n7 v4 r; ^# J6 }! S1 h2 a
    
14. # O: H' O. T( h\" [! l
    
15.   5.1251e+008
    2 K0 f+ g+ _' o4 g. M
16. 
      Y6 }4 [, C7 L9 q$ D
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];3 ?* O  Y* l' i. L$ Y8 v0 P& N7 w2 ?
   
4. f(x)=x+1;
   
5. ! Y( O, |# B8 {0 x$ t
6. mvar:
   
7. * E2 e% z$ E% _\\" n- w- [
8. t0=clock(),( c( A* P5 N( c$ U1 N' [4 S9 C' Q\\" ^
   
10. s=0,! h( C; N8 T/ W
    
12. k=1, while{k<=1000,
    
13. . |4 l0 X- c3 S/ J
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++},
    
15. * b& H5 M! F8 i  S
16.   k++
    
17. 3 _6 R7 D9 `1 }/ N( C( A; I
18. },
    
19. 4 H' `# ^* O: V4 f9 y1 J/ ]! k
20. s;, @  J; m4 X\\" K) b
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">3 ?* O  Y* l' i. L$ Y8 v0 P& N7 w2 ?</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">! Y( O, |# B8 {0 x$ t</span> mvar:<br /> <span style=\\"display:none\\">* E2 e% z$ E% _\\" n- w- [</span> t0=clock(),<font class=\\"jammer\\">( c( A* P5 N( c$ U1 N' [4 S9 C' Q\\" ^</font><br /> s=0,<font class=\\"jammer\\">! h( C; N8 T/ W</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">. |4 l0 X- c3 S/ J</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++},<br /> <span style=\\"display:none\\">* b& H5 M! F8 i  S</span>   k++<br /> <span style=\\"display:none\\">3 _6 R7 D9 `1 }/ N( C( A; I</span> },<br /> <span style=\\"display:none\\">4 H' `# ^* O: V4 f9 y1 J/ ]! k</span> s;<font class=\\"jammer\\">, @  J; m4 X\\" K) b</font><br /> [clock()-t0]/1000;

结果：
: x/ \& o/ w* q5 `' D512508953.6658993
* e  k& |0 g/ ?  M1.078
* A* g& l; \. M6 Y9 x: J6 N
可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
9 I" M; x: g# Z
/ O& Z: K( W) e) x& C. K7 B4 r) SMatlab代码：

1. >> clear all% j' z8 E9 w  T
   
2. a = linspace(1,1,1000);
   / W$ [  f5 h+ e\" E
3. tic( d- r7 P. U/ [
   
4. for k=1:100004 s8 F% O\" q' V3 u: d\" K! X
   
5.   for i=1:1000
   4 Q\" T6 t- ^\" M5 j; E( N
6.     a(i)=2;
   , V2 N, n$ A7 M
7.   end$ W2 s3 }+ x  q, q7 @
   
8. end6 ^3 q$ N. g3 ]$ S  r/ C; E3 b
   
9. toc
   9 J! S9 |2 X9 v
10. ' l3 K, O! A$ @* x% d! v
    
11. tic\" G: I% L% k6 p+ S6 I& T
    
12. s=0;# g1 v% z: D; }9 V* w: F6 `
    
13. for k=1:10000\" T  U; @4 ~1 C  v
    
14.   for i=1:10008 r, Z4 O/ x- [9 Z% X2 U* w
    
15.     s=s+a(i);
    . Y9 [7 g4 _. D* Y& O& H& I& A1 y
16.   end\" W5 a2 ^( q1 }$ D) e2 |1 h' r
    
17. end
    2 i8 l9 g$ P6 W, E\" R\" h
18. s
    # a( I1 s$ }7 I4 k, J4 _4 Z* `
19. toc
    ) i- V\" j\" U% o8 B
20. Elapsed time is 0.140658 seconds.2 g2 U- J) D+ c  ?9 [1 O5 x
    
21. 
    - l7 Y3 m+ b1 z+ S
22. s =
    9 ^7 c; A) n\" x* J
23. 
    0 V4 M) x( \6 K
24.     20000000
    \" a! u! ^& X4 A
25. 
    8 [( Z+ p7 o2 g7 W
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
4 S3 x6 j8 m: b5 S, C
Forcal代码：

2. !using["math","sys"];
   
3. ) U. x\\" Y5 ^5 I! ]: G
4. mvar:2 ?. h3 M; L& m
   
6. main0(::a)=
   
7. , @; q2 m6 o9 h' k' C
8. a = linspace(1,1,1000),
   
9. 3 B2 p' s/ T0 ?4 ]3 [4 @
10. t0=clock(),
    
11. : n; \# P3 }$ [+ G\\" S
12. k=0, while{k<10000,3 y3 J$ y' C9 \
    
14.   i=0, while{i<1000, a(i)=2, i++},\\" S( [- W5 S) m; U( i% Y
    
16.   k++% J' Z4 e' H# o4 h4 \
    
18. },
    
19. \\" P9 ?! Y$ K# h/ e
20. [clock()-t0]/1000;9 A/ I9 P& U- Q3 V. X\\" ^1 y3 |' R
    
22. 
    
23. 2 G0 B2 g0 q; |2 J$ X  f
24. main1(::a)=- O7 y2 E, I7 z/ c& S
    
26. t0=clock(), s=0,
    
27.   g; w8 e1 h' ~7 i  g\\" y( h
28. k=0, while{k<10000,
    
29. # H, r# a# s\\" }2 d5 f0 T
30.   i=0, while{i<1000, s=s+a(i), i++},. P/ R; \% T. s2 W+ e# C
    
32.   k++
    
33. % q- y9 C  t8 }$ C
34. },
    
35. ' d  Q) Z, f9 ^' W8 ~3 o- ~! A
36. s;0 b; }2 R% Z1 {
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">) U. x\\" Y5 ^5 I! ]: G</span> mvar:<font class=\\"jammer\\">2 ?. h3 M; L& m</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">, @; q2 m6 o9 h' k' C</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">3 B2 p' s/ T0 ?4 ]3 [4 @</span> t0=clock(),<br /> <span style=\\"display:none\\">: n; \# P3 }$ [+ G\\" S</span> k=0, while{k<10000,<font class=\\"jammer\\">3 y3 J$ y' C9 \</font><br />   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">\\" S( [- W5 S) m; U( i% Y</font><br />   k++<font class=\\"jammer\\">% J' Z4 e' H# o4 h4 \</font><br /> },<br /> <span style=\\"display:none\\">\\" P9 ?! Y$ K# h/ e</span> [clock()-t0]/1000;<font class=\\"jammer\\">9 A/ I9 P& U- Q3 V. X\\" ^1 y3 |' R</font><br /> <br /> <span style=\\"display:none\\">2 G0 B2 g0 q; |2 J$ X  f</span> main1(::a)=<font class=\\"jammer\\">- O7 y2 E, I7 z/ c& S</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">  g; w8 e1 h' ~7 i  g\\" y( h</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\"># H, r# a# s\\" }2 d5 f0 T</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">. P/ R; \% T. s2 W+ e# C</font><br />   k++<br /> <span style=\\"display:none\\">% q- y9 C  t8 }$ C</span> },<br /> <span style=\\"display:none\\">' d  Q) Z, f9 ^' W8 ~3 o- ~! A</span> s;<font class=\\"jammer\\">0 b; }2 R% Z1 {</font><br /> [clock()-t0]/1000;

结果：
. O; `. p* B* z15.484  //存数据所用时间，秒
20000000.
, h9 D3 L) j7 k" a8 Q5.531   //取数据所用时间，秒

－－－－－－－－－－
. {# j. Q% d5 T# {5 L/ V$ M
6 Y) M1 T3 H9 x用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. ) W: l0 O\\" U/ N8 F
4. mvar:
   
5. 4 K+ C1 _& I# `
6. main0(::a)=
   
7. / z/ t1 i$ p1 _- v4 T. ], {
8. a = linspace(1,1,1000),' N- d  u  T4 `) Y, A6 K$ p. o
   
10. t0=clock(),0 y! X. R% F5 w4 T  p\\" \6 Z
    
12. k=0, while{k<10000,
    
13. 6 V% e0 N0 |5 A\\" w& M% x
14.   i=0, while{i<1000, A(a,i)=2, i++},+ V$ i$ e* E7 H! ^
    
16.   k++
    
17. ) l  i$ M8 o  u7 q
18. },
    
19. ' d/ F8 Z4 F: N0 _
20. [clock()-t0]/1000;* x5 }. b5 }7 h
    
22. 
    
23. ( s/ v2 o/ ~& L- T\\" L  K
24. main1(::a)=# n* r+ M7 R5 f% F! s\\" e
    
26. t0=clock(), s=0,' P2 Z\\" H: ~0 T5 ]6 l
    
28. k=0, while{k<10000,
    
29. % C. p& z4 ?7 ]+ U
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. ' T- _* s+ [5 I/ o
32.   k++
    
33. 4 q/ ]% j% K: s  Y4 R  \9 S- R
34. },; B9 X8 x* s; p0 ?
    
36. s;/ D  `$ H2 I3 l  B% u3 G* a6 t
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">) W: l0 O\\" U/ N8 F</span> mvar:<br /> <span style=\\"display:none\\">4 K+ C1 _& I# `</span> main0(::a)=<br /> <span style=\\"display:none\\">/ z/ t1 i$ p1 _- v4 T. ], {</span> a = linspace(1,1,1000),<font class=\\"jammer\\">' N- d  u  T4 `) Y, A6 K$ p. o</font><br /> t0=clock(),<font class=\\"jammer\\">0 y! X. R% F5 w4 T  p\\" \6 Z</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">6 V% e0 N0 |5 A\\" w& M% x</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">+ V$ i$ e* E7 H! ^</font><br />   k++<br /> <span style=\\"display:none\\">) l  i$ M8 o  u7 q</span> },<br /> <span style=\\"display:none\\">' d/ F8 Z4 F: N0 _</span> [clock()-t0]/1000;<font class=\\"jammer\\">* x5 }. b5 }7 h</font><br /> <br /> <span style=\\"display:none\\">( s/ v2 o/ ~& L- T\\" L  K</span> main1(::a)=<font class=\\"jammer\\"># n* r+ M7 R5 f% F! s\\" e</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">' P2 Z\\" H: ~0 T5 ]6 l</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">% C. p& z4 ?7 ]+ U</span>   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">' T- _* s+ [5 I/ o</span>   k++<br /> <span style=\\"display:none\\">4 q/ ]% j% K: s  Y4 R  \9 S- R</span> },<font class=\\"jammer\\">; B9 X8 x* s; p0 ?</font><br /> s;<font class=\\"jammer\\">/ D  `$ H2 I3 l  B% u3 G* a6 t</font><br /> [clock()-t0]/1000;

结果：
! ]" p% J. f8 h2 g0 \10.469
" p2 l2 q& A3 r( \: h# l20000000.
" ~0 t* x- S) J" L$ G- s  |, e1 J3.781

－－－－－－－

在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
; p( Q4 c" D& v
matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
- e# {0 s5 {) q/ ^7 n
Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

4 r0 H5 |$ {3 H- R5 n5 `5、在循环中加点东西

( R. `% ^" a& j5 N8 ^: mmatlab代码：

1. >> clear all% _4 k0 d\" v, K! I: L) @' I
   
2. f=@(x)x+1;
   ' ^, |/ a/ F5 {( F
3. a = linspace(2,2,1000);
   + T' G. x2 H\" ?) X: d
4. tic/ N! b) w7 ^4 p1 k
   
5. s=0;
   : g\" W7 G8 B\" f+ r/ u* h( j
6. for k=1:1000# k3 s0 P% n) I  P( ]& i% C
   
7.   for i=1:10009 Z9 h# }1 ~# `: a2 j. g$ R9 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));6 B) K! E0 A; ^8 W
   
9.   end
   9 K% w% H0 P6 a9 G' r
10. end7 K% H- @  H, H. J% H
    
11. s* G( J0 R+ y! J\" J: c. Z' J\" w5 O
    
12. toc
    $ E  G\" \* ~4 D6 F4 B0 ?. v
13.   F! I+ d# N; S$ K& n  |
    
14. s =- m% k: g5 Z8 q. {& p& g: ?
    
15. 
    8 N* E! f1 ]1 n1 a
16.   2.4002e+007
    5 e  Q4 j/ W- }
17. ( Q: G* k2 S/ }3 e6 c/ A& c5 P' [
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];  x1 x7 [& D( |* R\\" H- e( O
   
4. mvar:
   
5. # t) W2 b3 T! J3 H  R$ y
6. f(x)=x+1;9 B- {8 r& H6 G$ v/ M
   
8. main0(::a)=' D  h5 s8 H: Z( L
   
10. a = linspace(2,2,1000),
    
11. 8 p1 q  p( |( r* m
12. t0=clock(), s=0,
    
13. / J: w8 ]+ ]* T- t7 F: f
14. k=0, while{k<1000,, y6 t5 S5 w+ v9 V\\" H\\" f6 j  w: ]
    
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++},0 J6 H$ L# y! X) N: e
    
18.   k++
    
19. \\" o! h9 X  z9 m) f; V( d
20. },2 B, O, E/ o: N# Z
    
22. s;1 u1 u, B6 I8 ~; b* e; d7 t. A
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">  x1 x7 [& D( |* R\\" H- e( O</font><br /> mvar:<br /> <span style=\\"display:none\\"># t) W2 b3 T! J3 H  R$ y</span> f(x)=x+1;<font class=\\"jammer\\">9 B- {8 r& H6 G$ v/ M</font><br /> main0(::a)=<font class=\\"jammer\\">' D  h5 s8 H: Z( L</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">8 p1 q  p( |( r* m</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">/ J: w8 ]+ ]* T- t7 F: f</span> k=0, while{k<1000,<font class=\\"jammer\\">, y6 t5 S5 w+ v9 V\\" H\\" f6 j  w: ]</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\\">0 J6 H$ L# y! X) N: e</font><br />   k++<br /> <span style=\\"display:none\\">\\" o! h9 X  z9 m) f; V( d</span> },<font class=\\"jammer\\">2 B, O, E/ o: N# Z</font><br /> s;<font class=\\"jammer\\">1 u1 u, B6 I8 ~; b* e; d7 t. A</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
8 ^% T0 L( x' @: q6 V$ P" W5.172

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

* ], I7 T' B  l* }－－－－－－－－－－
6 L+ u. p! w& P/ N
不过，显然matlab和Forcal的程序应该这样写：
+ d# Y# u% G/ ~- P! J1 S; ?, j1 {
matlab程序：

1. clear all
   * I. _8 K# {3 t
2. f=@(x)x+1;
   ; x3 Y; I& a, z$ F) s( @% b& d6 h
3. a = linspace(2,2,1000);. i7 T8 t0 Q3 L$ V
   
4. tic, e  r) l6 Z0 ^9 y' r/ a4 \
   
5. s=0;
   5 ^  Y( F+ c9 \7 v1 Y- o
6. for k=1:10009 i3 c7 t6 ]& d  {+ _
   
7.   for i=1:1000( j6 @) b( _1 e; {
   
8.     t=a(i);
   0 v; m7 X: F+ ^5 P- d
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);
   5 p1 |/ R5 V- |5 F; N& O
10.   end- A  d+ f2 i0 R) q- x0 w* g, b
    
11. end4 u4 M/ O: E9 j; O6 J
    
12. s
    * d$ U7 I8 |  Q, E
13. toc
    0 B! @5 N\" m1 J+ @3 h  u
14. 1 P8 ?; F! }\" J2 r
    
15. s =
    & \3 B) F- M' y# a\" I& X* ^8 _
16. 5 `0 d1 W+ y9 N1 S7 Z/ j
    
17.   2.4002e+007
    ) a. K, l! X$ n$ I% R
18. 4 o5 w/ ~- `3 F8 B
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];5 i3 s% C7 @3 h, e, r. ]
   
4. mvar:
   
5. ! @\\" j3 j! V1 C7 g+ W) A
6. f(x)=x+1;\\" k( b8 L& D* q: @  x7 e9 f
   
8. main0(::a)=
   
9. 9 c% Y* m7 r$ J1 ]  ~9 m
10. a = linspace(2,2,1000),( }' m5 f: X+ w5 w
    
12. t0=clock(), s=0,
    
13. 8 g% Y# z# c* Y0 M
14. k=0, while{k<1000,1 C1 K* Y6 |/ C3 w+ q' q8 u
    
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++},9 j+ I4 r& D5 R& S. E  U
    
18.   k++( w* M# D, x/ {& `3 D
    
20. },
    
21. # e5 x/ d, a+ v2 @
22. s;
    
23. + ~. ?\\" M5 W& e0 M$ p  F
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">5 i3 s% C7 @3 h, e, r. ]</font><br /> mvar:<br /> <span style=\\"display:none\\">! @\\" j3 j! V1 C7 g+ W) A</span> f(x)=x+1;<font class=\\"jammer\\">\\" k( b8 L& D* q: @  x7 e9 f</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">9 c% Y* m7 r$ J1 ]  ~9 m</span> a = linspace(2,2,1000),<font class=\\"jammer\\">( }' m5 f: X+ w5 w</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">8 g% Y# z# c* Y0 M</span> k=0, while{k<1000,<font class=\\"jammer\\">1 C1 K* Y6 |/ C3 w+ q' q8 u</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\\">9 j+ I4 r& D5 R& S. E  U</font><br />   k++<font class=\\"jammer\\">( w* M# D, x/ {& `3 D</font><br /> },<br /> <span style=\\"display:none\\"># e5 x/ d, a+ v2 @</span> s;<br /> <span style=\\"display:none\\">+ ~. ?\\" M5 W& e0 M$ p  F</span> [clock()-t0]/1000;

结果：
# |( P5 }4 N$ k- G24002271.69189384
: r7 P& ~7 L6 Q( p+ a5 }1.625

matlab效率已不及Forcal。

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

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
; G  x! W2 }: i5 j% Q
Forcal代码：

2. !using["sys"];
   
3. 3 B2 ~& @( a  J9 M
4. mvar:) ^8 V9 T' i2 t6 a  Q1 F
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组& @$ o5 V9 m  G+ c0 I7 e# w, s
   
8. t0=clock(),' b/ R6 U$ e& b5 Z3 o& q
   
10. k=0, while{k<10000,9 P\\" G9 j( R1 b6 d0 ^
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值- \' q9 y; b, Q. ?
    
14.   k++
    
15. ; P* K, \* _# i( [
16. },; n4 C1 F, P; h7 X1 q/ K
    
18. [clock()-t0]/1000;2 B8 l' K$ ]; H# Z3 G
    
20. 
    
21. 6 I! c2 ^& f! z+ ?; N
22. a="\&\[1000]"/4,    //用字符串模拟静态数组* S$ c. `3 `6 s6 Q3 t\\" s/ b8 B( A4 Q
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值9 W# e7 s' w) J6 c
    
26. t0=clock(),
    
27. \\" m* {' G% D7 L
28. s=0,
    
29. # P. B9 [+ N- u7 H2 d\\" U* H% a
30. k=0, while{k<10000,0 a1 u7 K0 {; _3 J7 `
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值' ?$ P- T8 p/ F9 q7 A5 B# m! P8 A
    
34.   k++
    
35. 1 q! x4 U  O; e4 V, E
36. },( l  u  z# W+ a* z2 \4 F4 W: n
    
38. s;) Y, e) y2 @, D- s# v3 q
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">3 B2 ~& @( a  J9 M</span> mvar:<font class=\\"jammer\\">) ^8 V9 T' i2 t6 a  Q1 F</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">& @$ o5 V9 m  G+ c0 I7 e# w, s</font><br /> t0=clock(),<font class=\\"jammer\\">' b/ R6 U$ e& b5 Z3 o& q</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">9 P\\" G9 j( R1 b6 d0 ^</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">- \' q9 y; b, Q. ?</font><br />   k++<br /> <span style=\\"display:none\\">; P* K, \* _# i( [</span> },<font class=\\"jammer\\">; n4 C1 F, P; h7 X1 q/ K</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">2 B8 l' K$ ]; H# Z3 G</font><br /> <br /> <span style=\\"display:none\\">6 I! c2 ^& f! z+ ?; N</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">* S$ c. `3 `6 s6 Q3 t\\" s/ b8 B( A4 Q</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">9 W# e7 s' w) J6 c</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">\\" m* {' G% D7 L</span> s=0,<br /> <span style=\\"display:none\\"># P. B9 [+ N- u7 H2 d\\" U* H% a</span> k=0, while{k<10000,<font class=\\"jammer\\">0 a1 u7 K0 {; _3 J7 `</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">' ?$ P- T8 p/ F9 q7 A5 B# m! P8 A</font><br />   k++<br /> <span style=\\"display:none\\">1 q! x4 U  O; e4 V, E</span> },<font class=\\"jammer\\">( l  u  z# W+ a* z2 \4 F4 W: n</font><br /> s;<font class=\\"jammer\\">) Y, e) y2 @, D- s# v3 q</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
* B, u: c! g! d20000000.
0.579      //取值用时间，秒