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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  l$ K8 F% P! `, F, f$ M4 b. V
! K0 R3 D. Z1 V! b( \, `: `＝＝＝＝＝＝＝＝＝＝＝＝

3 c4 y5 `( }7 u/ f* D; y- l; }Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   0 F: c2 x0 @. Y4 j9 j' U
2. tic;
   $ n6 E( {/ q. v' x' m
3. s=0;
   # [) c  Z; ]; ^
4. for x=0:10007 ]. [2 d# C% O
   
5.   for y=0:100; b; X3 W/ }9 N% ^
   
6.     for z=0:1000 e1 E! d% o' O: Q: I0 ~8 n
   
7.       s=s+f(x,y,z);
   8 B! W4 |/ m! a* v0 G
8.     end: L- L' ]! C- ]/ T/ f, s
   
9.   end
   ! _+ U' P; L6 _0 `: n, L
10. end
    & j/ a' @9 h; m) n' R
11. s
      K% T# J* \! m4 L- C
12. toc
    % A. h\" }4 Q) w0 p4 |6 q
13. 
    ! d# k$ ?5 X. y1 a/ }1 y- S) X
14. s =- x  v\" V. n$ d1 N: ?4 P
    
15.     6.126720600000000e+009
    : T) P# f( Y& g# {1 T* ~0 \0 T1 ]
16. + L; r6 F, k, |: e- F1 i/ }! J
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
7 V: O  G7 @6 t, ~* e) om函数：

1. %file xyz.m
   ! X! v4 j/ g\" ]- f. e# r
2. function c=xyz(x,y,z)$ u$ y# D& \' \6 ]$ h+ W' ~
   
3. c=x+y+z;* |' ]- B: g, b
   
4. end

复制代码

测试代码：

1. tic;
   3 c9 C, F9 N) b7 e* k1 @9 M4 W6 `7 n( @
2. s=0;
   / H( K: d; T) N/ R1 Y* ^; |9 p
3. for x=0:1000
   / z3 D# K. z. h* c7 @- G7 A
4.   for y=0:1004 j& a\" ]! ]7 G  Y% U7 U
   
5.     for z=0:100
   . \' |+ O\" J1 M
6.       s=s+xyz(x,y,z);' A: G4 f) ^$ J. y9 o9 R\" |
   
7.     end- R2 S! d\" a5 v7 P5 S
   
8.   end
   , Q3 w* G% s( n: T
9. end
   ( t4 h( v9 w, `
10. s
    7 B/ c+ W) K, I, B8 u2 N
11. toc
    ) k3 A1 r) b: [$ f$ h1 S9 J5 L4 i) T
12. , S6 P8 [6 P- R
    
13. s =
    ( m9 a, Z& b! D4 }- j
14.     6.126720600000000e+009+ P4 |; L' D5 [* F$ J; n# [
    
15. 
    ' x% A; \4 z9 d& s% }- g
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

' s' o$ ?( D: C  N( H) d* [Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   8 D& F& I' |. X
2. mvar:
   & l% R: l, s. x) x# c
3. t=sys::clock(),, e/ p) i; T' L4 V* e\" Z\" y, D( K
   
4. s=0,
   4 y0 H; M2 y& m( @$ D
5. x=0, while{x<=1000,/ Y9 d% C( x0 ?* X
   
6.   y=0, while{y<=100,# K+ u2 l9 V+ p3 h2 h  F# q4 W5 m
   
7.     z=0, while{z<=100,
   8 Y( V# g, t' P( }1 @
8.       s=s+f(x,y,z),1 D* J3 A2 b. y, u, j- N
   
9.       z++: o1 h* {# [; k  {; X- }\" U
   
10.     },
    # b& V& o8 K0 o
11.     y++: j4 b4 @2 |\" e, G
    
12.   },3 m. ]! R2 W( p9 b' |0 O
    
13.   x++- u8 S# v% C, W& [% ?
    
14. },, Y/ v1 @/ t5 G7 x  G\" l8 X# g. Z2 w
    
15. s;3 g5 P  T. ]( `
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

. Z* I: L1 L" {' m  m' G二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
! l" \; _7 m2 q3 Q* H7 B( Z

hxftntA

  N, M$ q; n( }$ Y

1 n3 _; |+ w  g8 {


6 H2 V2 J4 q. [0 Y& s% G$ u

  L* q# q* y! p3 f

4 s4 `1 ]% w& a6 R) u: K0 L, Z' R

, M, [4 Y* q: o; x# ^