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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
7 Y: B9 Z* g1 ^
+ u$ [: V: W+ ^3 i* J- NMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   * C/ ]0 \, J1 k+ @$ f9 j0 k$ h
2. tic;
   $ ?* w2 C2 e( _+ b3 l) e+ D+ K4 a
3. s=0;' L% a\" h+ D' f+ E+ n
   
4. for x=0:10009 }9 @4 d, B; Y& A& a' Z6 i2 Q, F
   
5.   for y=0:100' `) V# R. B0 X, ]6 A
   
6.     for z=0:100. ^8 v\" i3 Q\" L+ ^: c, v5 O
   
7.       s=s+f(x,y,z);$ w$ Q1 T2 e& H7 E% C# y4 F8 O- P
   
8.     end( A! D0 f. e+ a3 Y( b
   
9.   end6 y: _, U0 i  f' n
   
10. end
    * g8 c4 B6 T! H: p4 O: |/ F
11. s\" O\" [: K, v: `# s6 i$ |
    
12. toc+ o: s- v\" Y1 D\" e* s1 z
    
13. \" Q4 u: a% a3 m+ O# _% B2 F
    
14. s =
    3 X% o9 I; z! Z
15.     6.126720600000000e+009- j& h% X: J+ e& C+ P
    
16. 
    ; \/ b# ?2 K0 e2 g
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
  o( M* n6 y0 i6 ?+ gm函数：

1. %file xyz.m' w# N, B- ]0 s
   
2. function c=xyz(x,y,z)% I- `( w7 H5 @7 _2 ?
   
3. c=x+y+z;5 {$ r; b5 w, o9 S: @. l& b/ r
   
4. end

复制代码

测试代码：

1. tic;. R; k/ x. ]) \9 P/ h6 G8 z5 |+ w
   
2. s=0;  ^2 m% w% T' ^9 {6 U; K% n! }
   
3. for x=0:1000
     p$ c  L' u+ H) ^- y# u
4.   for y=0:100! W. b# A  ?. h, Y0 s9 n* r) v9 o9 K
   
5.     for z=0:100
   6 D9 Y\" |, |2 m, l
6.       s=s+xyz(x,y,z);
     z& B& h- d) a& m5 @
7.     end
   - }! q5 j0 j# `7 P! R, R# J0 V* g$ ^
8.   end
   ! H, c# B& d6 `; B# `9 e\" e
9. end6 {, |! T. m  Y, h( m2 Y  L) \' a
   
10. s7 A& }1 \5 w4 u+ f! ^* x\" E1 }+ D
    
11. toc  V$ n  o, ]\" U; G3 i
    
12. 
    + g' v; L: y. u* V6 c! ]; t( W, K
13. s =3 U7 _. K  \) P0 ~& L& u
    
14.     6.126720600000000e+009
    4 P' {; L( u8 T
15. 
    1 R* `) Q\" W, V* m0 `8 j: }& k
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

: E6 U! Z. x7 p" J, R. TForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;! M2 j: I, I$ e( w) O
   
2. mvar:3 q4 @. S! d9 B
   
3. t=sys::clock(),+ u2 G3 M* E* t4 m
   
4. s=0,
   ) \( w; H) U- u: T/ p
5. x=0, while{x<=1000,
   ' ]\" ]* E! }; a+ a$ B
6.   y=0, while{y<=100,
   * [) u$ V+ c' ?. T0 T, w* x
7.     z=0, while{z<=100,4 P! x/ L2 t# Q' y9 U% ?  q8 v! C5 }
   
8.       s=s+f(x,y,z),% h  n6 x! @% ]# C( I- e0 O' _
   
9.       z++
   * O\" |* U* W: r& {: f; k
10.     },
    - B  u\" L1 m! f# \\" q3 i
11.     y++6 `. J# @8 Z6 p3 n# ?9 l
    
12.   },$ Z& G4 a0 F  _5 c. Z\" S  T  }
    
13.   x++
      Q  p- L9 a9 C- |8 _1 P
14. },( L# v0 \) _3 a/ l. J/ f' t8 D
    
15. s;
    : H% p& g5 t$ G# v- |* K
16. [sys::clock()-t]/1000;

复制代码

结果：
; T" V0 g) x$ q/ N6126720600.
# H  ~' z1 a" F! C1.109  秒
) W- b7 o; f6 r3 ?- t$ S
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

" I' N# t6 ]+ e( u2 t5 p1 ]' W
8 l9 ^. d+ K$ d) e9 V
  P$ P" s5 m6 H( J3 S
6 p& ^3 ]7 f! G4 Y
8 h" m9 s4 ?  k- @& h


+ t8 x- c! j8 d# d; f
7 P( G, A1 c9 A! y) B$ P( B

7 \2 F" w- Q% {% ^1 ?
& U3 a; p# `5 X, r7 l3 w) x5 r7 a
% e9 Z% }( ?* n3 d* ?( _



" Z1 E' n7 _% A5 m" t" d6 t0 y