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

forcal

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

: _: Z; q2 V. P5 |1 Z/ R＝＝＝＝＝＝＝＝＝＝＝＝

( r+ K. D+ y9 ^  K3 ^8 MMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   - N8 E- Z4 Z2 v, v; M3 M
2. tic;( @$ h  C( Z1 U5 S
   
3. s=0;
   4 _8 m) d) f! S9 b5 r
4. for x=0:10000 w: q- k4 k+ M: A* O2 R6 F
   
5.   for y=0:100
   ) B\" D1 m3 H+ |: l: u; T8 |! T
6.     for z=0:1009 a) J' W4 O9 A. P7 [; \1 z\" u7 K
   
7.       s=s+f(x,y,z);' |5 i0 V8 |6 W8 P( W! d- [- N, `5 ^9 N8 }
   
8.     end
   : m6 @( ]8 o6 f
9.   end
   2 w0 d. @! w3 V- Z( @$ _- u
10. end1 g' P1 R! A! f% g- Z
    
11. s
    % h: x6 }5 t4 F( ^
12. toc; h+ Y, z/ n6 a
    
13. 
    \" Y) A( Z8 F4 r
14. s =
    0 C# K' A! ~* m; q, q
15.     6.126720600000000e+009
    2 ]6 p$ y# q: ?* w4 p0 G
16. 
    * k4 l/ _' Y- o) ?$ {2 ?/ I\" n5 z4 K
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
3 K( W+ P0 |, m+ m! _$ zm函数：

1. %file xyz.m
   8 {/ w7 D1 m* k( A
2. function c=xyz(x,y,z). r: \+ a/ s5 M4 k7 r! J1 F1 {
   
3. c=x+y+z;7 A. k% U! Q! d! Z, K
   
4. end

复制代码

测试代码：

1. tic;
   - o3 C; W# u\" ^7 d! A, ]7 R5 O
2. s=0;9 ?8 C2 k' m( r: Z6 F
   
3. for x=0:1000' C7 X. v: B' u7 j% ^; d
   
4.   for y=0:100
   ' G( }5 h& F' ~! C+ l/ f9 Z
5.     for z=0:100' a# F+ n* z6 _# R
   
6.       s=s+xyz(x,y,z);: U$ |3 H1 N1 k$ w\" S
   
7.     end( P! Q& k. ~3 z7 G
   
8.   end# @/ z$ m0 M1 }# b
   
9. end
   6 @1 ~( I/ [6 m\" G4 b/ x
10. s
    # c6 o, U% N6 w- m4 w
11. toc
    0 X  f# N+ m# j  V3 e
12. 
    0 j! z\" q/ l5 N\" Q0 p2 Q3 r
13. s =& K2 g, d% K1 W# j* g
    
14.     6.126720600000000e+0096 j) s3 \, k& g- W2 k( I9 @
    
15. 8 \: R6 A( w, |7 b
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

4 O0 I/ m4 ~; g; U+ k( f7 `4 \) B9 iForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;. H: p' C& I' h4 e$ L
   
2. mvar:, ]& ]8 b7 H5 }\" Y
   
3. t=sys::clock(),
   # p\" R% Z: S) h+ G
4. s=0,
   , N8 u: D1 ~4 @, N: c
5. x=0, while{x<=1000,
   ' T\" U! N5 h; g* l( M
6.   y=0, while{y<=100,
   ( {6 q3 J! E0 Q/ }/ ?7 H7 O
7.     z=0, while{z<=100,! S/ A+ z6 _- U& q& l9 l
   
8.       s=s+f(x,y,z),
   \" b, h( r: R# x\" _, {, T  C
9.       z++
   ( s1 _- z1 h: r; b! A$ B\" i) A, W
10.     },1 A2 d- p0 n- `5 W
    
11.     y++5 Z& a& z' Z1 I# [: j
    
12.   },( G/ H( i( U# X3 Z8 S
    
13.   x++' s2 b. x7 l, h& W6 Z# c
    
14. },
    # t2 g. z. A( l# o/ h
15. s;
    & j/ Y8 g3 F6 J' ]. x2 x2 C# q: g. b
16. [sys::clock()-t]/1000;

复制代码

结果：
1 R- g" p9 U) x$ c6126720600.
1.109  秒
6 I8 x5 c) P; ]* i" O; [
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

$ F6 K/ Y! B" z



* I- Z6 ~8 _! i- r8 {

3 B2 O7 y# U& W+ i$ }# |4 [0 J0 Y" K# K( ~
3 m5 Z' z0 A/ g# O& z


' `1 _: X2 j! b2 k
* Q7 k* v( L, k7 j




) s" I; h* A5 O; W  K, l