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

forcal

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

: S  y/ T* Y/ H9 z6 M＝＝＝＝＝＝＝＝＝＝＝＝
2 T$ q; h0 Q3 l# b' a  P
" f$ j, U  M4 K! G% f; XMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;: ]5 s1 U) m8 w( V% t3 Q3 H
   
2. tic;
   . L/ K# c9 r. }  L9 O9 f
3. s=0;. N* S, X: J* ]0 T3 N2 |
   
4. for x=0:1000
   9 ~  i! N# \  c6 y( f% L7 w
5.   for y=0:100; g3 h; }4 A% z9 F6 w2 ?- A
   
6.     for z=0:100
   + s! O9 |$ N0 L: D) L4 @5 x! l
7.       s=s+f(x,y,z);
   & E) `: H/ p' C% ^+ D\" q
8.     end8 b$ \2 ]% v2 C3 o( e2 d
   
9.   end' n2 W6 Y- U6 o0 T; J
   
10. end
    ; l6 s3 F3 i8 p$ |( y/ z/ R1 r
11. s
    / l& w. _+ ]4 D( Y2 }& Y4 p
12. toc
    / D3 w5 D- w2 j5 t. V- l7 J
13. 
    2 G# y\" v& p( G# h2 a
14. s =
    2 }$ v6 s* U+ e; Y2 v- B
15.     6.126720600000000e+009
    ' `' n7 Y* y) I2 P- ~! o
16. 
    3 J) l1 ^2 N) M- v4 D$ |2 k9 Y  u
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ( {- \5 L8 Z\" B, A4 v, I
2. function c=xyz(x,y,z)
   9 t9 C, k6 p3 W5 Z1 Y9 i
3. c=x+y+z;7 C+ b% _/ W+ f  V
   
4. end

复制代码

测试代码：

1. tic;
   ) e. L0 W. I+ O6 K. j: @
2. s=0;
   & f! Y% p! o( c  {\" o
3. for x=0:10001 F% E* b- r- u3 c4 i$ R; H- U& B
   
4.   for y=0:100& P9 t7 L3 _6 U/ X' q1 r
   
5.     for z=0:100$ {# z6 Q: W. q
   
6.       s=s+xyz(x,y,z);/ V! W  n: X+ Q- A! ]
   
7.     end& r) L) t2 i2 g) v- y1 e
   
8.   end
   ' ^9 s+ }' V! |5 o
9. end
   2 V3 S) U6 [* @% r' ]# h
10. s1 G- |* e( C  y8 r% S! U; N
    
11. toc
    \" A! `% d6 ^\" a* g) t
12. 
    0 D\" O: x( m1 h0 M5 w
13. s =
    1 m4 i( B: F& p6 @' d6 }
14.     6.126720600000000e+0099 l, U6 f' ?# r\" \1 h3 ~2 I: m5 f
    
15. 
    9 N2 x; `6 m; V. d/ G
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

  F# W/ F8 }  Y% ], lForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;. ^* T$ E* H+ Y% s
   
2. mvar:
   2 H* |1 L# m+ b\" H
3. t=sys::clock(),' `% Q6 ?\" R5 ^
   
4. s=0,0 f0 {( V& U' O0 W\" ?: D: `
   
5. x=0, while{x<=1000,: _1 d2 d( w, t* e
   
6.   y=0, while{y<=100,& T8 f  T6 L5 D( ?6 b9 k
   
7.     z=0, while{z<=100,' h) T# \' \/ _# v' A& B
   
8.       s=s+f(x,y,z),
   , V: j) e! T2 h7 A$ ~
9.       z++
   ! {' w. h+ H+ q5 n
10.     },8 d6 m& [7 |. i1 \/ M
    
11.     y++
    $ D* M  D& Q2 Z8 g: }+ [
12.   },' b. {( y$ _: Y  z: w; F+ L
    
13.   x++7 \) r. n% H\" x\" ]. ^4 E+ C
    
14. },
    0 u4 B$ m% ?! Y( {$ \
15. s;. V7 ?0 U; @9 _: }2 b
    
16. [sys::clock()-t]/1000;

复制代码

结果：
2 |- U2 B' c: d$ Z6 y3 i6126720600.
1 N/ F" W8 U0 M, Q1.109  秒

6 {7 W, n' r% k/ n5 Z$ z& u0 R二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
. d+ m+ m. x  n  ~# k2 ?- y

hxftntA

% J+ b/ Z0 u8 k# ^

8 O4 S* v, t* |. g2 Q
: l0 T. i4 g0 m2 T
# \% A! {% R9 J; Y+ B" ]
! B& N" k1 G! H' `) F5 M


/ V" R: w2 T8 \4 m/ E* F* P5 O) a
; M9 r- q/ B& ~) k* j
1 o5 E, ]$ Q2 ~% K+ B7 F. b


3 }/ I1 `0 k5 g. b7 O* l
2 ?) S: T! c" D1 h7 b; C- w1 X

/ ~' J" Q& @6 G; `/ S4 d/ \6 ~