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

forcal

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

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

# }# m7 s3 [4 ^' XMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   \" \* u! o' [3 T0 p$ L$ ?
2. tic;
   ( _5 m% r7 M! }/ q( D1 R$ q
3. s=0;! p. Z$ i1 H% H. }3 X5 f$ e0 |
   
4. for x=0:1000
   ; z2 q1 R\" O0 g/ U. o- }\" v
5.   for y=0:1006 o, _  M9 e$ G; N5 |
   
6.     for z=0:1007 N6 c& w5 H7 }0 O' S* q
   
7.       s=s+f(x,y,z);
   4 m2 n0 f- n, D4 b9 W2 S3 A0 Q\" `
8.     end
   3 \- ]% l2 _7 \
9.   end
   % Z( d+ K; ?7 N( n
10. end% d# f* s# ^7 @5 c- s# i- E: w* x: [
    
11. s
    ) ]7 [* t: T( n
12. toc3 [8 a' C( }! ]% U5 P6 b, H
    
13. 4 X% a4 U% n$ I5 Z5 M
    
14. s =- r3 s* B1 {' s$ j
    
15.     6.126720600000000e+009
    ! \6 Z1 A! Y7 I8 Q9 N
16. 
    4 t\" F/ \3 s) i8 g& i2 N3 `
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m& _3 ^% I3 c1 {. A8 u9 F
   
2. function c=xyz(x,y,z); i! B) e; B0 x* r
   
3. c=x+y+z;- o6 ?$ W7 F. e+ s/ R4 D+ z
   
4. end

复制代码

测试代码：

1. tic;\" F/ t! B4 j$ E+ g1 P9 b
   
2. s=0;. S4 B- G; @! M) C1 {% r
   
3. for x=0:1000
   % t/ m8 G% H6 |6 o
4.   for y=0:100
   ! b/ w! {2 ?# f0 w( L% p- S
5.     for z=0:100
   9 W* H5 p* F4 X* [0 w
6.       s=s+xyz(x,y,z);9 O& h- R# T2 U3 ^: G
   
7.     end
   \" `- d! I1 t# u# C: o0 n\" R
8.   end
   5 J' ^6 ^( }! b6 u# i/ I% N
9. end
   5 H7 p. E7 s3 P4 B
10. s\" I) n' G+ ~& e+ t( l1 d
    
11. toc7 Y: o\" N8 {; s, U% C4 T
    
12. 0 X  s5 s, n) C
    
13. s =8 _- F0 E% w. r; i3 B
    
14.     6.126720600000000e+0099 i9 r; ]8 C& h1 ]* I6 p
    
15. ; \' g+ A8 z% [0 V% N1 a5 f
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
/ U* Y( M  V5 h. w! v
& @, S( \+ H+ x2 p8 a* c# \: {Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;) b# a4 x8 w3 p; f
   
2. mvar:2 A! ^- C\" V2 F
   
3. t=sys::clock(),7 q5 K* G5 A0 a' C
   
4. s=0,
   / [' U\" U- Y8 m$ A4 K7 G3 Y; r
5. x=0, while{x<=1000,+ C& \* d4 _# O  h3 d; w
   
6.   y=0, while{y<=100,' ^\" p) a/ l& S5 m; V) n
   
7.     z=0, while{z<=100,$ w& ]5 ]2 h/ |& Y- v, E
   
8.       s=s+f(x,y,z),( C\" y+ x0 R& \4 c  ^
   
9.       z++) H& s+ B% ^+ r2 i+ X1 S+ c* X% F
   
10.     },
    * G\" [( q8 b' O* _8 z
11.     y++! ]# p8 g  l% ^) `7 l. `* e
    
12.   },1 v. J! E* L( w/ T5 c! d( Y; U# `
    
13.   x++, V3 \7 n9 w. B3 B/ z9 u9 u0 p
    
14. },
    ; V0 f# p6 d3 p+ }. k
15. s;
    1 y+ t* o% S' o0 z/ @- N
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
# o2 n& g8 U5 }) ^1.109  秒

  |; ?8 `0 J* z% S! a: h二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

1 b/ C, Y3 @) G7 \  ?( t* V8 k( d


7 Z" E+ @* G; Z, {9 R


3 N- v+ U4 y4 i2 E+ t
! `. ]( j* _4 V: f4 Q
2 g# }1 _1 e  q# @+ ?


  Q, }& X% c$ B* K' |1 @
* _$ Y- ~9 g# {  @) ]8 Z* R
% z. M7 n" W0 M0 E
7 r$ S4 h; B* o0 v# r9 y
8 n/ [" |* w5 V/ j