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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
" z5 ~! Y( m% q( k1 f
$ Y& |2 M3 E; a* D  m  f) g＝＝＝＝＝＝＝＝＝＝＝＝
: _3 {! H0 g4 T0 {4 R9 n# B6 x
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ; \6 w1 T7 ]$ @4 {
2. tic;' E  N6 ?* K3 l0 d$ }5 Q6 ?' P
   
3. s=0;
   8 a0 i9 l% T( c\" P9 P
4. for x=0:1000
   ) M, R& v5 e' G; H  X+ v4 @4 E
5.   for y=0:100
   - Q( j& `5 d& v  u
6.     for z=0:1001 u! d/ \3 P. J0 X# o
   
7.       s=s+f(x,y,z);
   - F7 h2 L# ~6 v* |) R! w0 k
8.     end
   & m& z+ U3 C7 H( x
9.   end: D$ c9 i! [: j' l- O& ~) g
   
10. end
    & b7 @- f& l, L- m) K3 M
11. s
    4 p# M6 }- S8 |
12. toc9 L* |8 x  d$ `! v( T
    
13. # q  j8 }4 E- f  ~\" v( m
    
14. s =
    4 O3 Q  z; @; C
15.     6.126720600000000e+009
    8 B6 o* T\" [( `& a# L' A
16. . F9 e' ^9 x2 ^; L
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
6 b- ^& E" F. T: xm函数：

1. %file xyz.m
   % y/ C8 U6 b2 l& m. _
2. function c=xyz(x,y,z)  X# ]2 L* I0 O9 p& m. Y
   
3. c=x+y+z;- g9 u. ?3 B1 [% b
   
4. end

复制代码

测试代码：

1. tic;6 l4 n/ }( j. L! ]& l
   
2. s=0;  i3 ^( x* O' n5 g9 F1 ]0 p
   
3. for x=0:10009 @* |0 f2 @' F; H4 {! t7 W& p4 c
   
4.   for y=0:1009 g' o# }$ w' P+ c+ E- r8 N
   
5.     for z=0:100! \9 g) j' E\" e
   
6.       s=s+xyz(x,y,z);
   % y# w- @; C0 s6 E4 p4 K& y* R
7.     end/ j5 y- t- F1 o! c* w1 ^' [
   
8.   end& k) E3 ^, N0 A( k
   
9. end; {$ p3 n. z: m) u' K- Z
   
10. s, P5 V$ i) k4 f/ x  n7 s
    
11. toc0 P) T# k4 ?$ H' G/ c( d9 z8 r8 b! c
    
12. ) Z; {- \9 N4 l! ]  Y\" t
    
13. s =
    ' ~) H% H* e2 A* T
14.     6.126720600000000e+009
      V. l! L  J9 d: Y4 d
15. 
    ; K6 w: o6 H4 W0 ~# k# W1 V
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;( Q% r3 |7 j; H& Y
   
2. mvar:
   2 H6 f# U- b8 \6 }5 T\" p: f5 w; K$ |
3. t=sys::clock(),
   , z8 r. a8 F# \( t; [1 q
4. s=0,
   \" ^! r\" G' Z( N) o$ g
5. x=0, while{x<=1000,
   7 D+ e, b% E$ ^$ P
6.   y=0, while{y<=100,8 W9 F9 e3 q0 ^' |( E# ~/ I
   
7.     z=0, while{z<=100,: I( i, p7 w4 x2 L\" l/ U
   
8.       s=s+f(x,y,z),0 ]. M, Y! D1 d+ i0 g& _
   
9.       z++  B7 p  P6 V' z2 }
   
10.     },
    5 K% `6 ~: Y* b; B\" ^4 C3 ]2 G
11.     y++( `2 Y4 q4 Y# ~5 H$ B8 `) e; b; _
    
12.   },! A8 }7 N3 Q  u1 \& f
    
13.   x++
    4 c6 s( f+ ^0 U
14. },6 O1 R$ @; \! P
    
15. s;. N; r' c- N0 O: A/ l: R8 ~
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
$ _, }- L' E" ?3 |1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

6 P" X( t$ c5 W0 c* s. u1 t4 M
* Q  `; d$ `. o3 E/ T6 }% u


/ Q  \  ^* w# K5 L
% o$ T. a! C4 v" O& c
+ a9 X$ S6 z# {" X

$ z/ N* `& W% O% W6 `( n

) \5 d( u6 i6 h) d, r2 {  o" _& m0 w9 D
1 @% t% @' [) v% y$ ~! O; h
2 K0 j$ x/ c0 d. C) W. U

, t( l' |5 C4 e) G1 U