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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
: [" J+ g$ s9 F9 R* X& ?
# ]) p+ C6 g  J! U8 G1 h＝＝＝＝＝＝＝＝＝＝＝＝

. w: H+ A+ d. bMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   & |2 y- e) v7 g4 W, K, s7 E1 R, K( j
2. tic;
   ' G+ @8 z8 ~& n& |+ }; a7 f
3. s=0;( M% @  a) G3 P$ J% J
   
4. for x=0:1000
   ( G  v  R: V. d$ B7 N
5.   for y=0:100' M5 ]; G4 q5 x- g
   
6.     for z=0:100
   \" m* _7 t  G+ o0 K9 H5 T  L( n+ X
7.       s=s+f(x,y,z);
   # W2 W6 R) P  U  a- Z' P
8.     end
   8 D6 \! u2 D: @
9.   end
   0 z6 P9 z4 K' m* G
10. end
    / w; v. p) H' b) F: d2 F9 b
11. s
    - T; q3 Q: {' |0 d7 y
12. toc! G* r5 P' F; j+ }4 N, F  x
    
13. : [! ?' a, A3 B1 ?8 {8 e
    
14. s =
    7 Q- ^) r2 N. b
15.     6.126720600000000e+009* o0 x2 N  A9 o$ `1 m) J
    
16. 1 ?( R. T# `' P/ U
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
) V' C/ [- e2 z  `! Z3 G4 jm函数：

1. %file xyz.m
   : ?' Q+ S* K& n
2. function c=xyz(x,y,z)6 G: k\" d7 ~8 m6 ~0 d9 [
   
3. c=x+y+z;
   % E& H( f+ x# ]3 g! |, ^: i
4. end

复制代码

测试代码：

1. tic;( w8 ]  G; V) B* r
   
2. s=0;$ t; K( m5 E% Y( z* [1 f
   
3. for x=0:1000! M# P' C  o; O
   
4.   for y=0:100. }* C/ F% a* _3 c2 x: F
   
5.     for z=0:100: \, \7 o  A$ E1 v& v0 I3 X
   
6.       s=s+xyz(x,y,z);
   + M; f. g9 K$ }. U6 I
7.     end$ ]0 ]5 f  e\" t0 m+ ]
   
8.   end& R/ C  I! j0 d' L
   
9. end. E/ T( S: E8 U/ a$ Q\" B5 M
   
10. s8 F8 J' k/ Z. ]0 ~4 q  L- L* F
    
11. toc9 W& S- M' V\" S9 D0 Q( l* V/ H4 I  t
    
12. $ z; M( O% y( U
    
13. s =
    . L8 @# t% Q/ m: B! i
14.     6.126720600000000e+009
    5 B6 S5 M3 B* e( S& U8 d
15. 
    ! ^9 U% ?: {0 J
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
) _! W0 s& ]' B; d0 c" n0 }8 C
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   8 ]( x6 i+ j3 _0 w! s; E! ]
2. mvar:
   7 k) c% R# d- [5 ?1 O
3. t=sys::clock(),7 P' V9 J- j0 e! S
   
4. s=0,/ B+ f6 o# z& q# e$ Q  ~/ }9 y
   
5. x=0, while{x<=1000,
   ' G\" V& R- s+ f2 x% B/ i  Z4 Q\" m
6.   y=0, while{y<=100,) [) d! s$ [+ Q9 ~; B* B
   
7.     z=0, while{z<=100,# }8 a/ v7 t\" Y% J9 o+ y0 u  p5 O: }
   
8.       s=s+f(x,y,z),
   , J- Z: B8 C/ @+ Q! [1 D* C: w! U
9.       z++
   ' w7 R7 h! u) |
10.     },  ?- w: p! ~5 y3 t/ c3 \: ^* z. g  Y
    
11.     y++- p! T5 _  f! S8 _2 |5 ?& C
    
12.   },' y4 h6 V( J# z/ h
    
13.   x++
    $ k; w* g7 {, ~& c# J8 M
14. },) w; p, s' p+ c( l
    
15. s;
    7 u/ _; c3 R* R0 G
16. [sys::clock()-t]/1000;

复制代码

结果：
- p1 K* J& y$ a' Z7 e7 l9 t6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
5 v# V' U( _$ c: X, b3 r' M  v- W  y

hxftntA

5 f6 a" [6 O, o* J0 e$ I6 k

% A+ L3 {+ U0 u* U- J3 m4 e6 G

: H; N$ r* c% \
0 K- r' i. Q6 F& F$ y( K  U6 h* Z! K

9 b& X0 {9 {: ]% k7 D/ X; r8 u

% n; z4 j! I5 K" ~) \
- [  b+ ^0 M2 p8 W3 c



9 h8 p: f; n( t