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

forcal

1、简单循环
' C! D0 p9 d* c8 s. x
  F4 b' e& _3 S1 c7 c: Y5 KMatlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

1. clear all
   : M0 R( g7 o3 q9 ]% I% w
2. tic  p3 |$ P1 h5 q- H$ t. A; {# A
   
3. s=0;
   \" A' O! H2 y: W
4. for k=1:10000# H4 u1 l/ m+ D2 |4 z( u
   
5.   for i=1:1000
   : o+ w* z0 |\" M
6.     s=s+1;
   4 I4 k5 ~; i; _# H
7.   end. }1 d& o+ e2 O) b* m\" }3 f
   
8. end
   0 H: P; Z5 t. |2 d; u2 T
9. s
   % }  X3 x$ s, T# M  k
10. toc
    ! l3 {6 a4 u+ k9 z0 U; o1 n
11. 
    * O/ B4 n$ Y1 Q2 ?7 s$ _
12. s =: z; o. P4 [8 I% ]6 s1 m
    
13. + H, ?1 S% I2 L/ C  }# A* h* W
    
14.     10000000
    & `+ r7 y9 t/ @  b9 l
15. ' r3 v9 i1 ]: m
    
16. Elapsed time is 0.092808 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
- [4 W: D$ S4 u+ X  I4 K
Forcal没有JIT加速器，代码：

2. !using["sys"];
   
3. 8 ^' H  L, q. ~
4. mvar:
   
5. 2 {+ g+ n' \! V! t7 J2 V
6. t0=clock(),# C, v: w: `- A  y4 u* w0 B/ e
   
8. s=0,- t, s9 o1 g( M
   
10. k=0, while{k<10000,
    
11. 6 J, q9 H$ L+ r$ y' ~# q1 q, H! R
12.   i=0, while{i<1000,s=s+1, i++},
    
13. ! C; w# S, p# ]4 U4 Z! ?& L' N
14.   k++) M2 W- o) m% H  L# u
    
16. },# n+ M- b0 ?$ S1 {) c
    
18. s;4 h\\" l6 [7 V/ @, c+ E
    
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">8 ^' H  L, q. ~</span> mvar:<br /> <span style=\\"display:none\\">2 {+ g+ n' \! V! t7 J2 V</span> t0=clock(),<font class=\\"jammer\\"># C, v: w: `- A  y4 u* w0 B/ e</font><br /> s=0,<font class=\\"jammer\\">- t, s9 o1 g( M</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">6 J, q9 H$ L+ r$ y' ~# q1 q, H! R</span>   i=0, while{i<1000,s=s+1, i++},<br /> <span style=\\"display:none\\">! C; w# S, p# ]4 U4 Z! ?& L' N</span>   k++<font class=\\"jammer\\">) M2 W- o) m% H  L# u</font><br /> },<font class=\\"jammer\\"># n+ M- b0 ?$ S1 {) c</font><br /> s;<font class=\\"jammer\\">4 h\\" l6 [7 V/ @, c+ E</font><br /> [clock()-t0]/1000;

结果：
10000000.
8 b# N& f, i" \; a: l0.704  秒
0 f# D4 L" V( t7 u$ a$ j
" m- n% Z* d& {二者运行效率相差7倍多。

2、循环中稍微加点东西
, B8 V! z1 S5 x. B$ I
Matlab代码：

1. clear all
   9 `  G! i+ I' Z9 s\" o: K
2. tic
   \" n4 x2 r8 Q& r) ]9 M  O' z, N! I& C3 d
3. s=0;
   , d3 b- M. m% U6 c6 q, ]& Q! N
4. for k=1:1000
   ( D$ v! c8 G% p  }& r1 A
5.   for i=1:1000! _3 O9 }) D- k
   
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);
   $ |* A7 ?& R/ n4 C+ i. q8 b) I
7.   end
   $ s/ p/ L; @8 D4 u/ F$ r
8. end
   ) }3 N% T# U& |7 y/ T- l0 ]8 x
9. s
   1 Y7 c6 e7 t  M  p  y
10. toc
    ( d5 Q4 t! F; @; I& v1 x
11. . }( i, I7 s: w9 U\" K5 N6 A4 U
    
12. s =$ d! d# n2 L% w/ c3 h
    
13. 6 [# T8 j+ I- y5 r& ?
    
14.   1.1009e+007
    & y, ~\" S4 @( Y6 [- ^
15. * Y2 z3 J- S2 _* u. c: l1 w2 c
    
16. Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. - F7 Z. u  z* s. i* E
4. mvar:
   
5. . m4 Y9 C9 k: |6 `4 y, N\\" p2 l) q
6. t0=clock(),
   
7. - K  Q1 [6 O* F
8. s=0,
   
9. * v' V# R5 }/ }
10. k=1, while{k<=1000,# P8 Z, N# w; G. N3 d# B
    
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. ' A% O\\" E0 b: |: ~- n
14.   k++
    
15. % i2 j5 y* Q# l; l# \. r2 S
16. },
    
17.   N7 C1 p! T6 ?* ]% p9 m
18. s;
    
19. 8 z( V) a6 X+ `$ l4 v
20. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">- F7 Z. u  z* s. i* E</span> mvar:<br /> <span style=\\"display:none\\">. m4 Y9 C9 k: |6 `4 y, N\\" p2 l) q</span> t0=clock(),<br /> <span style=\\"display:none\\">- K  Q1 [6 O* F</span> s=0,<br /> <span style=\\"display:none\\">* v' V# R5 }/ }</span> k=1, while{k<=1000,<font class=\\"jammer\\"># P8 Z, N# w; G. N3 d# B</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\\">' A% O\\" E0 b: |: ~- n</span>   k++<br /> <span style=\\"display:none\\">% i2 j5 y* Q# l; l# \. r2 S</span> },<br /> <span style=\\"display:none\\">  N7 C1 p! T6 ?* ]% p9 m</span> s;<br /> <span style=\\"display:none\\">8 z( V) a6 X+ `$ l4 v</span> [clock()-t0]/1000;

结果：
11008953.66602346
1.031
+ `; l% i+ q$ E  P7 e* U
7 S5 ^0 N4 T% E3 {* U二者效率相差已比较小。
) @1 y! _$ ?3 l+ U" [* e
3、再增加一个函数调用

Matlab代码：

1. clear all
   / Q, F; `( C7 K3 T6 p
2. f=@(x)x+1;
   , v! L$ V+ j1 B% K; \) a' g
3. tic
   $ e/ L\" L1 @3 r* c
4. s=0;
   0 p- y6 b+ l! O  w2 ^$ w: F
5. for k=1:1000
   0 W* c\" i) C/ b# H3 `
6.   for i=1:1000
   ; r$ W: ^5 \) |6 i) ^- x  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);
   . z3 s( |6 |% t4 M( W7 O: G
8.   end
   ! C7 U4 q$ R% _( U, `: O$ w6 g6 L/ Y
9. end
   $ Y9 G9 H9 m5 P! T9 Q
10. s/ {- H4 o# x; d) Q7 w
    
11. toc$ V- D\" h% }5 w: v4 g$ w
    
12. 
    8 H$ f9 t/ I- d3 V
13. s =
    5 L7 t' I$ E6 E% u
14. 
    7 V7 ]- K9 s+ n/ W# P
15.   5.1251e+008
      F7 R% |9 h6 q+ f5 l1 ?, @
16. 
    - X: ?0 f3 D9 j4 e% U
17. Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

2. !using["sys"];
   
3. ' `( i' M$ e0 p5 y4 I$ |
4. f(x)=x+1;
   
5. 9 y( J3 D$ A* i1 P6 K  F+ K* |4 m$ R
6. mvar:
   
7. ' V( T6 v% w\\" r$ O$ g5 v8 O
8. t0=clock(),
   
9. 7 A. J7 I! |% G2 x1 A
10. s=0,! n% l7 r1 t  J  o2 D
    
12. k=1, while{k<=1000,  O( M* L' c# x  }6 W
    
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. , |3 k* A8 L$ i4 M: Z
16.   k++
    
17. 8 `6 w# k! {' y- `* |* M
18. },
    
19. ' U4 b' H: o: Y+ O0 T
20. s;% e: T\\" l  [\\" E( r6 \0 I7 N
    
22. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">' `( i' M$ e0 p5 y4 I$ |</span> f(x)=x+1;<br /> <span style=\\"display:none\\">9 y( J3 D$ A* i1 P6 K  F+ K* |4 m$ R</span> mvar:<br /> <span style=\\"display:none\\">' V( T6 v% w\\" r$ O$ g5 v8 O</span> t0=clock(),<br /> <span style=\\"display:none\\">7 A. J7 I! |% G2 x1 A</span> s=0,<font class=\\"jammer\\">! n% l7 r1 t  J  o2 D</font><br /> k=1, while{k<=1000,<font class=\\"jammer\\">  O( M* L' c# x  }6 W</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\\">, |3 k* A8 L$ i4 M: Z</span>   k++<br /> <span style=\\"display:none\\">8 `6 w# k! {' y- `* |* M</span> },<br /> <span style=\\"display:none\\">' U4 b' H: o: Y+ O0 T</span> s;<font class=\\"jammer\\">% e: T\\" l  [\\" E( r6 \0 I7 N</font><br /> [clock()-t0]/1000;

结果：
- b7 w5 K& W, l0 b3 r512508953.6658993
1.078
+ }' O# b4 _3 T* p  j& }
可以看出，Matlab效率已不及Forcal。
4 \" Z8 @' ~' k% u: \故JIT加速器虽好，但效果似乎并不明显。

海水

forcal

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

9 {0 P7 u0 @# I5 J0 aMatlab代码：

1. >> clear all$ {* X: I  v* M; C8 `\" }! c
   
2. a = linspace(1,1,1000);
   , q8 I1 a) Q7 q0 o9 O' k
3. tic
   + |+ `. l$ W6 o  A9 E% H) [# d+ a9 \
4. for k=1:10000
     y6 U1 o+ v7 O' L' d
5.   for i=1:1000, t+ B. H6 u1 q4 K) v, ?. w( A
   
6.     a(i)=2;
   # Y4 d+ d$ E& K5 |\" ]$ e8 D\" Z
7.   end
   . F' ]0 B3 S0 a) t; H
8. end
   / K, w; Z; x  a$ G* [  }' n% f: Z
9. toc
   9 O\" G5 \\" A8 C  O$ {: D4 t& p
10. 
    0 c' X1 u7 B9 J( @0 c3 i
11. tic+ ?+ V) q0 f! P2 t1 a/ l2 x7 Z: s  J( Y
    
12. s=0;
    3 M# \, W: h; j* e4 N- H
13. for k=1:100001 X1 V( e' {* L\" m$ Q, Z7 v0 v5 o/ X
    
14.   for i=1:1000
    $ K  W1 e# C+ B& |. @1 f% o8 b
15.     s=s+a(i);' P. `1 F0 s9 N4 ~
    
16.   end
    7 C* b* ~6 E' w0 V  ~
17. end
    8 x5 @1 e/ ]1 ]8 V\" F
18. s2 E/ T, B, l7 n, y9 e' P0 A* W
    
19. toc6 K* K8 s& B. h- P6 x% d\" Q: m
    
20. Elapsed time is 0.140658 seconds.: G6 E. Y2 `, Q5 T% \  d5 N& D
    
21. 
      P- ^1 A# G7 S% s9 R) H4 j
22. s =
    8 U( {9 ?; u9 D' K8 o
23.   N3 i; {' U( V7 O  k8 j5 s6 g( h
    
24.     20000000, o0 w; x3 e- [
    
25. 3 D$ C  k/ D/ B( P& T# t, D
    
26. Elapsed time is 0.091794 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝
  z! b$ u5 ?1 Y
Forcal代码：

2. !using["math","sys"];# G% _$ E) A7 {- v9 P3 P
   
4. mvar:9 n. g4 \* K$ u% [  k
   
6. main0(::a)=\\" S% U+ w1 Q% A) w  P0 ?) L$ i% a
   
8. a = linspace(1,1,1000),
   
9. 2 e' g' l  S3 r
10. t0=clock(),: d  z. g+ @; u. p
    
12. k=0, while{k<10000,
    
13. 9 g$ N9 d8 ^5 r9 o
14.   i=0, while{i<1000, a(i)=2, i++},
    
15. ) W' k/ Q! ~; X2 v. I
16.   k++
    
17. ( c0 k4 H/ P\\" x0 |: Q9 j; [( e
18. },' C' q. h, T& T/ U1 T3 o
    
20. [clock()-t0]/1000;
    
21. 9 X* v7 n* T4 T2 k  x6 `8 b
22. 1 q) E; q0 B! D0 T+ B$ b
    
24. main1(::a)=
    
25. : `$ v+ v4 M- R\\" ^
26. t0=clock(), s=0,
    
27. , O! l$ r0 J1 u1 U! z( Q0 A5 o0 P
28. k=0, while{k<10000,; H& }4 r' E# c& i
    
30.   i=0, while{i<1000, s=s+a(i), i++},. L8 K9 Q3 M$ m8 }\\" T5 d3 p5 |9 ~1 A
    
32.   k++
    
33. ; H) G& ?\\" t# {. q9 f
34. },' B5 M9 z$ c1 C- L) p
    
36. s;, W\\" |2 t+ `8 l& G
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\"># G% _$ E) A7 {- v9 P3 P</font><br /> mvar:<font class=\\"jammer\\">9 n. g4 \* K$ u% [  k</font><br /> main0(::a)=<font class=\\"jammer\\">\\" S% U+ w1 Q% A) w  P0 ?) L$ i% a</font><br /> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">2 e' g' l  S3 r</span> t0=clock(),<font class=\\"jammer\\">: d  z. g+ @; u. p</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">9 g$ N9 d8 ^5 r9 o</span>   i=0, while{i<1000, a(i)=2, i++},<br /> <span style=\\"display:none\\">) W' k/ Q! ~; X2 v. I</span>   k++<br /> <span style=\\"display:none\\">( c0 k4 H/ P\\" x0 |: Q9 j; [( e</span> },<font class=\\"jammer\\">' C' q. h, T& T/ U1 T3 o</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">9 X* v7 n* T4 T2 k  x6 `8 b</span> <font class=\\"jammer\\">1 q) E; q0 B! D0 T+ B$ b</font><br /> main1(::a)=<br /> <span style=\\"display:none\\">: `$ v+ v4 M- R\\" ^</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">, O! l$ r0 J1 u1 U! z( Q0 A5 o0 P</span> k=0, while{k<10000,<font class=\\"jammer\\">; H& }4 r' E# c& i</font><br />   i=0, while{i<1000, s=s+a(i), i++},<font class=\\"jammer\\">. L8 K9 Q3 M$ m8 }\\" T5 d3 p5 |9 ~1 A</font><br />   k++<br /> <span style=\\"display:none\\">; H) G& ?\\" t# {. q9 f</span> },<font class=\\"jammer\\">' B5 M9 z$ c1 C- L) p</font><br /> s;<font class=\\"jammer\\">, W\\" |2 t+ `8 l& G</font><br /> [clock()-t0]/1000;

结果：
15.484  //存数据所用时间，秒
20000000.
5.531   //取数据所用时间，秒
; o# s; R; N. `! L+ @
; k: B  p$ W+ k( ~7 Q－－－－－－－－－－

" \. j9 \( t! a用函数sys::A()存取数组数据效率稍高些：

2. !using["math","sys"];; H. z0 Z! F1 I3 s1 N4 g
   
4. mvar:9 D7 p. H- [/ S. \4 u' \
   
6. main0(::a)=
   
7. / M6 Q+ Y. c% E( z7 P
8. a = linspace(1,1,1000),
   
9. . t/ Q( ~- m0 N3 H\\" o9 R& D0 S
10. t0=clock(),\\" v2 ?! v1 h1 d; I& L$ J- _
    
12. k=0, while{k<10000,
    
13. / X( ^8 N5 }7 t
14.   i=0, while{i<1000, A(a,i)=2, i++},! n% k: F$ c# ~8 M3 x$ z/ i  {' g# {
    
16.   k++0 h' i6 S' ~\\" e* A8 L' ?
    
18. },( P! G7 I) E: B# D. C  ]
    
20. [clock()-t0]/1000;
    
21. - c6 y$ E6 [5 y\\" q
22. / u' T$ z\\" b/ Q  E
    
24. main1(::a)=2 B; B, }# `( Q, A; C
    
26. t0=clock(), s=0,
    
27. 8 t/ [, U& C* F+ ?8 @' G
28. k=0, while{k<10000,
    
29. ! @3 t\\" C$ V* F3 n  |' n3 J
30.   i=0, while{i<1000, s=s+A(a,i), i++},\\" V3 W' N+ l/ B/ t. c
    
32.   k++2 k- E2 X. G7 i1 P1 h
    
34. },
    
35. + u. D. M5 l# U( B9 ~
36. s;0 E# T+ M8 e+ o; J3 V0 P
    
38. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">; H. z0 Z! F1 I3 s1 N4 g</font><br /> mvar:<font class=\\"jammer\\">9 D7 p. H- [/ S. \4 u' \</font><br /> main0(::a)=<br /> <span style=\\"display:none\\">/ M6 Q+ Y. c% E( z7 P</span> a = linspace(1,1,1000),<br /> <span style=\\"display:none\\">. t/ Q( ~- m0 N3 H\\" o9 R& D0 S</span> t0=clock(),<font class=\\"jammer\\">\\" v2 ?! v1 h1 d; I& L$ J- _</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">/ X( ^8 N5 }7 t</span>   i=0, while{i<1000, A(a,i)=2, i++},<font class=\\"jammer\\">! n% k: F$ c# ~8 M3 x$ z/ i  {' g# {</font><br />   k++<font class=\\"jammer\\">0 h' i6 S' ~\\" e* A8 L' ?</font><br /> },<font class=\\"jammer\\">( P! G7 I) E: B# D. C  ]</font><br /> [clock()-t0]/1000;<br /> <span style=\\"display:none\\">- c6 y$ E6 [5 y\\" q</span> <font class=\\"jammer\\">/ u' T$ z\\" b/ Q  E</font><br /> main1(::a)=<font class=\\"jammer\\">2 B; B, }# `( Q, A; C</font><br /> t0=clock(), s=0,<br /> <span style=\\"display:none\\">8 t/ [, U& C* F+ ?8 @' G</span> k=0, while{k<10000,<br /> <span style=\\"display:none\\">! @3 t\\" C$ V* F3 n  |' n3 J</span>   i=0, while{i<1000, s=s+A(a,i), i++},<font class=\\"jammer\\">\\" V3 W' N+ l/ B/ t. c</font><br />   k++<font class=\\"jammer\\">2 k- E2 X. G7 i1 P1 h</font><br /> },<br /> <span style=\\"display:none\\">+ u. D. M5 l# U( B9 ~</span> s;<font class=\\"jammer\\">0 E# T+ M8 e+ o; J3 V0 P</font><br /> [clock()-t0]/1000;

结果：
: `6 H0 U8 z* }  }) ~, ?4 A  A10.469
+ L! H: `: d# Q& a' n# z20000000.
3.781
8 w6 H7 S- `; n1 q5 U
7 Y% C" b* E* |! {2 B. X－－－－－－－
9 ^, g4 z0 t" @
9 E0 @) I8 V, X# z; i0 r在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
6 m0 q8 C, ^4 D/ N3 l0 J  _1 z在使用函数sys::A(  )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。
" c2 K+ ~# e. l& n$ C
+ w) X# q2 [! I5 v& I! _matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。
% a  u( a$ p) R) Z" g  t
Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。
/ m8 d0 W( b" [, L, g* e
5、在循环中加点东西

matlab代码：

1. >> clear all
   1 F9 r\" x. T. e0 W8 T) S  B7 `/ p! I
2. f=@(x)x+1;  {8 o' x5 [\" G  [
   
3. a = linspace(2,2,1000);
   7 P$ \4 F( M2 N2 H0 f# F8 G
4. tic
   - ?8 F+ E4 ]* p; u# P+ o* G. w
5. s=0;
   / a3 ~) g' B1 A3 G
6. for k=1:1000% {6 t+ H; W* {- @5 ?: \& V) Z1 `7 ^& Q
   
7.   for i=1:10005 ]7 O/ v, J- k8 ]# m* Q! ^1 o\" Y
   
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));4 n2 n\" C! a1 w5 X7 [0 q* ~
   
9.   end
   * v$ ^4 p2 D9 }5 F! J
10. end
    # \$ @* w/ r/ R8 N
11. s
    3 d9 F+ S9 y4 V9 g
12. toc
    & F2 n2 J+ e; c
13. 
    ; m\" Z& `0 w% `4 \
14. s =
    # J\" d# p! X2 a6 Q' T% k
15. , R7 e\" l* G# y  t4 d
    
16.   2.4002e+007- g+ M8 B# J! g: \
    
17. 
    # ]( _6 [& z8 ^) t- Q8 S1 d& t  ?
18. Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

2. !using["math","sys"];! ~4 h' _. U6 V5 }- T5 `+ |
   
4. mvar:/ c# S' w; m1 c5 n% z
   
6. f(x)=x+1;
   
7. 1 V* V! `0 q8 e4 ^/ H' {
8. main0(::a)=% j( i6 n7 j# N% h$ b: {7 H. t
   
10. a = linspace(2,2,1000),
    
11. + G$ n2 _8 S6 j/ V% |: g
12. t0=clock(), s=0,
    
13. / y7 Y. n1 U2 s4 p3 ^* p
14. k=0, while{k<1000,
    
15. % W9 g' o# I5 d% ]\\" d- 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++},$ f! b  ~, u! M& {\\" x( p3 A3 W8 E
    
18.   k++
    
19. 3 z% o% S0 q: K5 N) D- A2 A0 X
20. },
    
21. $ M3 T8 s8 o9 [
22. s;' o5 D0 i5 l/ g' L+ Q( T
    
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">! ~4 h' _. U6 V5 }- T5 `+ |</font><br /> mvar:<font class=\\"jammer\\">/ c# S' w; m1 c5 n% z</font><br /> f(x)=x+1;<br /> <span style=\\"display:none\\">1 V* V! `0 q8 e4 ^/ H' {</span> main0(::a)=<font class=\\"jammer\\">% j( i6 n7 j# N% h$ b: {7 H. t</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">+ G$ n2 _8 S6 j/ V% |: g</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">/ y7 Y. n1 U2 s4 p3 ^* p</span> k=0, while{k<1000,<br /> <span style=\\"display:none\\">% W9 g' o# I5 d% ]\\" d- n</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\\">$ f! b  ~, u! M& {\\" x( p3 A3 W8 E</font><br />   k++<br /> <span style=\\"display:none\\">3 z% o% S0 q: K5 N) D- A2 A0 X</span> },<br /> <span style=\\"display:none\\">$ M3 T8 s8 o9 [</span> s;<font class=\\"jammer\\">' o5 D0 i5 l/ g' L+ Q( T</font><br /> [clock()-t0]/1000;

结果：
5 ]  z6 a: j' l- `. B$ S24002271.69189384
  U; V7 u% n- U' }! B7 I9 b- V5.172

二者效率现在只有一倍差距。

－－－－－－－－－－

不过，显然matlab和Forcal的程序应该这样写：
% v2 d2 d4 U: t. c" C
# R8 d5 w( }4 \$ Z6 [& @; @matlab程序：

1. clear all\" c/ K) ?# u7 n# Z% Q- {) _
   
2. f=@(x)x+1;
   * Q& o- Y1 Q: d* v( P$ Z
3. a = linspace(2,2,1000);: d5 N2 M7 R3 p4 j# @0 c$ N
   
4. tic
   / m+ V( h$ c5 P1 y
5. s=0;
   + L1 K! D  I) p% m! S  h
6. for k=1:10007 T( R: N' ^' G4 [' ~* J$ v: o: _
   
7.   for i=1:1000
   & J\" U. [( A0 Y5 x+ ~, d5 C7 C
8.     t=a(i);% D1 `5 _3 p, g: ~0 t: v7 r
   
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);
   9 T: N8 o. y' k
10.   end) q' F$ ~+ }. l# n, d9 x& v
    
11. end' F$ c4 M( K$ \' C8 a- r  ?5 W
    
12. s
    + @5 i' Y& y3 R/ W) c, E* M# k
13. toc
    2 I! X% y  U; _3 s# k
14. , _% i  e9 Z0 `
    
15. s =
    - ^% A- N) x- [. j& M- b3 e
16. 
    ! s9 J1 q. Y5 a+ S  R7 J
17.   2.4002e+007% S) G7 r4 w0 B0 u8 s
    
18. 0 G. L$ i6 E/ \( W
    
19. Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

2. !using["math","sys"];1 m0 j- e) @* L0 r  _' t( }1 v6 P
   
4. mvar:
   
5. 0 q. Z- y4 V2 f+ q3 ]
6. f(x)=x+1;, G- j9 a7 b0 @
   
8. main0(::a)=2 A+ f/ J\\" z6 G$ X$ u
   
10. a = linspace(2,2,1000),
    
11. 4 B4 \  g) B' Q% X; ?
12. t0=clock(), s=0,
    
13. & @8 K2 D% u5 _# z2 D* r
14. k=0, while{k<1000,: D. ?* n: D* D8 O2 l. X0 S9 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. & p2 Y8 Z/ a$ I3 S
18.   k++
    
19. # Q; P, I  I, ?7 H4 f5 `
20. },; o; I+ v6 v8 T1 ^( _- ^' d+ ^
    
22. s;
    
23. , n& a  @' P. z4 ?9 y$ N
24. [clock()-t0]/1000;

!using["math","sys"];<font class=\\"jammer\\">1 m0 j- e) @* L0 r  _' t( }1 v6 P</font><br /> mvar:<br /> <span style=\\"display:none\\">0 q. Z- y4 V2 f+ q3 ]</span> f(x)=x+1;<font class=\\"jammer\\">, G- j9 a7 b0 @</font><br /> main0(::a)=<font class=\\"jammer\\">2 A+ f/ J\\" z6 G$ X$ u</font><br /> a = linspace(2,2,1000),<br /> <span style=\\"display:none\\">4 B4 \  g) B' Q% X; ?</span> t0=clock(), s=0,<br /> <span style=\\"display:none\\">& @8 K2 D% u5 _# z2 D* r</span> k=0, while{k<1000,<font class=\\"jammer\\">: D. ?* n: D* D8 O2 l. X0 S9 x</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++},<br /> <span style=\\"display:none\\">& p2 Y8 Z/ a$ I3 S</span>   k++<br /> <span style=\\"display:none\\"># Q; P, I  I, ?7 H4 f5 `</span> },<font class=\\"jammer\\">; o; I+ v6 v8 T1 ^( _- ^' d+ ^</font><br /> s;<br /> <span style=\\"display:none\\">, n& a  @' P. z4 ?9 y$ N</span> [clock()-t0]/1000;

结果：
5 Q4 r# Z3 K% j+ z* Y- P5 a24002271.69189384
1.625
8 f; T; _0 b$ x% z* a) X6 C
8 x. j3 ?& v7 x' Imatlab效率已不及Forcal。

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

forcal

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。
8 G1 J! i% R7 L# O  |% U1 W
$ L8 g+ `4 V' O* o3 A. H- P/ v2 YForcal代码：

2. !using["sys"];
   
3. 0 z. E7 t3 @6 S5 t5 D. [/ y9 @/ d, b
4. mvar:- ^. z9 h* S8 h( n
   
6. a="\&\[1000]"/4,    //用字符串模拟静态数组
   
7. $ a. G& o# L4 m0 R8 {
8. t0=clock(),7 N) `( ~( M- e9 S
   
10. k=0, while{k<10000,
    
11. 9 k! ~, b7 c! H9 a6 m$ `
12.   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值
    
13. 0 r4 k8 I* K# V4 z5 Y
14.   k++; k* Q  n! u, h2 W
    
16. },
    
17. , v, n+ T0 P\\" |. r, T7 a! P
18. [clock()-t0]/1000;+ s+ Z: r. W- n: S' Q) `
    
20. # R/ @% J8 w% _- Q7 D
    
22. a="\&\[1000]"/4,    //用字符串模拟静态数组
    
23. / {0 l. h& m( R6 ^! {
24. i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值
    
25. 0 W3 z( Q5 j3 r' Y
26. t0=clock(),
    
27. 0 Z, c  N1 g8 D
28. s=0,
    
29. ) n8 C6 I\\" X0 ^
30. k=0, while{k<10000,* Z5 H- t/ n, \, s/ f
    
32.   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值
    
33. 7 p2 V) O! {\\" Z  ~8 w; M' H7 ]2 V
34.   k++
    
35. ! e$ h6 v9 N3 {
36. },
    
37. ' L  V+ i! I: {  c/ R; K1 H
38. s;
    
39. 9 ?/ Y$ J* A\\" C4 \9 Q. w5 U3 l. M
40. [clock()-t0]/1000;

!using["sys"];<br /> <span style=\\"display:none\\">0 z. E7 t3 @6 S5 t5 D. [/ y9 @/ d, b</span> mvar:<font class=\\"jammer\\">- ^. z9 h* S8 h( n</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">$ a. G& o# L4 m0 R8 {</span> t0=clock(),<font class=\\"jammer\\">7 N) `( ~( M- e9 S</font><br /> k=0, while{k<10000,<br /> <span style=\\"display:none\\">9 k! ~, b7 c! H9 a6 m$ `</span>   i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">0 r4 k8 I* K# V4 z5 Y</span>   k++<font class=\\"jammer\\">; k* Q  n! u, h2 W</font><br /> },<br /> <span style=\\"display:none\\">, v, n+ T0 P\\" |. r, T7 a! P</span> [clock()-t0]/1000;<font class=\\"jammer\\">+ s+ Z: r. W- n: S' Q) `</font><br /> <font class=\\"jammer\\"># R/ @% J8 w% _- Q7 D</font><br /> a="\&\[1000]"/4,    //用字符串模拟静态数组<br /> <span style=\\"display:none\\">/ {0 l. h& m( R6 ^! {</span> i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值<br /> <span style=\\"display:none\\">0 W3 z( Q5 j3 r' Y</span> t0=clock(),<br /> <span style=\\"display:none\\">0 Z, c  N1 g8 D</span> s=0,<br /> <span style=\\"display:none\\">) n8 C6 I\\" X0 ^</span> k=0, while{k<10000,<font class=\\"jammer\\">* Z5 H- t/ n, \, s/ f</font><br />   i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值<br /> <span style=\\"display:none\\">7 p2 V) O! {\\" Z  ~8 w; M' H7 ]2 V</span>   k++<br /> <span style=\\"display:none\\">! e$ h6 v9 N3 {</span> },<br /> <span style=\\"display:none\\">' L  V+ i! I: {  c/ R; K1 H</span> s;<br /> <span style=\\"display:none\\">9 ?/ Y$ J* A\\" C4 \9 Q. w5 U3 l. M</span> [clock()-t0]/1000;

结果：
4 D; a" F% H2 S0 c0.5620000000000001  //赋值用时间，秒
# c$ n$ P. y& g0 Z# R20000000.
# T  i- R% p& f5 ^9 u2 L0.579      //取值用时间，秒