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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
# T9 w! \0 }* a# X) K; v9 B5 t
( H  h4 B9 }) w) X  g# zMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   6 m& _# J4 O+ N/ R\" \7 t
2. tic;
   # ?( i: E* O: c/ z' ?4 C' B0 k
3. s=0;
   \" z: d) \% A3 {8 s! R7 n! t% x
4. for x=0:1000
   : `3 S& y2 O( v$ I  E9 s! ?. D
5.   for y=0:100' S8 m8 V( I5 a/ X* ?7 ^( t
   
6.     for z=0:100
   ! O$ z) a' R6 T8 R! i( U5 j7 H# S\" l% D
7.       s=s+f(x,y,z);
   9 l! c! p) F6 @% x  R- V6 S4 |
8.     end  \1 A# |6 Q+ G9 M# d* G( ^# Y% G
   
9.   end3 Y; z& E/ O) Y: L* t
   
10. end4 W, {' _% U, f' y$ s7 n+ U
    
11. s2 H1 V- b' ~9 d7 n+ V+ }
    
12. toc
    ' o0 j) Q% e- V6 w
13. 
    ! v3 D5 b2 ?! x
14. s =) q' a6 H. O! l: b! _+ w& W+ [# u
    
15.     6.126720600000000e+009
    $ C6 }+ U/ a: H
16. 
    7 X3 f! p1 V7 y  N- k
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
# V% ]7 [( h+ m8 P/ t* k  hm函数：

1. %file xyz.m
   3 [: s' c9 o0 e# M6 _) _- X
2. function c=xyz(x,y,z)
   ) q9 Y: O) P2 ^& U\" f% ?4 a. o
3. c=x+y+z;4 F) L) D2 s8 i( Y% x: f) G
   
4. end

复制代码

测试代码：

1. tic;
   ; @; ~% D6 |6 U
2. s=0;
   8 Q2 x5 e  d' W5 x
3. for x=0:1000$ q9 S9 C2 ?; C% Y7 W4 z, G' @\" o
   
4.   for y=0:100
   $ I1 l8 O, m/ V3 m
5.     for z=0:1008 t0 V, p\" F; N# A; _
   
6.       s=s+xyz(x,y,z);( H8 X7 P\" Y, R3 |- E3 x. L
   
7.     end8 r2 C3 b8 Z; D8 P; Y/ |* Q8 f8 p% q
   
8.   end
   5 E+ B$ [! s4 g0 C, T\" m
9. end
   , X\" `* p\" O0 T+ C6 ^: L5 w  r
10. s
    ! x* [/ X* ]  l( s
11. toc2 j: Z% t, A+ J+ _1 Y5 J
    
12. 
    # W- H1 O* B) N3 I
13. s =
    4 A* f$ E' d0 K9 {
14.     6.126720600000000e+009
    - `1 b  r/ S* e, v
15. 
    4 ^$ u\" f9 ]  R/ S6 k+ O& l
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

) R0 f7 C7 C- V- [  s6 d1 ZForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ( i0 h7 N% p$ |8 W) {; N6 k1 p
2. mvar:; l9 p* _- F+ J$ _0 P6 h
   
3. t=sys::clock(),+ {8 B% @  i) g) x  b: E
   
4. s=0,
   : N/ l! h$ h/ ]; w. n% [
5. x=0, while{x<=1000,5 k& y' l/ d, r# Q/ n
   
6.   y=0, while{y<=100,( `! o( K* c5 M% I; b
   
7.     z=0, while{z<=100,
   ! d; s/ n- G* h: K
8.       s=s+f(x,y,z),
   7 D0 d) k# j+ D* q' w- c
9.       z++
   9 m8 f6 s5 P) }+ A$ B
10.     },5 i2 j\" a- m% C% z
    
11.     y++% ]' K2 \9 Q9 s* X: n; n& d
    
12.   },
    3 A/ N) @( s, W. o+ N6 {
13.   x++
    * u9 S0 L$ ~5 N0 U
14. },
    0 ~+ _3 p) _5 n\" i+ u  a* W$ W' q
15. s;
    \" `2 h! Z3 Z% {( j0 d( p/ Q\" \
16. [sys::clock()-t]/1000;

复制代码

结果：
5 {/ E- q* h( w7 l, K0 E6126720600.
* W3 w5 ~7 u9 P$ C0 {2 b1.109  秒

1 j3 H9 R& M# y$ S, h二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
9 @& O# }- D, r5 v

hxftntA

$ X- A; R9 P: A* m: \. P$ E* c
3 \$ H+ ?% U' @
; ~# D" L. p6 t$ R' S
+ B! |8 j. Y1 J) o4 s



1 b1 s1 N: K0 n
9 I& Y& a7 Z( j) Y' A

  I# b' x7 s9 {# A' `. t$ k3 A
$ B) @0 p7 R6 o9 V- ~
% a- z& m( d* h0 d! C6 i, u
8 k" E2 [! {" R: a
+ s3 I2 {/ w/ c5 w  J1 L! E