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

forcal

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

1 A- {0 m* D4 Z6 E＝＝＝＝＝＝＝＝＝＝＝＝

0 i! w5 Z% g3 i+ U% YMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;3 d5 c# N' Z! w; `; ?  r
   
2. tic;
   3 Y4 l7 Q% H) ?
3. s=0;5 W# L2 e8 G/ D1 c( ]: a% L
   
4. for x=0:1000
   . n' G2 D/ D2 @: e# v6 R) u+ ]- @
5.   for y=0:100
   8 U4 S5 X3 ]7 q8 G- N- U: c# ~$ c
6.     for z=0:100, q; R$ X& y3 G2 F$ Y5 l
   
7.       s=s+f(x,y,z);
   $ ?# \1 R% V. i: L
8.     end% n7 P- l& G$ Y8 }/ g
   
9.   end: E9 n( t/ Z0 F7 D: _7 j
   
10. end# G; K, j/ }: h  B
    
11. s
    \" k; Q( \$ j0 u\" r* P
12. toc% w; i) L$ B: }6 ?3 S, ]5 F
    
13. 
    & R$ k9 \* B- v% s\" R
14. s =3 t+ J\" W3 |4 s- _$ ~
    
15.     6.126720600000000e+0095 p- x/ I7 T1 l4 x. \
    
16. 
    # ?- m( P6 ?; U! {, K/ l
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
$ p+ ]( i5 ~6 r, M: B( H/ xm函数：

1. %file xyz.m5 c\" d9 V0 u( j7 ]- K
   
2. function c=xyz(x,y,z)
     ?# V' f( K* h9 i
3. c=x+y+z;2 A3 E7 q( G0 D
   
4. end

复制代码

测试代码：

1. tic;; B9 d- O5 c& T; O
   
2. s=0;
   ! q2 e  K# q8 u\" d
3. for x=0:1000$ P# g2 K9 J* y3 s6 k
   
4.   for y=0:100$ _1 c2 P% Z1 ]4 K' H5 s$ w
   
5.     for z=0:100
   , t, O: A+ {+ q, R# ^/ s
6.       s=s+xyz(x,y,z);) ~4 \9 V5 B* i5 x2 i; t\" X9 ]9 f9 X4 |
   
7.     end! j8 i3 c* e# b9 v6 u5 v- Z% V
   
8.   end- ^\" Q) P- R& D8 X
   
9. end
   * ?3 t# x9 R7 G0 k4 v
10. s
    5 }  X! G- @* d2 W! n
11. toc
    , e/ [8 j1 [! Z4 j8 x\" s) i% J1 w
12. 
    ( P( K/ P5 U+ h7 Z& l
13. s =- Y: q/ j; Y\" n. x6 \5 Z
    
14.     6.126720600000000e+009( u( N' k\" c3 c\" ]\" a. {. A
    
15. * [\" r: L5 Y$ J( h6 v\" c! J; L5 _
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;& e6 l\" t/ q' S3 Y% d) u4 h6 r' p
   
2. mvar:
   + b1 J+ [! M9 p0 i. O9 q
3. t=sys::clock(),
   \" Y( z5 Y( X+ _8 I5 ^$ v2 e- J0 {
4. s=0,
   ! O1 N/ j/ c) z7 @* c% e- Y/ ~
5. x=0, while{x<=1000,
   , M3 H* _\" G\" P
6.   y=0, while{y<=100,
   4 S% T- Q9 H7 u0 z
7.     z=0, while{z<=100,+ Q6 M: N; J/ s7 s
   
8.       s=s+f(x,y,z),
   \" N8 o7 Q3 O2 {8 G7 e+ d
9.       z++0 P/ ~. k: _5 m3 o+ x  W. A9 A
   
10.     },3 U$ q+ m- l0 f
    
11.     y++
    6 W0 E* A: W! Z, _8 R9 `
12.   },
    5 U# W; l- @; m# ?, R1 [5 w
13.   x++1 B7 X0 I5 P) A( }8 z$ s0 |. }0 s
    
14. },: V7 w( E5 Q- Z* ^- f6 X
    
15. s;
    1 I/ j  Q/ H. K# w2 i5 ^. t
16. [sys::clock()-t]/1000;

复制代码

结果：
; k. q+ ~7 f# A1 y' Q) n6126720600.
$ R" @, p& \2 f1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

- C9 d5 [4 B3 y8 S- D

  n6 V7 r9 N7 D! r- N( S# _

% f: @' N  d& O0 Y, U- W. L
: a3 C8 F2 X5 T
4 S. G6 K& v- N9 o0 G5 H7 X& [6 }

$ T8 o4 A4 q( ]1 S+ L
) j3 m8 t  h0 e) ?


) N' Q( d$ j3 p
' U% p9 J, g2 X$ w$ A5 d' I+ g

& m. W' b- c) R2 y$ x