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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。

＝＝＝＝＝＝＝＝＝＝＝＝
( h' [0 Y. {; s# o& G9 b( i( T
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   $ L5 ^+ L1 ]5 [, i\" P  ?
2. tic;( F  r5 p$ A& e2 t& p/ i, i
   
3. s=0;- y! T. j( H8 I$ z
   
4. for x=0:1000
   $ G- Y\" ]* B9 f3 T6 i
5.   for y=0:100
     W0 U7 k, R# S# C
6.     for z=0:1009 {7 o# |! H& I0 q7 i
   
7.       s=s+f(x,y,z);
   ; U# a/ J+ \( P5 [
8.     end
   % o. `) R- q& j2 \& O
9.   end
   2 j6 B* n: E' X$ `: n. f' |3 e
10. end) K% k1 j0 d: w6 |
    
11. s: Z2 @6 r( t( s2 C) X: ?. f8 |3 |
    
12. toc
    \" B' u5 y* ~0 D, K
13. 
      ]# S% m5 T5 E
14. s =\" G) y) ?$ e0 |/ L5 q; h
    
15.     6.126720600000000e+009
    . J$ {- w\" C  x/ K- O9 o; r7 |7 B
16. 
    * h: e; J( D* k& M
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m: o+ s* f1 x+ l* @3 j- j, }7 G4 k
   
2. function c=xyz(x,y,z)
   6 @5 y* i* o- p8 u6 E
3. c=x+y+z;' q/ e& @' _) y9 R! d
   
4. end

复制代码

测试代码：

1. tic;
   8 ]! D, }* l: x2 U& Y$ n4 H
2. s=0;+ m& @6 ~. T$ N, G$ w
   
3. for x=0:10009 C7 `2 p8 L$ P8 b
   
4.   for y=0:100
   # t6 h1 f( Y5 u$ a/ `6 w4 r7 u1 N
5.     for z=0:100) l: R6 `: D. t$ q' h\" @+ e
   
6.       s=s+xyz(x,y,z);
   , H' C, ?& ~( A1 }4 T/ w( j2 `! h
7.     end9 R4 I% B2 Q2 k/ c) a/ u5 K
   
8.   end
   4 q, f9 W8 q% b, {
9. end
   , W8 M3 Q4 [( k( |0 f\" E
10. s& V; U' |2 `( J
    
11. toc
    & p/ h: ?  P\" i' t5 w1 J, o
12. 
    & i3 P- {$ P. N4 ~  X
13. s =' {7 h: ?8 l; D: J
    
14.     6.126720600000000e+0096 J3 N9 U- p; L; a\" G
    
15. & t& C$ M1 S9 ]9 E0 z7 N* [
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

. Y9 f' \& D9 x! z3 C: I8 j5 t1 \/ Y  lForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   2 s1 \) p$ H, W8 @+ Z
2. mvar:
   * C9 M! W$ P) C$ V/ P
3. t=sys::clock(),
   0 @, {; z! X- }\" |
4. s=0,; A2 G  c5 x# Z. j, u
   
5. x=0, while{x<=1000,6 `* V! c6 E) N( `6 G
   
6.   y=0, while{y<=100,
   3 ]8 v$ }' }- }9 x/ F
7.     z=0, while{z<=100,; h/ Z. }% z1 H7 Q# k0 @
   
8.       s=s+f(x,y,z),
   0 E7 B\" {/ r( c# Q+ A8 c0 N
9.       z++
   1 }# X# S( h2 D) }7 {
10.     },7 w+ y\" {9 ~* y% u
    
11.     y++
    , p5 y: g& |; ?: ^* j7 J
12.   },
    7 n' V/ n9 g5 \! A6 E) c
13.   x++1 j, W2 }\" n2 g) ]  [, l
    
14. },
    \" `# J+ u( E- u- N4 U1 E
15. s;7 _! L9 m3 d4 @6 ]; {
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
+ p0 \6 U5 V- E5 e1.109  秒

3 N( n+ L7 `" k8 N3 R) L二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

0 n9 ?/ B0 f4 v
4 |7 X. d. P. k7 u( B6 c, I' S$ T
7 E% G2 N! w8 q# {( J9 O
1 {( p8 N$ N1 M6 |( f4 s



6 G( e( M) V  a
: l' g, Y$ P8 [
/ b0 I6 {( h& `* z7 X7 k! {



: a. w* K9 t, s% i
0 K+ p: A5 \& c4 N! ]/ B