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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
$ y( T" ^  U9 W0 s* ^/ d
＝＝＝＝＝＝＝＝＝＝＝＝
# ]3 P5 `, z& c4 C9 D
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;5 K3 \9 {\" c1 l
   
2. tic;: o\" v4 ?  g1 u3 ?  u+ A
   
3. s=0;
   : G& \* C1 U' L! M# v% s6 z
4. for x=0:1000
   1 l/ ]5 ]  M4 P  z9 L4 `
5.   for y=0:100
   # ]  q/ y+ t% u9 n6 O9 w
6.     for z=0:100
   6 t9 P' |\" V2 x6 ~  J6 c5 B0 q
7.       s=s+f(x,y,z);% J6 j+ C8 Z6 _  w/ P
   
8.     end; h9 }: z0 g\" `- \$ W  r+ W# p2 ^
   
9.   end
   ' P9 y& j. o7 g# `
10. end0 d\" s8 _: H2 f( U% k7 T
    
11. s\" Z, j3 F( Z: M% B$ E8 _- W' v
    
12. toc& G) f! s( C5 `\" Q
    
13. 
    9 Q* u3 ?$ B  I& N\" v5 d$ n7 Q0 f- g
14. s =: H& ~, Q- P2 R6 i, i: G
    
15.     6.126720600000000e+009
    9 [  @+ q\" o. J9 w7 e  b
16. 
    ( [: [9 U, T3 _2 K. x  b\" Y) y
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
4 S: C. z5 I( ]$ `. ?! H1 s- {m函数：

1. %file xyz.m2 y* F2 g  B8 Z! `- h8 i4 h+ |8 n
   
2. function c=xyz(x,y,z)
   $ d/ N\" h  ]. @. `6 v
3. c=x+y+z;2 Y1 b' ^. t# _% Z, M# V0 k\" f* U
   
4. end

复制代码

测试代码：

1. tic;
   ; K2 f! D# r& ~0 M. S3 D
2. s=0;
   ; q$ ~% k2 {: i
3. for x=0:10004 Y* t1 s( k7 e; n) O9 |. q
   
4.   for y=0:100! H* {; {, d( R! q, v
   
5.     for z=0:100
   , x, O9 r( l  E
6.       s=s+xyz(x,y,z);
   2 W6 y/ H- U7 b! b+ U6 C5 \
7.     end
   7 `& s6 R2 s3 t
8.   end7 O3 G0 @1 ^$ j/ T, o- J
   
9. end
   + b; _7 u0 z6 t! j- Q) B% _
10. s
    6 u  H5 P: s2 _1 ^$ B& s0 P( p
11. toc
    * I1 s! O! G- @( f& w6 j7 {8 e
12. 
    % r! b, F& w0 n% o/ j7 W
13. s =) l  z0 U0 i- k\" n. `. `+ |. @! D0 T
    
14.     6.126720600000000e+009
    / M: j. [/ P) T1 u$ |7 `5 X
15. 0 x7 U2 V7 A( w& B$ C
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
, }, o7 E8 g2 k! b( r1 h3 J( f
6 E" x( v% S2 W- r( H: ?; t1 OForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   5 {8 j4 ^- e5 `+ K3 `% l
2. mvar:! `' X# i% w) J2 s
   
3. t=sys::clock(),$ z4 f  I4 e# d+ |5 N# \, ?
   
4. s=0,% r, p* b- z5 s$ @
   
5. x=0, while{x<=1000,
   + L. {' C4 ^0 S$ J/ Y/ N
6.   y=0, while{y<=100,
   3 W/ }# x0 `6 C! r- v: G
7.     z=0, while{z<=100,
   9 E* ~0 W9 A2 e3 B, }$ t6 u
8.       s=s+f(x,y,z),- j0 a+ z( |: d
   
9.       z++
   5 e4 K; a- j% C8 e* G
10.     },/ e% r\" V! a0 D! l+ ]
    
11.     y++1 D, C! Q) w7 t\" I$ P* h. w' q
    
12.   },
    4 ?2 Z$ t& t& A1 \. T! u2 z
13.   x++$ o  q( c- _1 Y. k% T6 [& J3 c1 i
    
14. },
    & }; q2 [% j& j. _7 i
15. s;5 T5 R! X- g5 V) @* b) f* A9 l
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
0 _9 S/ ]/ f& R
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
! C1 T; A2 i- x

hxftntA

( Y4 B0 z, Q2 T$ P4 R' ~4 e1 `) n

' H( h+ v0 x/ t; z- `



  y, }( q) }/ |9 c7 O

9 }* M1 C& E9 J! J5 v# [, @
( U8 C/ I" h: [* P3 o$ w8 s
, o# \9 q/ a6 Z* d  B