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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
3 S4 i8 n! v8 f- M, U) E, \% {3 A
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   & Z8 {' c% @6 d- ~% a
2. tic;. z. x  q  i4 i; G  }1 f
   
3. s=0;
     w& f1 s) R' N9 E
4. for x=0:1000! w0 c! Y! X5 D3 v% N0 s1 O
   
5.   for y=0:100
   , n; q0 |  j+ K+ m\" l' I$ `
6.     for z=0:100  |# m# [/ |# z9 G. O
   
7.       s=s+f(x,y,z);
   * H! I) P3 [  K1 j! m, a
8.     end
   + D( C, o/ ?, I4 d8 f  h
9.   end' k& ^\" r\" T5 \( L
   
10. end
    ) F, k, t: D\" `\" F* _/ u/ h
11. s! b4 B, L5 H& Z3 x+ V1 Y
    
12. toc
    + w5 |& t5 e) o  m$ M
13. 
    % n! e. W& R* s1 ?
14. s =
    # B5 g  G3 N& q* h
15.     6.126720600000000e+009
    0 Z& p  A% \$ D! _
16. 
    ! `- H6 W9 C0 x5 u! }5 \) L9 m/ @. v
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
m函数：

1. %file xyz.m
   ( p8 d% b7 d; Q& f5 q# v+ Z
2. function c=xyz(x,y,z)1 j+ I7 _1 Y$ J
   
3. c=x+y+z;7 A5 `9 A& M6 W. J$ |% b
   
4. end

复制代码

测试代码：

1. tic;
   + i: E$ e8 Y1 @: z- K9 F
2. s=0;
   * O* A) [/ {3 c0 d
3. for x=0:1000
   6 z- N& \6 y* Z3 X1 _
4.   for y=0:1002 I0 K$ b6 g% W( s/ E. V2 c
   
5.     for z=0:100
   . I, }; a! R0 `
6.       s=s+xyz(x,y,z);
   5 N+ c% c+ w) V1 t6 i  |
7.     end( R2 N7 Q! [8 u% y
   
8.   end
   ' x8 B' S5 ^5 Q3 O3 ]2 U: G
9. end& r# ~) R3 m$ w; f- S
   
10. s
    , S2 r2 o9 }8 V! `/ Q! C
11. toc- X' h% y4 ]- [$ f) \2 T' t6 c
    
12. ' T+ Q, n4 N. d. ^
    
13. s =
    ! P/ @) B\" F3 Q, @# i\" }; ?' d) {% s
14.     6.126720600000000e+009' ^' a) |5 Y4 m- o- n0 C. @
    
15. 
    \" l$ [1 p; D( J% Z
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   % S, |3 f/ F5 z% @: c( U
2. mvar:; f9 W! g8 Q9 @8 Y8 Z
   
3. t=sys::clock(),1 n; U8 |\" s4 j& n
   
4. s=0,
   8 j0 C' L+ \$ ]& @\" m
5. x=0, while{x<=1000,6 v7 R' g5 H% W/ ~7 w) k  ^) S8 w
   
6.   y=0, while{y<=100,9 L7 `: W7 m0 G: x; a
   
7.     z=0, while{z<=100,  c) u& f+ k2 a  J1 P
   
8.       s=s+f(x,y,z),
   ; F* Y. b  u; I\" H9 K) X6 a9 X
9.       z++& ~3 D6 n5 h$ K8 o' {
   
10.     },/ N( e, x/ N. v  I  d* R- ~
    
11.     y++% b5 g8 o8 x4 M4 r1 V6 O. K) s
    
12.   },
    # J0 X- _; L' W9 ^
13.   x++: N1 A2 T9 D+ ^. ^3 ]+ d
    
14. },
    % q9 R, f* i6 W6 S
15. s;
    ' T$ u0 k: W) N
16. [sys::clock()-t]/1000;

复制代码

结果：
' W+ m( ^& V* G* I8 T6126720600.
6 ~6 R8 P1 G8 @% i6 l6 j; {1.109  秒
: R( W' `1 O5 v
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

$ C2 a) B' g; G% J; Y. G



3 f& I: J3 _% K8 L$ q! @% e- |3 [

+ a! e! k* \' {




2 q. P% g! e7 w
  s5 h- o) S* i( b; T% l' I


/ K' Q6 {' F4 Q) O- S2 Z3 x