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

forcal

1、简单循环
- S" d8 I3 N- I  u0 k2 k* o2 n; K
1 G: j2 y( `0 nMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   ; V# T* M, i2 t9 L4 E
2. tic
   # F  n# e8 G% c6 J, ~1 H) M3 J
3. s=0;
   ' @6 B- Q: _1 V8 x' k6 j, b
4. for k=1:10000
   . d: g- f4 s\" o
5.   for i=1:1000
   & ]9 j( [1 f! u8 v
6.     s=s+1;
   + _0 ~# h. [+ {$ q% k\" Q. _
7.   end; c( g* G' W4 r$ O
   
8. end
   6 x' }8 r' k  @* h! v0 N5 V
9. s
   2 p: n+ U- B+ n\" d, d
10. toc
    : y; W/ ^: ?: @. g
11. 2 `8 O3 l2 D( z  s- P& S0 _! N1 O
    
12. s =0 P0 L& X7 X  Q7 m. }! w; T
    
13. 
    ! U- M, |& q: A% W; @
14.     10000000
    2 v( T5 X, n+ |( [6 G. A0 @
15. 8 g, x1 G* s( J, S\" C! U. G  y4 f8 @
    
16. Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. ) j) X' u% l/ v- ~\\" p
4. mvar:* A/ Z! |$ d+ l& R
   
6. t0=clock(),
   
7. 7 t& v+ L! t  _/ f: f7 v% \
8. s=0,
   
9. : w1 e* l. l\\" e( O% Y7 P* P
10. k=0, while{k<10000,
    
11. , f6 D/ R4 `, K1 P
12.   i=0, while{i<1000,s=s+1, i++},
    
13. 0 l: t6 g. [& j4 ]  g
14.   k++\\" t' l* e) s6 I: m6 z3 D  x- w
    
16. },. T: t) G0 \1 s1 \! g/ O
    
18. s;  }: k% O; w+ a- u, }/ G+ |
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">) j) X' u% l/ v- ~\\" p</span> mvar:<font class=\\"jammer\\">* A/ Z! |$ d+ l& R</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">7 t& v+ L! t  _/ f: f7 v% \</span> s=0,<br /> <span style=\\"display:none\\">: w1 e* l. l\\" e( O% Y7 P* P</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">, f6 D/ R4 `, K1 P</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">0 l: t6 g. [& j4 ]  g</span>   k++<font class=\\"jammer\\">\\" t' l* e) s6 I: m6 z3 D  x- w</font><br /> },<font class=\\"jammer\\">. T: t) G0 \1 s1 \! g/ O</font><br /> s;<font class=\\"jammer\\">  }: k% O; w+ a- u, }/ G+ |</font><br /> [clock()-t0]/1000;

结果：
* X- d4 n1 N! i) _9 H$ s10000000.
0.704  秒
4 V6 z% p2 }% k' t
二者运行效率相差7倍多。
/ {3 ]/ C% `0 P% ~0 r! O
( V; ^( |% t4 _$ R6 C2、循环中稍微加点东西
$ y4 w$ {! e, [" {, A  H* R
Matlab代码：

1. clear all5 l\" z6 F' p: l; X9 A, L- B
   
2. tic' _+ H6 G/ C; W  K. ?+ Z& S6 F& O
   
3. s=0;
   & X( b2 d/ Z9 S3 Q8 v
4. for k=1:1000
   * X( y& U1 \8 f* U6 s
5.   for i=1:1000
   5 {7 m5 v9 B/ q7 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);
   # D5 K( ?+ m\" O9 j+ w! I) ?
7.   end
   ( b) l3 r7 M$ i- P# C5 [0 v
8. end; z& n7 y) P+ C9 v5 E9 C- n8 i
   
9. s  w( E6 d7 l# n3 Q
   
10. toc# i1 I8 t1 W* ~7 H. z6 U
    
11. 
    ' x7 ?% i  }: ~2 a  J6 h3 p
12. s =3 t) G; b  N( ~' L- A
    
13. 
    0 G4 F+ [/ m/ O; g% b
14.   1.1009e+007* v! A& r. i' e$ O8 Y6 I
    
15. 
    7 D( y- z) p. _+ k- C- m
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];& X1 M, l3 |* _1 Y. E- A
   
4. mvar:
   
5.   y6 `3 n) i0 E7 M
6. t0=clock(),
   
7. * _( ?\\" d! Z: y3 J9 c4 |5 E: \
8. s=0,
   
9. 5 s* u* m5 M\\" p
10. k=1, while{k<=1000,
    
11. + ?8 w# c2 _9 Q; ?8 D) @; N, ]
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. , ?, ~) u9 u; c$ S
14.   k+++ X! |0 g. m# o# z. e+ n
    
16. },  P% t& G, }  Q7 i\\" Q
    
18. s;1 c6 {8 S, Q  Q$ {. L: R
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">& X1 M, l3 |* _1 Y. E- A</font><br /> mvar:<br /> <span style=\\"display:none\\">  y6 `3 n) i0 E7 M</span> t0=clock(),<br /> <span style=\\"display:none\\">* _( ?\\" d! Z: y3 J9 c4 |5 E: \</span> s=0,<br /> <span style=\\"display:none\\">5 s* u* m5 M\\" p</span> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">+ ?8 w# c2 _9 Q; ?8 D) @; N, ]</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\\">, ?, ~) u9 u; c$ S</span>   k++<font class=\\"jammer\\">+ X! |0 g. m# o# z. e+ n</font><br /> },<font class=\\"jammer\\">  P% t& G, }  Q7 i\\" Q</font><br /> s;<font class=\\"jammer\\">1 c6 {8 S, Q  Q$ {. L: R</font><br /> [clock()-t0]/1000;

结果：
! _" z" \8 C" ^& [; f7 k# ^11008953.66602346
% g5 E7 \# {  K, N4 p1.031
+ C: g% T" h% L( o# Y
& d) ?) s* y6 x, ?. M5 c二者效率相差已比较小。

3、再增加一个函数调用
8 N& a0 X$ o  q8 [( z1 w
8 D8 Y  |6 b8 n' o' vMatlab代码：

1. clear all: y, ^  n0 U8 q. E: v
   
2. f=@(x)x+1;
   3 f8 u; a* X, e  m* P! C
3. tic! p) s. y3 l5 j0 h2 u9 V* f6 C, v
   
4. s=0;
   ) q) p' K- \7 r# o
5. for k=1:1000
   $ t: E; }5 G; j$ N0 _- J. b$ \
6.   for i=1:10001 `  C9 d: Y! K7 \
   
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);/ e; ]( ?6 m. @  J\" i
   
8.   end
   6 D\" ]1 \7 j# O8 u3 |  Q; a
9. end
   : g6 \0 j+ V  R
10. s9 q, Y& s- f/ ]
    
11. toc% [4 X  _) I3 x$ a. x* |8 d
    
12. 
    / `4 w& b6 E5 _+ d4 M2 m! }$ e
13. s =! o0 f8 `: |$ u: z
    
14. : X, N9 }8 k. @
    
15.   5.1251e+008
    * @* x/ O7 A* I& `
16. ( t' E6 }0 V9 D( s' Q& y
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 0 }: E9 r& O; @/ m& l  Z
4. f(x)=x+1;
   
5. $ G$ S* Z' n: U% C7 W4 a5 c
6. mvar:
   
7. % j9 i! n9 O\\" u6 }. Y
8. t0=clock(),
   
9. ) I$ D; @6 A\\" J  L. }/ r4 X
10. s=0,
    
11. / Z$ t; V9 O+ Z\\" L) i! v$ V( F- q% Y
12. k=1, while{k<=1000,/ G3 q# J6 e. b% K1 A! ?
    
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++},# R2 ^  B* T' E0 I
    
16.   k++
    
17. ) {2 \0 r+ L2 ]6 J8 E
18. },& A5 u/ n1 X, j
    
20. s;# K( B# R) `  h; K9 F
    
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">0 }: E9 r& O; @/ m& l  Z</span> f(x)=x+1;<br /> <span style=\\"display:none\\">$ G$ S* Z' n: U% C7 W4 a5 c</span> mvar:<br /> <span style=\\"display:none\\">% j9 i! n9 O\\" u6 }. Y</span> t0=clock(),<br /> <span style=\\"display:none\\">) I$ D; @6 A\\" J  L. }/ r4 X</span> s=0,<br /> <span style=\\"display:none\\">/ Z$ t; V9 O+ Z\\" L) i! v$ V( F- q% Y</span> k=1, while{k<=1000,<font class=\\"jammer\\">/ G3 q# J6 e. b% K1 A! ?</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\\"># R2 ^  B* T' E0 I</font><br />   k++<br /> <span style=\\"display:none\\">) {2 \0 r+ L2 ]6 J8 E</span> },<font class=\\"jammer\\">& A5 u/ n1 X, j</font><br /> s;<font class=\\"jammer\\"># K( B# R) `  h; K9 F</font><br /> [clock()-t0]/1000;

结果：
1 u/ h- o* ^& _512508953.6658993
1.078

可以看出，Matlab效率已不及Forcal。
$ R) ^4 a- [6 L' |5 p% l# M故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
" s' q% p% I' |" j2 M
Matlab代码：

1. >> clear all
   6 m0 f/ ?, ]2 \5 v
2. a = linspace(1,1,1000);
   - }) [# P' ?; g* T! f0 n
3. tic
   . _* s8 S\" z% B6 v4 s+ J0 {
4. for k=1:10000
   . i7 f4 Q8 V! T1 l6 D4 q* D0 Z9 l
5.   for i=1:1000\" e3 n1 s/ s& }* t( R0 q% ?
   
6.     a(i)=2;
   ' M/ R  p& K' P# q
7.   end
   9 a1 `2 r1 z7 b% g1 c) p) ?: @
8. end* y. p% R, s5 }& t5 k4 G: }! D\" W
   
9. toc2 z$ A\" }6 ]. e3 B' g& ?+ K$ G
   
10. ( b5 U% s8 ?; G# B2 }1 D
    
11. tic
    9 n\" d2 p' T, ?* K% S
12. s=0;3 g' y- {* E* i! E% s/ T$ m
    
13. for k=1:10000) r1 p* r9 n8 G, L% V3 f2 J' c8 T
    
14.   for i=1:1000( u. B$ H# g\" s3 B' J4 K
    
15.     s=s+a(i);# v2 S1 }% ]$ j\" X. R5 ?1 Q
    
16.   end
      z% Y# J( c; ~8 a
17. end
    2 y3 d( `/ f! a\" [# G
18. s( y0 f3 O3 O' T$ E
    
19. toc3 ~  y* v/ R4 X! r5 K+ X
    
20. Elapsed time is 0.140658 seconds./ v& N/ M8 N( z+ ^7 r
    
21. 
    \" [7 H- v2 z7 P5 t6 h
22. s =/ N% x  H8 L5 h3 _/ ?! o$ P5 |5 N
    
23. 
    3 E* q5 ?\" o0 E1 E5 X& \6 P& N7 g
24.     20000000+ j6 ~# A  x* P1 G/ e5 F
    
25. . d& ^! w1 i$ \  s, M
    
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
9 g; B% j. I) j0 j7 G+ |  b
/ o, s8 Q: j& `- N9 rForcal代码：

2. !using["math","sys"];- \) q3 T) C- {; H; B+ P9 r' U
   
4. mvar:
   
5. 4 ~5 Z( q6 v/ L  c- }% D0 F
6. main0(::a)=6 d5 U3 \5 U5 j7 i1 r- J
   
8. a = linspace(1,1,1000),+ t3 v- B; [& m4 D* J4 \/ L
   
10. t0=clock(),
    
11. ) x/ }, q7 d; z; i0 x( T* F
12. k=0, while{k<10000,
    
13. \\" F, O\\" A; N( V1 h
14.   i=0, while{i<1000, a(i)=2, i++},1 T; k4 }\\" G* u8 G% C, ~
    
16.   k++5 m: l/ w. P) l0 ?7 m. ]4 f
    
18. },. S# b+ N4 T\\" X, o1 d! J' q$ k
    
20. [clock()-t0]/1000;, _, U, G3 z  D7 X\\" d* t) v
    
22. , ^; Q' q' S) h
    
24. main1(::a)=2 |8 P4 c8 d/ T9 x# E9 w( |
    
26. t0=clock(), s=0,
    
27.   N4 D; s1 K, ^  W: n( Y6 P
28. k=0, while{k<10000,
    
29. $ S, }% r  W, j0 _1 K, w
30.   i=0, while{i<1000, s=s+a(i), i++},( [& i' V: B6 a
    
32.   k++5 c7 x  U9 E. G- d6 W) s4 D
    
34. },6 v; x' {3 p1 D0 e) \; w
    
36. s;! G% I$ f0 c3 m\\" ?
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">- \) q3 T) C- {; H; B+ P9 r' U</font><br /> mvar:<br /> <span style=\\"display:none\\">4 ~5 Z( q6 v/ L  c- }% D0 F</span> main0(::a)=<font class=\\"jammer\\">6 d5 U3 \5 U5 j7 i1 r- J</font><br /> a = linspace(1,1,1000),<font class=\\"jammer\\">+ t3 v- B; [& m4 D* J4 \/ L</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">) x/ }, q7 d; z; i0 x( T* F</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">\\" F, O\\" A; N( V1 h</span>   i=0, while{i<1000, a(i)=2, i++},<font class=\\"jammer\\">1 T; k4 }\\" G* u8 G% C, ~</font><br />   k++<font class=\\"jammer\\">5 m: l/ w. P) l0 ?7 m. ]4 f</font><br /> },<font class=\\"jammer\\">. S# b+ N4 T\\" X, o1 d! J' q$ k</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">, _, U, G3 z  D7 X\\" d* t) v</font><br /> <font class=\\"jammer\\">, ^; Q' q' S) h</font><br /> main1(::a)=<font class=\\"jammer\\">2 |8 P4 c8 d/ T9 x# E9 w( |</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">  N4 D; s1 K, ^  W: n( Y6 P</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">$ S, }% r  W, j0 _1 K, w</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">( [& i' V: B6 a</font><br />   k++<font class=\\"jammer\\">5 c7 x  U9 E. G- d6 W) s4 D</font><br /> },<font class=\\"jammer\\">6 v; x' {3 p1 D0 e) \; w</font><br /> s;<font class=\\"jammer\\">! G% I$ f0 c3 m\\" ?</font><br /> [clock()-t0]/1000;

结果：
! |. I6 E& f" Y; i7 A" Y15.484  //存数据所用时间，秒
6 I4 p2 T1 Q( @! x, m20000000.
  n) A" i8 ?/ V2 o+ A5.531   //取数据所用时间，秒

－－－－－－－－－－
9 W: y- @# T; g5 u* {+ k  \$ D
用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];
   
3. , u. v8 m' ^. r. y# p! P
4. mvar:# h  `' {/ f7 D* J  z0 ^) g: ?
   
6. main0(::a)=3 H4 B+ ~/ ]- i5 T! y
   
8. a = linspace(1,1,1000),
   
9. 2 Y6 K( ^5 x4 t3 V* m
10. t0=clock(),
    
11. * m3 K. h. A$ A& t/ M; D
12. k=0, while{k<10000,
    
13. 5 \, c2 f% i+ ~# l4 p1 |: F
14.   i=0, while{i<1000, A(a,i)=2, i++},
    
15. 3 Q& q\\" e  C* G( C\\" ]. c3 s7 O
16.   k++
    
17. 2 r# t% t4 `! j\\" U! Y( F6 E4 H
18. },6 s) v( U) U( e
    
20. [clock()-t0]/1000;
    
21. 1 R8 Z- o# e9 v. L1 l8 ~
22. ! d0 u/ I  Z\\" P3 R; O7 h
    
24. main1(::a)=, t. t. X( l5 G; l7 l( k& d
    
26. t0=clock(), s=0,
    
27. 5 h' u& g\\" d# R4 [( T
28. k=0, while{k<10000,
    
29. . B; w* G( E7 {. \1 v# B! c- \: N7 S
30.   i=0, while{i<1000, s=s+A(a,i), i++},( U: Q: Y7 }- }% f
    
32.   k++
    
33. , |6 g) `! N$ H- @/ O7 Y
34. },/ ^\\" ]) O' X- Y8 m
    
36. s;
    
37. 2 ?2 X7 w0 g7 M9 K4 e: D\\" j
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">, u. v8 m' ^. r. y# p! P</span> mvar:<font class=\\"jammer\\"># h  `' {/ f7 D* J  z0 ^) g: ?</font><br /> main0(::a)=<font class=\\"jammer\\">3 H4 B+ ~/ ]- i5 T! y</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">2 Y6 K( ^5 x4 t3 V* m</span> t0=clock(),<br /> <span style=\\"display:none\\">* m3 K. h. A$ A& t/ M; D</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">5 \, c2 f% i+ ~# l4 p1 |: F</span>   i=0, while{i<1000, A(a,i)=2, i++},<br /> <span style=\\"display:none\\">3 Q& q\\" e  C* G( C\\" ]. c3 s7 O</span>   k++<br /> <span style=\\"display:none\\">2 r# t% t4 `! j\\" U! Y( F6 E4 H</span> },<font class=\\"jammer\\">6 s) v( U) U( e</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">1 R8 Z- o# e9 v. L1 l8 ~</span> <font class=\\"jammer\\">! d0 u/ I  Z\\" P3 R; O7 h</font><br /> main1(::a)=<font class=\\"jammer\\">, t. t. X( l5 G; l7 l( k& d</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">5 h' u& g\\" d# R4 [( T</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">. B; w* G( E7 {. \1 v# B! c- \: N7 S</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">( U: Q: Y7 }- }% f</font><br />   k++<br /> <span style=\\"display:none\\">, |6 g) `! N$ H- @/ O7 Y</span> },<font class=\\"jammer\\">/ ^\\" ]) O' X- Y8 m</font><br /> s;<br /> <span style=\\"display:none\\">2 ?2 X7 w0 g7 M9 K4 e: D\\" j</span> [clock()-t0]/1000;

结果：
$ Z& D0 h4 c$ b6 W) i/ a1 O10.469
8 W) n0 ^; t7 v6 w- @20000000.
; D9 p: |' W0 \# k3.781
! n7 h$ |, S- z. g: t$ u# W! |( ], S5 ~
－－－－－－－

在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
6 j5 e/ n) \3 x' P( C6 i7 F6 r在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

/ {5 D5 g3 f8 c1 XForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

9 `9 i3 ~* P2 {, k6 W& N5、在循环中加点东西
( a1 t( i9 }7 P7 C
0 ^7 \" A/ g- X$ u5 u6 m) Vmatlab代码：

1. >> clear all\" I. A3 ?  T2 G! v5 O9 b
   
2. f=@(x)x+1;  o$ v; A$ C0 X2 L, v
   
3. a = linspace(2,2,1000);
   * |3 [5 p/ ~$ {  ^
4. tic+ v/ E9 r! B* D$ c( {
   
5. s=0;7 P$ O. S6 U  m% O\" d  l+ J
   
6. for k=1:1000
   3 B5 ]' y6 Y# i4 c5 ]* T
7.   for i=1:10001 M' W8 T$ s\" C- V3 I4 y0 m! 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));% I2 }$ p0 U/ a: w: b0 Y
   
9.   end
   ; I' b1 i7 d0 D
10. end) B4 o\" Z* b% @/ e7 G
    
11. s* i3 \$ g+ Y0 i( }* O
    
12. toc
    % b9 n: w! O2 G3 I
13. - |' q% S+ E) O# d# Y+ L
    
14. s =& ]* {! b2 u: a9 G! z+ g
    
15. 9 a' P3 q6 y7 G6 a# \3 b1 z1 y
    
16.   2.4002e+007
    : c' v9 d\" K0 \% E/ S
17. 
    ; K9 `2 O. \0 i! L, _+ S
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];* s) y, ~/ i+ ]1 o
   
4. mvar:
   
5. 5 X- }4 F6 Z4 u
6. f(x)=x+1;6 ~' j8 W9 L; t; G2 e& d1 K
   
8. main0(::a)=
   
9. * @  }6 \9 V3 h2 t
10. a = linspace(2,2,1000),
    
11.   E; L$ i- |* w% C
12. t0=clock(), s=0,5 W8 v2 J. d5 C
    
14. k=0, while{k<1000,8 O) y! N; p. A' N
    
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++},6 I3 C, y: T5 ]
    
18.   k++3 I2 ~& D1 I2 m  d
    
20. },* n! V2 v! q( [; u. B/ }1 t' G) K) _
    
22. s;
    
23. ' c% K3 X) |- D( e
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">* s) y, ~/ i+ ]1 o</font><br /> mvar:<br /> <span style=\\"display:none\\">5 X- }4 F6 Z4 u</span> f(x)=x+1;<font class=\\"jammer\\">6 ~' j8 W9 L; t; G2 e& d1 K</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">* @  }6 \9 V3 h2 t</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">  E; L$ i- |* w% C</span> t0=clock(), s=0,<font class=\\"jammer\\">5 W8 v2 J. d5 C</font><br /> k=0, while{k<1000,<font class=\\"jammer\\">8 O) y! N; p. A' N</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\\">6 I3 C, y: T5 ]</font><br />   k++<font class=\\"jammer\\">3 I2 ~& D1 I2 m  d</font><br /> },<font class=\\"jammer\\">* n! V2 v! q( [; u. B/ }1 t' G) K) _</font><br /> s;<br /> <span style=\\"display:none\\">' c% K3 X) |- D( e</span> [clock()-t0]/1000;

结果：
24002271.69189384
/ d" K  A# [8 F1 c$ d8 ^3 i" U5.172
% U. X5 S( B0 S* v1 C4 `
二者效率现在只有一倍差距。
" ^& P+ W* J( u6 x
3 J$ J4 n- r% s－－－－－－－－－－

不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

1. clear all; G9 d* m! f% L) l5 m# M
   
2. f=@(x)x+1;
   9 y  J7 P0 g% z) T
3. a = linspace(2,2,1000);& l1 u: `4 _# \. C) b7 h) D
   
4. tic
   9 ~) X  w; X# k2 {
5. s=0;
   , v# V5 G2 z+ g% U
6. for k=1:1000
   ) e# K0 t, x. c( c$ ~
7.   for i=1:10001 N0 {$ _3 J! a9 b: r3 v4 j  v
   
8.     t=a(i);# F* K) e) \8 l# u% D- n; H+ C
   
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);
   ( u8 G' q0 U4 }8 K0 L
10.   end
    0 h% I- f. |/ J1 I* x4 B
11. end
    + G$ `7 t5 ~: W  z/ a
12. s
    $ r) Y- A7 E, ]* v: o5 Z& L
13. toc; Y* q, y) N! q7 W8 d5 \
    
14. 
    4 G1 [$ D# m+ w( G
15. s =
    % F* x$ ?  t1 d+ f( m% R% d
16. ; D/ z, C  R9 k
    
17.   2.4002e+007
    8 }) f7 o2 r4 S( T9 S7 L) G
18. 
    1 `9 a, e/ t% a3 ~! p+ P* \' I1 k
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];
   
3. 3 t8 V3 Z$ ?2 i' L
4. mvar:
   
5. 6 Q3 B  z+ Z: \* @% |6 q+ p
6. f(x)=x+1;
   
7. 6 N' |! b/ m7 \. h! I( _( m
8. main0(::a)=
   
9. ' t$ G, n( k5 f2 X, g1 D  p
10. a = linspace(2,2,1000),2 {: F5 D1 P9 K( ?/ R
    
12. t0=clock(), s=0,- S; v\\" D8 }0 j& k
    
14. k=0, while{k<1000,
    
15. 9 _; E0 H# _6 m( o! j# x
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. # X9 M5 }4 `) p. ]2 x
18.   k++
    
19. 7 Y5 _1 ?& o' V% E
20. },+ g7 i5 t* b; ~! t\\" m5 |% }
    
22. s;
    
23. , ^6 q9 R\\" Y/ j/ m; ?; p6 ~
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">3 t8 V3 Z$ ?2 i' L</span> mvar:<br /> <span style=\\"display:none\\">6 Q3 B  z+ Z: \* @% |6 q+ p</span> f(x)=x+1;<br /> <span style=\\"display:none\\">6 N' |! b/ m7 \. h! I( _( m</span> main0(::a)=<br /> <span style=\\"display:none\\">' t$ G, n( k5 f2 X, g1 D  p</span> a = linspace(2,2,1000),<font class=\\"jammer\\">2 {: F5 D1 P9 K( ?/ R</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">- S; v\\" D8 }0 j& k</font><br /> k=0, while{k<1000,<br /> <span style=\\"display:none\\">9 _; E0 H# _6 m( o! j# x</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\\"># X9 M5 }4 `) p. ]2 x</span>   k++<br /> <span style=\\"display:none\\">7 Y5 _1 ?& o' V% E</span> },<font class=\\"jammer\\">+ g7 i5 t* b; ~! t\\" m5 |% }</font><br /> s;<br /> <span style=\\"display:none\\">, ^6 q9 R\\" Y/ j/ m; ?; p6 ~</span> [clock()-t0]/1000;

结果：
24002271.69189384
8 `( ^8 {4 f  ~9 N1.625
, a2 H: }, t+ R) _% l3 e
matlab效率已不及Forcal。
# O) B$ q5 t3 v1 m! S  j$ B, A" w
提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
: P6 Q: \* J5 y, b6 I3 Z7 b
3 U! I5 \# \5 N) y$ uForcal代码：

2. !using["sys"];; k/ i8 s9 s( u' S8 a: \
   
4. mvar:
   
5. # H, B: V; u' o6 M% V. E
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. $ c- Z' G4 I. G9 d: ]
8. t0=clock(),# I; P5 s8 T4 Z& f# D! }4 |
   
10. k=0, while{k<10000,
    
11. * ~7 J* w/ S4 a. O* J: ^
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. 2 c8 @. o# h  I( g9 t6 v) ]' k
14.   k++  |2 R6 i* M& @( s6 O
    
16. },% c  u# g$ e; |& S# ^8 e
    
18. [clock()-t0]/1000;
    
19. 0 R2 ^8 v& ^1 n/ |1 P1 ^+ R* |
20. 
    
21. / y# C1 f5 X+ [; s
22. a="\&\[1000]"/4,    //用字符串模拟静态数组# i' p5 j- J0 g\\" c$ v2 Q
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值6 s6 e$ N8 L1 b5 _
    
26. t0=clock(),
    
27. 8 f6 }\\" J+ n# h  {, D. j. n' `/ g- k
28. s=0,
    
29. \\" U& f5 v$ p7 r' d/ d
30. k=0, while{k<10000,
    
31. 8 @\\" ^* h0 F) H& y\\" h. _) _
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. # f& ^$ w% v5 P: J
34.   k++: c\\" O# {1 W. C# Q
    
36. },\\" T% e% z- i0 Q1 X4 _
    
38. s;! D7 }5 o7 A0 v\\" a) h. o! m( I
    
40. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">; k/ i8 s9 s( u' S8 a: \</font><br /> mvar:<br /> <span style=\\"display:none\\"># H, B: V; u' o6 M% V. E</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">$ c- Z' G4 I. G9 d: ]</span> t0=clock(),<font class=\\"jammer\\"># I; P5 s8 T4 Z& f# D! }4 |</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">* ~7 J* w/ S4 a. O* J: ^</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">2 c8 @. o# h  I( g9 t6 v) ]' k</span>   k++<font class=\\"jammer\\">  |2 R6 i* M& @( s6 O</font><br /> },<font class=\\"jammer\\">% c  u# g$ e; |& S# ^8 e</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">0 R2 ^8 v& ^1 n/ |1 P1 ^+ R* |</span> <br /> <span style=\\"display:none\\">/ y# C1 f5 X+ [; s</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\"># i' p5 j- J0 g\\" c$ v2 Q</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<font class=\\"jammer\\">6 s6 e$ N8 L1 b5 _</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">8 f6 }\\" J+ n# h  {, D. j. n' `/ g- k</span> s=0,<br /> <span style=\\"display:none\\">\\" U& f5 v$ p7 r' d/ d</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">8 @\\" ^* h0 F) H& y\\" h. _) _</span>   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\"># f& ^$ w% v5 P: J</span>   k++<font class=\\"jammer\\">: c\\" O# {1 W. C# Q</font><br /> },<font class=\\"jammer\\">\\" T% e% z- i0 Q1 X4 _</font><br /> s;<font class=\\"jammer\\">! D7 }5 o7 A0 v\\" a) h. o! m( I</font><br /> [clock()-t0]/1000;

结果：
0.5620000000000001  //赋值用时间，秒
20000000.
5 Y: m6 }9 Q  u- w  m0.579      //取值用时间，秒