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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

$ M: M% i2 D, W% \; DMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;( x  V8 d! Z/ E% ?1 h& x3 ~/ s' z
   
2. tic;3 S, l- u# O  V# y% j3 S( `
   
3. s=0;
   - {1 F; C. S\" C4 S; ^. V* y
4. for x=0:1000- X& _3 n5 H& S
   
5.   for y=0:100
   7 ?5 `& B2 y$ M; P$ [- P, q
6.     for z=0:100
   ' x\" C0 X) |- ^% `5 T
7.       s=s+f(x,y,z);
   & j/ t% w: m/ M1 r
8.     end: R+ Y- s; B! V0 d7 p
   
9.   end* K4 ~7 b, Y& |- C5 R
   
10. end  z+ t$ B6 S1 N! I) Z7 t* n  \
    
11. s' F6 z- B6 T4 i, e
    
12. toc
    * Z- T- v3 J. O5 Q\" G( u% ^- c
13. ' t, O$ o: X8 p. _- R& G
    
14. s =1 w$ _\" q  R& Z, S& D\" W
    
15.     6.126720600000000e+009. ?5 {2 d. c8 ~0 W7 r
    
16. 0 s& F# }2 z, @, _5 w* x3 n
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m+ E0 s; K' T( X8 |+ P3 F
   
2. function c=xyz(x,y,z)
   ( T! y5 J/ n2 R% o  _
3. c=x+y+z;0 R- O% O- z# @0 a
   
4. end

复制代码

测试代码：

1. tic;
   ( K& P3 M: Y) ]  U- U\" q5 }
2. s=0;
   0 k5 ?3 o$ k% |0 u$ n1 C) G
3. for x=0:1000
   : }' E\" s\" w  P0 ~
4.   for y=0:100
   . T0 j& M+ Z# J. f$ n( J
5.     for z=0:100
   ! y\" c% O/ Z7 U* ~, c
6.       s=s+xyz(x,y,z);: b+ K+ H1 x& K0 H
   
7.     end
   + H8 F- s$ n7 E\" [; O+ e) l
8.   end8 |% k( |( u% ]3 z$ E
   
9. end
   0 N1 w, }6 m7 @4 j
10. s5 X\" |5 V! D0 q\" |) n9 M3 G6 x
    
11. toc% S9 A  }& ^  T& X9 o
    
12. 4 l& Z9 S\" b& ~' V- W+ }* }) f3 O
    
13. s =) ]- T2 X* c* }$ }( h% L\" @
    
14.     6.126720600000000e+009! Y$ _& y. Y! c$ n0 a# R
    
15. 
    3 J/ f  N- g. f% m7 Q
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
( v7 n$ J9 {! z
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;( g: e  Z0 z/ Y$ i  {
   
2. mvar:
   0 D/ ?; f8 h) z  Y
3. t=sys::clock(),
   / f' l$ v0 l' Q9 x$ u  }
4. s=0,
   + }( t% y1 x' N; b
5. x=0, while{x<=1000,
   * d( q) h' `+ X. h
6.   y=0, while{y<=100,& {1 ^# d6 K' e1 a/ Y
   
7.     z=0, while{z<=100,1 ^0 n' ^\" J\" G7 M- b7 g
   
8.       s=s+f(x,y,z),( J- y7 q: C3 @3 q
   
9.       z++
   9 d0 G! o3 V4 j, S* Q2 }- `6 D7 N
10.     },! g: ^+ V3 N7 H5 M
    
11.     y++1 N' b! l\" o+ e: G, J# Z  ^( ~
    
12.   },, ?4 }/ p0 s& {* |8 L7 W
    
13.   x++/ Y) t\" `! I  m4 `. U3 p: J
    
14. },+ @7 a4 X! Z# S' _- l! r. x# R
    
15. s;
    5 Y) ~0 m$ }& x- l7 z% A6 y* `6 s
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
! r+ C; u/ q* M' O8 Z1 i1.109  秒

二者效率有几倍的差距。

hxftntA

$ L1 a% h' H0 s7 e5 e$ ?# ^
9 v0 k+ x, m( l# D
$ X9 J# S+ {: Y+ w* R/ W
) A. `  ~" \" D  ?1 N" ?



, a& [) C) Y  v8 p0 u: R

  f; C) `, R) B5 M9 c- o7 O9 f
( V# t4 W0 d8 e

ehi28

嗯,不错,支持一下.
# H2 o; G: n0 t) C