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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
0 x' B, M* q. Y  A$ S0 I
＝＝＝＝＝＝＝＝＝＝＝＝
6 X+ \$ n: X; V0 E0 K9 N6 ]1 B, I
1 y# [3 w) W' R& Y$ J+ O" oMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;* T9 I3 g$ @+ a) E9 P
   
2. tic;: _, [$ @\" e2 k' I
   
3. s=0;( p2 C4 b+ j* [& |/ D
   
4. for x=0:1000+ k! B, W' [, [7 H8 @
   
5.   for y=0:100
   : X: x5 q6 O# ~
6.     for z=0:100
   * J\" _, X5 x, P
7.       s=s+f(x,y,z);
   2 Q+ k) K! M& [0 X- R  i
8.     end+ U. B, E9 A( R9 @
   
9.   end% i( j, e4 ~  O2 _- P\" s
   
10. end
    $ p# p8 Z* t\" @
11. s
    5 I. r  A6 z% s# t8 S
12. toc
    9 r5 i! ]\" V. S\" V. j9 ]! \  P
13. 
    + B8 g# M$ I! C. t( K' Y
14. s =
    0 E1 B% O; A( z: Q
15.     6.126720600000000e+009( X. A% s5 X) w8 H1 x2 Z
    
16. , L% b7 y  N5 i; q1 C( A
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
' g8 q, S# l4 V1 pm函数：

1. %file xyz.m3 T$ b3 D\" [+ B\" r
   
2. function c=xyz(x,y,z)2 M0 Z\" i8 e% i& X. r! \
   
3. c=x+y+z;
   + @4 {  S8 {1 Z5 O) T* O
4. end

复制代码

测试代码：

1. tic;: W) F1 N6 V9 V+ [7 A! h. @  Q
   
2. s=0;7 b, r0 ~1 {: d) \- [
   
3. for x=0:1000
   0 |& T; _2 \4 y4 x
4.   for y=0:100  z# `; p1 Q\" T3 T
   
5.     for z=0:1005 A3 v/ V& k/ X7 j/ q
   
6.       s=s+xyz(x,y,z);
   ) R) d- [0 R7 \1 F. e7 m
7.     end! p* D9 @4 O# e3 s! o\" h
   
8.   end, L4 s3 ^5 n  T
   
9. end) v6 v4 P. I6 s' i) u
   
10. s
    5 _# k1 P# X' W+ K) |
11. toc
    & F% ?/ G5 q7 O: G\" a3 w( {
12. / R% l6 b0 m  Q; g$ j7 I2 m* R. ]: }
    
13. s =
    : e& u3 ~- u8 l( o. e
14.     6.126720600000000e+009
    - S% d6 s# B) n4 ~- s4 b7 N
15. 
    ) e& M1 @. C! j' _+ M$ {
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

, {$ n6 H9 H9 U; ]; tForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;% L+ ]! l) D/ ]; i5 U* e
   
2. mvar:  P; t1 X  t, [  R$ K( K6 H
   
3. t=sys::clock(),
   & |5 f, X* \3 u0 P& H& G1 V
4. s=0,8 i* Q, ?7 N* h, V* X
   
5. x=0, while{x<=1000,1 }: C1 j. |0 R, m6 t& w
   
6.   y=0, while{y<=100,+ k9 K- j; H9 m+ o  L
   
7.     z=0, while{z<=100,\" ~7 P( G) n! _
   
8.       s=s+f(x,y,z),
   ; i$ K+ k2 F8 o* X
9.       z++6 M; z$ B( i: q
   
10.     },
    1 l+ G9 b\" g\" m6 Y( G0 z
11.     y++
    ! m- t' d4 U4 }6 W3 }6 m
12.   },
    3 A% `  a; l4 X2 S9 t
13.   x++- S1 D, k4 N' m! @* z& e
    
14. },7 \8 U* U# u/ i3 }
    
15. s;
      C: B$ D- V5 p1 ~
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
- C3 W. |. `/ }4 |# M1.109  秒
7 P" T8 C" P) e7 C( |8 f
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
! r3 W- l, M% i. p- J; J

hxftntA

' [4 [3 b9 r5 {. r1 t3 O+ a
1 @- d& a3 N( d1 k" E

! X  f7 w# b) ?6 i
$ t3 k. Z  M  s! Y+ t( Y
- D4 u/ M7 o" e6 w. X2 G6 X& N
; ~' U" _, P  W1 N% J: ?
7 ]1 G$ U( t- F. }% ^

, Y( p6 J- [/ L" r- ]* y( f- ]


1 d' |( G- Q9 t" T5 e+ N
" R' u) C& P7 ~) T