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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
% S% Y  E3 `# x" Y
＝＝＝＝＝＝＝＝＝＝＝＝
6 j- T8 Z$ U7 L
) C; |# ^5 K; m0 B' I/ A, {Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;, l5 H- v- g+ J; _* ]( L% o% y$ y( V
   
2. tic;
   6 t9 a3 Q, T3 t5 X( e+ }
3. s=0;
   - |! r- \5 J$ @6 d; |
4. for x=0:1000
   ) r\" `( s' i3 R1 z. }+ r- X. y
5.   for y=0:100+ p- O/ x+ N! o  u
   
6.     for z=0:1003 g% F9 W' j( B  {; D: M
   
7.       s=s+f(x,y,z);: O( u\" N5 z2 v; T/ I& i/ G
   
8.     end. [) o. `6 j( F! g\" U: y7 E1 w
   
9.   end) `4 D3 F# j* V0 Z. V' F9 l
   
10. end
    1 E1 A6 W4 t2 A) c- R1 Y
11. s5 i, h( B4 m! g3 S! Q
    
12. toc
    - R/ f( i) K0 l; v
13. 3 E- W' `: J( x5 }  I6 i6 h
    
14. s =, e$ ?3 h5 ]% i1 r% x
    
15.     6.126720600000000e+009
    0 U0 E; s5 u9 G# G+ W
16. $ l/ f2 `- A0 m4 H
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m  f$ d' C  R) O- O
   
2. function c=xyz(x,y,z)
   8 c: J+ y+ s* E& }  ~# [
3. c=x+y+z;, [) p1 R0 q4 ?3 q
   
4. end

复制代码

测试代码：

1. tic;2 g. T7 W( f% J! G& i
   
2. s=0;
   : F3 K1 K; {( h
3. for x=0:1000
   , s' S& U  z' @' h0 O7 t
4.   for y=0:100. j2 W2 f  x6 @' |. j) N
   
5.     for z=0:100. Q  K8 ~7 o  ]# B
   
6.       s=s+xyz(x,y,z);
   + v& O& Y3 o0 c, D& F( C. w
7.     end
   / A' J$ Z+ L: B8 `, B, |
8.   end
   0 [5 L& e# H: y% e
9. end2 Z7 p! v8 X  K; N
   
10. s7 g9 E$ R9 N  A) m7 A, [# _
    
11. toc
    % Y- {# q7 N8 H5 @7 J
12. 1 p* q: |+ i: W! |0 W; h' B
    
13. s =
    9 H; g& m4 i3 _1 t
14.     6.126720600000000e+0097 R( w4 I, ^6 |# t& I
    
15. 
    ) w7 D! n6 ~* G
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ) J0 ~- B( d) y, n/ }\" ]
2. mvar:0 g, f& S8 m1 k8 ?' V6 y' V, a
   
3. t=sys::clock(),; L5 l2 k1 [6 |5 s) O5 u% Z
   
4. s=0,
   6 |5 g, b2 d; c! |5 j2 Y) t
5. x=0, while{x<=1000,0 A2 J4 h% P# s
   
6.   y=0, while{y<=100,
   ! T& G3 O* u+ m% m2 n+ H
7.     z=0, while{z<=100,2 I: K- i- p3 R% q4 E
   
8.       s=s+f(x,y,z),: O\" ~  H8 j$ |) L  E$ j  Z# w: i
   
9.       z++
   8 [% W% k$ {( B+ E
10.     },
    \" m# |+ ]. d  X8 R! K0 s$ R% C
11.     y++2 g- \6 \* m6 N$ Z: I: [; n' E
    
12.   },
    8 [- g3 X$ I8 Q- _/ U0 H% ]
13.   x++/ O9 w) G3 D4 o
    
14. },6 }! G( E- R) h; V
    
15. s;
    . O) z  M7 f& c; I. u
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
& [3 u3 I, ~- U2 R8 _1.109  秒
' {6 M0 U# y* J: o9 V7 U+ W
( R: ~" X  l. {) g  P: |二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
* j" {5 ?2 e: j. C9 x' B6 s$ I1 l

hxftntA

% f# C+ g7 \5 K9 m' i
4 ]* }7 b+ ]. `+ h) B9 W6 [- W
7 T1 Y2 I% d( d. Z( ~
8 p0 _9 @7 [- z1 q' i/ x- m

8 n: [: ?5 N3 E. \+ U2 o
4 i2 K, I6 u. f) \; o
& f6 {# @4 u' I  M* g" @
3 x, d/ o9 s% B% Z: w0 R; ]
  O- ?5 \3 i! s& g

& N0 [! ^) f/ J3 [% ?
7 {3 n3 Y: b1 j2 ?