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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
" C2 W9 E9 \( b9 n% a
＝＝＝＝＝＝＝＝＝＝＝＝
- ?7 R, m. A7 U$ M
" h# M& p( P5 \' I$ j$ J" ?$ bMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;8 }, U* u2 a- L9 }5 i7 Z
   
2. tic;
   5 W4 H. Q% ]: l9 _7 `
3. s=0;& ~/ K* g) n- U, ~( }3 K& ^
   
4. for x=0:1000; t$ }, W+ J8 r& _
   
5.   for y=0:100
   , A2 @\" a3 ]' m5 O5 t& L- |7 h
6.     for z=0:100, ?$ }, D0 _. ^4 r/ p  ?, l
   
7.       s=s+f(x,y,z);
     n. T' {; Z6 g* R# x  t
8.     end1 O2 r* ^' b4 B' S
   
9.   end
   4 k5 f\" a\" K, T, x5 p3 d1 C1 J
10. end
    , O$ E/ s! l8 R+ F6 N' F$ W' L
11. s) x6 a\" N/ Y! L4 j
    
12. toc
    # v( E8 d  b& S
13. 
    ) ~1 T( l9 h, M- H' K
14. s =5 [\" _* c( w& K$ w' _1 s3 M\" b
    
15.     6.126720600000000e+009) C8 A6 N9 x! @, t+ j& h* _
    
16. / W: S1 X' d3 F: ]
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m$ h  v! R# C9 z7 O3 s
   
2. function c=xyz(x,y,z)
   5 V# u, @, }8 u
3. c=x+y+z;' ?. E3 y* D8 f- r& g( L
   
4. end

复制代码

测试代码：

1. tic;$ l2 o( {- P( j; o( ?+ K1 g
   
2. s=0;, K, o; n\" T! T$ V! C4 v) A
   
3. for x=0:1000
   + f: j- }$ U( @& J; _
4.   for y=0:100' ?* t6 b' ~- p
   
5.     for z=0:100
   / ?, x5 g4 R2 C* }. @
6.       s=s+xyz(x,y,z);
   $ R. q5 N\" G8 C, P6 ~  `7 Y
7.     end
   0 U0 g7 D# a3 i
8.   end& r! X2 m& T# J. ]
   
9. end' |6 T5 J( d7 Z# O* D2 Z9 C
   
10. s6 P/ g! J2 q1 C6 s\" O. P
    
11. toc4 G3 H1 F# T$ |5 G3 E( X% L/ z
    
12. 9 @/ g' Z0 h5 S, g, }
    
13. s =
    ; t# t4 O2 h% c\" R! K
14.     6.126720600000000e+009
    ! d! U+ @' S8 h- \4 ?; |- Y/ g\" l
15. 8 r: |. J0 Y: _8 ?
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;1 I6 {9 ~/ Y3 U/ O8 ~' f* U' i
   
2. mvar:
   : M' G& e' B7 Z' Z0 f
3. t=sys::clock(),; e, {8 w- Z9 g9 p3 N  A\" p
   
4. s=0,
   1 G9 ?. v+ q\" ?& F: |+ d8 _
5. x=0, while{x<=1000,
   : j9 N  F2 S9 Q9 J/ d/ Q/ z% J1 P! t) ^
6.   y=0, while{y<=100,
   9 b) N) D  n. ~2 s  _) [' W
7.     z=0, while{z<=100,
   9 J% v/ O! [# g8 O5 t4 ]
8.       s=s+f(x,y,z),* P- w- n/ \  [% [3 M9 v
   
9.       z++
   ! ~+ b! f  k6 A  ~) x
10.     },
    - S1 U# v+ Q- ~+ L5 i
11.     y++
    6 M% Y  w; I* Q' d
12.   },
    8 L+ Y- P! a  J. J\" _$ V
13.   x++: L' D# F  \1 T) g\" _
    
14. },
    8 t! T& y: D2 M6 F
15. s;
      \. X. {* N& E$ u$ Y, @
16. [sys::clock()-t]/1000;

复制代码

结果：
. D9 w. }: M0 B; a. g# n' L9 t& E6126720600.
1.109  秒

* ]( D/ B1 W; v2 J5 _  x二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
* f1 Y. c: o, M. Z. |7 {7 [8 f( H

hxftntA

/ s+ T: Y0 l! q) H) ^: T

) J! |6 ]$ m( [4 L
8 C- _+ o: U9 W& P' u
  x/ r! n. Z5 v" R4 e5 W! Q3 o  K
" t( p1 i1 z! _5 e+ H( b8 l0 ~' L+ w


* i' W) }+ K% S
+ h2 T# f1 \0 b+ O$ p

, i* E1 i4 J/ C7 X4 N& j


0 a! `- d) E1 n4 `