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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   3 @0 U1 ?2 N% b
2. tic;
   ' g6 b, P4 R( j% d. N: B- _- D\" B
3. s=0;
   & W6 X& m( p7 }/ d$ j9 M
4. for x=0:1000\" [3 ]* q& w8 Y3 _
   
5.   for y=0:100
   % g  ]% C3 a) T2 B# {
6.     for z=0:100' S- ~* I: h7 n\" `; [  j5 N
   
7.       s=s+f(x,y,z);
   + p. A, }+ l  Q5 h& \/ o& i
8.     end, ~% F1 [6 N* h
   
9.   end
   6 J! v# f  z2 r, j- R
10. end
    9 K7 R' {. K9 d; `
11. s; Y0 |' R$ m. w. m8 G% E( ?* [
    
12. toc
    2 H' H0 ~  ?6 F& S; `& n2 [\" R1 n
13. 
    7 v. B# s4 R# k9 d! K% X4 t
14. s =
    & `4 t* \4 b8 s& w: ?
15.     6.126720600000000e+009- t# F7 p: Y* Z
    
16. 
    + f# P& g! [7 w! }- E
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
9 K8 W  b( m* a, K6 [2 vm函数：

1. %file xyz.m8 y- q3 Y  H: k2 i
   
2. function c=xyz(x,y,z)( _0 N5 r/ k- n, x# B$ Z' }! m2 `
   
3. c=x+y+z;) x4 \/ p% i1 u
   
4. end

复制代码

测试代码：

1. tic;
   4 p* a- @7 y- T) ?! C; w, `2 C
2. s=0;
   - O0 w; d\" b8 R% G. G
3. for x=0:1000
   3 X  C: C: u8 {1 W; o* i
4.   for y=0:1005 K3 C( Z' h. g& ]& n9 N6 e
   
5.     for z=0:1006 X# N2 I; G* r: j
   
6.       s=s+xyz(x,y,z);
   * L) z+ s! p' V\" D7 I1 V& k4 Y
7.     end) Q0 \! D1 _0 e4 K1 j
   
8.   end
   ( x4 I8 A) h\" P7 T1 ]4 m3 c
9. end9 \8 A: w& ~% m) d
   
10. s
    + Z3 z7 D& T: P% ^& H
11. toc+ a( f: {, v( r3 e& _  M6 w
    
12. 
    ( G\" J& r; }& J& W2 @
13. s =
    \" O, `; G1 t, g
14.     6.126720600000000e+0091 p& Z0 K6 n! k5 u+ K
    
15. 
    4 X6 l& |4 T4 X) G  _6 r
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
  [/ i" r6 p. M. Y
, T3 t5 @: F7 Y; _7 u, p+ @6 |$ PForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;3 r( L& l7 \$ v+ G
   
2. mvar:& o: E3 B1 D6 O3 H
   
3. t=sys::clock(),
   4 L# [* ~$ T3 j. o- W. I
4. s=0,
   , b\" n* c% q  ~9 r
5. x=0, while{x<=1000,
   ( T0 ^$ I/ R% o- x% F/ A
6.   y=0, while{y<=100,
   0 x' g# z/ v' }2 x. ?# v2 {
7.     z=0, while{z<=100,\" \  D) U* G3 o
   
8.       s=s+f(x,y,z),
   1 b& y+ g& O2 p4 X\" Y2 G
9.       z++
   - P4 S$ l1 L( U6 L
10.     },' \9 ~7 E- }8 L) R! w
    
11.     y++
    + L9 q/ ?+ ?, I6 m! y5 e6 K
12.   },# k, A9 W( W4 }, Q/ O
    
13.   x++
    $ V( o5 W3 W# M2 A! {
14. },  d% P5 b  ^. h; f) G2 y8 K
    
15. s;
    # F7 M, {$ G\" {) x6 _4 S; E2 L
16. [sys::clock()-t]/1000;

复制代码

结果：
! {% G0 C6 u2 P# o9 `' |6126720600.
1.109  秒
9 W4 H) @9 w! q6 n% c7 @
* U  S, T2 Z; B# M0 f0 W' k+ u二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
$ \9 ?; Q6 D! w. m1 N

hxftntA

& X& Z+ J! h' {/ s4 u8 b7 ]
. }6 }* j" @4 h' x5 F% e7 E
: D3 O( F$ {  l; J6 `& D4 G- _: n+ w
0 o: U4 a& D; B8 X  E

, p- s5 R* ^1 j: |  D


( M+ Y7 G5 j& ^8 J8 B

( N/ l2 T! K# D- W9 A: [. w
/ O* m+ z. {" N0 s; q( N. K& X


+ f1 |3 s: X) f) ^) v
' G( Y& ~: g/ Y4 T  y2 N0 p3 @