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

forcal

1、简单循环
" a- f- p2 t9 n' K# k
Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all1 t3 [; ~- J3 [8 E- P
   
2. tic
   ' _0 s; v! v7 m+ T- A9 J
3. s=0;
   # u! i. F: M# Z3 G& ~9 V, {) C
4. for k=1:10000) p* m5 i9 q: r2 L9 g1 c7 X9 F
   
5.   for i=1:10001 U4 T\" z/ |0 B- f  ~
   
6.     s=s+1;7 D, c) W& G& R% ^3 a  P5 K9 g
   
7.   end
   ) `% _. z5 j) L& h- I9 {
8. end
   3 m8 o+ a$ v* I& W
9. s\" M2 ?3 Y1 ?6 G5 l( M) f- T\" W
   
10. toc
    4 H* {$ ?. A! C, K7 p9 k( j/ ~8 G
11. * B* p# _* A: c: W3 {+ x, M. L, w; g
    
12. s =
    1 s  h7 ^( S$ ~  j+ O0 b$ ^4 R
13. 4 p$ w+ H$ n8 Q8 z1 j! ?7 i
    
14.     10000000
    ' V- i5 m+ o9 B8 I
15. ) t* n( P  r7 y9 Y# V! o( Z
    
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
# n  J& `& \9 k( R2 \3 W. N
Forcal没有JIT加速器，代码：

2. !using["sys"];7 _9 J  G! ~2 x8 T- n9 Y
   
4. mvar:
   
5. - h( a7 K\\" l/ ]. d2 C3 M
6. t0=clock(),
   
7. ; ^( _) V& p! f- a4 F1 o8 c
8. s=0,
   
9. 5 k+ x+ d8 K( e, n' e4 ?: T2 T2 l
10. k=0, while{k<10000,2 t6 r# {2 F) }9 d
    
12.   i=0, while{i<1000,s=s+1, i++},. G3 d' Z% K2 z5 _3 f& N! [
    
14.   k++# m$ V! l$ N, c! T4 @0 [$ G! r
    
16. },
    
17. 5 K3 f/ a; y; P% z
18. s;
    
19. ( x5 C( U0 t\\" c6 f0 Q+ N7 m  A+ t+ A
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">7 _9 J  G! ~2 x8 T- n9 Y</font><br /> mvar:<br /> <span style=\\"display:none\\">- h( a7 K\\" l/ ]. d2 C3 M</span> t0=clock(),<br /> <span style=\\"display:none\\">; ^( _) V& p! f- a4 F1 o8 c</span> s=0,<br /> <span style=\\"display:none\\">5 k+ x+ d8 K( e, n' e4 ?: T2 T2 l</span> k=0, while{k<10000,<font class=\\"jammer\\">2 t6 r# {2 F) }9 d</font><br />   i=0, while{i<1000,s=s+1, i++},<font class=\\"jammer\\">. G3 d' Z% K2 z5 _3 f& N! [</font><br />   k++<font class=\\"jammer\\"># m$ V! l$ N, c! T4 @0 [$ G! r</font><br /> },<br /> <span style=\\"display:none\\">5 K3 f/ a; y; P% z</span> s;<br /> <span style=\\"display:none\\">( x5 C( U0 t\\" c6 f0 Q+ N7 m  A+ t+ A</span> [clock()-t0]/1000;

结果：
* _7 H4 l+ J3 G0 G10000000.
: Q" u1 ]% H2 h/ E6 l8 y5 s9 L0.704  秒

二者运行效率相差7倍多。
6 s. L" \" r. Y3 v
+ j9 Q$ k' E1 j1 y2、循环中稍微加点东西
5 j; }# S' `$ @2 `! q* |! R: m3 K
* Q$ ?* |5 Y1 ^7 E( V& g, mMatlab代码：

1. clear all
   1 X1 @\" C1 c# C
2. tic
   / g2 p: N% @  G. D' y$ P, l$ x
3. s=0;) n' J1 U! q$ S1 F
   
4. for k=1:1000- P2 v! B\" m3 N: C. Q
   
5.   for i=1:10000 k9 R, _) W% j
   
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);. c$ O1 }  _8 O- l# z/ ?
   
7.   end  p1 A4 {+ _8 O& q- [, N0 _
   
8. end; c% s8 ~5 e\" s- P
   
9. s\" n! ]# z3 n; v5 r3 l/ i1 `
   
10. toc
    0 G- F$ y4 L\" B% H( @! a7 h
11. 
    + Z% a& ^- U) \5 S# H% g
12. s =7 q$ S% u- t! x
    
13. 
    - m3 g- u9 C0 m& \% }# r
14.   1.1009e+007
    2 L* F) B* X# n$ {
15.   _0 H+ [9 D! D
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];5 s1 X# }4 o  `* [: h
   
4. mvar:/ ]4 a! m+ w5 A* |, E) ?$ C
   
6. t0=clock(),
   
7. ' `, ?3 f: P& i* U5 u# ~7 ?
8. s=0,
   
9. 5 ]+ H+ S' ]* j+ D0 D6 J
10. k=1, while{k<=1000,; F7 g) g* Z  v, r2 M* Z1 f
    
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. 2 k* I7 ^+ v4 B8 K$ L3 }. i: H
14.   k++
    
15.   C- m1 m- t2 ~+ k\\" R9 h
16. },, s5 A/ V% `' V. y7 V! P- Z5 m
    
18. s;; {% X3 a  p0 r\\" |
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">5 s1 X# }4 o  `* [: h</font><br /> mvar:<font class=\\"jammer\\">/ ]4 a! m+ w5 A* |, E) ?$ C</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">' `, ?3 f: P& i* U5 u# ~7 ?</span> s=0,<br /> <span style=\\"display:none\\">5 ]+ H+ S' ]* j+ D0 D6 J</span> k=1, while{k<=1000,<font class=\\"jammer\\">; F7 g) g* Z  v, r2 M* Z1 f</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\\">2 k* I7 ^+ v4 B8 K$ L3 }. i: H</span>   k++<br /> <span style=\\"display:none\\">  C- m1 m- t2 ~+ k\\" R9 h</span> },<font class=\\"jammer\\">, s5 A/ V% `' V. y7 V! P- Z5 m</font><br /> s;<font class=\\"jammer\\">; {% X3 a  p0 r\\" |</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
9 |6 v( \3 Q4 U! b1.031

二者效率相差已比较小。

# p  d. ~6 r1 P3 X6 g9 i3、再增加一个函数调用
8 X6 _' K/ J) c
Matlab代码：

1. clear all
   , v6 Z8 g# }2 L
2. f=@(x)x+1;2 E8 j! ~4 @) |+ Z
   
3. tic7 K  ~4 H/ \9 X\" W# K) K: @3 C
   
4. s=0;
   8 K, |& b! I; p
5. for k=1:10006 x4 z9 F$ G0 [! n( D4 t
   
6.   for i=1:1000
   8 B+ G/ p5 L2 J
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);6 O9 |  @0 L$ C# q5 T8 `! c
   
8.   end8 y4 K3 N' u5 ]
   
9. end6 T0 d' Q9 v! r0 ^* h) t7 Q. z5 C4 q
   
10. s\" Z$ ~2 ^' S; z5 Z
    
11. toc& C6 _\" h3 w; ?. N  W! w# M
    
12. 
    7 [4 Z- W3 w2 J/ r
13. s =8 q# x\" Z4 P# N4 @9 ]5 I/ F
    
14. 
      v- n( v2 x; R! A% x$ b
15.   5.1251e+008
    / y8 `+ g' ]; i0 I+ U
16. 
    7 m! p) U$ Z) {\" m/ n5 ?& y3 ?
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];1 g7 r! H+ O  U\\" c
   
4. f(x)=x+1;; P: {' G' Q. D
   
6. mvar:
   
7. ) \( j\\" O/ x' I0 I  r/ t
8. t0=clock(),
   
9. ' _! C; b8 y9 b
10. s=0,
    
11. . @4 H8 A* L9 v1 U5 S' r( g
12. k=1, while{k<=1000,1 L7 G1 s( W6 M
    
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++},; t) K$ F/ _9 Y- C
    
16.   k++$ t3 \3 p6 T\\" I* H; g% o
    
18. },
    
19. & E! t/ m# G! E2 T& r; |$ \+ B
20. s;  V, S* _& P& U\\" [5 m/ z! F
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">1 g7 r! H+ O  U\\" c</font><br /> f(x)=x+1;<font class=\\"jammer\\">; P: {' G' Q. D</font><br /> mvar:<br /> <span style=\\"display:none\\">) \( j\\" O/ x' I0 I  r/ t</span> t0=clock(),<br /> <span style=\\"display:none\\">' _! C; b8 y9 b</span> s=0,<br /> <span style=\\"display:none\\">. @4 H8 A* L9 v1 U5 S' r( g</span> k=1, while{k<=1000,<font class=\\"jammer\\">1 L7 G1 s( W6 M</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++},<font class=\\"jammer\\">; t) K$ F/ _9 Y- C</font><br />   k++<font class=\\"jammer\\">$ t3 \3 p6 T\\" I* H; g% o</font><br /> },<br /> <span style=\\"display:none\\">& E! t/ m# G! E2 T& r; |$ \+ B</span> s;<font class=\\"jammer\\">  V, S* _& P& U\\" [5 m/ z! F</font><br /> [clock()-t0]/1000;

结果：
/ G6 @9 R9 H$ M/ d512508953.6658993
1 R0 s" z( p5 X1.078
: u) r; f  H2 t) f# ~. z) p( n
可以看出，Matlab效率已不及Forcal。
( V) T, }" e2 f, r故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
! T; x* D7 |1 d& ?( f: ?& t
/ a( ~- J- [# T: t9 F# e2 B7 OMatlab代码：

1. >> clear all+ x: _; q% I( m; B- n+ b7 b$ O' u
   
2. a = linspace(1,1,1000);
   9 q2 W2 V% E* W
3. tic
     L- x- f: i( q0 F
4. for k=1:10000
   3 d3 ^6 E  `8 L# w& a
5.   for i=1:10003 B9 ]/ `1 y) q, c# l# V
   
6.     a(i)=2;
   7 \: g4 ^! M9 N! x  \. c
7.   end5 o' J, A6 n1 z5 a+ y
   
8. end
   ; q+ m$ B5 c$ J+ f' v
9. toc3 s  v4 c0 E+ T) s* Q# ]
   
10. ( P\" g6 B6 H( m7 }6 Y/ G  b9 ]
    
11. tic3 m& _7 y* _, _+ y' n, C
    
12. s=0;
      |! }: M' c4 k\" F8 {/ G
13. for k=1:10000
    1 V( `6 F4 C. R% c+ Y! t' D. w
14.   for i=1:10008 [: o3 T7 }% n5 s- D9 |4 ~  d
    
15.     s=s+a(i);/ s2 T& y\" B$ n9 r& N' a) `; M
    
16.   end
    3 q2 l7 M: [2 z, Y
17. end
    ! R% Z, y8 }% T+ J7 ?( G* o/ B
18. s8 i2 V) @4 S; g! u+ ~$ N( }\" i+ \
    
19. toc' ], Z! P& t+ g# }. `+ ]. M
    
20. Elapsed time is 0.140658 seconds.
    2 F) H\" t% k+ `\" d7 U
21. + A9 N7 Y! D+ X
    
22. s =
    * o% u2 c, |( X* z
23. \" z, J8 y* C$ s9 F3 w0 ]
    
24.     200000000 L) d5 ]3 h, z( l$ j
    
25. , ], g) x# ]- q& w\" j6 u  G
    
26. Elapsed time is 0.091794 seconds.

复制代码

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

! I' Z$ N$ M( t2 e3 r8 aForcal代码：

2. !using["math","sys"];
   
3. 5 I* j+ l4 S# V: {
4. mvar:
   
5. 8 o5 l8 ~- W/ Q- n+ [. Y% l* @' X
6. main0(::a)=
   
7. . O: u, C/ x4 @; i
8. a = linspace(1,1,1000),) e2 R* y# J6 ~\\" e8 b8 R1 F2 ?
   
10. t0=clock(),\\" C; b' s% g, G( m/ Q# o
    
12. k=0, while{k<10000,( d8 j9 ~- t3 k$ a
    
14.   i=0, while{i<1000, a(i)=2, i++},) u. v6 v+ A5 T0 Y
    
16.   k++
    
17. ( j\\" @4 u* c& W7 a9 z2 Z3 b
18. },
    
19. % G, q/ B* d3 T8 {
20. [clock()-t0]/1000;0 o# y) o+ Q1 C% N
    
22. 
    
23. 3 ]# Z/ z4 v% q, `- w& r
24. main1(::a)=
    
25. ( @5 ?9 T3 w1 m. X/ y9 b, ]+ g
26. t0=clock(), s=0,
    
27. & V+ }2 b; F7 E, [6 C8 F
28. k=0, while{k<10000,7 Q+ x( D! Z1 b8 \0 M* k
    
30.   i=0, while{i<1000, s=s+a(i), i++},\\" a& @: V/ U5 Q& U* ~9 o6 E  p- [# u
    
32.   k++
    
33. ( Y6 d. W; {9 S8 L: T, y4 S7 {5 B
34. },
    
35. 1 S) t6 w$ S2 ^+ W2 p
36. s;
    
37. 2 j+ H/ Y/ D& N8 ~& }) E- j\\" T1 _: T
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">5 I* j+ l4 S# V: {</span> mvar:<br /> <span style=\\"display:none\\">8 o5 l8 ~- W/ Q- n+ [. Y% l* @' X</span> main0(::a)=<br /> <span style=\\"display:none\\">. O: u, C/ x4 @; i</span> a = linspace(1,1,1000),<font class=\\"jammer\\">) e2 R* y# J6 ~\\" e8 b8 R1 F2 ?</font><br /> t0=clock(),<font class=\\"jammer\\">\\" C; b' s% g, G( m/ Q# o</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">( d8 j9 ~- t3 k$ a</font><br />   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">) u. v6 v+ A5 T0 Y</font><br />   k++<br /> <span style=\\"display:none\\">( j\\" @4 u* c& W7 a9 z2 Z3 b</span> },<br /> <span style=\\"display:none\\">% G, q/ B* d3 T8 {</span> [clock()-t0]/1000;<font class=\\"jammer\\">0 o# y) o+ Q1 C% N</font><br /> <br /> <span style=\\"display:none\\">3 ]# Z/ z4 v% q, `- w& r</span> main1(::a)=<br /> <span style=\\"display:none\\">( @5 ?9 T3 w1 m. X/ y9 b, ]+ g</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">& V+ }2 b; F7 E, [6 C8 F</span> k=0, while{k<10000,<font class=\\"jammer\\">7 Q+ x( D! Z1 b8 \0 M* k</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">\\" a& @: V/ U5 Q& U* ~9 o6 E  p- [# u</font><br />   k++<br /> <span style=\\"display:none\\">( Y6 d. W; {9 S8 L: T, y4 S7 {5 B</span> },<br /> <span style=\\"display:none\\">1 S) t6 w$ S2 ^+ W2 p</span> s;<br /> <span style=\\"display:none\\">2 j+ H/ Y/ D& N8 ~& }) E- j\\" T1 _: T</span> [clock()-t0]/1000;

结果：
& f# S  o# n" e15.484  //存数据所用时间，秒
20000000.
) _7 Q  L1 ~7 s" n- M; i5.531   //取数据所用时间，秒

0 t5 c* a* I$ b- o. _" x$ o－－－－－－－－－－
9 R, ]$ r) ^! Y) O
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];; u9 f+ \\\" f# n6 t9 B\\" ^5 z
   
4. mvar:
   
5. 0 e. Y3 `- m) @  a
6. main0(::a)=
   
7. ) ?5 h. R9 b# B: B% t/ {+ Y
8. a = linspace(1,1,1000),8 P' G/ J9 B# I+ S/ Z# A( k
   
10. t0=clock(),3 R3 Q3 o5 G( E+ k. S# k/ O9 @\\" i
    
12. k=0, while{k<10000,
    
13. 4 }3 S! l$ F3 f  T
14.   i=0, while{i<1000, A(a,i)=2, i++},+ a7 L3 i% i# h1 C' r( I
    
16.   k++; A5 n5 o  H6 L) E: `# w
    
18. },
    
19. * c; V; Y) J0 }: c
20. [clock()-t0]/1000;
    
21. * }, {. H6 _8 J& ^! ]
22. 
    
23. * p% V2 i% R8 [( Z; m. ]
24. main1(::a)=  C0 |3 Y7 n. [) I
    
26. t0=clock(), s=0,* G% G9 v% `( A# j3 T
    
28. k=0, while{k<10000,7 v! |& {' q1 ?
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},& ?5 b; v: i, e- e
    
32.   k++) k1 @4 f# G5 {; c2 @& ^
    
34. },, u- @* P! g% b+ B2 D
    
36. s;
    
37. 6 S' j( X' P* `: ^2 G2 n. Y
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">; u9 f+ \\\" f# n6 t9 B\\" ^5 z</font><br /> mvar:<br /> <span style=\\"display:none\\">0 e. Y3 `- m) @  a</span> main0(::a)=<br /> <span style=\\"display:none\\">) ?5 h. R9 b# B: B% t/ {+ Y</span> a = linspace(1,1,1000),<font class=\\"jammer\\">8 P' G/ J9 B# I+ S/ Z# A( k</font><br /> t0=clock(),<font class=\\"jammer\\">3 R3 Q3 o5 G( E+ k. S# k/ O9 @\\" i</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">4 }3 S! l$ F3 f  T</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">+ a7 L3 i% i# h1 C' r( I</font><br />   k++<font class=\\"jammer\\">; A5 n5 o  H6 L) E: `# w</font><br /> },<br /> <span style=\\"display:none\\">* c; V; Y) J0 }: c</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">* }, {. H6 _8 J& ^! ]</span> <br /> <span style=\\"display:none\\">* p% V2 i% R8 [( Z; m. ]</span> main1(::a)=<font class=\\"jammer\\">  C0 |3 Y7 n. [) I</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">* G% G9 v% `( A# j3 T</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">7 v! |& {' q1 ?</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">& ?5 b; v: i, e- e</font><br />   k++<font class=\\"jammer\\">) k1 @4 f# G5 {; c2 @& ^</font><br /> },<font class=\\"jammer\\">, u- @* P! g% b+ B2 D</font><br /> s;<br /> <span style=\\"display:none\\">6 S' j( X' P* `: ^2 G2 n. Y</span> [clock()-t0]/1000;

结果：
- ?5 ]; B  l- ?# r( z10.469
20000000.
8 q" I& q0 @5 ?8 U" n3.781

7 E1 C# c& K6 {( ~% ~8 _－－－－－－－

  l2 X! Z1 I7 H* q在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
# n  ]. x% i. S在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

9 E3 k+ @) u0 ^6 a! k2 i; Mmatlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
1 ]2 p7 ^3 Z5 `% {7 r& h. Q
9 @# @( U2 i5 `$ T, v5 ^* c5、在循环中加点东西

matlab代码：

1. >> clear all
   # J\" q: z) t7 S$ I7 ?! l
2. f=@(x)x+1;/ T+ j\" c) t\" F* w1 Q
   
3. a = linspace(2,2,1000);8 D1 l! b$ \7 {1 u- t8 ~3 z4 k
   
4. tic
   4 T0 y  l* W4 ?* p! m
5. s=0;, t# v/ d9 ]5 x! f6 W
   
6. for k=1:1000. w5 `7 D* t5 S, d
   
7.   for i=1:1000
   , |7 l* |: T1 p: 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));
   # `8 c+ U\" ^/ [5 Y6 S' ^
9.   end
   ! \: d: ~\" R6 F4 T$ b
10. end6 @# M3 `% d. A0 m( }' R
    
11. s
    / o7 y3 k1 B# D% p, q0 b
12. toc
    ' e: n. M& _2 B6 y3 J0 a( T  n
13. 
    3 b4 ]0 W+ s3 k% x8 T/ l( P$ z
14. s =% F- t' j% ~8 W- U. P) ^. g
    
15. / h% h$ ~5 x0 v3 M
    
16.   2.4002e+007
    / u! j+ L: P4 ?; C. i/ o6 e$ ~
17. & ~, z3 L; W! ~8 U+ c, S. w& {
    
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];% ?/ {$ v  B8 J# g* s! r. Z8 T: u  b5 Z
   
4. mvar:
   
5. 4 A' `0 M1 }4 E0 H5 W3 ~
6. f(x)=x+1;
   
7. # V- {, j; c4 k0 @/ B* c. J  t
8. main0(::a)=
   
9. * g: b: \5 _0 ]0 C0 `
10. a = linspace(2,2,1000),
    
11. $ ?3 T$ Y3 x1 ?) e
12. t0=clock(), s=0,
    
13. / c5 c' n8 ^7 i- M
14. k=0, while{k<1000,
    
15. % T# }4 Y7 K5 V9 b- {
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.   @% o( u' G# A/ H\\" z- H1 J
18.   k++
    
19. - I/ e, ?+ P- [8 g: S
20. },
    
21. . {! r- B  U' q! \; G% f
22. s;
    
23. # L4 T; k: V' e% Q( n
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">% ?/ {$ v  B8 J# g* s! r. Z8 T: u  b5 Z</font><br /> mvar:<br /> <span style=\\"display:none\\">4 A' `0 M1 }4 E0 H5 W3 ~</span> f(x)=x+1;<br /> <span style=\\"display:none\\"># V- {, j; c4 k0 @/ B* c. J  t</span> main0(::a)=<br /> <span style=\\"display:none\\">* g: b: \5 _0 ]0 C0 `</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">$ ?3 T$ Y3 x1 ?) e</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">/ c5 c' n8 ^7 i- M</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">% T# }4 Y7 K5 V9 b- {</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\\">  @% o( u' G# A/ H\\" z- H1 J</span>   k++<br /> <span style=\\"display:none\\">- I/ e, ?+ P- [8 g: S</span> },<br /> <span style=\\"display:none\\">. {! r- B  U' q! \; G% f</span> s;<br /> <span style=\\"display:none\\"># L4 T; k: V' e% Q( n</span> [clock()-t0]/1000;

结果：
24002271.69189384
5.172
2 W) q$ w: L/ y7 a2 ]! u
二者效率现在只有一倍差距。

8 ~6 g  e7 s6 c! }: _－－－－－－－－－－
: E) Y3 [6 i! B( x6 D: p
不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

1. clear all/ J9 G$ s$ O* S6 u/ Z8 p
   
2. f=@(x)x+1;
   + ?3 d. s# E* n
3. a = linspace(2,2,1000);\" S& L$ u' r  s; e\" w1 t
   
4. tic  l' }: v9 `: L' J3 P
   
5. s=0;
   \" _  T2 g: e! t3 x
6. for k=1:1000
   ( ]' {- Z3 k# C8 m6 n1 a4 A, ^
7.   for i=1:10000 _; L) S4 l9 O  t5 j
   
8.     t=a(i);
   2 Z, r( B. U! V( c: f* ?: u2 B! N
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);
   ) B1 X, L0 o7 w  ^
10.   end+ e5 i2 \( e\" y9 J
    
11. end
    / b& E# `) j0 W/ f# Y, b. g4 w
12. s
    : u7 _' Z1 Q) M: O
13. toc
    / w\" Y8 _\" P# \/ i) G
14. 
    % }6 `! L* @  _6 ?
15. s =
    \" i9 ~% r  `\" ]8 X9 e& I( d  N+ x
16. 
    ! l& D( @# k3 G' f8 P: \' b
17.   2.4002e+007: b9 U: n/ c5 P3 z. o3 K6 @) Q4 E
    
18. 
    1 S; P$ Q& p. W5 s3 }
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];3 u5 |% K\\" X: h& v- v3 t. x
   
4. mvar:
   
5. + m0 R\\" N1 K' y
6. f(x)=x+1;
   
7. : W/ A0 Q$ t\\" G0 o
8. main0(::a)=5 r: b4 K2 |) ?
   
10. a = linspace(2,2,1000),
    
11. % i; J8 ^; I5 @: U9 [8 E& R4 A
12. t0=clock(), s=0,  `3 `$ i9 Z5 Q* T! f
    
14. k=0, while{k<1000,
    
15. 2 D4 g# q; z* S  k9 u+ f' g+ o  O
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++},. l0 {9 [' S6 ]
    
18.   k++1 f; {( P# Y7 C9 d* z- P) |# y
    
20. },+ E5 I8 }, T7 J& {  F2 C0 t
    
22. s;
    
23. , q! Y* R1 h/ Q% m! J
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">3 u5 |% K\\" X: h& v- v3 t. x</font><br /> mvar:<br /> <span style=\\"display:none\\">+ m0 R\\" N1 K' y</span> f(x)=x+1;<br /> <span style=\\"display:none\\">: W/ A0 Q$ t\\" G0 o</span> main0(::a)=<font class=\\"jammer\\">5 r: b4 K2 |) ?</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">% i; J8 ^; I5 @: U9 [8 E& R4 A</span> t0=clock(), s=0,<font class=\\"jammer\\">  `3 `$ i9 Z5 Q* T! f</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">2 D4 g# q; z* S  k9 u+ f' g+ o  O</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++},<font class=\\"jammer\\">. l0 {9 [' S6 ]</font><br />   k++<font class=\\"jammer\\">1 f; {( P# Y7 C9 d* z- P) |# y</font><br /> },<font class=\\"jammer\\">+ E5 I8 }, T7 J& {  F2 C0 t</font><br /> s;<br /> <span style=\\"display:none\\">, q! Y* R1 h/ Q% m! J</span> [clock()-t0]/1000;

结果：
7 s8 Y0 S# ?! F9 u5 O24002271.69189384
. g/ R) o5 u" G+ D  Q1.625
( V4 o4 V6 n( `' G/ @7 ^
matlab效率已不及Forcal。

  \  ~# a0 @2 G8 U. ?, z  ~+ N) s( k; [提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

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

7 |/ M( [$ x+ H5 C( ~+ G# ?. `Forcal代码：

2. !using["sys"];
   
3. # D1 W/ {1 A+ H# j8 q3 [; `0 @
4. mvar:
   
5. ; a0 I( `: ]' Y
6. a="\&\[1000]"/4,    //用字符串模拟静态数组0 R/ F: R8 k6 u* E. T4 E7 @! f9 V$ _
   
8. t0=clock(),
   
9. - @* d- _/ i( W! b) p; k& C
10. k=0, while{k<10000,: G' C' I: B0 u  Q( @/ ]0 C
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值  T' D/ M, B$ A2 P( }6 Y4 L! k
    
14.   k++
    
15. ) v+ b- t) E. d# D\\" c, _
16. },
    
17.   A3 Z3 F; R7 ]: x, ?# V
18. [clock()-t0]/1000;
    
19. . V: I/ v# m* f4 Z7 L- p
20. 
    
21. ! X  G; c: C1 ]6 Z; Y4 N  `$ c
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. , r  H' q' w' Y' l  v# T$ X
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值8 W\\" l. E, T0 X5 ~: E/ ~4 B: g$ P
    
26. t0=clock(),3 P; o- N+ @' O& V0 I5 B' L1 `
    
28. s=0,
    
29. 8 O3 W- E3 |' y* w
30. k=0, while{k<10000,  U' J( B  X\\" l- C
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值! G+ o, h5 S( v9 j0 ~7 A! K
    
34.   k++! D8 G2 J9 R1 ^* Z
    
36. },\\" y+ n% h3 ^1 M' e
    
38. s;
    
39. ) q  d; }* N% E' ~- P
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\"># D1 W/ {1 A+ H# j8 q3 [; `0 @</span> mvar:<br /> <span style=\\"display:none\\">; a0 I( `: ]' Y</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">0 R/ F: R8 k6 u* E. T4 E7 @! f9 V$ _</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">- @* d- _/ i( W! b) p; k& C</span> k=0, while{k<10000,<font class=\\"jammer\\">: G' C' I: B0 u  Q( @/ ]0 C</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">  T' D/ M, B$ A2 P( }6 Y4 L! k</font><br />   k++<br /> <span style=\\"display:none\\">) v+ b- t) E. d# D\\" c, _</span> },<br /> <span style=\\"display:none\\">  A3 Z3 F; R7 ]: x, ?# V</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">. V: I/ v# m* f4 Z7 L- p</span> <br /> <span style=\\"display:none\\">! X  G; c: C1 ]6 Z; Y4 N  `$ c</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">, r  H' q' w' Y' l  v# T$ X</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">8 W\\" l. E, T0 X5 ~: E/ ~4 B: g$ P</font><br /> t0=clock(),<font class=\\"jammer\\">3 P; o- N+ @' O& V0 I5 B' L1 `</font><br /> s=0,<br /> <span style=\\"display:none\\">8 O3 W- E3 |' y* w</span> k=0, while{k<10000,<font class=\\"jammer\\">  U' J( B  X\\" l- C</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">! G+ o, h5 S( v9 j0 ~7 A! K</font><br />   k++<font class=\\"jammer\\">! D8 G2 J9 R1 ^* Z</font><br /> },<font class=\\"jammer\\">\\" y+ n% h3 ^1 M' e</font><br /> s;<br /> <span style=\\"display:none\\">) q  d; }* N% E' ~- P</span> [clock()-t0]/1000;

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