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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

  ?0 b4 N4 t' i0 RMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;1 j% ?% d% y/ s
   
2. tic;7 |) e' \6 d% ~3 Q. r
   
3. s=0;: M- G3 L+ d; T( {; A  m. H
   
4. for x=0:1000
   3 Y4 ^4 M( b4 ^# u
5.   for y=0:100( Y\" p9 g! |% t# L, M- q
   
6.     for z=0:100
   3 u* g3 g) h; y\" s' I2 V* A/ u
7.       s=s+f(x,y,z);
   * b# w9 u$ h9 S) m1 z( u
8.     end+ Z& K' k& W5 m- G$ Y6 E0 j
   
9.   end! [* v5 P+ v& _7 I+ |* M
   
10. end7 ~* G; }\" V& S( s: S
    
11. s  X3 ]3 a& M4 K\" Y
    
12. toc
    4 H9 L\" @9 k0 R& y' o
13. 
    , N! F# ~: f) H* X\" s9 J+ q7 ~
14. s =: ?1 c5 P1 \8 x: B\" x( w4 J3 h8 ]/ C
    
15.     6.126720600000000e+009/ g7 \1 U& G( I+ r' \% q, S
    
16. 
      \0 e& m' Y- d$ v8 z9 G4 o& D7 ^
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
% j; g5 c5 l1 z" a4 S% nm函数：

1. %file xyz.m, f1 [% M; k\" a( ~4 m
   
2. function c=xyz(x,y,z)6 y\" O: i\" R/ f. t8 Y
   
3. c=x+y+z;( P) g) f0 j) U3 D
   
4. end

复制代码

测试代码：

1. tic;
   ' Q3 \3 }8 d6 Q! W9 t( j0 e
2. s=0;4 t% I; J+ _, h3 z! ?7 _
   
3. for x=0:1000
   ' Y\" h8 p: K6 E
4.   for y=0:100
   8 t3 U2 I5 O$ j
5.     for z=0:100
   ) R7 `+ a5 r* {- e4 V7 @4 x. t7 r, e
6.       s=s+xyz(x,y,z);. Y7 D. s+ l9 ^$ @; G6 Z# j
   
7.     end
   ) k, `8 `. ?3 A; K
8.   end/ |5 w( |) r% j. o4 m
   
9. end$ k8 u4 p' m! p( @2 w
   
10. s
    # Z: U! G' V# Q3 N\" |) |2 t
11. toc. E) }- q* @0 {7 K3 [- m
    
12. 3 i' p6 F& f$ S' m- e5 ^
    
13. s =- l\" M# V, F0 S$ v9 R6 k5 F/ Q8 t0 ]
    
14.     6.126720600000000e+009+ E0 d' {  Q' y1 H
    
15. 
    ! r/ H4 ~3 a  i
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
+ P2 |+ f$ G6 n$ F
1 K" i: |9 O* f$ ]Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;. h8 W8 t! g3 r' \3 }  k
   
2. mvar:
   $ L/ t$ W% @+ W; W
3. t=sys::clock(),
   % ]% [5 s1 [- J1 R
4. s=0,( D& E3 a1 R1 P6 ~* _. u, z+ S
   
5. x=0, while{x<=1000,
   $ u$ n% B, M* Z: }6 |\" K
6.   y=0, while{y<=100,
   ; j- @3 g7 x( i
7.     z=0, while{z<=100,; g- M( C8 s6 P; e
   
8.       s=s+f(x,y,z),
   $ m2 P# \8 `0 n( {# h, n! h2 P
9.       z++$ G0 A8 S# W6 {
   
10.     },
      |: ~- \0 e- W' ~$ v
11.     y++' O) s' [) v5 V( P8 {$ x4 U4 Z6 K
    
12.   },% J4 X& l4 S& N6 o6 z( p
    
13.   x++* O; q, G( R2 O9 `
    
14. },
    * w! a1 W5 I, |& |  p
15. s;
    . y4 M- q! f( V6 }/ \
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
3 Q5 x: K% G% I8 W; i
% @  b5 K. M- Y二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
" r" S: r. j# i. X" C4 r5 G

hxftntA

9 A; i4 t' H% ~, X9 @# c0 x- T7 O, z
3 H0 E1 a6 g& ^# h  K3 I
5 Z( |0 J9 A. F+ F

; @0 C! O9 f; t4 d& G/ f





' ]4 v1 ]) k3 M8 y% y
! o1 V" t4 P9 e9 r



5 |& H- D- c  n, t7 G& p% x3 Y