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

forcal

1、简单循环

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

1. clear all' n! E+ R8 c# K* p, t( T
   
2. tic8 Y) x* h/ W! F
   
3. s=0;  o' X- P# ]; R; \, P$ j
   
4. for k=1:10000& F  T, U- ?  s8 v6 M
   
5.   for i=1:1000& [! e# {9 ^2 |
   
6.     s=s+1;
   . G( n7 C- i! {. S4 X9 f/ Q! m
7.   end
     g' V2 S& J  _& e( v0 W+ W
8. end5 c# }! c3 S3 _0 Q3 P
   
9. s9 I5 A. {  H3 R3 j0 H
   
10. toc6 v# k3 q; P0 P2 c  z' `7 U\" }- y
    
11. 
    ; }+ R: m! I$ E1 R! k% T+ \% X
12. s =
    + g6 b$ v, d\" h
13. 
    ( f* R  `8 w\" G6 a! N' n
14.     10000000
    7 Q$ e) A; x8 q% Z# Z
15. 
    2 L2 x8 b7 Q' X0 i( B( T7 I
16. Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. ; l6 P+ q! f- J\\" b- O* e+ [$ R
4. mvar:4 s$ w' A$ D; D
   
6. t0=clock(),
   
7. + L0 N' {% \; d; z+ r% ]
8. s=0,- \4 i6 [$ S6 o) t/ \! J9 ^
   
10. k=0, while{k<10000,
    
11. # H+ f* M. l$ M\\" [2 G7 b4 K! z. ]
12.   i=0, while{i<1000,s=s+1, i++},
    
13. # N- D. a' {, r
14.   k++
    
15. 5 t7 U' _\\" g5 R: ?
16. },
    
17. ) K+ S\\" v- _* \+ N9 N) g. W$ U
18. s;
    
19. : _- H6 Z7 B* b' Y8 J
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">; l6 P+ q! f- J\\" b- O* e+ [$ R</span> mvar:<font class=\\"jammer\\">4 s$ w' A$ D; D</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">+ L0 N' {% \; d; z+ r% ]</span> s=0,<font class=\\"jammer\\">- \4 i6 [$ S6 o) t/ \! J9 ^</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\"># H+ f* M. l$ M\\" [2 G7 b4 K! z. ]</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\"># N- D. a' {, r</span>   k++<br /> <span style=\\"display:none\\">5 t7 U' _\\" g5 R: ?</span> },<br /> <span style=\\"display:none\\">) K+ S\\" v- _* \+ N9 N) g. W$ U</span> s;<br /> <span style=\\"display:none\\">: _- H6 Z7 B* b' Y8 J</span> [clock()-t0]/1000;

结果：
10000000.
0 f# f0 p: }" s& _# g( t0.704  秒
# q! F: W3 l9 ?, ?
; S# a! ]# _) y' b7 `二者运行效率相差7倍多。
3 I1 }; I7 f# d  Y2 v( n
+ C+ H9 c7 @) |+ H8 r0 M% r7 @2、循环中稍微加点东西

( y" s: X6 i1 q6 X: s/ uMatlab代码：

1. clear all
   : P1 N; ], Q' R\" U8 J: x* t* h
2. tic
   * }/ S1 Z9 M4 A. d2 o; f
3. s=0;! L  ?: F$ F3 F4 Z
   
4. for k=1:1000
   & P. ~4 r: g' ^2 r$ l/ W# T. |* y4 [8 R
5.   for i=1:1000
   0 x, l. F: z! \# e$ t
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);
   # g, _/ I) i9 L. I3 z
7.   end* I' S8 C$ r- x9 l* j
   
8. end
   ) {5 C& \% z5 \, {! j
9. s& }( J2 D8 |. b1 m7 b, o) I( L
   
10. toc; i! U+ ~4 o) K! r
    
11. 6 H' M6 A* k% m1 V. E
    
12. s =- I; z2 A) ?) q
    
13. 
    + Y3 @: w* g/ H! j0 k' J6 Q
14.   1.1009e+007
    ( s  z( I4 W+ }( y
15. 8 t9 I5 y( X0 F
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 4 b$ t, z! f- }+ F
4. mvar:
   
5. % D- E5 `9 Z) F5 u$ @- u\\" q
6. t0=clock(),
   
7. % z2 ?5 }* t9 h* d; x2 S
8. s=0,# r1 k  E* y, k' X* g1 g' O0 E
   
10. k=1, while{k<=1000,6 M) X0 K( a& s( n$ t1 F0 D: H: o
    
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. 1 W7 Z$ l. B  A0 v# o) e  x
14.   k++; D! H, x2 A7 Z3 \; W\\" C
    
16. },; Q1 F3 s) d8 o) Z\\" J, [
    
18. s;
    
19. . Q1 }9 e4 m+ f: l' x# z9 T
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">4 b$ t, z! f- }+ F</span> mvar:<br /> <span style=\\"display:none\\">% D- E5 `9 Z) F5 u$ @- u\\" q</span> t0=clock(),<br /> <span style=\\"display:none\\">% z2 ?5 }* t9 h* d; x2 S</span> s=0,<font class=\\"jammer\\"># r1 k  E* y, k' X* g1 g' O0 E</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">6 M) X0 K( a& s( n$ t1 F0 D: H: o</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++},<br /> <span style=\\"display:none\\">1 W7 Z$ l. B  A0 v# o) e  x</span>   k++<font class=\\"jammer\\">; D! H, x2 A7 Z3 \; W\\" C</font><br /> },<font class=\\"jammer\\">; Q1 F3 s) d8 o) Z\\" J, [</font><br /> s;<br /> <span style=\\"display:none\\">. Q1 }9 e4 m+ f: l' x# z9 T</span> [clock()-t0]/1000;

结果：
11008953.66602346
1.031
( n+ G) k: J+ v% X
. b, X5 j$ s" {3 {2 C  Z' a二者效率相差已比较小。
  i- G/ A& F2 ]% W3 @( z
3 E1 F  h( O& q* h3、再增加一个函数调用

6 H- O# V: R( U; YMatlab代码：

1. clear all- L' [! J. \. p4 I5 _$ G! C
   
2. f=@(x)x+1;
   0 \: w9 W; [# c# p- B3 V% s5 J) F\" a. M
3. tic
   . A) J6 j4 `% W# ~  b
4. s=0;6 i/ U4 {+ d- z  g  H  h
   
5. for k=1:10004 X+ b: L7 U0 {3 _8 ?% ?
   
6.   for i=1:1000
   4 Z* E0 z, n* d* l9 }0 Z8 H6 r
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);; I6 V% ?- k& ^5 ~/ j
   
8.   end8 {% W* H; ?8 C9 ~
   
9. end
   ! M2 D$ l) e0 Q: _+ R
10. s5 H! |  R' Q$ |  U' W
    
11. toc
    8 X; K+ h2 |$ u/ d  P
12. 
    - Y& Y7 R) C$ }& w5 @# @! B
13. s =
    ' t4 }5 V1 [: n
14. ) b8 h, g. O( `
    
15.   5.1251e+008
    , f( v9 \) j. v) c! H
16. 7 Q: t  z) m1 D
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. ) ~7 j& w& f5 ]' [7 M5 k
4. f(x)=x+1;& |. b1 C\\" p6 j\\" S
   
6. mvar:
   
7. $ y. P/ K9 R; [4 ]\\" Q3 A! E
8. t0=clock(),
   
9. 8 ]; v5 g' r( Y. o9 ~! E1 H6 E
10. s=0,
    
11. 8 ^3 [1 q+ D' z% _
12. k=1, while{k<=1000,
    
13. 7 n7 @. f, z0 ~2 v+ v
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 z% j\\" u4 N# V, L) G+ X
    
16.   k++2 {. l9 B4 @( Z8 Z! R5 b/ ^
    
18. },
    
19. 5 `1 ]0 j) [3 c0 @6 E2 G
20. s;2 O8 Y5 v\\" f9 x( u
    
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">) ~7 j& w& f5 ]' [7 M5 k</span> f(x)=x+1;<font class=\\"jammer\\">& |. b1 C\\" p6 j\\" S</font><br /> mvar:<br /> <span style=\\"display:none\\">$ y. P/ K9 R; [4 ]\\" Q3 A! E</span> t0=clock(),<br /> <span style=\\"display:none\\">8 ]; v5 g' r( Y. o9 ~! E1 H6 E</span> s=0,<br /> <span style=\\"display:none\\">8 ^3 [1 q+ D' z% _</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">7 n7 @. f, z0 ~2 v+ v</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 z% j\\" u4 N# V, L) G+ X</font><br />   k++<font class=\\"jammer\\">2 {. l9 B4 @( Z8 Z! R5 b/ ^</font><br /> },<br /> <span style=\\"display:none\\">5 `1 ]0 j) [3 c0 @6 E2 G</span> s;<font class=\\"jammer\\">2 O8 Y5 v\\" f9 x( u</font><br /> [clock()-t0]/1000;

结果：
& I6 }3 a9 [& z% k) t; S512508953.6658993
1.078
3 }+ M$ ]( d: Q, B( m7 P; c# e
可以看出，Matlab效率已不及Forcal。
( Y% `3 |  D$ F+ p- e故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

Matlab代码：

1. >> clear all
   # I; `0 q; b/ L
2. a = linspace(1,1,1000);3 g& O  y  b, g& ~! m% o7 }, L
   
3. tic
   7 K+ ?0 K0 R# ^3 \) G. z
4. for k=1:10000  @\" O1 Q( X: m2 \( W3 |, O
   
5.   for i=1:1000
   5 K0 G8 o+ Z! z# ~
6.     a(i)=2;! J% |\" [: |: [) m. I
   
7.   end
   5 n0 i7 S% [( u/ H, D. s5 ]* p
8. end
   ' Y% o+ R$ Y5 ~5 M  _: _
9. toc, Z, Q: Y# L: D4 m* Q\" a
   
10. 
    \" f2 w) q& |: v( U+ z: ]/ }$ V, d
11. tic; \1 S% ]' i6 X+ x' S, W/ z+ L9 s
    
12. s=0;
    7 x& c( S( b( n2 G0 T\" G
13. for k=1:100005 b- U6 L2 p3 a  h
    
14.   for i=1:1000
    0 ~$ D$ U1 w& L\" O1 D
15.     s=s+a(i);
    . t' b+ a  x, G5 e
16.   end\" b. U$ A, B* z! L7 ]
    
17. end
    \" Y/ x& C3 e/ k! ?/ y( T
18. s
    6 q0 I) q0 v4 |8 {, t\" s! C7 E* H6 w
19. toc8 ?$ ]2 u! a& y& S! e$ T
    
20. Elapsed time is 0.140658 seconds.
    , J. M' \5 L; t8 W# x& _
21. ' l) p4 b1 J; B8 o8 t% W' _
    
22. s =
    7 g6 |, ?! w9 h8 J* F) M5 s
23. $ N( Q) G5 X0 w8 X+ C. u
    
24.     20000000. O! c$ c\" m9 K5 K; O
    
25. 
    4 T- [+ _+ P5 E
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
& ~: ]9 L+ E6 H& k* M
Forcal代码：

2. !using["math","sys"];/ G4 v* [$ y$ `' N\\" \: k+ A5 D
   
4. mvar:
   
5. : L$ D9 t; h3 a
6. main0(::a)=\\" @8 B' t! S9 q
   
8. a = linspace(1,1,1000),
   
9. 3 F2 x1 O% u' p7 o3 r7 {' ~6 X
10. t0=clock(),
    
11. ( T\\" ^\\" }# N' ^& Y
12. k=0, while{k<10000,4 }- f2 N  H% w4 q
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. * e9 R7 l1 u4 A4 P! ?
16.   k++
    
17. # `( G0 Q# f, e4 X! V+ o: p  W
18. },
    
19. , U; I: a- u3 R: W- v
20. [clock()-t0]/1000;
    
21. % \$ B0 f- }( w3 J* i
22. 
    
23. ! H. V( t; n+ z! z
24. main1(::a)=/ l& W' r5 o, Q9 Q2 J$ v$ u# L
    
26. t0=clock(), s=0,
    
27.   W8 p7 I# L4 t8 |) w
28. k=0, while{k<10000,
    
29. $ f# S' E. |9 b1 \- E
30.   i=0, while{i<1000, s=s+a(i), i++},8 k0 N, E3 Y$ I4 z! c  o
    
32.   k++$ r+ ?! l' p, K# ?1 a
    
34. },
    
35. 0 q! b! L0 v: x, x6 g\\" {7 e
36. s;1 B# g7 g\\" b$ i' I\\" I
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">/ G4 v* [$ y$ `' N\\" \: k+ A5 D</font><br /> mvar:<br /> <span style=\\"display:none\\">: L$ D9 t; h3 a</span> main0(::a)=<font class=\\"jammer\\">\\" @8 B' t! S9 q</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">3 F2 x1 O% u' p7 o3 r7 {' ~6 X</span> t0=clock(),<br /> <span style=\\"display:none\\">( T\\" ^\\" }# N' ^& Y</span> k=0, while{k<10000,<font class=\\"jammer\\">4 }- f2 N  H% w4 q</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">* e9 R7 l1 u4 A4 P! ?</span>   k++<br /> <span style=\\"display:none\\"># `( G0 Q# f, e4 X! V+ o: p  W</span> },<br /> <span style=\\"display:none\\">, U; I: a- u3 R: W- v</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">% \$ B0 f- }( w3 J* i</span> <br /> <span style=\\"display:none\\">! H. V( t; n+ z! z</span> main1(::a)=<font class=\\"jammer\\">/ l& W' r5 o, Q9 Q2 J$ v$ u# L</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">  W8 p7 I# L4 t8 |) w</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">$ f# S' E. |9 b1 \- E</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">8 k0 N, E3 Y$ I4 z! c  o</font><br />   k++<font class=\\"jammer\\">$ r+ ?! l' p, K# ?1 a</font><br /> },<br /> <span style=\\"display:none\\">0 q! b! L0 v: x, x6 g\\" {7 e</span> s;<font class=\\"jammer\\">1 B# g7 g\\" b$ i' I\\" I</font><br /> [clock()-t0]/1000;

结果：
# _+ p0 T. S' M  t* ]15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒

－－－－－－－－－－

用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];: M& |: H. w$ O- ~. i
   
4. mvar:
   
5. 5 P: J  [, z: H5 E
6. main0(::a)=
   
7. 9 a& u/ [& j8 i\\" o$ m
8. a = linspace(1,1,1000),# N5 T- Y! H. [3 l, h) {6 J
   
10. t0=clock(),4 s: J: f1 C& X. c  J( u4 D
    
12. k=0, while{k<10000,
    
13. - Q) n# {, o5 ?% j) L9 o
14.   i=0, while{i<1000, A(a,i)=2, i++},0 r; z* n2 k+ N8 E
    
16.   k++1 u\\" ?; S, H) U* `
    
18. },
    
19. ) ?8 p' B: O( _! c
20. [clock()-t0]/1000;: ~5 O3 j2 _9 p; T\\" v& s
    
22. 
    
23. 2 i2 u$ R2 `4 k! \# J# ~% s
24. main1(::a)=
    
25. $ t/ A  X) n) V4 }7 y- e7 Q
26. t0=clock(), s=0,5 I* p; G! q: n6 V$ U: L
    
28. k=0, while{k<10000,* z* S( p1 u\\" e\\" I4 ]
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},* h/ \2 {4 l, v7 ^( A
    
32.   k++- j# S- f* N9 _\\" a$ e9 B+ J
    
34. },
    
35. 1 a8 y' M9 o6 P% o( e1 L8 M
36. s;
    
37. ' X* n! ?8 Y, {. J
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">: M& |: H. w$ O- ~. i</font><br /> mvar:<br /> <span style=\\"display:none\\">5 P: J  [, z: H5 E</span> main0(::a)=<br /> <span style=\\"display:none\\">9 a& u/ [& j8 i\\" o$ m</span> a = linspace(1,1,1000),<font class=\\"jammer\\"># N5 T- Y! H. [3 l, h) {6 J</font><br /> t0=clock(),<font class=\\"jammer\\">4 s: J: f1 C& X. c  J( u4 D</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">- Q) n# {, o5 ?% j) L9 o</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">0 r; z* n2 k+ N8 E</font><br />   k++<font class=\\"jammer\\">1 u\\" ?; S, H) U* `</font><br /> },<br /> <span style=\\"display:none\\">) ?8 p' B: O( _! c</span> [clock()-t0]/1000;<font class=\\"jammer\\">: ~5 O3 j2 _9 p; T\\" v& s</font><br /> <br /> <span style=\\"display:none\\">2 i2 u$ R2 `4 k! \# J# ~% s</span> main1(::a)=<br /> <span style=\\"display:none\\">$ t/ A  X) n) V4 }7 y- e7 Q</span> t0=clock(), s=0,<font class=\\"jammer\\">5 I* p; G! q: n6 V$ U: L</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">* z* S( p1 u\\" e\\" I4 ]</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">* h/ \2 {4 l, v7 ^( A</font><br />   k++<font class=\\"jammer\\">- j# S- f* N9 _\\" a$ e9 B+ J</font><br /> },<br /> <span style=\\"display:none\\">1 a8 y' M9 o6 P% o( e1 L8 M</span> s;<br /> <span style=\\"display:none\\">' X* n! ?8 Y, {. J</span> [clock()-t0]/1000;

结果：
10.469
20000000.
; H9 z  {4 A! k& z4 ?3.781

1 F4 ?  r7 K0 X; Z" r7 j: c# o－－－－－－－
* G" |! y- I+ l4 K) m1 `6 K) X! o
在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
7 `4 e7 U4 U, ^4 q. M! t在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

* S% p- T6 `% l' m9 Y! |matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
% d) y$ u+ s" _! u7 P. c7 {
Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

& z1 |! J8 S5 O( Z( x! O+ A) O# ?& q5、在循环中加点东西

matlab代码：

1. >> clear all
   \" n2 d\" ?% |& ^6 z2 v+ e
2. f=@(x)x+1;2 c- v\" Y, @$ b5 N* g6 Q9 i4 h
   
3. a = linspace(2,2,1000);/ ^& A. t' b; o
   
4. tic
   / ]\" G' b) n3 [' p
5. s=0;' u' }! Z6 t, w9 A  ~0 |
   
6. for k=1:10000 b4 {& ~) {; D7 e
   
7.   for i=1:1000/ N1 s3 x* Y, R( }' h0 F
   
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));7 h2 k; m9 c: q! L( o- G' }/ v5 n
   
9.   end
   4 V- l- H* B! H! t
10. end
    & F* L! p8 ~' C
11. s
    / t5 z! o4 d/ Z& T) `
12. toc7 x( n! O# n$ f6 ]+ W
    
13. ; s2 I6 r7 V. [; `0 J9 B
    
14. s =
    : k  @$ x# P\" j1 p
15. 
    + \7 |  @! M. y1 i! X  [) K0 m
16.   2.4002e+007+ |4 ?- Z1 w' u$ x+ A
    
17. & x9 ?7 o7 j! O/ y8 k1 Z* V
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];+ s: c- H9 [\\" K$ ^& ?1 O0 [
   
4. mvar:
   
5. 6 L' z7 a7 i7 }, B! T+ k5 o$ o
6. f(x)=x+1;: v& Q2 B! D2 E& S& I
   
8. main0(::a)=7 [0 Q+ z- }- z& ~
   
10. a = linspace(2,2,1000),# v* C* h$ H\\" `\\" W# G5 v
    
12. t0=clock(), s=0,
    
13. * I# `4 V\\" `: E+ p( @+ H5 _- A$ |* [
14. k=0, while{k<1000,
    
15. 7 _' U1 _6 ^* f; U/ S
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* F5 j! n\\" M1 K! c( u
18.   k++
    
19. 8 T6 h# K- r& T1 N: T7 P
20. },
    
21. 1 X- j6 L% o1 h5 w5 x
22. s;5 ^\\" a( }1 q0 g$ Q. D+ I7 F6 e
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">+ s: c- H9 [\\" K$ ^& ?1 O0 [</font><br /> mvar:<br /> <span style=\\"display:none\\">6 L' z7 a7 i7 }, B! T+ k5 o$ o</span> f(x)=x+1;<font class=\\"jammer\\">: v& Q2 B! D2 E& S& I</font><br /> main0(::a)=<font class=\\"jammer\\">7 [0 Q+ z- }- z& ~</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\"># v* C* h$ H\\" `\\" W# G5 v</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">* I# `4 V\\" `: E+ p( @+ H5 _- A$ |* [</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">7 _' U1 _6 ^* f; U/ S</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++},<br /> <span style=\\"display:none\\">2 L* F5 j! n\\" M1 K! c( u</span>   k++<br /> <span style=\\"display:none\\">8 T6 h# K- r& T1 N: T7 P</span> },<br /> <span style=\\"display:none\\">1 X- j6 L% o1 h5 w5 x</span> s;<font class=\\"jammer\\">5 ^\\" a( }1 q0 g$ Q. D+ I7 F6 e</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
8 [) s* l0 L6 D3 w5.172
1 o7 f# h% `; \+ {, S% w- \& t) @2 b, }
( y3 e, j, [1 a; p) b二者效率现在只有一倍差距。

－－－－－－－－－－

$ G/ U3 n( m, E2 w; Q! l6 h( i不过，显然matlab和Forcal的程序应该这样写：
' R" f  t) b/ }8 {$ X) T! Q
# `5 X7 n: s8 f9 H; T! M7 E/ Gmatlab程序：

1. clear all1 |# r) C: k& Q+ Q
   
2. f=@(x)x+1;
   & }  B& Z2 C8 s+ m2 y
3. a = linspace(2,2,1000);
   ; k! D% V' D\" N0 X' E
4. tic5 g& |& [) b/ a9 h
   
5. s=0;( q\" w' }5 O2 U  @& M8 Q! a5 h# y4 L8 z
   
6. for k=1:1000, n1 b, g7 r4 J# E2 N5 ]; n- _3 g
   
7.   for i=1:1000
   ; r( R; p. E% _# o) ~/ `
8.     t=a(i);2 k. `2 s6 g3 b+ ?- e4 f
   
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);* h) |2 G; o/ e' o; a
   
10.   end- b) k, K! q8 k( N1 Q6 {0 G3 S' e
    
11. end$ [6 ^, N& l& U+ _& `; E- b8 {
    
12. s
    ; N! V2 C5 h4 H2 g' Y# s
13. toc
    # b8 q5 O  ~4 c) z% S: B
14. 
    : S7 p$ a& P) ^
15. s =, m9 [* m7 ~2 ~% @! P; x5 e( _/ Q
    
16. 
    3 Y! o) _& ]4 y* K! G% o
17.   2.4002e+007' Z\" n) ]7 {! \, o
    
18. & y5 f( I1 w8 J\" r
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];# g1 O7 A( |* b: T& v. h7 w( x
   
4. mvar:
   
5. ' C1 f( p# H( K8 a& \
6. f(x)=x+1;
   
7. . I+ h, p/ {' e- t\\" f3 X
8. main0(::a)=\\" p/ T: e/ w- Z) }
   
10. a = linspace(2,2,1000),' {1 u\\" _7 \' Z4 ]\\" k1 a1 P8 N
    
12. t0=clock(), s=0,
    
13. $ ]0 Y, W7 h+ v$ \- t
14. k=0, while{k<1000,8 e6 B. H% e- t  }; ]2 t3 [8 {
    
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. 6 Z$ i( ~1 F* c, \0 T: N, J
18.   k++# O& n: h9 N$ y/ |  J
    
20. },
    
21. ) `$ S# ~% @( L
22. s;( E4 p  N+ L+ }$ y! ]
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\"># g1 O7 A( |* b: T& v. h7 w( x</font><br /> mvar:<br /> <span style=\\"display:none\\">' C1 f( p# H( K8 a& \</span> f(x)=x+1;<br /> <span style=\\"display:none\\">. I+ h, p/ {' e- t\\" f3 X</span> main0(::a)=<font class=\\"jammer\\">\\" p/ T: e/ w- Z) }</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">' {1 u\\" _7 \' Z4 ]\\" k1 a1 P8 N</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">$ ]0 Y, W7 h+ v$ \- t</span> k=0, while{k<1000,<font class=\\"jammer\\">8 e6 B. H% e- t  }; ]2 t3 [8 {</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\\">6 Z$ i( ~1 F* c, \0 T: N, J</span>   k++<font class=\\"jammer\\"># O& n: h9 N$ y/ |  J</font><br /> },<br /> <span style=\\"display:none\\">) `$ S# ~% @( L</span> s;<font class=\\"jammer\\">( E4 p  N+ L+ }$ y! ]</font><br /> [clock()-t0]/1000;

结果：
1 E' _: x( W2 z6 `24002271.69189384
1.625

/ ^: D  P  _: }! G8 E- ematlab效率已不及Forcal。

# h3 D6 e  w! P7 C1 {5 M2 }4 h" E提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

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

! O" y% d& e3 k  pForcal代码：

2. !using["sys"];
   
3. : h& a  f/ E7 `* s& s$ V1 F) b
4. mvar:
   
5. 4 V5 T8 w( M) b% z3 p2 c
6. a="\&\[1000]"/4,    //用字符串模拟静态数组3 a1 R  D5 N& e* T: r; p
   
8. t0=clock(),3 T  |8 b6 x$ g; G* T. E0 W
   
10. k=0, while{k<10000,
    
11. , y/ g3 l) ^9 x2 z+ L0 C+ u2 r( X
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值9 ^1 Q* ^1 Y; W1 J( i6 C
    
14.   k++; d/ |. @4 f4 F+ P2 K1 P
    
16. },  I+ h/ H, n1 D
    
18. [clock()-t0]/1000;
    
19. - Q. L6 h$ I+ e1 u7 q
20. % j0 E+ M1 Y0 r\\" G, O9 G& W2 I
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. ( D7 |) f0 c5 s$ f
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. , U: }+ V9 N( L9 w8 I9 @3 b
26. t0=clock(),- T$ o' Q& q+ H; }; {3 q! M
    
28. s=0,; m; R9 ?5 w; y5 G2 A5 ^
    
30. k=0, while{k<10000,$ ?7 w$ G& o\\" y; a+ ~\\" c% d  I0 v% E
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值' ?  Z7 X2 T& c$ n9 ?2 D
    
34.   k++
    
35. * F* C4 [; h9 N0 S8 x1 i1 A# i
36. },
    
37. 3 G8 o$ d* H0 U: D( x7 M9 ^
38. s;# e, n; ^% B0 A& T
    
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">: h& a  f/ E7 `* s& s$ V1 F) b</span> mvar:<br /> <span style=\\"display:none\\">4 V5 T8 w( M) b% z3 p2 c</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">3 a1 R  D5 N& e* T: r; p</font><br /> t0=clock(),<font class=\\"jammer\\">3 T  |8 b6 x$ g; G* T. E0 W</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">, y/ g3 l) ^9 x2 z+ L0 C+ u2 r( X</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">9 ^1 Q* ^1 Y; W1 J( i6 C</font><br />   k++<font class=\\"jammer\\">; d/ |. @4 f4 F+ P2 K1 P</font><br /> },<font class=\\"jammer\\">  I+ h/ H, n1 D</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">- Q. L6 h$ I+ e1 u7 q</span> <font class=\\"jammer\\">% j0 E+ M1 Y0 r\\" G, O9 G& W2 I</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">( D7 |) f0 c5 s$ f</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">, U: }+ V9 N( L9 w8 I9 @3 b</span> t0=clock(),<font class=\\"jammer\\">- T$ o' Q& q+ H; }; {3 q! M</font><br /> s=0,<font class=\\"jammer\\">; m; R9 ?5 w; y5 G2 A5 ^</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">$ ?7 w$ G& o\\" y; a+ ~\\" c% d  I0 v% E</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">' ?  Z7 X2 T& c$ n9 ?2 D</font><br />   k++<br /> <span style=\\"display:none\\">* F* C4 [; h9 N0 S8 x1 i1 A# i</span> },<br /> <span style=\\"display:none\\">3 G8 o$ d* H0 U: D( x7 M9 ^</span> s;<font class=\\"jammer\\"># e, n; ^% B0 A& T</font><br /> [clock()-t0]/1000;

结果：
7 x9 P8 [# w7 ~+ T  U! _( f0.5620000000000001  //赋值用时间，秒
20000000.
0.579      //取值用时间，秒