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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
% H8 `- G) E' |' e
＝＝＝＝＝＝＝＝＝＝＝＝

4 V3 ?3 \( U9 [9 p* j( Y+ \4 P! H7 l; TMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ) k& [1 @# u* }5 T& Q1 m' y
2. tic;& E( ?1 P' H8 ?1 U/ F' H, g' n
   
3. s=0;
   9 Y4 N) K( S/ {1 {
4. for x=0:1000: k0 s* L/ H4 [2 f3 u
   
5.   for y=0:100
   - c1 f$ J1 o. P; P$ u( z
6.     for z=0:1007 D& ~- E& v3 |2 ?7 v
   
7.       s=s+f(x,y,z);
   ( h' C4 w! e. M& U8 Y9 S
8.     end
   0 Z1 q4 T# m3 T\" B9 j/ q
9.   end5 t4 J3 R. J7 l0 H! L; Q& Z: q7 ?
   
10. end! h- W4 P* p+ ?9 }. `1 F\" w
    
11. s
    / q# P# i* {8 u\" a, }
12. toc: u2 h; b/ D. i. N9 O
    
13. 
    2 f2 p& `9 x8 ]% C. M& Y- [
14. s =
    5 }& H$ s- g/ p\" @  W
15.     6.126720600000000e+009
    $ O  U' N. E. v( L\" p  {5 N
16. 6 g\" p, \( ~. z. |# _3 M! E# S7 Y5 v
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
5 @- z9 ?( ^, ~6 `: n* rm函数：

1. %file xyz.m
   2 D/ [6 \7 L. t% c
2. function c=xyz(x,y,z)! B( M6 P$ q$ v9 O5 z0 F! y& b3 l. ]1 v
   
3. c=x+y+z;' n) ?* V# g! Z1 n
   
4. end

复制代码

测试代码：

1. tic;
   + d( i% i6 c! w  y9 ]
2. s=0;
   3 ^6 y6 f: d9 e! V: [3 d( Y- k: P  Z
3. for x=0:10001 S. `+ Y# ~6 T- v9 E, r$ P
   
4.   for y=0:100
   3 p( p* v9 a1 `% r* N  K) r
5.     for z=0:100/ C# M1 M1 E/ y6 o
   
6.       s=s+xyz(x,y,z);- U7 e# F- ]* _& m
   
7.     end
   5 F0 {- l; s% B! i  P! J6 u  r
8.   end
   + U8 C4 W. M9 p( f8 t  f, t/ E\" I
9. end
   # q) |3 S' F9 N7 S2 k8 t
10. s
    1 I* c8 f: ~' E, y( c& P* }' W
11. toc
    . _& E5 Z( b9 U2 j) k
12. 
    6 O- Q& k: ~( i( E1 E
13. s =
    $ w0 A\" i* y/ X, ]+ `! i- h
14.     6.126720600000000e+009
    ) R2 `* P; N+ t
15. # W. u; H9 [$ o% j; Y2 N
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
: ~1 N: K  G7 m
9 k/ l$ z' p' T0 @9 W* Z% ?Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;& w* e  u  z* ~0 S. I
   
2. mvar:
   . b( h  U0 ^/ n- u. g
3. t=sys::clock(),
   - d! e5 s$ M8 E) M
4. s=0,
   * w3 Z2 B. k4 F$ g2 E4 X9 O6 f6 D8 C
5. x=0, while{x<=1000,+ z3 q- S: S9 |5 K
   
6.   y=0, while{y<=100,
   * u9 Q* z  Z, ?3 I
7.     z=0, while{z<=100,2 v2 o% G$ ]% Q  d
   
8.       s=s+f(x,y,z),/ B& N/ B# X/ w- T- s
   
9.       z++- x2 i6 `$ O3 V* }$ k
   
10.     },1 X; R) O) Q\" [. ]! \\" l
    
11.     y++
    # Z; n0 E2 P0 V2 o% F1 s
12.   },
    ( p  v: x0 l! Z0 `% {- B5 q2 q
13.   x++
    9 ~7 d6 k7 _, V% x
14. },6 _- U  U4 I' s; M\" ~; m' ]3 Y' O
    
15. s;
    , k( n\" X! Z% E+ Y; j! b
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
0 I9 H/ r( g/ \1.109  秒
' m8 {4 b6 E! F' H( s/ j
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
* ?' w' H1 Y7 T. n0 k

hxftntA

' Z) G0 u. D  _% H% J

$ r, u4 h4 Q: F
, ~' n" S, ~0 Q% n# }/ _: S
3 H" ~- C' B4 W1 l1 o
  M  Y! J1 E# ~! X1 r
6 j% j6 G5 n5 [% O2 o4 L( C3 {

+ q7 f5 }) O' w5 b1 f5 h
3 L4 g# H# {3 q- P4 m9 y% O: S
1 c9 W) l  S' `$ J. H

5 T, ]5 m# T6 P$ R& u% |9 `


- a9 _; k5 X) e+ t# j
5 U. E; Y( e8 P5 ?# _6 {+ A