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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
. Z& J/ m! g: o, _8 Y. m
8 _" h" h% y  g4 \# _. O* e6 N＝＝＝＝＝＝＝＝＝＝＝＝
' x7 F6 F4 V2 }
0 P& l0 ?2 C  o2 H$ iMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   0 U0 V1 k# \; s4 q: D/ C* x! p
2. tic;
   ! X\" z' L$ `! a4 g
3. s=0;
   0 u1 k7 z; Y' h8 ]
4. for x=0:1000
   5 D) O* z  q( Y
5.   for y=0:100
   9 @% k+ o+ i& M1 R' F
6.     for z=0:100
     `- a( K( B  a, `
7.       s=s+f(x,y,z);7 u: v% \; N& o/ l$ J
   
8.     end: T% o: a& r% F
   
9.   end
   1 V\" l( g* i9 {) t/ c3 M
10. end
    ( k4 e6 c' v* g. s
11. s& p\" y6 b8 w# z( p
    
12. toc
    \" ]- U' {: e8 ^5 q5 N7 a9 ]
13. 
    ' b% j5 ~+ k& Y* W3 z% k6 g
14. s =& R- e\" y3 h! r2 P7 j. I
    
15.     6.126720600000000e+009
    5 n% j1 A+ Z* N- @3 E3 H\" ?
16. 
    4 e$ K% [0 L3 N( ^/ d5 }. @% e
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
! Q' u9 f& p& V3 Q$ rm函数：

1. %file xyz.m
   1 M( T( ^3 o) Z9 t' U8 {
2. function c=xyz(x,y,z)6 a- u' H* |5 r
   
3. c=x+y+z;/ M+ L9 G3 j3 K; A& z% s
   
4. end

复制代码

测试代码：

1. tic;  U( [% P2 ]' a& A0 y% a
   
2. s=0;# W+ b$ S; \! T; B
   
3. for x=0:1000
   3 v$ Z, m; x6 _8 m
4.   for y=0:100
   $ o3 g, e& r; @: y9 k
5.     for z=0:100; s3 `4 }+ T0 g# z# V+ E
   
6.       s=s+xyz(x,y,z);
   ( J1 W! K$ h& y5 Y- G# E
7.     end3 A. q: E& g3 v2 j5 H
   
8.   end& R\" `/ x( Q! w\" Z  Z. {
   
9. end
   % V( q5 s$ @6 @, |; j1 l
10. s
    % G* \, a* B% v( u( @- f$ S- o4 p8 K
11. toc
    8 N* X+ J6 i; d- b
12. 
    ) L; O* }9 c0 v  x\" A# v
13. s =6 ^% Q* I9 D* w8 F4 D& ]1 c: b5 l
    
14.     6.126720600000000e+0094 ]1 H) w0 T3 q6 m: {. }
    
15. , z2 f+ G& I1 S: {' i7 ~
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
" }' `/ `2 f! ^+ x" V
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ! S0 @, u' s2 n# e- u
2. mvar:, Y# i) x7 D, \7 G2 T
   
3. t=sys::clock(),
   , f7 B6 h' ^5 Z4 ~
4. s=0,
   ( C\" E% ~& O8 r: J. [8 [- ?
5. x=0, while{x<=1000,, ~7 q. `) C! U7 j6 }
   
6.   y=0, while{y<=100,8 A: f/ o6 {' Q
   
7.     z=0, while{z<=100,: l+ ^& ]- A1 X5 o' n\" h
   
8.       s=s+f(x,y,z),
   \" k$ F1 u, r% r3 X. W
9.       z++% w; q% _( i2 }0 b$ B) {9 e
   
10.     },
    8 q/ Z8 n  t8 P4 X) q. H2 C% `
11.     y++* u  u6 r7 |; X( Z7 Q4 t6 C0 j
    
12.   },
    \" L+ }, l# N( H5 o0 V8 K
13.   x++
    1 u8 e2 ]) l. O3 n; o, n$ U
14. },
    ; A4 {. ?9 Y0 U8 h6 o2 ^9 C$ x
15. s;
    ) ?, n# i4 T; \5 ~) q$ B
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
1 ~8 V2 j4 ?$ {% J4 _/ b& P! S

hxftntA

! S3 w* R2 t5 l


% h1 e, k0 i3 \6 o/ N
4 L) q/ R& L# F



, {9 Z1 I2 K. C: A$ M
( o0 O8 R: y: m5 L7 d9 ~8 I
( G0 ?! _! [  z, F+ a
5 X# Z/ H; c) B$ b0 }
# F5 `# g  H' s- n) J5 Y
; Z2 d4 p! u" |9 j
5 u# @) \7 N5 {& I8 a
0 R# y5 B0 V. M! K