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

forcal

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

. Y: z1 i6 l( s% L; ^& H＝＝＝＝＝＝＝＝＝＝＝＝
( v3 r/ |6 d9 N/ ]; ~+ U: e$ o2 g
4 L9 N  [, A3 W( ]" CMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   , h\" \* H0 s, O* D
2. tic;9 {# x5 A0 `' Y\" q
   
3. s=0;: M' X0 @% K/ X+ Y3 B
   
4. for x=0:10004 \1 f- @; v9 e8 D. O% A
   
5.   for y=0:100& ]! h( Y3 l3 l5 ^  g
   
6.     for z=0:100
   ; M, f* w/ p, T7 y/ b0 ]
7.       s=s+f(x,y,z);, a6 Z* k) A+ d( |2 d
   
8.     end
   7 `\" c* t; J& D, G
9.   end
   $ x$ |. V' }+ `\" o* b- W( ]( B# ]
10. end
    . O- J- I' i& y  I# y& G$ ~- W
11. s
    ( l* G# z! X+ u  t; E
12. toc
    1 |  Z2 p1 Y6 }! r3 k* [2 W( u
13. 
    - y8 ^5 w1 i7 e2 d
14. s =; F% S) M  Y7 [8 T, }; f
    
15.     6.126720600000000e+009
    , \# a8 _, L. E2 f, r$ J! Q\" F& U$ J2 r( Q
16. 
    5 Y/ a/ m; |' h  V! {8 x: u0 W( a! G
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m; ~% r/ r) S4 N0 g0 d( r: ^
   
2. function c=xyz(x,y,z)
   : @/ m( i( X  l4 P& B
3. c=x+y+z;, L: S: z/ |7 n$ f+ b
   
4. end

复制代码

测试代码：

1. tic;
   7 o& V6 e8 I, P
2. s=0;9 M( g* ~* {' L2 y7 d\" M( T: N
   
3. for x=0:1000
   ; e6 |% e; ~$ ?
4.   for y=0:100; l% D8 b: `/ |9 h6 n
   
5.     for z=0:100: S  |\" W; k: [7 @5 J1 T$ F
   
6.       s=s+xyz(x,y,z);\" c: d4 A& b# X7 ~- ]
   
7.     end1 ?1 ~- @9 ^/ T- P7 g+ J7 r; x; O
   
8.   end
   ; |. v, a' s. N  n9 W
9. end+ u\" X\" b6 I$ A3 X  f( z, _
   
10. s
    \" e; J. Z) _, v0 y+ D# a: |+ I\" V
11. toc
    & j0 w( `# o4 a3 Z% V  V
12. 1 w0 B4 g& F  B4 r
    
13. s =
    7 E$ |& x\" d$ i; w7 W
14.     6.126720600000000e+009
    , G% b# R! D% L. R
15. ! b% D: }- D' Q* v$ u
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;0 m  o\" _2 R3 l\" V, i/ \# V4 {2 Y
   
2. mvar:/ q* b3 q( y$ L+ q
   
3. t=sys::clock(),
   & {9 j9 H; ?! L: X1 D
4. s=0,% S6 Q. L1 d) L2 K- a% a
   
5. x=0, while{x<=1000,, H7 N5 D$ }0 q8 O
   
6.   y=0, while{y<=100,
   ' H: F, [+ u9 t% V
7.     z=0, while{z<=100,/ p' u# b- }- O7 g
   
8.       s=s+f(x,y,z),0 Y3 W3 \& G' r' g
   
9.       z++) Z' D9 C\" }9 w. S& t1 Z
   
10.     },
    ) Z5 D  j, x0 ?% l0 `5 N- _\" L
11.     y++1 T4 J6 `* U) h) M6 m
    
12.   },: \- \* e* K2 S' l4 P$ V& A% q; l7 L: N
    
13.   x++6 I/ V: u7 i$ W9 F
    
14. },5 }! b) u. C, V# L9 V
    
15. s;
    3 Z& ]8 l6 A0 i& B2 `( f# ?
16. [sys::clock()-t]/1000;

复制代码

结果：
" w5 t# a5 @+ w& w* b1 G) z5 G. J6126720600.
1 K( L2 V. c% {- j5 k( E: q- u1.109  秒
! }0 n) h0 ]- j7 O1 Z: t+ W
4 D; C" w6 a) |! w# I9 {二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
& B: g! C% X- Q8 H. D2 E& S, j

hxftntA

- {" I0 J7 Z( m: I$ g7 m1 `# N
, H4 b6 [. ~9 ]! J( M+ ?



, A6 n' B5 _! J4 S  T/ f# B) t

1 p7 h; c4 O5 A/ d" x

3 o+ I  e/ B: x7 N4 C% e9 N+ Q) b
" I" _# d" c! ~8 a$ i3 f: |" y

& H' @2 Q7 Q3 h' e2 M


" A% b8 D/ k. l7 Y8 |