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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
9 l% r/ F  A9 x4 S; w
% I9 `3 ?' w5 p& H. y% TMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;% R, U& Z+ ?\" c! M# {( l1 R
   
2. tic;! \! ~$ _5 \3 v4 f; ]\" N! H
   
3. s=0;
   8 @0 n, y: s2 N% o; u+ |! c# X
4. for x=0:1000  U. A\" D/ ?9 C
   
5.   for y=0:100: D\" j6 \9 |$ V
   
6.     for z=0:100
   7 X% f' s# _/ t7 u
7.       s=s+f(x,y,z);
   2 U) F0 W3 L% O6 {
8.     end
   4 O, e& R1 }( ^6 {( c
9.   end
   : f\" c; g4 J* N2 |- U2 g
10. end7 R( ]% k\" Q& H9 a) h8 i+ h8 d
    
11. s
    . \\" |, K9 f' Y9 ?, `3 c$ X4 c
12. toc
    2 T  Y, u; L( m; j$ K- ?
13. 
    1 m+ P9 `. v( z1 O0 K: o% }( d
14. s =
    1 K$ j- A0 L2 U( X
15.     6.126720600000000e+0095 M6 ~) w2 K9 i- T
    
16. % s% N9 y: j* v
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m( [\" v2 _3 w3 n& Q* v2 }* l: c
   
2. function c=xyz(x,y,z)% ]: s: s% ]  ?
   
3. c=x+y+z;( q* E6 R. H/ B
   
4. end

复制代码

测试代码：

1. tic;
   ; V* }/ J2 D0 W
2. s=0;# @0 \: @1 d: v3 Y' t
   
3. for x=0:1000
   \" Q. K9 H7 D* {2 Y% M
4.   for y=0:100& ~4 u  O; N4 }1 b
   
5.     for z=0:100' W% ?2 y$ a9 J) h9 |
   
6.       s=s+xyz(x,y,z);( t, p+ i* S, C% z/ c6 H; e
   
7.     end* Q% |& R! Q- U
   
8.   end4 w; b3 l( v% g\" O' L# m
   
9. end
   * e. s! g8 W$ C  j0 i$ c
10. s
    - l) l& M9 J, m* N1 E
11. toc  e. [# S* U+ ~' H! C
    
12. 9 h( p: k. z( B  i' e' R0 Z* v/ I
    
13. s =& \\" j  t% q\" Z3 H2 a# l& Q
    
14.     6.126720600000000e+0091 I- n\" p0 M. b; |, ^: T
    
15. 
    4 \% _1 Q0 O+ s* {
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;. r9 n( J$ Y1 I: m1 `5 B$ m+ c
   
2. mvar:
   * T0 b* {/ W* a6 x8 P6 }+ j
3. t=sys::clock(),
   ' @0 Y* \# \0 V! D$ [
4. s=0,
   ! j  w( X' i, F) x) K( n- n
5. x=0, while{x<=1000,
   ' I) y) _\" e; I: y\" q: p
6.   y=0, while{y<=100,
   : p, t' Y- Y: A
7.     z=0, while{z<=100,1 u  Y9 V$ O\" Z8 e5 y: f. U+ p
   
8.       s=s+f(x,y,z),$ S0 o% K% }2 r  e
   
9.       z++/ m7 I- j* ~( c, |) p0 y
   
10.     },
    $ f2 U2 I! q# N6 R- O* O
11.     y++
      B& ^9 r; q4 \& F  S+ @8 s
12.   },
    6 D) ^2 s5 k* m- t# P
13.   x++7 z3 a# |4 Q8 C* V1 n3 d
    
14. },
    , d$ m8 B% I! I6 D( y\" d9 A
15. s;# s% L5 c! Y. u' Q; J2 E
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
% x& {( x" B; Z8 h8 G
, A3 I( n  h0 u1 T& r# _7 P/ R二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

1 b. u" b4 F' f: v0 W7 Q( [$ f
6 |4 \) C- a+ f$ ^) y
0 Z2 U  i1 u& W/ M$ _





9 r: {6 N  k2 {7 @3 x+ b  ~# @