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

forcal

1、简单循环
# a! {7 L: z8 X7 i
1 V" \4 S9 z/ zMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   5 Z( S+ M2 c3 [6 s
2. tic5 \, J, y2 W, u/ _0 H- b
   
3. s=0;, J& L+ Y8 \3 u+ L' f
   
4. for k=1:10000
   1 M! Q9 f+ @. y5 i$ \
5.   for i=1:10001 f9 F6 e! j  {
   
6.     s=s+1;
   8 j! \; I3 U9 A) F
7.   end\" I' ]! a- L& _& C& Q
   
8. end
   + V( m\" d; n- k& L4 y& Y
9. s
     }* {  F4 V7 O0 w
10. toc5 R6 S. a8 r) K* \9 `+ U1 @
    
11. 8 H' j+ k: D8 e; k* q
    
12. s =
    8 [3 X6 G) K7 A5 o8 D
13. 0 q- h8 g5 e$ i+ D7 ?2 g3 Y# |
    
14.     10000000
    ) o% K6 r, F, z0 N) W2 L
15. 
    , x3 M; U\" I& R
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
. H2 v; d" @6 i% _1 N5 b
Forcal没有JIT加速器，代码：

2. !using["sys"];  ~! x, {& F8 v; j4 O. Q3 G. @
   
4. mvar:
   
5. ) ~: s6 B- u6 i; G
6. t0=clock(),
   
7. # G1 M: z\\" }\\" K5 ]
8. s=0,
   
9. + l# ]8 C7 `7 Z, z
10. k=0, while{k<10000,0 ?' l, }$ A\\" v/ _' O# ?. g
    
12.   i=0, while{i<1000,s=s+1, i++},6 {6 U7 v8 R0 H: K6 J
    
14.   k++$ A* K/ J5 i/ C; T4 E
    
16. },
    
17. \\" ~: v! k0 |) t0 X+ _6 w6 W
18. s;4 V, e+ ~8 F( }8 a) `3 P: k
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">  ~! x, {& F8 v; j4 O. Q3 G. @</font><br /> mvar:<br /> <span style=\\"display:none\\">) ~: s6 B- u6 i; G</span> t0=clock(),<br /> <span style=\\"display:none\\"># G1 M: z\\" }\\" K5 ]</span> s=0,<br /> <span style=\\"display:none\\">+ l# ]8 C7 `7 Z, z</span> k=0, while{k<10000,<font class=\\"jammer\\">0 ?' l, }$ A\\" v/ _' O# ?. g</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">6 {6 U7 v8 R0 H: K6 J</font><br />   k++<font class=\\"jammer\\">$ A* K/ J5 i/ C; T4 E</font><br /> },<br /> <span style=\\"display:none\\">\\" ~: v! k0 |) t0 X+ _6 w6 W</span> s;<font class=\\"jammer\\">4 V, e+ ~8 F( }8 a) `3 P: k</font><br /> [clock()-t0]/1000;

结果：
0 U( }# Q9 U7 v2 L  d10000000.
, `0 E  f4 ^4 g0.704  秒
" f+ D+ I' k3 _1 [0 l
, u7 N; g: \- ~) k' f二者运行效率相差7倍多。

2、循环中稍微加点东西
5 e: `  i. s' e9 V
1 Y" m' Y( z# F, Y) ~Matlab代码：

1. clear all
   . O! P1 x  R8 C4 o0 X5 d
2. tic
   % q! f0 M- e' Z6 Z4 ^2 D1 f$ W
3. s=0;- B0 |  g) }% d9 E& L6 f; @' C
   
4. for k=1:1000
   % u  Q0 A* E% t2 e3 ~3 ~
5.   for i=1:10001 g6 j! R7 B6 `4 |( q6 |2 O- \. j) 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);; U1 d- f/ j1 G; b
   
7.   end* u/ ]* j% `  N. u, z
   
8. end
     x4 p, M2 w/ k
9. s
   : u! p8 X- y4 {- _: S7 r% }& p
10. toc
    3 w4 G2 j\" c\" Y# }1 p+ z
11. 6 E' c\" a) q- K' t! v$ r! [* ^& H& H
    
12. s =
    2 L& I; j$ Z! x9 H' J% }& h; G2 W
13. 
    \" s& i9 G% C' z! l+ h1 h$ D/ E& s
14.   1.1009e+0074 F9 q0 r9 T& i4 F: Q
    
15. 
    : _0 {4 S! W# S5 Q. O5 m9 u' o( Z
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. ; ~% {2 Q) u2 T
4. mvar:- q9 R. s. G/ ~: m, f6 r
   
6. t0=clock(),
   
7. 0 F1 y8 w9 n0 L+ c' K
8. s=0,
   
9. / `: j. j7 O4 w0 \& C  k
10. k=1, while{k<=1000,1 I+ ^* f$ f: S( P1 W+ T
    
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++},. X& v  d( m% b% H3 b0 x  I  h
    
14.   k++0 e& l8 E; r# u9 I\\" `. Q% Q. M
    
16. },
    
17. - E& \0 I: ~7 o6 {# x$ N) Z
18. s;
    
19.   U8 t2 E4 F+ T1 e: V
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">; ~% {2 Q) u2 T</span> mvar:<font class=\\"jammer\\">- q9 R. s. G/ ~: m, f6 r</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">0 F1 y8 w9 n0 L+ c' K</span> s=0,<br /> <span style=\\"display:none\\">/ `: j. j7 O4 w0 \& C  k</span> k=1, while{k<=1000,<font class=\\"jammer\\">1 I+ ^* f$ f: S( P1 W+ T</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\\">. X& v  d( m% b% H3 b0 x  I  h</font><br />   k++<font class=\\"jammer\\">0 e& l8 E; r# u9 I\\" `. Q% Q. M</font><br /> },<br /> <span style=\\"display:none\\">- E& \0 I: ~7 o6 {# x$ N) Z</span> s;<br /> <span style=\\"display:none\\">  U8 t2 E4 F+ T1 e: V</span> [clock()-t0]/1000;

结果：
' w2 o% v* K( ]# ]) L11008953.66602346
1.031

二者效率相差已比较小。

3、再增加一个函数调用
7 T+ I+ c, i( b2 j: m8 X
Matlab代码：

1. clear all
   9 G8 P+ c* Q, v  g6 v8 w8 w
2. f=@(x)x+1;\" l! a& t$ r* }; D\" A6 f
   
3. tic
   8 G\" a  A, G3 i5 V% g# u& C
4. s=0;
   2 r9 C1 g; i9 t7 Y* ~; r/ t
5. for k=1:1000
   . u' c& B+ i& R
6.   for i=1:1000
   6 N9 [) ^- N- A% c& O5 J5 E\" m
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);
   9 h! p/ _7 Z\" m
8.   end4 _# m6 H+ w3 H$ a
   
9. end
   : N$ y4 e  a4 q
10. s  x) Y3 y, ?& m' x& P/ L, K
    
11. toc9 H& Z2 j: s) ^8 [& O1 U0 g. c+ c\" }
    
12. 
    \" H: j; [, M! H* @+ ^3 Y0 \0 o
13. s =9 F- p& q* b$ `8 B- M
    
14. % B& |, x3 C\" {) w( v: ~
    
15.   5.1251e+008\" z' m9 |$ A7 _. ?
    
16. 
    9 v' }' D& o+ Q; Y7 [
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 9 G. B9 U\\" H5 O. c' Z. k
4. f(x)=x+1;: t2 I& o( q3 {& S4 p& [$ r- c
   
6. mvar:  ^2 `3 I* ?: S* T% h$ h* x
   
8. t0=clock(),
   
9. ' {& Q2 Q. P9 \
10. s=0,9 O0 k1 E+ P% z/ s
    
12. k=1, while{k<=1000,
    
13. % W* q3 e* ~1 E& V; b2 z7 b
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. 4 v/ a$ ]- w3 E1 l
16.   k++
    
17. 0 E/ `1 U, M- w- I
18. },
    
19. ) C9 }0 _, o9 w; L0 k
20. s;8 f9 d) X) v, `* x: U& U, f
    
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">9 G. B9 U\\" H5 O. c' Z. k</span> f(x)=x+1;<font class=\\"jammer\\">: t2 I& o( q3 {& S4 p& [$ r- c</font><br /> mvar:<font class=\\"jammer\\">  ^2 `3 I* ?: S* T% h$ h* x</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">' {& Q2 Q. P9 \</span> s=0,<font class=\\"jammer\\">9 O0 k1 E+ P% z/ s</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">% W* q3 e* ~1 E& V; b2 z7 b</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\\">4 v/ a$ ]- w3 E1 l</span>   k++<br /> <span style=\\"display:none\\">0 E/ `1 U, M- w- I</span> },<br /> <span style=\\"display:none\\">) C9 }0 _, o9 w; L0 k</span> s;<font class=\\"jammer\\">8 f9 d) X) v, `* x: U& U, f</font><br /> [clock()-t0]/1000;

结果：
512508953.6658993
1.078

( I" a0 q' x, O8 X) S可以看出，Matlab效率已不及Forcal。
8 M( h" ~1 u+ ~故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

Matlab代码：

1. >> clear all
   : U\" a7 |+ l$ t( D& A. ?' g; ?
2. a = linspace(1,1,1000);5 L\" `$ e2 O- O, G6 D
   
3. tic' R5 g8 U\" G$ `* y) R\" P5 P/ R% j1 U
   
4. for k=1:10000
     T5 W\" m$ X  I- b8 g
5.   for i=1:1000* u6 x) O- E4 m6 R, h8 C3 R
   
6.     a(i)=2;- q8 {0 E4 X2 @9 i) K
   
7.   end6 w% ?  V; f% W9 l! [
   
8. end. G9 h3 |5 I1 [6 `0 ~6 e# f
   
9. toc
   / g. J4 }/ B8 C3 Y7 A! S
10. 
    + Y0 ~: N: _$ H6 O
11. tic
    + [5 r, k1 W; `  M! L1 a; `
12. s=0;' n, k3 a7 D# s  T4 w$ N4 r\" h; B9 ]
    
13. for k=1:100002 l\" R1 X0 x$ o0 A5 ^, r
    
14.   for i=1:10000 R5 x( `/ j7 U/ d9 I
    
15.     s=s+a(i);
    8 Y  }. \8 f1 B( K
16.   end
    ! o2 R0 C7 F: P2 s  _) ^3 a
17. end7 y: I8 ^/ R( P3 r% a; B4 f\" _
    
18. s6 h. Z1 U+ a  ?1 l- A: t9 C% G, J
    
19. toc
    $ x& Z3 `3 u- ?( b0 U8 d
20. Elapsed time is 0.140658 seconds.9 A1 B0 t6 Z$ F1 ]; r
    
21. 
      w) l% i6 y0 H2 }
22. s =
    $ s* U& o' i6 [) T* Q% M5 S4 I7 U
23. 
    * ~+ y5 f5 {1 l2 K5 _
24.     200000007 C\" P8 K( O7 |# u+ Z( F
    
25. 
      D# N% G\" h4 i% @6 u$ ~
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
  B; M) E) n/ s) @
Forcal代码：

2. !using["math","sys"];7 [8 A/ Y+ t6 w5 @8 j
   
4. mvar:
   
5. ' I% Y1 V& y7 D. W' p8 Q1 c/ P
6. main0(::a)=8 B' O8 N- ]/ K' h# A' a
   
8. a = linspace(1,1,1000),4 G' z/ w; b9 q: I
   
10. t0=clock(),7 ?\\" X  w: D/ h+ r5 v- J5 C4 o
    
12. k=0, while{k<10000,
    
13. 4 l) M/ X# F/ g3 C% y\\" e& C1 e
14.   i=0, while{i<1000, a(i)=2, i++},# u7 c0 O( u0 G  J; i' ^% g* J
    
16.   k++7 `( I  O9 Q% u8 X$ [
    
18. },
    
19. * P; L' `: T) w$ d& U  Q
20. [clock()-t0]/1000;
    
21. 1 z4 r8 ^4 D! F0 N! [& Q* M$ Q
22. 
    
23. ; Q: S+ k# p7 ~\\" ^
24. main1(::a)=
    
25. - s: ]7 \! E\\" a5 j\\" ]
26. t0=clock(), s=0,- K: K- j/ g- g\\" m2 o
    
28. k=0, while{k<10000,! {# F3 R( v+ ~8 {
    
30.   i=0, while{i<1000, s=s+a(i), i++},\\" z3 _\\" u, K8 F9 _' D\\" j
    
32.   k++$ }) M% Y9 |* _6 q+ A
    
34. },
    
35. 9 v& h) ]- B/ v
36. s;
    
37. ( C' t( x6 X% d! R; z- h, }
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">7 [8 A/ Y+ t6 w5 @8 j</font><br /> mvar:<br /> <span style=\\"display:none\\">' I% Y1 V& y7 D. W' p8 Q1 c/ P</span> main0(::a)=<font class=\\"jammer\\">8 B' O8 N- ]/ K' h# A' a</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">4 G' z/ w; b9 q: I</font><br /> t0=clock(),<font class=\\"jammer\\">7 ?\\" X  w: D/ h+ r5 v- J5 C4 o</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">4 l) M/ X# F/ g3 C% y\\" e& C1 e</span>   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\"># u7 c0 O( u0 G  J; i' ^% g* J</font><br />   k++<font class=\\"jammer\\">7 `( I  O9 Q% u8 X$ [</font><br /> },<br /> <span style=\\"display:none\\">* P; L' `: T) w$ d& U  Q</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">1 z4 r8 ^4 D! F0 N! [& Q* M$ Q</span> <br /> <span style=\\"display:none\\">; Q: S+ k# p7 ~\\" ^</span> main1(::a)=<br /> <span style=\\"display:none\\">- s: ]7 \! E\\" a5 j\\" ]</span> t0=clock(), s=0,<font class=\\"jammer\\">- K: K- j/ g- g\\" m2 o</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">! {# F3 R( v+ ~8 {</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">\\" z3 _\\" u, K8 F9 _' D\\" j</font><br />   k++<font class=\\"jammer\\">$ }) M% Y9 |* _6 q+ A</font><br /> },<br /> <span style=\\"display:none\\">9 v& h) ]- B/ v</span> s;<br /> <span style=\\"display:none\\">( C' t( x6 X% d! R; z- h, }</span> [clock()-t0]/1000;

结果：
$ j" \, z' l3 O% \% K& j7 r15.484  //存数据所用时间，秒
/ i* }# f; [; ]4 b20000000.
5.531   //取数据所用时间，秒

－－－－－－－－－－
0 C3 B+ q7 [5 {& E' M, @
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];2 Z  G2 x* j6 s\\" |* u
   
4. mvar:\\" M; ~. N7 e& D& C( X: c8 _) v
   
6. main0(::a)=
   
7. 2 n9 F( S9 `1 U+ G* T/ i
8. a = linspace(1,1,1000),
   
9. * v7 j3 F1 v4 r) w$ |
10. t0=clock(),
    
11. ; ]3 z0 b# E. O5 h7 r7 m: h+ d
12. k=0, while{k<10000,
    
13. 4 ]8 E# N* g4 a
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. ( l9 X6 i) d, l
16.   k++
    
17. ( [. V9 j% O! L
18. },
    
19. : Y' |# u0 x) ~6 z) }
20. [clock()-t0]/1000;0 y  j: I5 r7 b& ]: d
    
22. 
    
23. \\" H9 h; ~1 U2 m. x0 i( U3 D
24. main1(::a)=6 \8 }  s( u! n( d' J/ }4 h! w: H: ]
    
26. t0=clock(), s=0,
    
27. * Q- h3 d9 F& I+ D
28. k=0, while{k<10000,# ]; M1 \7 k8 `6 e/ D' E3 |  _; _
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. * N' \& ?- S) |9 D
32.   k++
    
33. # `+ v7 |6 {7 x
34. },& R2 C( H2 i6 e7 k6 @! a
    
36. s;
    
37. + T/ B$ ^% V2 C9 H& J9 ^3 c
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">2 Z  G2 x* j6 s\\" |* u</font><br /> mvar:<font class=\\"jammer\\">\\" M; ~. N7 e& D& C( X: c8 _) v</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">2 n9 F( S9 `1 U+ G* T/ i</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">* v7 j3 F1 v4 r) w$ |</span> t0=clock(),<br /> <span style=\\"display:none\\">; ]3 z0 b# E. O5 h7 r7 m: h+ d</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">4 ]8 E# N* g4 a</span>   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">( l9 X6 i) d, l</span>   k++<br /> <span style=\\"display:none\\">( [. V9 j% O! L</span> },<br /> <span style=\\"display:none\\">: Y' |# u0 x) ~6 z) }</span> [clock()-t0]/1000;<font class=\\"jammer\\">0 y  j: I5 r7 b& ]: d</font><br /> <br /> <span style=\\"display:none\\">\\" H9 h; ~1 U2 m. x0 i( U3 D</span> main1(::a)=<font class=\\"jammer\\">6 \8 }  s( u! n( d' J/ }4 h! w: H: ]</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">* Q- h3 d9 F& I+ D</span> k=0, while{k<10000,<font class=\\"jammer\\"># ]; M1 \7 k8 `6 e/ D' E3 |  _; _</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">* N' \& ?- S) |9 D</span>   k++<br /> <span style=\\"display:none\\"># `+ v7 |6 {7 x</span> },<font class=\\"jammer\\">& R2 C( H2 i6 e7 k6 @! a</font><br /> s;<br /> <span style=\\"display:none\\">+ T/ B$ ^% V2 C9 H& J9 ^3 c</span> [clock()-t0]/1000;

结果：
10.469
- H7 k6 g& g) g' X8 F/ P20000000.
* w2 X: E) i7 c: A0 e* x1 H9 T' @3.781
: n/ y9 l2 p2 B3 m6 {5 M
－－－－－－－

1 u, t( ]" N* p9 [' E- Q1 {/ M7 \在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
# j+ K* [/ s/ Z8 W) q* y) N在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
. f8 D/ q$ I0 o5 |' M! Y) c# g
  g$ d! ~" n4 r$ Y2 u. Imatlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
% [6 b/ H( u) H2 v
! ?; A: ]7 ^% F' a  a+ K. @Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
) j$ l3 ~4 H" X
" e# L+ Q% |5 F6 d! b$ J  z' B5、在循环中加点东西

) j3 v* |0 a$ @1 @0 ~matlab代码：

1. >> clear all
   3 `\" z8 U' \3 ^  b6 [
2. f=@(x)x+1;
   9 T8 Z( I( u4 y- N: ^3 @$ C
3. a = linspace(2,2,1000);
   ! K9 s4 W2 N7 z, |
4. tic% o! P- i9 d2 ]& G: Z( ?! Y/ }
   
5. s=0;. S5 i\" [$ Y6 e; ?, `9 O1 P. Q0 d
   
6. for k=1:10009 {  K\" B( V: c3 i4 {. V4 O
   
7.   for i=1:1000
   5 N% ~4 u9 p\" X
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));! C4 E( n: a9 x$ S  h
   
9.   end5 g: o: O( V2 Q2 ~6 ^  p
   
10. end% ~4 R. W, z3 r4 c8 X& H
    
11. s7 b5 o, G  }! `! f$ N3 Q
    
12. toc0 c) H5 ]: v+ _8 P2 C  k/ X
    
13. 0 n) W. I5 `# R9 B$ f
    
14. s =& l  N& N5 @7 e7 _9 }# \
    
15. 
    2 M% f, ?7 W$ \! _, l6 j
16.   2.4002e+007
    1 _* W7 B1 x' X0 t! W
17. 0 G4 ~- E3 V% |$ U1 T+ j; B. n
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. \\" g' ~. k\\" [. T9 w. B
4. mvar:
   
5. . f% l1 m( J\\" M
6. f(x)=x+1;
   
7. ' [7 m7 D% {( M
8. main0(::a)=, T. u- K( |6 K- O\\" t
   
10. a = linspace(2,2,1000),
    
11. 6 \  `6 ]: s4 e. K/ p. i
12. t0=clock(), s=0,$ h: H% O) B\\" d0 ?, b* O) \
    
14. k=0, while{k<1000,
    
15. : `/ D- r\\" g# [, O6 f# ]4 o* Q
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++},8 [  x# t# t- B4 ~2 A
    
18.   k++
    
19. * G  c+ {+ T8 z8 e# v
20. },
    
21. \\" K2 U* L. h6 n$ P
22. s;, ]\\" n( |1 Z# C6 @- u* A; {
    
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">\\" g' ~. k\\" [. T9 w. B</span> mvar:<br /> <span style=\\"display:none\\">. f% l1 m( J\\" M</span> f(x)=x+1;<br /> <span style=\\"display:none\\">' [7 m7 D% {( M</span> main0(::a)=<font class=\\"jammer\\">, T. u- K( |6 K- O\\" t</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">6 \  `6 ]: s4 e. K/ p. i</span> t0=clock(), s=0,<font class=\\"jammer\\">$ h: H% O) B\\" d0 ?, b* O) \</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">: `/ D- r\\" g# [, O6 f# ]4 o* Q</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\\">8 [  x# t# t- B4 ~2 A</font><br />   k++<br /> <span style=\\"display:none\\">* G  c+ {+ T8 z8 e# v</span> },<br /> <span style=\\"display:none\\">\\" K2 U* L. h6 n$ P</span> s;<font class=\\"jammer\\">, ]\\" n( |1 Z# C6 @- u* A; {</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
9 [* r$ Q1 Z5 {. s/ x9 m5.172

6 o9 U/ V/ l# Z$ t+ p6 {4 \+ n二者效率现在只有一倍差距。

－－－－－－－－－－

6 d+ r0 X& h& v+ H( U不过，显然matlab和Forcal的程序应该这样写：
/ M$ `7 X; q3 l  n- F% Y8 g: L
matlab程序：

1. clear all9 ^6 K% v2 {/ b/ t' f+ {2 L8 L
   
2. f=@(x)x+1;
   : v- }  q7 [) b* R6 G! n# P/ a
3. a = linspace(2,2,1000);
   6 Q! V  _9 q2 c3 X0 t
4. tic
   & H6 N1 m! O$ {! o+ b
5. s=0;! `$ l; E+ j% ~' I% {* C8 B5 v
   
6. for k=1:10000 ^2 B3 A* W9 u, \. |
   
7.   for i=1:10009 B$ \. W4 c9 A, F/ x\" \' g
   
8.     t=a(i);* O& z& v# O- d( Z- k\" Y4 l
   
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);& S1 c$ u* J9 W0 E: W3 }
   
10.   end
    ' |0 V2 |1 E, E, e8 @) {, i% y6 U
11. end3 w+ B! x9 a6 J$ @% C) `/ `
    
12. s
    $ ~  C! y9 ~+ T/ c$ |5 G9 B/ r
13. toc
    \" Q9 w7 @7 z; k+ @
14. 
    : ^. F8 `5 G( D* b2 [, o
15. s =+ w  {5 @7 V  M2 f' W( S; x
    
16. 
    / {) g  p8 v: U! I9 ?
17.   2.4002e+007
    0 ]: `( h% T0 ~
18. 0 W; D9 X2 w9 ?; m6 `; i
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];# o& q! N6 q4 l4 ~
   
4. mvar:' ~' j1 U9 F8 l6 A! n; }0 }
   
6. f(x)=x+1;. P8 u1 z/ w  [4 P9 e
   
8. main0(::a)=/ r* w& h) R* m) J
   
10. a = linspace(2,2,1000),) p; r4 @$ D) ^# t/ l1 F
    
12. t0=clock(), s=0,( m' f! [& h9 _0 r! [
    
14. k=0, while{k<1000,- C: m+ T1 j% f$ c4 S8 n; p4 c
    
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++},  A* p5 P) v! D( m% d% c
    
18.   k++' _% j/ u) ?1 f- V1 K
    
20. },
    
21. 8 M- `+ y  x9 N) p5 c; o
22. s;) O- @! d( p7 y# \
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\"># o& q! N6 q4 l4 ~</font><br /> mvar:<font class=\\"jammer\\">' ~' j1 U9 F8 l6 A! n; }0 }</font><br /> f(x)=x+1;<font class=\\"jammer\\">. P8 u1 z/ w  [4 P9 e</font><br /> main0(::a)=<font class=\\"jammer\\">/ r* w& h) R* m) J</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">) p; r4 @$ D) ^# t/ l1 F</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">( m' f! [& h9 _0 r! [</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">- C: m+ T1 j% f$ c4 S8 n; p4 c</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\\">  A* p5 P) v! D( m% d% c</font><br />   k++<font class=\\"jammer\\">' _% j/ u) ?1 f- V1 K</font><br /> },<br /> <span style=\\"display:none\\">8 M- `+ y  x9 N) p5 c; o</span> s;<font class=\\"jammer\\">) O- @! d( p7 y# \</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
* }& d; A9 j3 W- z* P" C& J; [1.625
% \7 i+ f# }+ i+ Q$ Z
matlab效率已不及Forcal。

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

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
- M! ^0 J% ?- \# A
  {4 e& i7 w7 N0 EForcal代码：

2. !using["sys"];+ q3 _; M& o5 h( \, c& M
   
4. mvar:
   
5. \\" Z& Q/ a$ O# c
6. a="\&\[1000]"/4,    //用字符串模拟静态数组1 H7 U2 `/ S% M; Q! |- J
   
8. t0=clock(),
   
9. 8 r/ C' [0 z3 p- ~  D) K
10. k=0, while{k<10000,
    
11. ! `: m8 Y3 r$ k+ s0 e4 k3 ^
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值\\" Z, E3 b# z0 I  r
    
14.   k++
    
15. # [) g, z8 U3 {  f9 t* U3 N
16. },
    
17. 5 ]8 d( _6 {% K7 q% u8 S\\" b
18. [clock()-t0]/1000;3 a$ M  x  P, e- k/ o! H+ n8 f
    
20. - A$ I+ Y; t7 R1 P) S
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组7 y9 u5 T% c5 y# p( Z0 M0 G9 n
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. 7 L& C+ U% s/ X, L$ z/ q
26. t0=clock(),
    
27. . j1 J9 `% c& l  s1 `
28. s=0,
    
29. : |! h7 G( i# O' }' [, R) l$ W
30. k=0, while{k<10000,& x( U2 i) P% O. k0 [# H
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. 5 a1 }\\" q' Q/ q8 X8 e6 [
34.   k++$ `1 C' y& y, ^# q( u, B$ g* ^; X
    
36. },
    
37. $ ]4 Q4 c8 j1 s/ ?: G  `3 ?7 i
38. s;
    
39. 9 ^5 Q5 {* V6 ]* s/ a1 s5 z
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">+ q3 _; M& o5 h( \, c& M</font><br /> mvar:<br /> <span style=\\"display:none\\">\\" Z& Q/ a$ O# c</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">1 H7 U2 `/ S% M; Q! |- J</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">8 r/ C' [0 z3 p- ~  D) K</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">! `: m8 Y3 r$ k+ s0 e4 k3 ^</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">\\" Z, E3 b# z0 I  r</font><br />   k++<br /> <span style=\\"display:none\\"># [) g, z8 U3 {  f9 t* U3 N</span> },<br /> <span style=\\"display:none\\">5 ]8 d( _6 {% K7 q% u8 S\\" b</span> [clock()-t0]/1000;<font class=\\"jammer\\">3 a$ M  x  P, e- k/ o! H+ n8 f</font><br /> <font class=\\"jammer\\">- A$ I+ Y; t7 R1 P) S</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">7 y9 u5 T% c5 y# p( Z0 M0 G9 n</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">7 L& C+ U% s/ X, L$ z/ q</span> t0=clock(),<br /> <span style=\\"display:none\\">. j1 J9 `% c& l  s1 `</span> s=0,<br /> <span style=\\"display:none\\">: |! h7 G( i# O' }' [, R) l$ W</span> k=0, while{k<10000,<font class=\\"jammer\\">& x( U2 i) P% O. k0 [# H</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">5 a1 }\\" q' Q/ q8 X8 e6 [</span>   k++<font class=\\"jammer\\">$ `1 C' y& y, ^# q( u, B$ g* ^; X</font><br /> },<br /> <span style=\\"display:none\\">$ ]4 Q4 c8 j1 s/ ?: G  `3 ?7 i</span> s;<br /> <span style=\\"display:none\\">9 ^5 Q5 {* V6 ]* s/ a1 s5 z</span> [clock()-t0]/1000;

结果：
- C2 ^! @& y8 G) Z6 h  X0.5620000000000001  //赋值用时间，秒
20000000.
6 T! w4 B- p0 g0.579      //取值用时间，秒