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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
' [: d2 i. o7 D& e. {: I# f+ I
# j' l# Q/ r! P" s1 [＝＝＝＝＝＝＝＝＝＝＝＝
" _5 R  B* Y4 e. j
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   2 E  p- y3 y+ ^1 u+ W6 ^# G
2. tic;
   . C( F8 L: x1 H$ O- o. q. _
3. s=0;
   ( q% W4 r: w$ o' e\" t6 b  j2 n
4. for x=0:10004 W6 r; {% C5 ^' i  ]8 V+ H5 P
   
5.   for y=0:100
   , V0 R\" |$ V\" }+ U  c% `
6.     for z=0:100
   8 f, ^# s, ^8 q8 n' u$ R, M
7.       s=s+f(x,y,z);+ L2 e  o6 @. L5 g6 k% N1 I\" g
   
8.     end
   6 f# |4 N& B% Z3 v% B% i
9.   end# k- H/ a) S% c6 d* d
   
10. end
    \" K6 G' \8 h6 v
11. s
    5 O5 ?  Q, W  e
12. toc6 \. W( M5 q- Z4 P* C* @
    
13. 
    ) X) x# b4 ~1 j9 o5 C8 z/ S, _8 [+ l
14. s =3 ^% y* D$ M# p; N  M- U
    
15.     6.126720600000000e+009) }5 B5 y5 a& p( ]
    
16. ( i7 I8 M* s3 R. H, c( F
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   # E) p# h) m& W9 K
2. function c=xyz(x,y,z)
   4 L8 [# I- B\" w' X* u
3. c=x+y+z;9 V- R  |; R2 a
   
4. end

复制代码

测试代码：

1. tic;
   . f! W% m/ e0 u4 C& K1 `# Y/ R; W
2. s=0;
   8 I2 M0 a( K$ y/ [6 {3 A  E
3. for x=0:1000
   # D  I% y- i  W# m
4.   for y=0:100. z+ Q: U9 r( }1 @
   
5.     for z=0:1000 [2 m/ G( h  [2 F3 z& i
   
6.       s=s+xyz(x,y,z);0 o0 @$ X- r$ |9 ?. M
   
7.     end
   , m) A+ c+ D4 Y% b+ I& S8 F
8.   end
   - q& U7 J# j6 r7 x2 U
9. end
   * T8 R* ]: C  U* I* X\" s3 A
10. s
    - z  }* g5 r, E$ {4 N
11. toc
    + [; W2 ]2 `7 U2 i6 {4 |
12. 
    ) C2 K; _: c9 u2 |7 V
13. s =
    / U/ W, Q; a0 {( r' S. z9 c' I
14.     6.126720600000000e+009) C# g5 ?: m1 F8 J$ q
    
15. 
    $ R$ h6 s7 ]' n
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

# V0 J1 c* r  bForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & s, n/ L, I( |. `' x
2. mvar:
   4 I/ v  k& O6 u  O& p
3. t=sys::clock(),
   2 d% F8 y0 F1 _# A2 a- J# x) q\" W
4. s=0,
   : z; l' ^1 q) X  A0 P7 ]# o
5. x=0, while{x<=1000,/ Z) Z' r7 }- p) _% V5 W
   
6.   y=0, while{y<=100,' D  X5 z\" [+ P2 T5 T4 ?
   
7.     z=0, while{z<=100,# g! D9 t! f( ~: z: Q6 s! a
   
8.       s=s+f(x,y,z),. z7 E( H$ m. o# P* C1 C) E; A
   
9.       z++\" j2 ~7 Y$ r* K1 x/ z6 n  M7 w
   
10.     },
    # I- E1 `! P3 w6 X0 b
11.     y++1 r& R. ~5 _+ t; V
    
12.   },
    % |$ w- R( z/ b, q& S1 R
13.   x++* t7 M2 B7 g! S  W( G
    
14. },
    5 w/ t2 F: J( t4 H9 S' y3 O
15. s;
    6 i\" z9 H- f; m$ F' ]
16. [sys::clock()-t]/1000;

复制代码

结果：
5 h+ r7 n4 n+ f+ B8 }6126720600.
# l& d% n, P/ j+ J4 R) {+ b1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

$ E3 W" M* `& m" T' `
. f: P; B9 ~& s8 u8 [/ s) N8 a- v
3 c* H* A/ o; O' E$ ?% @7 b

1 n- Q. h! n+ r2 ~
# L. z4 ^& Z6 u0 Z: O


1 _' d1 {- `" r& [, @* k# S1 @, n! e1 \' o
; d! K5 b$ D( q- i$ T6 z9 x
4 M8 }% `, I" ~' g1 G. V  S
  W3 G, C4 \4 C* v
4 ^! @; u! }% k1 Q8 d

- l. B# _! @+ a5 Z% o

9 h5 [+ u7 F3 w3 A' B+ G: U5 @