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

forcal

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

- K  V+ }# }2 E  p. K) H+ V＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;$ i\" W5 [$ X/ A. Z* k
   
2. tic;: a( y9 p: d- [2 y1 U% w$ v
   
3. s=0;
   - v& E& M3 y! K
4. for x=0:1000
   5 \8 t% X2 n\" z) e0 |5 R4 Y
5.   for y=0:1001 O. e' T4 m- |4 g: S8 N. p
   
6.     for z=0:100
   % h# A, T* B9 z  k0 P& j
7.       s=s+f(x,y,z);, B+ S9 a' _' w3 H0 u7 x* a
   
8.     end4 A/ `! L9 x0 _, n% J
   
9.   end8 {! \7 i- M5 @; V4 k
   
10. end+ _$ \7 R* ?\" e1 b; W
    
11. s
    + e) }: Z, Q+ ^9 d/ `5 j1 V
12. toc; p( X1 h; g\" N) s! R\" }; Y' }/ r( p
    
13. 5 n6 V, e, D( `7 m
    
14. s =
    ; q* q5 b/ B) C
15.     6.126720600000000e+009
    ( }4 _: g9 Y* s# f+ `
16. 7 {  F: H7 _0 w/ x* K+ t
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
% S6 f3 d% {5 s1 u& [5 L: em函数：

1. %file xyz.m
   5 X8 \9 l7 F5 }7 I( F$ w4 i
2. function c=xyz(x,y,z)- y! W( f! r6 y7 p3 Q5 a1 w0 `( c
   
3. c=x+y+z;6 M/ b, y* h) i) f8 b
   
4. end

复制代码

测试代码：

1. tic;& U4 Z. d/ F+ @
   
2. s=0;
   + I* C* l9 D6 f4 E  z+ b; K
3. for x=0:1000
   6 J4 {4 y  x+ F* I& {
4.   for y=0:1005 c! ?1 s# {6 P: |
   
5.     for z=0:100
   0 {0 k1 |+ x6 p5 C+ y3 b\" ~
6.       s=s+xyz(x,y,z);) n4 b* a! f6 T0 y: i
   
7.     end
   9 A- d' k  L( R( ~. w2 z5 L\" o1 A
8.   end6 p- c; y$ l. `/ S: ~- F
   
9. end+ `, V8 q( m) Z$ p  f0 `5 H' y
   
10. s& H1 i0 k3 R7 k4 T
    
11. toc) Z& X3 B, N$ R) g0 F0 }, E
    
12. 
    8 ]! V, G5 m6 e6 F; y
13. s =
    ! p. K% h3 j  l' L  D4 \
14.     6.126720600000000e+009( s\" ]* {# q! L/ [- Q9 R
    
15. 
    8 `$ T) w+ _% b3 x4 l# A
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

! f; N+ M0 K$ t0 }9 _% H+ ZForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;6 b) w2 s$ r$ m
   
2. mvar:
   4 v' Q/ \5 ?# s' r: \
3. t=sys::clock(),& g0 K4 c7 J  L2 O* ]& F/ I) x
   
4. s=0,
   ( P+ G6 o2 ?0 k- i
5. x=0, while{x<=1000,
   % W, _. i+ G3 t7 q9 y! L6 |
6.   y=0, while{y<=100,& C- [+ \6 a( a% L6 A\" l
   
7.     z=0, while{z<=100,
   + ~7 @% M, {: Y* ^
8.       s=s+f(x,y,z),
   # p/ R1 }\" ^2 f  h0 d, R2 \+ n  C
9.       z++  X7 O( E* r$ b$ U8 G  K
   
10.     },
    ! P+ o3 t, f# |6 m: @8 w
11.     y++
    9 m- l4 v3 Q0 J5 V. `( _7 L5 m
12.   },% W: Q5 m5 X) e3 s+ N7 q
    
13.   x++! ?  q& X, l% L0 m) O
    
14. },3 ~, M8 g4 s1 e; T* |
    
15. s;/ S+ ]0 v' C- O7 b5 B7 E% \
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
# z9 F) v2 N4 A- n2 V2 S$ {6 j1.109  秒
/ Z  H& _; b% R3 @& l' d
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

3 X" f- U/ H9 o; ^/ J% Y8 Y, w. C
/ @5 J  o* _1 H: j& b' N

+ n( j, ~  b: w, }. v& E- S  p5 E



( ?' ^! x, L' r/ D& G


0 K3 M! n+ f% N8 D4 R0 d  ]
) o1 ?1 U& p. T

" D! V+ X& ?& W" C
2 [  j7 s* C% O3 e7 W! n- X3 H

9 `: u* u4 D+ v  V: `