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

forcal

1、简单循环

+ S) d" Y+ p  \. HMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all6 n/ u  L  w- `3 F: J
   
2. tic9 u+ e8 \& {. U( S
   
3. s=0;4 q6 K4 `- v9 T8 U7 ~
   
4. for k=1:10000
   / z& y; E0 h5 M/ c/ f* i8 @
5.   for i=1:1000: O2 M% P+ @5 O$ G& K
   
6.     s=s+1;
   . \. [, y& F\" k
7.   end
   7 l/ m( F* ?% ?8 a, M/ `7 l
8. end  h, q) m9 b9 s  n* Y# ~9 ?# l
   
9. s/ L. {/ S) Y8 q9 n' B( T
   
10. toc
    + s  z- Z, r8 o0 `/ E* d* S2 B
11. \" \; N  b- p3 F9 K
    
12. s =
    2 i0 @& s3 ~' L! R
13. 
    + d! |0 f4 Z. T\" p# t
14.     100000008 B: ~\" ]' s7 G7 U
    
15. - A% `  l/ H# w\" v9 F, U
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

; B. F; E/ C, ^- a& o8 M% AForcal没有JIT加速器，代码：

2. !using["sys"];
   
3. + P6 c/ q9 w: Y5 g* S2 v& q
4. mvar:. f& d; g9 z, B
   
6. t0=clock(),
   
7. ( ~: E( ]8 J' J
8. s=0,3 G# T$ y6 U' g3 P! B$ e
   
10. k=0, while{k<10000,2 Q% d: ~( n6 G7 J! `- l/ g- q
    
12.   i=0, while{i<1000,s=s+1, i++},
    
13. - U, a7 T9 ^$ q+ X+ U
14.   k++6 H, h5 Q& @* q4 C7 W
    
16. },4 W7 k5 ?  l0 d
    
18. s;
    
19. 1 A1 N  L\\" I& ?9 q
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">+ P6 c/ q9 w: Y5 g* S2 v& q</span> mvar:<font class=\\"jammer\\">. f& d; g9 z, B</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">( ~: E( ]8 J' J</span> s=0,<font class=\\"jammer\\">3 G# T$ y6 U' g3 P! B$ e</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">2 Q% d: ~( n6 G7 J! `- l/ g- q</font><br />   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">- U, a7 T9 ^$ q+ X+ U</span>   k++<font class=\\"jammer\\">6 H, h5 Q& @* q4 C7 W</font><br /> },<font class=\\"jammer\\">4 W7 k5 ?  l0 d</font><br /> s;<br /> <span style=\\"display:none\\">1 A1 N  L\\" I& ?9 q</span> [clock()-t0]/1000;

结果：
6 A7 I) N- Z" S1 S) u+ [) g10000000.
0.704  秒

& @2 W' i8 ^8 z6 d5 C二者运行效率相差7倍多。

2、循环中稍微加点东西

Matlab代码：

1. clear all: g  L9 a- L% f( N
   
2. tic- ?8 c! r4 G2 c6 O
   
3. s=0;: {1 ~% H\" {; [! R
   
4. for k=1:1000; D8 ~3 Z\" c& r; {\" r5 Q, X
   
5.   for i=1:1000' }& k/ s+ J3 n! q
   
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);, u( |. t. A3 H- F9 A
   
7.   end
   ! x* L8 K  e+ ^
8. end
   # _# Z' |3 q: n
9. s
   # K) s9 r5 M! G$ V1 H\" S
10. toc$ L% v; ~) R- _, s0 W, V' W
    
11. 
    ( r, D* Y6 B0 B- o
12. s =
    \" e) B+ j9 o\" b( j: e8 \
13. ' U) ^3 g- h8 T- E0 p# c& _
    
14.   1.1009e+007
    $ E/ U* O. Y# |. j2 m1 n
15. 
    * s9 H8 d0 x\" D/ Q% {
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. : @2 M/ V4 f9 x3 U' `
4. mvar:
   
5. ' _) q. z- J\\" K% c: b
6. t0=clock(),
   
7. : c: F2 a$ u- e- E/ }
8. s=0,
   
9. 1 j) y* E$ P( F1 u9 |
10. k=1, while{k<=1000,
    
11. ) [  G6 y8 _* D) |9 i$ ?
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++},0 g% J' Z' T6 O0 M0 m3 k
    
14.   k++' R( e$ o% g$ F
    
16. },( S6 D$ m( [  Z( g
    
18. s;
    
19. ! f4 v5 g7 l7 K) o- {5 d
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">: @2 M/ V4 f9 x3 U' `</span> mvar:<br /> <span style=\\"display:none\\">' _) q. z- J\\" K% c: b</span> t0=clock(),<br /> <span style=\\"display:none\\">: c: F2 a$ u- e- E/ }</span> s=0,<br /> <span style=\\"display:none\\">1 j) y* E$ P( F1 u9 |</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">) [  G6 y8 _* D) |9 i$ ?</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\\">0 g% J' Z' T6 O0 M0 m3 k</font><br />   k++<font class=\\"jammer\\">' R( e$ o% g$ F</font><br /> },<font class=\\"jammer\\">( S6 D$ m( [  Z( g</font><br /> s;<br /> <span style=\\"display:none\\">! f4 v5 g7 l7 K) o- {5 d</span> [clock()-t0]/1000;

结果：
11008953.66602346
1.031
$ U% h/ M0 G$ x0 v6 ~! y
. m; @- Y3 Q% g二者效率相差已比较小。
( G+ K3 e' n$ e/ l9 a3 X
3、再增加一个函数调用
+ ?0 W0 T9 S( {, c! T8 ?
Matlab代码：

1. clear all1 W/ `2 l! N9 U4 R
   
2. f=@(x)x+1;7 j$ Y8 p3 {) k, a9 c8 h
   
3. tic# O6 [5 Q! Q6 G6 N
   
4. s=0;
   ) ?2 b1 \; k, e7 @
5. for k=1:1000/ E# Q. e, u3 Q1 W3 O/ ~1 v5 t
   
6.   for i=1:1000, W, [- u+ ~' U
   
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);3 F) M( C6 Y/ M$ M, u$ Y
   
8.   end
   - a4 L! O' r+ a% b
9. end
   / L; u& I* y0 d( h+ e& x
10. s
    ( o3 G/ A( N; Y9 L, {
11. toc
    5 f) v: i/ P( i$ F/ J* k; k* \
12. 
    ' Q+ K! g6 Q4 X' p
13. s =% J\" v& _  ^6 O; z' M! D\" B2 d
    
14. 
    4 g* S! s) v% c\" ], A; |
15.   5.1251e+008
    9 D6 k- T. u5 q! ]: q
16. & s4 X( A8 c1 A* s4 i
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];* C+ k7 I( S' H+ W* t' \
   
4. f(x)=x+1;3 _* A% y- V6 i& x1 P# g2 z
   
6. mvar:; O6 }$ i9 J7 Y% ~
   
8. t0=clock(),2 V  O% ]2 h, A7 h) x# D2 D% F
   
10. s=0,+ Q' v+ ?7 u$ d' l# }7 B& n
    
12. k=1, while{k<=1000,\\" ?: ?) k! J/ j* L5 Y) 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++},
    
15. * e/ e; q3 f# H6 [4 D1 ^& M
16.   k++
    
17. 0 @$ `. [1 Z) L
18. },
    
19. 5 Y/ P1 `\\" w$ p+ P0 ^
20. s;' C1 r0 ~+ \. Y% j8 A
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">* C+ k7 I( S' H+ W* t' \</font><br /> f(x)=x+1;<font class=\\"jammer\\">3 _* A% y- V6 i& x1 P# g2 z</font><br /> mvar:<font class=\\"jammer\\">; O6 }$ i9 J7 Y% ~</font><br /> t0=clock(),<font class=\\"jammer\\">2 V  O% ]2 h, A7 h) x# D2 D% F</font><br /> s=0,<font class=\\"jammer\\">+ Q' v+ ?7 u$ d' l# }7 B& n</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">\\" ?: ?) k! J/ j* L5 Y) E</font><br />   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\\">* e/ e; q3 f# H6 [4 D1 ^& M</span>   k++<br /> <span style=\\"display:none\\">0 @$ `. [1 Z) L</span> },<br /> <span style=\\"display:none\\">5 Y/ P1 `\\" w$ p+ P0 ^</span> s;<font class=\\"jammer\\">' C1 r0 ~+ \. Y% j8 A</font><br /> [clock()-t0]/1000;

结果：
6 U& l" N5 w3 s0 C512508953.6658993
1.078
& m" d6 o! y  J$ g+ C, B
( c" b2 `3 V, |6 L+ P) O可以看出，Matlab效率已不及Forcal。
1 ]$ p4 p2 h0 n6 R/ b4 z故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
4 b2 m  c/ B3 z+ ~+ M
* r5 E  T+ t3 `/ k6 kMatlab代码：

1. >> clear all
   * W0 N% w2 j# D: V3 x- a
2. a = linspace(1,1,1000);
   ' [+ a  [2 e! ~* P3 `/ V% b
3. tic$ F& K& P2 C& D0 G
   
4. for k=1:10000& g, [3 r' I& ~- [$ c
   
5.   for i=1:1000
   + G7 A4 E: s; v, l
6.     a(i)=2;$ ?5 H$ F8 c4 Q0 \$ N& n
   
7.   end5 e# o& i- l  X/ N6 }) _
   
8. end* e$ f# Z! U5 X% b2 A' z/ f
   
9. toc
   , w, {; L* W: j8 P9 A- k\" M
10. ) C: f9 X# Z' }. L* Q- \! l# Z
    
11. tic
    0 u4 m3 T* Q# C% C
12. s=0;0 z6 s% _. l. P  Q. J8 E0 O& R
    
13. for k=1:10000
    ; q2 z; _5 \2 R9 V0 n6 Y. l
14.   for i=1:1000
    ' h. E, G4 S, ]: z
15.     s=s+a(i);
    ! {; g4 V  F0 @
16.   end
    1 |, n5 n9 R# r! P& S- p0 i: s
17. end
    , w1 J0 {  g( M7 Y$ x9 A! V
18. s
    - B& `. T% H: B: D; X+ b' A
19. toc9 c( Y) j0 h! G9 y
    
20. Elapsed time is 0.140658 seconds.
    3 Z\" ], w& M) I* T1 N- Q
21. 9 ]5 A4 s3 E& a# B$ b' b' S4 R/ |* t
    
22. s =0 C/ k\" v: F1 x* L\" @5 _& }* v
    
23. 5 h2 {; U7 J. m* O2 x
    
24.     20000000
    4 h% c2 b\" n( z+ u' c  L$ R1 l
25. , \0 V, A% i, ?! G
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

2. !using["math","sys"];
   
3. - j4 Y8 X  d6 |# v: j% r  D
4. mvar:0 s\\" g- o& C( N% [# {
   
6. main0(::a)=
   
7.   Y5 j% X4 N+ O7 j3 @( `1 c/ M
8. a = linspace(1,1,1000),
   
9. , A2 h  E/ L$ p! i! @- R
10. t0=clock(),
    
11. # N2 j3 A2 |/ T2 h; ~
12. k=0, while{k<10000,
    
13. 9 n  R. R0 u! q5 p! h
14.   i=0, while{i<1000, a(i)=2, i++},& i! i/ ?$ {/ m  g3 R
    
16.   k++
    
17. % Y5 P# z3 T  k
18. },) z! ^- P\\" s) ^! v  q
    
20. [clock()-t0]/1000;
    
21. 7 M# j; S/ ~5 l, u1 _8 k; a' a- Q
22. $ q9 u7 M8 Q7 s$ J# o1 [
    
24. main1(::a)=' J+ w* f+ U. _\\" P( C
    
26. t0=clock(), s=0,
    
27. . C# w! B1 K1 z2 ?
28. k=0, while{k<10000,
    
29. + F9 P6 A2 |6 x0 Y# o0 A4 Z
30.   i=0, while{i<1000, s=s+a(i), i++},
    
31. 2 U) a  X- J+ p, |
32.   k++
    
33. , b7 ~/ f7 w5 \
34. },/ ^# z  m8 n0 M! u% }) ]
    
36. s;: O1 a0 i. s0 d5 x! d$ G
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">- j4 Y8 X  d6 |# v: j% r  D</span> mvar:<font class=\\"jammer\\">0 s\\" g- o& C( N% [# {</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">  Y5 j% X4 N+ O7 j3 @( `1 c/ M</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">, A2 h  E/ L$ p! i! @- R</span> t0=clock(),<br /> <span style=\\"display:none\\"># N2 j3 A2 |/ T2 h; ~</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">9 n  R. R0 u! q5 p! h</span>   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">& i! i/ ?$ {/ m  g3 R</font><br />   k++<br /> <span style=\\"display:none\\">% Y5 P# z3 T  k</span> },<font class=\\"jammer\\">) z! ^- P\\" s) ^! v  q</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">7 M# j; S/ ~5 l, u1 _8 k; a' a- Q</span> <font class=\\"jammer\\">$ q9 u7 M8 Q7 s$ J# o1 [</font><br /> main1(::a)=<font class=\\"jammer\\">' J+ w* f+ U. _\\" P( C</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">. C# w! B1 K1 z2 ?</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">+ F9 P6 A2 |6 x0 Y# o0 A4 Z</span>   i=0, while{i<1000, s=s+a(i), i++},<br /> <span style=\\"display:none\\">2 U) a  X- J+ p, |</span>   k++<br /> <span style=\\"display:none\\">, b7 ~/ f7 w5 \</span> },<font class=\\"jammer\\">/ ^# z  m8 n0 M! u% }) ]</font><br /> s;<font class=\\"jammer\\">: O1 a0 i. s0 d5 x! d$ G</font><br /> [clock()-t0]/1000;

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

; N% X" c% W: ?) ?  G: e－－－－－－－－－－

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

2. !using["math","sys"];
   
3.   C! }8 ]\\" x5 Z& ~\\" N
4. mvar:
   
5. 3 r/ A$ R: F2 K
6. main0(::a)=5 {- G0 w% N3 d8 v9 ]; Q
   
8. a = linspace(1,1,1000),2 z  s9 W: S$ C% x9 y, _8 _
   
10. t0=clock(),
    
11. $ i* U1 I! P9 d6 i( ]
12. k=0, while{k<10000,
    
13. 2 x3 M* g6 `; @* _
14.   i=0, while{i<1000, A(a,i)=2, i++},2 Z$ [5 o; f, g
    
16.   k++
    
17. # Q\\" V8 i/ {. O$ s8 f& Z
18. },
    
19. ; T1 \) N% M1 u
20. [clock()-t0]/1000;
    
21. + A- t3 H1 K1 g+ z) a
22. 1 x; x5 e2 n1 A( i3 X* d- N6 i; N
    
24. main1(::a)=/ k$ C0 Q& m) n: U: v( z  D- q& S
    
26. t0=clock(), s=0,
    
27. ! L( V7 ?\\" l$ F& h
28. k=0, while{k<10000,1 n2 _5 S- K% v5 P2 n. e. r/ u
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},
    
31. * M. F7 u/ s7 _. g' i
32.   k++
    
33. , X; U0 t/ o: x+ o8 ]/ V
34. },
    
35. 0 i/ O- `( [1 R4 K6 E% u* W
36. s;
    
37. 3 m3 }; l) P8 a0 N! l
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">  C! }8 ]\\" x5 Z& ~\\" N</span> mvar:<br /> <span style=\\"display:none\\">3 r/ A$ R: F2 K</span> main0(::a)=<font class=\\"jammer\\">5 {- G0 w% N3 d8 v9 ]; Q</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">2 z  s9 W: S$ C% x9 y, _8 _</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">$ i* U1 I! P9 d6 i( ]</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">2 x3 M* g6 `; @* _</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">2 Z$ [5 o; f, g</font><br />   k++<br /> <span style=\\"display:none\\"># Q\\" V8 i/ {. O$ s8 f& Z</span> },<br /> <span style=\\"display:none\\">; T1 \) N% M1 u</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">+ A- t3 H1 K1 g+ z) a</span> <font class=\\"jammer\\">1 x; x5 e2 n1 A( i3 X* d- N6 i; N</font><br /> main1(::a)=<font class=\\"jammer\\">/ k$ C0 Q& m) n: U: v( z  D- q& S</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">! L( V7 ?\\" l$ F& h</span> k=0, while{k<10000,<font class=\\"jammer\\">1 n2 _5 S- K% v5 P2 n. e. r/ u</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<br /> <span style=\\"display:none\\">* M. F7 u/ s7 _. g' i</span>   k++<br /> <span style=\\"display:none\\">, X; U0 t/ o: x+ o8 ]/ V</span> },<br /> <span style=\\"display:none\\">0 i/ O- `( [1 R4 K6 E% u* W</span> s;<br /> <span style=\\"display:none\\">3 m3 }; l) P8 a0 N! l</span> [clock()-t0]/1000;

结果：
10.469
20000000.
3.781
9 m/ i7 G) t7 j! W
, I  V: {/ q8 X3 R－－－－－－－

( s- W3 n% q- U* z在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
: Z+ z: G' Y/ C; X4 p# m
5 f; _5 p) A3 `' {$ m9 _4 N' }6 k5 ^' Jmatlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
: l& H( Z/ g; _; [8 X- Q0 w$ f: h
5、在循环中加点东西

matlab代码：

1. >> clear all
   / [+ B# G' |+ c- y$ p
2. f=@(x)x+1;; Y5 x8 n' X+ j/ C$ U  X
   
3. a = linspace(2,2,1000);) I# i+ C& R9 p! d0 b4 S8 s\" J: b& m
   
4. tic9 W# D7 l  Y  H8 |& l; r% y+ f
   
5. s=0;
   , k6 y4 f8 F' _4 z' u
6. for k=1:1000, T8 F! F4 O9 }' Y$ w
   
7.   for i=1:1000  m3 J- k; A1 k- x$ N
   
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));
   + g1 g. ~& w% o! O3 x8 o
9.   end
   / t\" `2 F/ j( j8 v5 S7 v
10. end
    ) d0 G\" {3 e# y4 J
11. s3 U9 a! D( i5 s4 P7 ~! J7 X
    
12. toc1 U' S! p  F% k2 ^8 y' D
    
13. 6 X, D. T2 z* }3 d8 k  u\" [
    
14. s =$ Q! v( V+ ?* m3 L# m
    
15. 4 R! w) h3 P0 k, b3 n; U; S* H! z8 K
    
16.   2.4002e+007
    8 p1 D' p. @* O7 [- P: J( A
17. 8 Z) ^: [7 y7 h: C
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];  u6 U) f6 o% S# b3 ?
   
4. mvar:. W8 H5 L+ t& ?* a
   
6. f(x)=x+1;1 d/ W! p1 t. l, ?2 n3 k  x4 u
   
8. main0(::a)=5 X8 `2 I% u6 m8 ], Z5 u3 i6 L$ Q
   
10. a = linspace(2,2,1000),. c( e/ B5 `5 p1 ^* ?
    
12. t0=clock(), s=0,1 o$ u: i1 F! _' k$ h
    
14. k=0, while{k<1000,
    
15. 9 t, f2 J+ A2 ?; Q  D# P
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. - o9 h\\" @2 D/ B' X0 g( Q
18.   k++
    
19. / D- j2 L! I8 I4 P3 Y( J7 b
20. },2 q/ n% v3 b$ {( c, v3 q. q8 `- H' l# |
    
22. s;7 Q% `. Z1 j! T4 C( |0 B
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">  u6 U) f6 o% S# b3 ?</font><br /> mvar:<font class=\\"jammer\\">. W8 H5 L+ t& ?* a</font><br /> f(x)=x+1;<font class=\\"jammer\\">1 d/ W! p1 t. l, ?2 n3 k  x4 u</font><br /> main0(::a)=<font class=\\"jammer\\">5 X8 `2 I% u6 m8 ], Z5 u3 i6 L$ Q</font><br /> a = linspace(2,2,1000),<font class=\\"jammer\\">. c( e/ B5 `5 p1 ^* ?</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">1 o$ u: i1 F! _' k$ h</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">9 t, f2 J+ A2 ?; Q  D# P</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\\">- o9 h\\" @2 D/ B' X0 g( Q</span>   k++<br /> <span style=\\"display:none\\">/ D- j2 L! I8 I4 P3 Y( J7 b</span> },<font class=\\"jammer\\">2 q/ n% v3 b$ {( c, v3 q. q8 `- H' l# |</font><br /> s;<font class=\\"jammer\\">7 Q% `. Z1 j! T4 C( |0 B</font><br /> [clock()-t0]/1000;

结果：
24002271.69189384
% [# B, Q* ?# G# G+ ~/ y5.172
. A" x! y: P" `# [: |) ?' s
二者效率现在只有一倍差距。

' P6 O  u2 ^) @  w3 l－－－－－－－－－－
6 F7 r6 O, G8 z9 s* G, x
9 m/ M% o5 ]7 ?1 z" o! S不过，显然matlab和Forcal的程序应该这样写：
; a+ \2 C' J0 u2 }9 q
matlab程序：

1. clear all
   % u: H  P0 g5 m* w3 f5 x5 x* F9 a. b
2. f=@(x)x+1;
   4 ?, U7 ]& t$ X3 T
3. a = linspace(2,2,1000);- c+ ?# K\" D! r
   
4. tic0 d! R' D- \$ c1 \
   
5. s=0;% S# W4 ?* E1 K\" m/ W$ \7 Q4 ~\" V
   
6. for k=1:1000
   % Z) o4 A7 J% }. P5 a3 _$ d$ r5 p
7.   for i=1:1000
   6 D% q  w2 T# u6 f. K, I
8.     t=a(i);7 d; l3 E6 z: J9 ]& K& P, P. g
   
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);
   / m9 I* g: O, F3 D: C  X6 ~! x
10.   end
    - u; P# X) x/ _5 O2 A4 l
11. end* }2 n; [3 x' m, d* n
    
12. s
    7 m; b3 U+ |5 ^) X; C
13. toc2 R  p% l, Z# m& |5 T7 I0 o
    
14. # ^, f- Q8 V+ F
    
15. s =/ u1 r% |$ y' c! l2 ~: X
    
16. 
    ) Y5 `* p6 G  n
17.   2.4002e+0074 j' E! v1 P  X  ?
    
18. 
    6 Z% m. }7 L4 D0 q# y5 x. O1 f  a
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. 1 Q( N: Q) ~6 m  \% u' Y
4. mvar:
   
5. # Z& a\\" Z; `6 n/ u* c# l: I, ?
6. f(x)=x+1;6 N! G, i$ ]3 U! ~
   
8. main0(::a)=
   
9. 7 w; j: p6 i8 ^% n) E: E
10. a = linspace(2,2,1000),
    
11. - D+ ^: k4 c7 v3 j! ~; Z  v* [
12. t0=clock(), s=0,  M3 r# e- n7 U: Q' C2 ]* ?/ P+ S
    
14. k=0, while{k<1000,
    
15. ( z% V! b) k  C\\" E2 n
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. 4 t8 G! M0 A; ^
18.   k++
    
19. ) I- s) N' a+ K/ \1 I- o
20. },
    
21. - p\\" A8 p& n, d( O2 `
22. s;
    
23. , [$ Y$ G, W  \4 ~0 f
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">1 Q( N: Q) ~6 m  \% u' Y</span> mvar:<br /> <span style=\\"display:none\\"># Z& a\\" Z; `6 n/ u* c# l: I, ?</span> f(x)=x+1;<font class=\\"jammer\\">6 N! G, i$ ]3 U! ~</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">7 w; j: p6 i8 ^% n) E: E</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">- D+ ^: k4 c7 v3 j! ~; Z  v* [</span> t0=clock(), s=0,<font class=\\"jammer\\">  M3 r# e- n7 U: Q' C2 ]* ?/ P+ S</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">( z% V! b) k  C\\" E2 n</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++},<br /> <span style=\\"display:none\\">4 t8 G! M0 A; ^</span>   k++<br /> <span style=\\"display:none\\">) I- s) N' a+ K/ \1 I- o</span> },<br /> <span style=\\"display:none\\">- p\\" A8 p& n, d( O2 `</span> s;<br /> <span style=\\"display:none\\">, [$ Y$ G, W  \4 ~0 f</span> [clock()-t0]/1000;

结果：
4 ~$ {1 \9 X& y/ R5 \24002271.69189384
1.625
0 r# E* H# J- }9 K7 M' [
% l. h" p; h8 A  Umatlab效率已不及Forcal。
7 C# [/ N' u5 V+ y, C; P
提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
' a6 ^1 ~5 i' A1 m! I: p/ U
Forcal代码：

2. !using["sys"];
   
3. 3 q1 ^' M& L; i/ G3 B0 s/ Q
4. mvar:\\" ?; }( ]' ~/ j
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7.   p) S# ?2 S' j# _  c
8. t0=clock(),
   
9. * q  A* G( v8 w3 M. p. S1 U+ Y
10. k=0, while{k<10000,5 z! E- M6 X- W/ n; Z3 v
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值( M) J( H( _5 i+ W* C0 Q7 d
    
14.   k++
    
15. + v6 q( Q4 D! n
16. },
    
17. # d( e% v9 {4 d6 Z; h\\" m1 o: K
18. [clock()-t0]/1000;  c# ~# s  w. d% x* c9 e) u. x1 A5 l
    
20. 
    
21. - I$ E( |( G2 B5 K, I. }: @
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23.   I$ [9 G! ~% g& m
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值: Z/ X: M9 P! ^2 L, y  q+ m$ e
    
26. t0=clock(),\\" O2 n7 h4 y  |
    
28. s=0,
    
29. / t+ T5 |7 ~* _8 h' m
30. k=0, while{k<10000,
    
31. 2 ?5 o/ R8 p# [) U1 l
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值0 A$ P# S* b9 }. m2 y9 m
    
34.   k++
    
35. . o& X# |) E7 Q$ v* f
36. },
    
37. $ Z, @* C+ z+ N' ^! X# n0 `
38. s;
    
39. $ I8 v7 {/ k: h0 Y4 ]3 H& [
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">3 q1 ^' M& L; i/ G3 B0 s/ Q</span> mvar:<font class=\\"jammer\\">\\" ?; }( ]' ~/ j</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">  p) S# ?2 S' j# _  c</span> t0=clock(),<br /> <span style=\\"display:none\\">* q  A* G( v8 w3 M. p. S1 U+ Y</span> k=0, while{k<10000,<font class=\\"jammer\\">5 z! E- M6 X- W/ n; Z3 v</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">( M) J( H( _5 i+ W* C0 Q7 d</font><br />   k++<br /> <span style=\\"display:none\\">+ v6 q( Q4 D! n</span> },<br /> <span style=\\"display:none\\"># d( e% v9 {4 d6 Z; h\\" m1 o: K</span> [clock()-t0]/1000;<font class=\\"jammer\\">  c# ~# s  w. d% x* c9 e) u. x1 A5 l</font><br /> <br /> <span style=\\"display:none\\">- I$ E( |( G2 B5 K, I. }: @</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">  I$ [9 G! ~% g& m</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">: Z/ X: M9 P! ^2 L, y  q+ m$ e</font><br /> t0=clock(),<font class=\\"jammer\\">\\" O2 n7 h4 y  |</font><br /> s=0,<br /> <span style=\\"display:none\\">/ t+ T5 |7 ~* _8 h' m</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">2 ?5 o/ R8 p# [) U1 l</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">0 A$ P# S* b9 }. m2 y9 m</font><br />   k++<br /> <span style=\\"display:none\\">. o& X# |) E7 Q$ v* f</span> },<br /> <span style=\\"display:none\\">$ Z, @* C+ z+ N' ^! X# n0 `</span> s;<br /> <span style=\\"display:none\\">$ I8 v7 {/ k: h0 Y4 ]3 H& [</span> [clock()-t0]/1000;

结果：
' e) Y+ l3 c$ y5 }0.5620000000000001  //赋值用时间，秒
+ Q6 i/ j+ H! ^! K+ Y- l  A20000000.
0.579      //取值用时间，秒