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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
5 Q* z( z2 ]$ u
＝＝＝＝＝＝＝＝＝＝＝＝
) s' O  E4 Z& Y; H/ r
( k: f# s0 W! |0 W' vMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;+ F2 {* F! J- |- d6 i; K
   
2. tic;
   % ~1 d  C! g6 u8 H9 Y8 T! D
3. s=0;
   / L9 W8 u7 ]3 _/ i
4. for x=0:1000; j( b5 C$ _/ H( k  i9 P0 O
   
5.   for y=0:100
   4 o+ O7 m7 w3 o: q1 k' Y' y
6.     for z=0:100+ f/ I# o& r) a, @
   
7.       s=s+f(x,y,z);
   4 y6 w6 i; |\" k7 |6 T+ a
8.     end6 F8 U: e2 I4 u1 R
   
9.   end+ O9 P/ L2 H/ L9 l! P9 p3 E. v( U
   
10. end  T+ ?) e, \, @; u: v* T
    
11. s
      N4 C1 @% \6 i
12. toc/ B( ~, `( l% R5 ?  X
    
13. 
    5 I; b$ l- F9 H& U$ S* P4 Y
14. s =) p# w9 q/ W- D( o% u
    
15.     6.126720600000000e+009
    ' ]) {& R( k6 _% i
16. 
    9 q* C$ Z/ D* J
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
) k  K# k( o9 M8 b) J) R0 ]# tm函数：

1. %file xyz.m+ j; M$ X\" R; [: Q  ^. U
   
2. function c=xyz(x,y,z)
   - n8 M9 d4 `: j1 U& q
3. c=x+y+z;# T% `! q\" t4 r4 g4 w
   
4. end

复制代码

测试代码：

1. tic;! G) K\" H) C  v  \\" K4 {
   
2. s=0;
   2 k0 P7 x3 e9 y\" V4 n$ n
3. for x=0:10001 W: @5 y6 j2 M
   
4.   for y=0:100
   0 w0 \8 \, C7 H9 \9 @( Z
5.     for z=0:100
   $ q# o! h; z$ M0 t5 \) r3 x
6.       s=s+xyz(x,y,z);: n6 c7 m) b1 ?0 T( I
   
7.     end* u, m\" e! D+ F2 q0 S\" Z. v' ~
   
8.   end
   - b/ @; ~' J- b% R
9. end0 ~' u) @$ ~  C( c. {) _9 e
   
10. s2 S+ |8 f9 z- C& t' ~- ?
    
11. toc
    # M7 D3 \2 O4 H* m9 B; p: z
12. 
    3 _1 U( A! N- _# B$ f
13. s =# e8 T1 P# |$ t% P& X
    
14.     6.126720600000000e+009  ^% |/ [, s. z: d; J/ q9 J$ t5 @
    
15. 
    ! N- _) |# q* \; F- g& X, c
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

0 o, @" f& b. S' e( O2 uForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;' @( A2 n\" Z# f5 v
   
2. mvar:  D\" y2 Z6 V) j# l  \% \
   
3. t=sys::clock(),1 G; o! A' [% ^- u. m
   
4. s=0,
   + j$ j9 x\" A% e7 X# w
5. x=0, while{x<=1000,
   2 b  `# `* _0 v
6.   y=0, while{y<=100,
   8 [' H, ^# p4 }% W8 Z2 O
7.     z=0, while{z<=100,0 F( H0 q* N( s
   
8.       s=s+f(x,y,z),
   4 J) d% z4 X/ S# W' |6 t+ }9 n
9.       z++0 }5 P9 O3 r\" E5 b7 g
   
10.     },! @% a$ B$ [4 @- q& d
    
11.     y++) z: e0 t# H' @* P
    
12.   },0 C0 j! h  Y  X% Y! z/ b2 l5 J6 C
    
13.   x++
    ( D$ n+ B; g& D5 Z* [5 e
14. },! \$ Q( A  a+ n* U* [1 q
    
15. s;$ f! g/ I, m( I( {$ h\" ^  A& {! ~7 h4 Q
    
16. [sys::clock()-t]/1000;

复制代码

结果：
9 G4 `& d. d' S8 r' v  M6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

% s- [7 V: B# r* ?5 G

8 [* d9 L* B! N6 G5 }1 A



% T8 U. Z/ {7 l: M' s4 J$ G8 w
% U9 k2 b4 @. K


3 a5 Q( U/ P4 \) L
4 ~& t( f4 y' h$ y3 T
- g: O; |2 g% L* x- I# X
5 Q; y& g+ ]+ q3 e5 }