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

forcal

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

5 _1 p" P8 J. L6 u: l% I! E! M" m＝＝＝＝＝＝＝＝＝＝＝＝

0 M3 }& e/ ~5 F% x: ^% [- t5 U. z1 oMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;3 ?/ g; W8 t4 f1 d- \
   
2. tic;
   ! I3 Y* C8 ~  @
3. s=0;: i6 L: ~9 c6 y8 Y, D6 ?
   
4. for x=0:1000; b4 a# P2 s& F9 \
   
5.   for y=0:100
   1 o8 M* x7 E3 ]& E3 l
6.     for z=0:100
   9 o! `& T7 u! ?; ^4 g/ t
7.       s=s+f(x,y,z);
   - B1 m\" k$ ]  z+ {5 m9 o' W
8.     end: J& b/ o% t9 d# [; R  s1 k
   
9.   end* d# n8 m( q4 e: P0 Y/ C
   
10. end
    6 i; j6 A' C: N# m# f  N
11. s  @5 t, ]7 |. b$ S& F4 y
    
12. toc
    ' h- v8 H; f\" U) i
13. 
    % M& q9 D: h, y% u% L. J
14. s =, N: s4 }/ \% N# g
    
15.     6.126720600000000e+009
    . h, X3 W9 I# |' W; K1 G
16. 
    3 O( e4 K( j) ?7 r1 G; N
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
$ T6 v/ s8 p" x/ s% e$ Q4 s9 Hm函数：

1. %file xyz.m
   ! E\" ^\" }: H4 N$ v
2. function c=xyz(x,y,z)7 R\" F$ z. B5 J* j+ h\" W
   
3. c=x+y+z;
   : R, C7 [- E2 y& g; w/ Z. ?
4. end

复制代码

测试代码：

1. tic;
   $ h7 a( H: `2 J5 G# v
2. s=0;# u8 u! V2 X( u' K+ A+ f* o
   
3. for x=0:10001 e- {4 L  q- u
   
4.   for y=0:1006 v+ l: J6 Z% W2 P( @
   
5.     for z=0:100, ]6 }' T, b* Q% |2 Q
   
6.       s=s+xyz(x,y,z);
   4 F# w$ j8 A' ]$ b8 ~8 _1 k) v
7.     end
   2 x1 H+ t& z; ^7 A
8.   end* h9 n\" D. G; R4 V% h% h% e, V7 ]
   
9. end# F( z; o# Q- r5 }- c
   
10. s
    4 w7 D5 J' G# @9 O9 z; L- h\" _3 W
11. toc
    9 v+ s  ~6 F$ M0 C+ n
12. 
    $ p  ]& ^* j1 L& a6 k
13. s =
    * _- M/ w; Q, g* \# \6 G* K
14.     6.126720600000000e+009
    . b2 E  I. ^9 v
15. $ [7 @. T5 @. [* ?0 R
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
/ @. ]( z/ U0 y) C+ s) z7 A
# d. ?3 |4 L& l5 j- N2 \Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;- ?4 p8 R5 o\" A; i
   
2. mvar:
   ! |: c8 i4 B  v- ?  l
3. t=sys::clock(),6 L; Y& X' A. F& X\" O5 q0 l: t
   
4. s=0,4 \, V3 p) L* W3 ^
   
5. x=0, while{x<=1000,, D! Q. x/ S3 \
   
6.   y=0, while{y<=100,
   4 u* w  O& M' A\" }5 `! b9 l
7.     z=0, while{z<=100,  v5 ?9 R  i! g\" {: }5 X
   
8.       s=s+f(x,y,z),
     I  ?5 W' S) N5 g5 g' I
9.       z++
   $ N0 l% V0 s  k$ x' p4 t6 J
10.     },
    2 W\" q/ t: }' Z$ v6 [  _3 y1 r
11.     y++& E- W5 L% K. {. Z
    
12.   },$ C; q. F! I; B, Q# c: ^
    
13.   x++, |  A% b8 D+ z( F3 ~) t
    
14. },. n' K: v) I1 y9 k& H& I6 J/ P
    
15. s;
    5 @& E+ O: ?- B# z7 `
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
7 t4 ~4 y) s  h# H/ v2 ?" y
( ^" u; W0 g0 m+ h/ Q二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

" `, L2 ]3 I' [8 o! F


. n) l, C( {, }: I0 g

& X( e2 B/ |/ l, m
6 h, R+ `0 Y, c! s( V: Y

" c" {" d6 W: |+ w* y8 |+ [1 |
2 N5 U: B5 d% z
4 ]1 g" ^. O; J. E) d6 M
3 v6 ^2 b  k4 f; W7 L: l+ R+ ^
( u% v  Z, \- {+ g+ O. V
) p- L: ]2 [: }3 X
: [  A# l& z/ \$ B( w6 s. Y