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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
; Z2 {8 t7 g7 L- }8 _
＝＝＝＝＝＝＝＝＝＝＝＝
, W2 K. [+ Z4 p6 R& S" @' P+ V
: m" ^+ }3 {( ~. m# p3 y/ O4 \' w( \Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   , N. _0 w5 K) }6 i$ Q
2. tic;
   # Y1 F5 P6 M3 i4 [' }) T
3. s=0;8 n- `5 O; _/ o( O4 x* p
   
4. for x=0:1000
   7 k$ C' n1 O2 }, _, u& I& K
5.   for y=0:100
   2 P0 {; w  ]8 Z\" o
6.     for z=0:100
   . _1 V4 p# v+ \( y8 e4 S
7.       s=s+f(x,y,z);
   5 r0 ~. u& T! B' O
8.     end* L4 F* M- h2 h. X. n& {7 ]
   
9.   end5 E0 ^6 ?1 G7 \5 E& G/ ^  ]\" V
   
10. end
      N+ q' P6 _; i3 ?& b
11. s' m\" ]\" T* t- I; e. r6 R* L1 I- u
    
12. toc
    & u4 M2 r  u1 e# ^9 `8 ^
13. ) z+ X( k7 @. n. }; \3 g
    
14. s =) j. Z3 S' F; L! m( z- Z
    
15.     6.126720600000000e+009
    2 i! p( c5 ]# S2 a, W* f
16. 6 A3 A( P\" ~: z) O
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   & {1 r+ {' \9 j
2. function c=xyz(x,y,z)
   . i# A3 T! w+ P: o0 T6 x
3. c=x+y+z;
   6 F! H6 Z# m) a  n% V
4. end

复制代码

测试代码：

1. tic;
   * C- Q4 |7 |: l9 M& Y) D
2. s=0;- I( I; N, }  W3 b
   
3. for x=0:10008 ?$ x' _  D) r0 p# s
   
4.   for y=0:100
   $ x) a\" w; F( h
5.     for z=0:1000 F) G! ~) p5 R! X& p
   
6.       s=s+xyz(x,y,z);
   # l' W( \, ^. q! Q6 V
7.     end# S; \2 h/ w9 R' a3 D- @
   
8.   end6 R2 I( t  D8 [  R
   
9. end
   $ |3 _; i6 E/ \) z: ^8 E
10. s3 z0 }: i3 K* _- P/ I\" _( @/ V. S* f
    
11. toc
    + g: _\" p8 C+ ]$ t$ K: r
12. 3 K/ p& l4 ~. t  c/ o
    
13. s =0 ?2 ~; \: \0 i# Z+ I
    
14.     6.126720600000000e+009
    7 C! }0 h# Q\" f2 t\" n7 }8 ]
15. $ H  x9 P\" ]8 U
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
' n3 e! _) S4 q  {) C1 H
8 }1 P- F' F. z: y. aForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & A3 |+ }/ v) o9 B5 @# J1 H* O( g: r
2. mvar:
   4 c/ K; k( s4 i, a& j3 `4 S/ y
3. t=sys::clock(),1 Y8 G9 w& {7 i
   
4. s=0,: Z2 I: e: H5 [1 d+ z1 c' L  t8 c
   
5. x=0, while{x<=1000,
   - l1 s6 X' ^# O3 X, }
6.   y=0, while{y<=100,3 e8 g. P3 V. A4 K/ Z
   
7.     z=0, while{z<=100,% }) d% G* a' D; q
   
8.       s=s+f(x,y,z),* v. Z, j' O2 S$ L\" l
   
9.       z++1 g; D$ M& K6 @' q$ W6 `
   
10.     },+ g' r9 c\" p) v1 d# t7 {& [1 \
    
11.     y++
    + d& d( H\" y; R- m3 }  ~( \
12.   },
      L* R; Z. R  a\" h  t' j) q5 t
13.   x++
    3 k, H1 f6 Y+ B7 r% y
14. },  m9 ^0 Y# U5 g2 h/ n
    
15. s;  m5 A9 g$ Z\" w
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
8 P" }% E* Y% `- v+ I1.109  秒
6 m& q/ Z0 q! ]7 K
- T9 Q! Y& a# D/ V% _$ w二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
* t5 a* {) T, g7 [& q% N* @

hxftntA

5 Q0 G6 `: f! {! \. |. E8 q0 i
( m1 Q& y1 N* _6 g3 ?( u
3 H( T0 F. _& I' o
5 S) I, ^- D% [' b) ]

! {4 J, G* c7 t+ W  E0 n- L

9 y& J; n5 K7 |: T

; t* J# G# @' K7 r4 e( A
! M4 Z6 V5 P8 p

1 o: A" `4 P  P( ~) h- Y+ t
+ {/ e$ \' h2 P1 A* n: q