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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
$ r% B' o5 c; J( P4 a
! ]; a- i) Q3 I9 B: y) r4 T! m6 q＝＝＝＝＝＝＝＝＝＝＝＝

- G% n( q  ?  V( P3 q/ Y# }- f( gMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   % X: g& |! B( u. ]' N
2. tic;
   ! E\" A3 s3 R: k3 g
3. s=0;
   - R# [( u7 g1 U
4. for x=0:10006 s* I2 w) O- s' L6 P9 k
   
5.   for y=0:100
   & d3 L' b( k% x1 x2 z  v, J( k
6.     for z=0:100: V4 l! x  O3 `3 x* [5 g, i) |
   
7.       s=s+f(x,y,z);1 U; C\" A1 q% K0 o* D
   
8.     end8 {& p: f0 O* e9 n/ M' H
   
9.   end
   0 w* }- x/ p% \$ n; J) u
10. end3 K3 z8 S9 F: L8 y) r5 v, U
    
11. s# Q5 K5 y\" c; f* V  d& ?
    
12. toc2 A+ z, c6 |, S7 o7 h\" S6 O2 ^! r
    
13. 
    & y( J6 c  W6 o( e\" W! o
14. s =
    1 w& B9 y/ C& `: B- a/ _
15.     6.126720600000000e+0091 e& {3 z) X6 j+ J1 ~9 D7 H& F5 S
    
16. 6 E9 A/ u9 V9 [, |
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m  G2 Q& ]$ {. n8 L9 W+ |
   
2. function c=xyz(x,y,z)7 I+ H/ ?) b8 }5 s& s; y
   
3. c=x+y+z;! f0 P+ H1 }: \. c8 c  P
   
4. end

复制代码

测试代码：

1. tic;+ x; ~, i7 H* J% H, `/ r& h! J* V6 N
   
2. s=0;2 Y0 T7 m- R5 S( r+ I5 p
   
3. for x=0:1000- A8 K. E& \8 s6 l
   
4.   for y=0:100
   8 \$ H  e8 |7 O% j% A
5.     for z=0:1007 f! c1 m- z* q4 ~6 p\" T) U
   
6.       s=s+xyz(x,y,z);
   , l+ h0 w  w3 J2 L! t
7.     end3 \( S/ s0 T9 D4 B* M  ]
   
8.   end
   3 J: P' ?! L2 R9 e3 J
9. end8 j* F) X, ^0 I0 ^& a  j* u
   
10. s; n$ ^/ y1 D6 v9 f
    
11. toc; f# I! t* V/ w\" X  [) w
    
12. \" M$ r- A( B9 [
    
13. s =
    , @: W$ s0 k6 ?) M' F4 c
14.     6.126720600000000e+009
    \" K  g) g- N, r7 v+ f1 c, B
15. 
    ' z* b\" C# r7 D& E
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
' D- [9 i4 S5 x- h" U9 W2 ^7 V. i' T
$ D$ w. v( }' @6 H7 }  n  z, U2 NForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   - g4 f: M9 I' M( \3 i+ D9 }
2. mvar:3 d& C8 H/ E8 n  A8 u: }\" {
   
3. t=sys::clock(),0 |\" v. ~) f0 g
   
4. s=0,
   / ^0 T: X* _1 M+ |3 R
5. x=0, while{x<=1000,
     `' W* d9 V6 ~; }% Q0 U
6.   y=0, while{y<=100,
   2 Z( i! }9 K- w4 {2 r, n
7.     z=0, while{z<=100,
   9 q* c! P; }- u+ j- \+ a5 J! ?% ?/ E
8.       s=s+f(x,y,z),$ ?( i( D. G- M; g
   
9.       z++
   # V9 {  W' [\" l
10.     },
    . a: h' Y8 [% E$ k5 r\" }
11.     y++
    7 {9 f; K5 ~7 d/ c) n$ Y2 H
12.   },/ B/ o$ s, S4 |
    
13.   x++. `7 p$ O; J+ c+ @! R) p
    
14. },+ d' Q6 i/ y* T, ]( W! S& v/ v
    
15. s;+ s' ?8 ~; D$ D2 r2 G
    
16. [sys::clock()-t]/1000;

复制代码

结果：
# H0 f( G* X0 J1 w/ V6126720600.
) @# [8 B# S9 }4 |1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
. V" g) {/ L- [8 v

hxftntA

0 x" r; t3 R# W  Q: o( h/ A
$ q; Q/ z0 i4 U, e1 w% p, z+ N6 i


& o! J2 |9 {- U9 m

, e4 i6 I: ?1 ?: J
- o. p  c+ x9 Y# e9 |/ W4 s



( S5 g8 N& T$ O8 A0 f
4 n' h, r( J+ s$ @% K