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

forcal

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

! {& u! d0 I1 U& m0 {: _: K＝＝＝＝＝＝＝＝＝＝＝＝

/ m9 \4 j+ d. Y! BMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;1 H; F' j7 s8 W
   
2. tic;% c2 |* n\" A, R6 I
   
3. s=0;# i/ R; `* z# \1 u' x
   
4. for x=0:1000
   : i% W5 T1 I2 G
5.   for y=0:100
     i) w4 L, R9 n9 n, Q- _  G& ^
6.     for z=0:100: H, F. j& V\" L3 X7 g8 b
   
7.       s=s+f(x,y,z);
   \" O2 x# ]! o  e! R' Q\" l* ]
8.     end# ~- g$ z' i+ D2 \; E\" q
   
9.   end
   0 c5 o% `% M8 k: h! t\" e5 i
10. end
    ' g- O$ ]( x; X. i' R* P/ k
11. s1 w. X+ b! L' X
    
12. toc
    $ b4 \8 b\" R! g9 p* ?+ i: v\" ~) L
13. / N, q2 M. u. _3 q9 v
    
14. s =\" m+ l8 L% A! ?0 M: g0 w
    
15.     6.126720600000000e+0093 M# M+ E3 g: b
    
16. 7 H1 Z7 k3 x# C+ r4 h7 N8 g: u
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
' _+ \7 y3 e( Z- G: f' H0 rm函数：

1. %file xyz.m
   5 E! U7 u; p+ g2 Q\" \
2. function c=xyz(x,y,z)  t$ t, V. n6 g/ ]
   
3. c=x+y+z;\" s: R* W4 S& I6 W
   
4. end

复制代码

测试代码：

1. tic;
   5 B! L- H% I& n+ W. x: J
2. s=0;; ?0 H* f: Z3 P4 x
   
3. for x=0:10004 T9 C+ s* @. i1 C) V/ D* X# v
   
4.   for y=0:100
   : l, a. h' Q0 R9 w# S& @
5.     for z=0:100& S0 u\" _5 S# @! N/ x
   
6.       s=s+xyz(x,y,z);
     i4 r6 A' {8 D# g$ C( v% d
7.     end- ?. K\" c- k7 ?3 z$ f
   
8.   end
   ) A5 ?( e& p# K\" J8 |9 F\" `* j
9. end
   / k2 M$ P) ^: l. o, g& z+ D3 y+ L
10. s
    8 ]8 E9 |\" r\" L! Y9 g
11. toc3 ~5 U; l  [% V/ `\" O
    
12. 
    . z/ t5 f\" |5 Q7 l\" ?, ]9 j
13. s =  E( n: S5 F! P; Y% p4 B
    
14.     6.126720600000000e+0095 s) q\" o: Z& U2 L2 w, I
    
15. 3 u! V' b! l! @
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
$ B7 R; X" ?" |! P6 x
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;' g6 y: x' m7 s0 v- f  ?) X* }
   
2. mvar:
   ( I- I4 K9 m, D3 Q! h
3. t=sys::clock(),% s! c# w0 |4 O  b8 p
   
4. s=0,' T( q% y9 I4 u7 B
   
5. x=0, while{x<=1000,3 Y0 b, F! f& k* x
   
6.   y=0, while{y<=100,
   ( @( [& J- M2 H. Q! y0 r
7.     z=0, while{z<=100,
   ) r9 S( o0 ]- Z\" o! G0 c' G\" T9 a
8.       s=s+f(x,y,z),
   ) X7 ~+ T, X/ R. p: Q
9.       z++
     i4 v/ h2 D8 o4 u+ f0 }
10.     },
    # f5 F  N\" D* Q4 Q
11.     y++
      u, K\" k' L. q2 v7 f4 G
12.   },7 K0 t, t7 P% M1 o# R3 y
    
13.   x++/ H, u3 n1 f3 G. d2 z% \2 p6 W
    
14. },) ^. f$ M. L: f
    
15. s;& }8 M* K) O6 X  U3 W6 p: p
    
16. [sys::clock()-t]/1000;

复制代码

结果：
- K! _9 P- u5 r% x3 l$ U' F6 N6126720600.
1.109  秒
! u) M8 i( S8 X& E# K" Q
5 M2 x( n* E! ~) |1 M% S* s* u3 V二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
' Y" `" R+ G7 `! u- q6 I

hxftntA

0 R% U5 h& n' M* v, q
) O( R/ ~* r6 F7 O7 _; ?  G
) Z; V; [  p) a
9 h, Y8 Z! m3 L2 i

" B% i* N# _# L6 o, ^2 |* U$ S

- D5 y) A% `) v