极限测试之Matlab与Forcal普通函数调用效率

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。

& B  y  ?& Z' Y% j) y9 W( |5 ?＝＝＝＝＝＝＝＝＝＝＝＝
! O2 k0 s5 D1 v$ i
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;( ^  }& i: G$ G9 r$ e2 h% h* a
   
2. tic;
   ' [8 M8 d8 }: y2 i0 e
3. s=0;0 N2 G: s1 ]( j) U
   
4. for x=0:1000
   : x9 D/ y+ d6 d9 z! d/ E& ?
5.   for y=0:100
   + X8 Q9 Q, y0 j( X
6.     for z=0:100
   6 K7 Y5 j7 |( s& ?
7.       s=s+f(x,y,z);, z; K6 p2 k& h5 Y0 y* p5 C
   
8.     end
   & E3 X& i( L# b% n- f8 D  x
9.   end
   8 X6 m; G  p4 R7 O9 X
10. end% I2 L  y  e\" ^! E
    
11. s6 u! c/ }- |/ W
    
12. toc; E' y  z% ]# L
    
13. 
    9 O# F# {0 |6 q. m6 P! c2 B5 j
14. s =  d; A. q! Q( S0 ?; B
    
15.     6.126720600000000e+009
    ' u; \4 G7 C. L; o; b# b
16. 
    5 g( ]  g9 p\" \; F$ W9 f$ q% R
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
m函数：

1. %file xyz.m
   0 c) d0 a  Q& m$ j4 ]1 s, u
2. function c=xyz(x,y,z)
   , n+ ^, m, @( ?9 q  n9 g+ c
3. c=x+y+z;& N2 H\" t( T! t' ~
   
4. end

复制代码

测试代码：

1. tic;
   8 {; n4 ~  F; q6 O4 V
2. s=0;: l8 |+ Z, H! ?9 _8 f( t! f
   
3. for x=0:1000
   ; T9 k$ h) d' n# c; L! Z& E
4.   for y=0:100- Y! |: ^; Q. O- j% z
   
5.     for z=0:100% r! @  U5 h5 s
   
6.       s=s+xyz(x,y,z);% f0 G! R4 M' }$ F0 m4 z0 d1 P2 ]
   
7.     end
   $ Y+ A1 n\" K: e! L% n
8.   end
   : P0 W7 d, O+ r8 c\" E
9. end  H* x5 H' k* y
   
10. s
    \" i; I# j8 U. _* h% I; |  C
11. toc* ~, x1 W! o4 F$ `0 O
    
12. 
    + c8 y6 T: s1 s& P' m- h\" O
13. s =
      W, Z. n  l+ X* h2 H; g( Y
14.     6.126720600000000e+009
    & ^( ?' l\" }) H\" c
15. % G+ U$ z- d# u$ U
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

6 L5 c  b& p$ G  x  N) Z; }0 M6 \Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;7 b& B7 _+ r1 g( a6 ?& `! t
   
2. mvar:! Y4 L) _6 C$ I8 R
   
3. t=sys::clock(),
   9 I. O/ K) ~/ A6 O
4. s=0,+ X1 v) i0 [  [6 j; j/ J+ y% ]0 O
   
5. x=0, while{x<=1000,8 Z8 c5 Z' Z0 ?+ N! O5 h
   
6.   y=0, while{y<=100,
   2 G7 G/ w: l, J5 S
7.     z=0, while{z<=100,+ M) H. a) Y) f& F; D; _/ _1 {
   
8.       s=s+f(x,y,z),& Y/ _; \' M4 J  `) a: h  Z0 O( L
   
9.       z++
     E3 G' X; ~! p( |1 U\" f
10.     },
    $ ?1 b+ }  U$ k/ F, e* q
11.     y++9 {' [* ~6 c5 S6 ?# V+ b9 |6 ]
    
12.   },
    3 f! L) `& H) N4 l3 G
13.   x++
    $ p6 v0 i\" b( G# P* N1 V9 ?
14. },- \( N6 q' e; |
    
15. s;
    9 u/ B( D) z( t$ Q9 L5 n' S/ t+ |
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

' o0 X7 D. \7 |+ P, W; W二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
9 j4 Q/ W! H: R2 W/ x; ~

hxftntA

, b( c. t2 w! }- _2 y  j4 ]
" o/ S7 y& s* f, e6 S$ g



9 N6 x. ?4 P: i1 y  J( N

/ J- K2 m6 Y; A" s5 @* U
* p% ?7 g) T* C+ D$ {

. k5 T! \- S9 }6 B- A



4 {4 M7 `; K5 w0 d$ Q
7 @/ |- ~+ L: H- {