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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
1 D! S5 p7 |5 s* H8 E' k
＝＝＝＝＝＝＝＝＝＝＝＝

" h4 s" d& g* c% rMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;: q* H: e' w\" ?
   
2. tic;
   + N/ D4 O6 Z* I$ u# z
3. s=0;
   ! }6 i. k0 y4 P/ ]
4. for x=0:10005 E* E- O7 b\" Q6 S3 n7 G
   
5.   for y=0:1009 M  ^) ?\" k: n4 v% I- N0 G# ^
   
6.     for z=0:100
   ! M0 ^1 S  C# n# Z( t, ]
7.       s=s+f(x,y,z);+ F4 k+ A. E8 A$ X# n# F
   
8.     end7 ~: v! v4 S8 P- p3 V: D2 R5 d
   
9.   end
   , W\" z) ^2 T2 e8 P
10. end
    % ^) f# G0 f- O# h* F) g& g$ K9 S
11. s
    % i) X- d, V' a- q\" @9 Q5 ?
12. toc. A, j1 `; r& \' E+ C
    
13. . U% i, H9 w( o( a3 Q  f( H/ i7 L0 V
    
14. s =' H/ R- c/ X& z- U
    
15.     6.126720600000000e+009
    5 S4 k: \' I  ^. B+ m
16. ; E1 \9 g7 ^; z
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m- u  c0 l9 s# ^' h# c% g$ P
   
2. function c=xyz(x,y,z)
     `3 N6 \/ x: R; n9 `/ c. l
3. c=x+y+z;: c, v1 y5 b/ u+ O6 k+ E9 j% ^
   
4. end

复制代码

测试代码：

1. tic;
   7 |& \! d% q6 I4 b* X, r  d/ \
2. s=0;. m' c5 _5 h/ l- \2 v! A% M0 ?0 G
   
3. for x=0:1000$ t7 F# D! Q. g8 k: p9 Z: V6 L+ p
   
4.   for y=0:100
   + ^) {  j\" q0 i; a
5.     for z=0:100
   . [, }& ?% x8 l
6.       s=s+xyz(x,y,z);& u% H. q3 E: W& F5 i
   
7.     end
   3 v0 p( c) a. _5 e1 x& P
8.   end, a$ |7 W# [6 o2 V
   
9. end2 ^9 d5 c5 @# u' z. P
   
10. s
    $ C0 K# }! S3 X8 Q* n$ T' A+ X# q
11. toc
    0 u% Y; d2 N! c: i! E/ H; y
12. 
    % Z, k! L  C: A* x$ I* \
13. s =- f6 t3 h: |3 k9 z
    
14.     6.126720600000000e+009
    , {7 X1 h- ?8 C4 D5 I! U
15. 6 x- u) q9 m4 x2 B) L/ |$ S! i- A, l
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   6 Y$ [  f: Q2 m4 e6 n
2. mvar:7 H\" I/ u2 h, v, Z; K
   
3. t=sys::clock(),
   1 |: [4 L; g' M3 w( q; A, y- D1 J* k# y
4. s=0,
   ; ?6 |2 e5 h* p2 k. J4 w  X
5. x=0, while{x<=1000,5 f& q/ ]2 W$ A' \
   
6.   y=0, while{y<=100,
   5 ?( n! ]# w\" j$ C* W
7.     z=0, while{z<=100,6 l9 ]\" T# Q2 ]  h: Q5 F5 i) v/ d. o
   
8.       s=s+f(x,y,z),\" Z5 ]+ A6 D& [$ ^. ?( R
   
9.       z++
   + `& Q$ P, U# D3 E: N( ?6 z
10.     },
    4 ~4 V- w2 M2 L) Q6 X
11.     y++
      o. h7 n6 g# A( b$ h9 m7 v. X. s' |
12.   },) P9 l& q9 x: t
    
13.   x++% I4 }( c, S; @4 {* s! h0 [
    
14. },
    8 T/ B' _2 i0 ~( y
15. s;1 F! ^& [5 `& e. H* X( S
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
+ d8 v8 V  u& O  D2 a' ^' r1.109  秒
+ [  O# E  a# \8 Y& L! ~* o& H
& O* M" Y7 V3 F4 s二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
2 B7 {% c5 F; V5 t/ m5 W. i

hxftntA

: h2 `: T0 c+ k& C5 j( P

' V) T1 ?. B! t- K

( p: r4 U0 N" A  l
/ J' r3 C* A$ q7 z: f9 O0 X5 ^
  J0 @; g  ]0 |' F
' X0 ?. @9 D8 U. x) t9 f



. ^0 `- n* Q+ f4 D" G
# f' D. m$ v- c. g


% v  k$ M6 I( |. y4 n' X6 d; g5 P
) \4 d3 t+ j* Y8 }. p" U. ?% n