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

forcal

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

* b+ l: y) O- M% j: V" n＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;/ [/ J1 ~3 Y3 p- l' V& c\" z. N+ x
   
2. tic;! p3 c! q+ A! r9 g: B
   
3. s=0;
   - [2 p8 q9 N; x, P/ i! W
4. for x=0:1000$ k: M+ r\" \: Y# I, q
   
5.   for y=0:100
   5 {8 \. i0 h, M# m3 \\" k
6.     for z=0:1009 `+ J0 R% Y4 @' D
   
7.       s=s+f(x,y,z);9 W$ W6 j3 E9 T+ ?. b6 \: U
   
8.     end- m3 D, B7 ]7 p' F
   
9.   end
   7 A1 T* o- |7 l5 R, H! L* q6 C
10. end) @& L& _& ]  r# {) z
    
11. s
    8 X- x/ C0 ?8 U! S+ q  M$ `\" D1 E
12. toc\" V$ |) U2 x( f8 U
    
13. 
    ) k9 U: G9 ~: H7 R: `1 x4 ?, f
14. s =% q. d% F! ]% V$ X! o
    
15.     6.126720600000000e+009
    % H\" F1 H9 t- w4 X! |7 f2 Q
16. % k! q- H# O! E, c3 B; L) a0 A
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m1 J\" U5 [8 u+ x) r! C- g
   
2. function c=xyz(x,y,z)5 v+ J+ P/ H. C2 L
   
3. c=x+y+z;
   ' T4 g+ N+ s( A$ F7 f: Z, X, p/ C
4. end

复制代码

测试代码：

1. tic;
   * y. b5 I& v/ a! d6 D
2. s=0;
   - t$ o0 o\" |0 W4 J  G* ^\" j
3. for x=0:1000
   # ?8 J- N1 x( x5 f2 a( P6 {! y6 v' U! K
4.   for y=0:100
   ( ?* k4 r/ q) {+ T
5.     for z=0:100) L3 \& I( @: P8 m- Y
   
6.       s=s+xyz(x,y,z);
   ! w6 T: Z7 A$ H7 l
7.     end
   4 R3 u7 q* h$ B9 Y' G
8.   end
   * u1 V; x- q$ c% |\" Z
9. end( x: ^! \8 T2 I2 t; Z: g8 y
   
10. s
    $ e% U( Z+ x% ~. ?2 g! K
11. toc3 Y4 T2 [7 D\" [3 a5 e/ D. S
    
12. 
    7 m& f5 {- u* m\" H
13. s =
    * P) |4 A% T! X% @. y
14.     6.126720600000000e+009: v; I2 F( q. _5 z
    
15. 
    + H2 g& M% g1 c# l! Y0 m  _
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

; G) w/ J$ D( z0 q/ }& }; XForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;: t& d7 X8 G4 {5 ^! n& a8 P$ D; q
   
2. mvar:
   / c/ {0 u' @3 q  _, A$ O! Q0 c
3. t=sys::clock(),
   2 o( ?2 p! K' F( `5 f5 U) Q
4. s=0,, O2 e\" V5 I' R
   
5. x=0, while{x<=1000,1 X' N* b& ~4 k% |4 \% B
   
6.   y=0, while{y<=100,' g2 i4 r6 v$ A
   
7.     z=0, while{z<=100,# ~5 b* Q+ u  k  ]$ ?1 P+ W\" i: B
   
8.       s=s+f(x,y,z),( b: s' ?& |7 `/ r
   
9.       z++- f( b6 Q; S4 X: O1 j0 r7 V, N8 Y
   
10.     },
    - d7 O. h: U0 y- p2 \  ?' m& ~9 I5 H
11.     y++
    : i- {. O# n$ f$ W. v3 [6 T
12.   },2 m9 i$ @% u& L5 b! t+ E
    
13.   x++* u. `\" }' [/ j' U+ r
    
14. },* l0 B& X; p! ]# a) T& F' g- _
    
15. s;4 o) s2 i8 B( A' t* [$ l0 e! G* L  X
    
16. [sys::clock()-t]/1000;

复制代码

结果：
- ~1 M( Y. A% J/ o5 u; I  X  G6126720600.
1 B: K- l! }  j7 R; ?3 Z1.109  秒
1 A% s3 \/ p  n: ]6 z
% e4 s; ~: o9 t7 w二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

* D5 ]8 A8 A& b0 N2 k) `" l1 X$ z" @





) h7 |2 q( Z- P% W0 t
$ s# ?6 O% B. S3 ~' v. y
5 {2 O  |& w1 H( G4 \0 h; y
1 Z7 Y2 [0 a2 D3 U1 d4 s& Y
& H$ F: X* }1 U# @; w


4 w/ H0 Y: ^( @& s3 s! n
8 C1 F1 P( e6 y

3 w. n1 ~# t! H