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

forcal

1、简单循环

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

1. clear all# b- A  g/ x$ A0 G6 E8 q# V
   
2. tic
   5 r\" p/ h) t/ H5 p, O; _: H9 l/ M8 D
3. s=0;) Y3 f4 r4 m0 w* N7 W
   
4. for k=1:10000
   / J- f5 [5 z& ^1 R8 L
5.   for i=1:1000
   5 I7 j! H+ D' e' ?\" U* C! x+ v
6.     s=s+1;
   - t/ b* T* B% X  ?* e
7.   end: T2 K8 B7 @$ t- H9 `, Q6 D7 M$ N
   
8. end
   3 M2 O+ ^3 O1 A* ~
9. s
   1 ^  N, d2 k5 T' ]2 [6 ^% m
10. toc
    ( x\" V0 X\" j) K$ P, a; E- o
11. 
    # B6 Q( l' n! V/ i  I
12. s =
    \" B. \9 f0 l8 R\" t& A
13. 
    2 A; _% K: s0 t* l: d8 R6 G
14.     100000007 {5 c0 H) H- q2 I3 {
    
15. 
    0 J$ ^  s& J+ i9 z0 ?- {\" q- z% ^
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
# u! x0 f0 J! t" i' k0 A
2 c+ y- p! @% c# ~Forcal没有JIT加速器，代码：

2. !using["sys"];  j5 g5 r; J9 z5 J8 H0 j7 I
   
4. mvar:6 Z' [4 f4 {2 ?0 D
   
6. t0=clock(),
   
7. 8 Y2 B$ A6 M# N, \$ g  t
8. s=0,
   
9. 2 m# i+ V$ v! p
10. k=0, while{k<10000,) M7 P; d- \. h, A- G
    
12.   i=0, while{i<1000,s=s+1, i++},
    
13. 2 I\\" x$ F, @: {- Z  h
14.   k++
    
15. ; c$ z4 Y\\" G' L2 R1 F) M- ?6 N
16. },, j' S9 y  f9 u' d\\" g2 d
    
18. s;! J% U) Y1 F) I% M) ~
    
20. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">  j5 g5 r; J9 z5 J8 H0 j7 I</font><br /> mvar:<font class=\\"jammer\\">6 Z' [4 f4 {2 ?0 D</font><br /> t0=clock(),<br /> <span style=\\"display:none\\">8 Y2 B$ A6 M# N, \$ g  t</span> s=0,<br /> <span style=\\"display:none\\">2 m# i+ V$ v! p</span> k=0, while{k<10000,<font class=\\"jammer\\">) M7 P; d- \. h, A- G</font><br />   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">2 I\\" x$ F, @: {- Z  h</span>   k++<br /> <span style=\\"display:none\\">; c$ z4 Y\\" G' L2 R1 F) M- ?6 N</span> },<font class=\\"jammer\\">, j' S9 y  f9 u' d\\" g2 d</font><br /> s;<font class=\\"jammer\\">! J% U) Y1 F) I% M) ~</font><br /> [clock()-t0]/1000;

结果：
' u- Q  e& P4 l" \" {/ `( R10000000.
0.704  秒

二者运行效率相差7倍多。
9 u) Z8 k# y% ]8 u3 R6 O
& f% G; P3 C0 p1 V$ I! v2、循环中稍微加点东西
8 P5 J( p( L5 [3 y/ E2 `# G
5 Y; H" k. m/ |Matlab代码：

1. clear all
   8 ]; T* f4 C( n+ E) D; U  b5 M2 ^
2. tic, V; U! K9 Y  S! }8 z
   
3. s=0;
   % N5 s2 `4 S1 g\" M3 w2 p* X
4. for k=1:1000
   * I% M4 h% L4 V/ U4 ]% e+ E# L
5.   for i=1:1000
   % U- ?. M, l3 K0 e
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);' k. K. x6 ?  y2 t+ Q: U
   
7.   end+ R, G$ i: y9 @, C/ o* t8 \
   
8. end; _  K# S5 a' v' f
   
9. s
   ( G5 {  i. j1 i! F+ T
10. toc% Y9 @' S% a8 `6 p
    
11. $ N. D6 F; k( F6 i4 w! {
    
12. s =
    ! P/ V8 [, Q) w- P, b& L( q& T7 O/ L
13. 
    ' [  a/ i- E0 @* R
14.   1.1009e+007
    7 A( ?  Q5 M, p8 g) \% |\" x# T
15. ) w( w% l7 C- c  y
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. 6 l  R  f% J/ d* ^4 l
4. mvar:- S, Y$ g( }+ u\\" ]
   
6. t0=clock(),/ l. ?0 H) w8 D+ Q' X
   
8. s=0,- Y& y$ }3 z2 B% B; K
   
10. k=1, while{k<=1000,
    
11.   c4 o\\" \\\" X5 Z+ K$ ^$ c
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. $ X\\" i) U5 V$ c* `! t+ e
14.   k++1 N+ z' s0 `  a) u3 x# Q) W7 f3 k
    
16. },\\" Z; }& V! r% c+ }
    
18. s;0 [1 j/ n- ]& h, w* m
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">6 l  R  f% J/ d* ^4 l</span> mvar:<font class=\\"jammer\\">- S, Y$ g( }+ u\\" ]</font><br /> t0=clock(),<font class=\\"jammer\\">/ l. ?0 H) w8 D+ Q' X</font><br /> s=0,<font class=\\"jammer\\">- Y& y$ }3 z2 B% B; K</font><br /> k=1, while{k<=1000,<br /> <span style=\\"display:none\\">  c4 o\\" \\\" X5 Z+ K$ ^$ c</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\\">$ X\\" i) U5 V$ c* `! t+ e</span>   k++<font class=\\"jammer\\">1 N+ z' s0 `  a) u3 x# Q) W7 f3 k</font><br /> },<font class=\\"jammer\\">\\" Z; }& V! r% c+ }</font><br /> s;<font class=\\"jammer\\">0 [1 j/ n- ]& h, w* m</font><br /> [clock()-t0]/1000;

结果：
11008953.66602346
- I$ @* \7 p  i3 v( e; ?8 x1.031

: y8 N3 w) w* D' k, \$ T% p二者效率相差已比较小。
+ {" `+ q- ?" b3 s/ I
3、再增加一个函数调用

Matlab代码：

1. clear all& V8 W2 v# k5 b( s, @3 J- N5 g( S
   
2. f=@(x)x+1;' N) P: L# ~2 s8 N: k
   
3. tic
   ) V2 b8 f8 b* c7 ], E1 s& o0 W\" Y5 n; x
4. s=0;  G9 V- ?0 O: N/ ?* Z, ^
   
5. for k=1:1000/ N+ Q& y4 z1 S' l/ x0 U& B
   
6.   for i=1:10000 _' m2 Y' g: w1 \: _
   
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 f/ O! ?! f9 S+ P7 Q. g9 ^. X
8.   end
   ( \\" Q+ W. V# a' ]5 a2 \
9. end' D  n2 ], j3 A- X2 D& Y+ X; B
   
10. s& G) m# w- X4 S& G) D4 d: F
    
11. toc. C; E- ~) Q5 C4 E# \9 N
    
12. 
    ( O6 v( l- Q$ j; r/ e: Z: g
13. s =
    $ m! `4 q/ T\" O, Q$ T3 x
14. 
    & n/ M9 Z6 n; Q- R2 P% I# Z( _. [# U
15.   5.1251e+008
    2 Y% S  v$ j4 U
16. 9 v0 V1 w: L  E
    
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];7 t& ?6 X& J6 c3 b, A
   
4. f(x)=x+1;+ v% l) {2 m+ s$ ^' ?  ~/ N. J
   
6. mvar:
   
7. ( S: c. D1 y  Q  V- ^2 p4 P5 C
8. t0=clock(),
   
9. 9 X: P8 a( D( K: Y% \0 r2 I
10. s=0,
    
11. ! J- r( A: q3 |/ z( o
12. k=1, while{k<=1000,6 Z* }, v6 O- o5 M' @* u( 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++},6 Y1 C6 h0 s  U
    
16.   k++
    
17. \\" k1 \! [. D5 o. L7 ^% T+ Y& X8 v3 |
18. },; Z4 f3 m. H) t1 D' |* z
    
20. s;' z/ u1 p$ s. K
    
22. [clock()-t0]/1000;

!using["sys"];<font class=\\"jammer\\">7 t& ?6 X& J6 c3 b, A</font><br /> f(x)=x+1;<font class=\\"jammer\\">+ v% l) {2 m+ s$ ^' ?  ~/ N. J</font><br /> mvar:<br /> <span style=\\"display:none\\">( S: c. D1 y  Q  V- ^2 p4 P5 C</span> t0=clock(),<br /> <span style=\\"display:none\\">9 X: P8 a( D( K: Y% \0 r2 I</span> s=0,<br /> <span style=\\"display:none\\">! J- r( A: q3 |/ z( o</span> k=1, while{k<=1000,<font class=\\"jammer\\">6 Z* }, v6 O- o5 M' @* u( b</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\\">6 Y1 C6 h0 s  U</font><br />   k++<br /> <span style=\\"display:none\\">\\" k1 \! [. D5 o. L7 ^% T+ Y& X8 v3 |</span> },<font class=\\"jammer\\">; Z4 f3 m. H) t1 D' |* z</font><br /> s;<font class=\\"jammer\\">' z/ u1 p$ s. K</font><br /> [clock()-t0]/1000;

结果：
+ B' j, W3 c% k- K; B512508953.6658993
1.078
$ X+ [6 C5 l6 G; Y
可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

4、在JIT加速器下，Matlab的数组存取效率非常高
% P1 }# B! J) g# W
Matlab代码：

1. >> clear all7 i8 W2 Y+ Z' W' B3 ~
   
2. a = linspace(1,1,1000);5 V1 r0 c* l4 |7 ~\" H) t
   
3. tic
   ( i7 r+ `& n( A0 _6 S
4. for k=1:10000. U\" n1 A3 P% g9 S/ c. h, F, A3 s: r
   
5.   for i=1:10009 ^. p  m\" v/ F
   
6.     a(i)=2;
   # y: L8 V% B  f4 n& N6 H! o
7.   end
   * p$ {. J* P2 i$ {: G3 z- S
8. end5 G, y1 p1 ?: \2 q  i3 \8 F8 Y
   
9. toc
   # e$ H2 o9 _9 p8 C) k
10. 
    ; `6 a8 P6 b, x9 ~# y
11. tic# o- ]9 x9 `2 a* j9 E
    
12. s=0;
    1 I& j' z7 s* G; E; |$ X
13. for k=1:10000
    * d- |/ H1 k0 i) F* |' {  [$ U) }
14.   for i=1:1000
    . p7 B1 s9 y# C  z* q- p
15.     s=s+a(i);% y5 e( Q\" v9 a+ b0 W
    
16.   end. B8 k1 K' T, x2 n; ]' [
    
17. end
    * v  D0 [' ?5 Z4 W6 H- m) v
18. s: ^2 G* h7 j0 \
    
19. toc6 ~( O0 P2 b6 Z, C# A/ H. Y  y
    
20. Elapsed time is 0.140658 seconds.
    5 m\" j& ]: ~# d4 k
21. 
    ; B4 o  n& ^+ \
22. s =
    1 m% ]0 \0 a! F0 D# _
23. 
    , e$ B: ?8 t+ C5 `3 c
24.     20000000* a+ G9 I( d2 {. c/ m
    
25. 8 y/ q; N4 H2 R' v  @7 Z
    
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
$ ?. N# e  h% k2 F, y! a; }" `  g
9 l+ I- p. Z7 c- l# L2 H8 dForcal代码：

2. !using["math","sys"];
   
3. \\" C$ P0 @1 I. `' C' R
4. mvar:7 b+ k7 d0 i$ o0 D\\" t* `
   
6. main0(::a)=2 M4 T( S5 j, O$ T4 y
   
8. a = linspace(1,1,1000),
   
9. # k- d9 A7 u\\" I. Y* _. J! R; j6 U
10. t0=clock(),6 ?# P; w/ W6 q+ h: z
    
12. k=0, while{k<10000,$ I- e* N. _6 B; y; T$ Z
    
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. % ^0 I1 X: \7 p8 |. J3 z
16.   k++
    
17. 7 l- S9 r4 m9 w, q# A9 w
18. },, V2 T' G2 t: \4 Z* a  Q! o+ T% @/ F
    
20. [clock()-t0]/1000;
    
21. \\" G8 Q( N\\" {8 ]2 L$ R
22. % x7 b' {$ h! u5 l7 b
    
24. main1(::a)=* F% q! j8 e8 ]  t& X( r
    
26. t0=clock(), s=0,, S5 H& `/ j! S: `+ u9 F' \
    
28. k=0, while{k<10000,
    
29. 2 n5 E( ]* d4 B% e
30.   i=0, while{i<1000, s=s+a(i), i++},: X; p+ {* Y# n2 O: x1 }
    
32.   k++8 Y6 i7 d( h8 b+ o
    
34. },
    
35. , s+ n$ h, x. _! f1 d4 K
36. s;4 f9 s* K9 O& q& @# n! Z/ I% y6 F
    
38. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">\\" C$ P0 @1 I. `' C' R</span> mvar:<font class=\\"jammer\\">7 b+ k7 d0 i$ o0 D\\" t* `</font><br /> main0(::a)=<font class=\\"jammer\\">2 M4 T( S5 j, O$ T4 y</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\"># k- d9 A7 u\\" I. Y* _. J! R; j6 U</span> t0=clock(),<font class=\\"jammer\\">6 ?# P; w/ W6 q+ h: z</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">$ I- e* N. _6 B; y; T$ Z</font><br />   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">% ^0 I1 X: \7 p8 |. J3 z</span>   k++<br /> <span style=\\"display:none\\">7 l- S9 r4 m9 w, q# A9 w</span> },<font class=\\"jammer\\">, V2 T' G2 t: \4 Z* a  Q! o+ T% @/ F</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">\\" G8 Q( N\\" {8 ]2 L$ R</span> <font class=\\"jammer\\">% x7 b' {$ h! u5 l7 b</font><br /> main1(::a)=<font class=\\"jammer\\">* F% q! j8 e8 ]  t& X( r</font><br /> t0=clock(), s=0,<font class=\\"jammer\\">, S5 H& `/ j! S: `+ u9 F' \</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">2 n5 E( ]* d4 B% e</span>   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">: X; p+ {* Y# n2 O: x1 }</font><br />   k++<font class=\\"jammer\\">8 Y6 i7 d( h8 b+ o</font><br /> },<br /> <span style=\\"display:none\\">, s+ n$ h, x. _! f1 d4 K</span> s;<font class=\\"jammer\\">4 f9 s* K9 O& q& @# n! Z/ I% y6 F</font><br /> [clock()-t0]/1000;

结果：
+ A* ~# N2 w# u! m7 A# s15.484  //存数据所用时间，秒
" A; t! P- j" M/ C- N7 U: G20000000.
5.531   //取数据所用时间，秒

7 E5 p8 ^7 S" E. N2 {－－－－－－－－－－
1 r/ S: c. k- B1 Q, p
: `* \5 U, j! |9 q8 y5 s* O用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];! G( C6 E2 x* x% z8 Y  i
   
4. mvar:
   
5. 5 k: b% I\\" u* @2 D/ H$ H
6. main0(::a)=
   
7. 4 T/ ^# O1 F, F8 u3 T/ \+ v9 q
8. a = linspace(1,1,1000),4 x6 [8 Z$ x$ l4 k$ I
   
10. t0=clock(),9 d+ I\\" f& V) R- G7 B
    
12. k=0, while{k<10000,
    
13. 2 S! J& ~4 T1 W' @/ u
14.   i=0, while{i<1000, A(a,i)=2, i++},4 |# j; o/ U  S, L% v
    
16.   k++
    
17. ! i* V8 M- k) G/ D
18. },
    
19. ) P: r0 f8 ^, T1 R
20. [clock()-t0]/1000;
    
21. 9 \* t% i: l! \! w' ]. V1 u6 w. b: i1 s. P
22. 6 g% Q; q\\" a4 K6 ]6 n4 n! i\\" }% v5 s
    
24. main1(::a)=: N7 r& x: Y9 b( T7 w  h
    
26. t0=clock(), s=0,
    
27. ' }1 Y4 y. v; }9 g$ e' V
28. k=0, while{k<10000,/ R) ^# `. f! c/ ]- ?
    
30.   i=0, while{i<1000, s=s+A(a,i), i++},+ t0 U! r3 s& w4 D4 [
    
32.   k++- H' O1 M7 r9 j5 ^
    
34. },( d3 a! [# `; Z3 Q9 m
    
36. s;
    
37. # H* a# f' ^  s( Y: d4 q9 X
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">! G( C6 E2 x* x% z8 Y  i</font><br /> mvar:<br /> <span style=\\"display:none\\">5 k: b% I\\" u* @2 D/ H$ H</span> main0(::a)=<br /> <span style=\\"display:none\\">4 T/ ^# O1 F, F8 u3 T/ \+ v9 q</span> a = linspace(1,1,1000),<font class=\\"jammer\\">4 x6 [8 Z$ x$ l4 k$ I</font><br /> t0=clock(),<font class=\\"jammer\\">9 d+ I\\" f& V) R- G7 B</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">2 S! J& ~4 T1 W' @/ u</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">4 |# j; o/ U  S, L% v</font><br />   k++<br /> <span style=\\"display:none\\">! i* V8 M- k) G/ D</span> },<br /> <span style=\\"display:none\\">) P: r0 f8 ^, T1 R</span> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">9 \* t% i: l! \! w' ]. V1 u6 w. b: i1 s. P</span> <font class=\\"jammer\\">6 g% Q; q\\" a4 K6 ]6 n4 n! i\\" }% v5 s</font><br /> main1(::a)=<font class=\\"jammer\\">: N7 r& x: Y9 b( T7 w  h</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">' }1 Y4 y. v; }9 g$ e' V</span> k=0, while{k<10000,<font class=\\"jammer\\">/ R) ^# `. f! c/ ]- ?</font><br />   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">+ t0 U! r3 s& w4 D4 [</font><br />   k++<font class=\\"jammer\\">- H' O1 M7 r9 j5 ^</font><br /> },<font class=\\"jammer\\">( d3 a! [# `; Z3 Q9 m</font><br /> s;<br /> <span style=\\"display:none\\"># H* a# f' ^  s( Y: d4 q9 X</span> [clock()-t0]/1000;

结果：
10.469
20000000.
  t3 P, |" {6 |9 _  Q( ~  q$ n3.781
6 f0 ~5 Q: |1 w9 L1 @' r
－－－－－－－
5 ^- B. j3 q1 M% H0 ~& z7 N
在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
; f) J% |/ E+ a, K& D4 u2 [在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

! R  |/ ~, ?$ s* |$ D; ?matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
. a2 b/ F: o$ k6 ?
  W( d: I* O! o: X6 lForcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
3 c6 ^# i1 n4 U2 T9 o1 {
0 K9 f+ u8 x! g# l% V  Q5、在循环中加点东西

matlab代码：

1. >> clear all
   7 ]% u# X  z  l. v' _! y/ t
2. f=@(x)x+1;0 C3 \7 i\" \6 q; |+ G2 G) w2 f0 v  R$ r
   
3. a = linspace(2,2,1000);
   & k1 a\" c6 O; Z' @1 `/ ?
4. tic
   # K5 g5 k3 V* r2 R' A8 n' R; W
5. s=0;
   1 G1 g% P/ ~5 |* Z/ x6 r
6. for k=1:10001 e& _3 B; H\" i! e6 P
   
7.   for i=1:1000
   6 B3 m! i/ `' t0 q6 r8 [
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));
   6 q+ o( |; }# U\" E
9.   end
   , k6 _+ M* ?; `4 p# m3 J
10. end5 u5 t; ^, I  s6 f) H7 p8 Z
    
11. s8 {. c7 x. l1 }  I/ T9 m
    
12. toc
    / Z# }2 I! [- w. E% T
13. 
    9 v+ @$ A9 T% f7 F7 Z' G
14. s =
    ( K: C0 M/ a) \) d
15. 
    0 x- f/ Q9 L\" n
16.   2.4002e+007
    \" N; N4 H+ w3 f& Q
17. 
    8 J8 t5 E  }# K2 w
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];
   
3. + O1 g* |5 s- i+ t
4. mvar:7 D; Z' h) d5 |5 x
   
6. f(x)=x+1;
   
7. , p. Y% {$ u; G7 m: I( W6 M+ k
8. main0(::a)=
   
9. * \1 E4 R' q4 l& ~
10. a = linspace(2,2,1000),
    
11. 7 s, W$ H% i. ~0 K* O: r% n# ^
12. t0=clock(), s=0,
    
13. ' R$ v% h7 P& a7 b2 O  T% u6 p
14. k=0, while{k<1000,# s3 c. v% g3 ~* l7 v! 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++},
    
17. 6 w; |\\" @1 y- ^0 x2 z
18.   k++
    
19. 0 m% V$ H* D# D. L/ {8 j
20. },: e6 o& s) s& y/ I7 b7 t; N5 ~; H# k
    
22. s;
    
23. 6 t! k1 i3 ?6 e
24. [clock()-t0]/1000;

!using["math","sys"];<br /> <span style=\\"display:none\\">+ O1 g* |5 s- i+ t</span> mvar:<font class=\\"jammer\\">7 D; Z' h) d5 |5 x</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">, p. Y% {$ u; G7 m: I( W6 M+ k</span> main0(::a)=<br /> <span style=\\"display:none\\">* \1 E4 R' q4 l& ~</span> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">7 s, W$ H% i. ~0 K* O: r% n# ^</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">' R$ v% h7 P& a7 b2 O  T% u6 p</span> k=0, while{k<1000,<font class=\\"jammer\\"># s3 c. v% g3 ~* l7 v! 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++},<br /> <span style=\\"display:none\\">6 w; |\\" @1 y- ^0 x2 z</span>   k++<br /> <span style=\\"display:none\\">0 m% V$ H* D# D. L/ {8 j</span> },<font class=\\"jammer\\">: e6 o& s) s& y/ I7 b7 t; N5 ~; H# k</font><br /> s;<br /> <span style=\\"display:none\\">6 t! k1 i3 ?6 e</span> [clock()-t0]/1000;

结果：
24002271.69189384
5.172
. p" t2 a$ ~) A/ s$ q
二者效率现在只有一倍差距。
: g: T& B. ^4 x* x1 T
1 r% Z" u: }% e5 L8 ?7 a2 v5 R－－－－－－－－－－

不过，显然matlab和Forcal的程序应该这样写：
% ~) A) c! K6 `& x7 K2 x
matlab程序：

1. clear all$ S, u3 c3 k; U) [
   
2. f=@(x)x+1;9 O8 k2 Z) S* k- p* _
   
3. a = linspace(2,2,1000);# k5 `* ]. \8 o, C# Y- K+ B$ N
   
4. tic
   ( l6 v2 A5 T5 X# v  J5 D5 t  X
5. s=0;
   + i, p3 v' x\" Q
6. for k=1:1000
   8 J4 v+ I+ Q( B
7.   for i=1:1000
   & K- N* Y2 ]: T  `/ y' i0 P! A
8.     t=a(i);) a. h+ L' Q5 [% Z0 P; l! y# s1 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);
   ( a% q: U  n, T1 O/ d( D
10.   end
    8 ]\" `, a8 q5 p- I' Y  z
11. end  I( H( ]- a! Z- X2 z0 Z$ C& j  t' M
    
12. s8 ~2 P( E5 Z7 ]4 W1 B0 c4 \5 u
    
13. toc
    6 _# _# c1 w/ ^: H( T$ W' n: }
14. ' u$ H) [) s3 D% C
    
15. s =: `9 F8 J9 k$ ?\" U
    
16. 
    7 X5 k0 @0 e; a/ \1 M
17.   2.4002e+0078 {& e3 I8 g+ O) s) E
    
18. 
    7 a, V, J- c( ~0 X; D8 ]
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];, P% n* p5 B  P* _/ f$ f* M
   
4. mvar:
   
5. : z8 [7 t# A7 S7 Z3 Q# M
6. f(x)=x+1;
   
7. 0 X5 _  ]1 M1 Y3 n
8. main0(::a)=* k% G9 C) F- t1 G7 Q
   
10. a = linspace(2,2,1000),
    
11. 4 |5 ^& c' [+ m# e/ l  Q
12. t0=clock(), s=0,
    
13. ( S. Z! [. j! a\\" q
14. k=0, while{k<1000,
    
15. \\" C) e- x' i1 v# Y
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++},. K6 c! q6 ?$ k: \
    
18.   k++  d( u4 D/ ?; D1 F
    
20. },
    
21. 9 s' h( ]1 v( W* a, O
22. s;
    
23. ) y% N1 k- x( ?  o5 G$ W4 a
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">, P% n* p5 B  P* _/ f$ f* M</font><br /> mvar:<br /> <span style=\\"display:none\\">: z8 [7 t# A7 S7 Z3 Q# M</span> f(x)=x+1;<br /> <span style=\\"display:none\\">0 X5 _  ]1 M1 Y3 n</span> main0(::a)=<font class=\\"jammer\\">* k% G9 C) F- t1 G7 Q</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">4 |5 ^& c' [+ m# e/ l  Q</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">( S. Z! [. j! a\\" q</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">\\" C) e- x' i1 v# Y</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\\">. K6 c! q6 ?$ k: \</font><br />   k++<font class=\\"jammer\\">  d( u4 D/ ?; D1 F</font><br /> },<br /> <span style=\\"display:none\\">9 s' h( ]1 v( W* a, O</span> s;<br /> <span style=\\"display:none\\">) y% N1 k- x( ?  o5 G$ W4 a</span> [clock()-t0]/1000;

结果：
24002271.69189384
1.625

8 D& a5 E* B  `/ W! y* G, b- Jmatlab效率已不及Forcal。
$ t* W$ R' c! e7 G- P
: w  F$ E2 ?1 u提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

forcal

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

& C, x2 m  c1 O! _4 Q6 t: v, _Forcal代码：

2. !using["sys"];
   
3. : l( [4 r, S+ G( j/ G6 E% Z
4. mvar:
   
5. $ {0 \. ~- c/ p0 o8 ~) K% d. \# o
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. / {4 b$ ^1 g5 x4 Z9 H
8. t0=clock(),
   
9. + q& d\\" ]: V0 @/ W4 s
10. k=0, while{k<10000,# R: ]3 p, V: K/ g7 R/ V9 j3 c
    
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. . l( D. r* j  T
14.   k++/ F; y/ W. L! I+ y\\" g
    
16. },* B# C( w  u2 ]& V
    
18. [clock()-t0]/1000;( p, Z6 O/ X7 I0 ^\\" C$ Y0 c2 E
    
20. ' d& N/ D7 J: f3 P9 c: ?; _
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组) _* @+ b0 j; T- T0 \
    
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. ) o! v' j% C6 s5 Q1 U
26. t0=clock(),
    
27. / K1 u, Z. f8 L, g
28. s=0,) M5 R2 j/ X1 }+ \
    
30. k=0, while{k<10000,0 f5 B2 Y( J4 \; x
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值  ~2 [  s9 c' \7 w  _' X9 O- d
    
34.   k++
    
35. $ J: K0 k0 ~2 w5 `
36. },
    
37. 3 T\\" ~' i0 G% [: j* Y8 [5 w
38. s;
    
39.   ~+ F; d( Z0 Z5 J
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">: l( [4 r, S+ G( j/ G6 E% Z</span> mvar:<br /> <span style=\\"display:none\\">$ {0 \. ~- c/ p0 o8 ~) K% d. \# o</span> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">/ {4 b$ ^1 g5 x4 Z9 H</span> t0=clock(),<br /> <span style=\\"display:none\\">+ q& d\\" ]: V0 @/ W4 s</span> k=0, while{k<10000,<font class=\\"jammer\\"># R: ]3 p, V: K/ g7 R/ V9 j3 c</font><br />   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">. l( D. r* j  T</span>   k++<font class=\\"jammer\\">/ F; y/ W. L! I+ y\\" g</font><br /> },<font class=\\"jammer\\">* B# C( w  u2 ]& V</font><br /> [clock()-t0]/1000;<font class=\\"jammer\\">( p, Z6 O/ X7 I0 ^\\" C$ Y0 c2 E</font><br /> <font class=\\"jammer\\">' d& N/ D7 J: f3 P9 c: ?; _</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<font class=\\"jammer\\">) _* @+ b0 j; T- T0 \</font><br /> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">) o! v' j% C6 s5 Q1 U</span> t0=clock(),<br /> <span style=\\"display:none\\">/ K1 u, Z. f8 L, g</span> s=0,<font class=\\"jammer\\">) M5 R2 j/ X1 }+ \</font><br /> k=0, while{k<10000,<font class=\\"jammer\\">0 f5 B2 Y( J4 \; x</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<font class=\\"jammer\\">  ~2 [  s9 c' \7 w  _' X9 O- d</font><br />   k++<br /> <span style=\\"display:none\\">$ J: K0 k0 ~2 w5 `</span> },<br /> <span style=\\"display:none\\">3 T\\" ~' i0 G% [: j* Y8 [5 w</span> s;<br /> <span style=\\"display:none\\">  ~+ F; d( Z0 Z5 J</span> [clock()-t0]/1000;

结果：
, ^' T& O! X5 }+ e! H! A) U2 [0.5620000000000001  //赋值用时间，秒
& H: q5 B/ g. V3 k20000000.
! `) c/ {7 L& G. j$ V/ K- o1 v3 O' z0.579      //取值用时间，秒