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

forcal

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

( t' i6 E1 f3 {' Q. u+ l; U/ [* ]5 ~＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;- L4 I; Y& k: r
   
2. tic;: U3 d\" O, S' \, k# b: r* i' d
   
3. s=0;
   - u/ F, B9 r8 b6 A2 r3 i& W
4. for x=0:10003 o6 c( b3 n! k, d/ \- _  G
   
5.   for y=0:100
     l1 H# e. z  M) y; t* [  @7 f4 k
6.     for z=0:1004 L4 W% q) p* G) q, ^
   
7.       s=s+f(x,y,z);
   \" `9 Z# V- G- m4 s/ [& n3 t
8.     end1 d$ r/ z1 s% Y% G4 o* I6 k
   
9.   end
   - X6 K! D' S' s' `4 o, y
10. end
    / X3 E5 ?( O) z+ o+ E. i# H1 ^' y
11. s
    # s' o* D$ f7 }
12. toc# p  p& O7 m/ b7 G' n% Q2 v+ k
    
13. ; f0 y6 p. ?! }6 C. ^+ y( a
    
14. s =
    / m! u; H. [( }# t
15.     6.126720600000000e+009
    \" r5 \+ H3 Z7 T/ ^( i2 A4 I( {$ O
16. 
    5 A  X\" ^+ u( ]
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
! ~: _+ Z0 H0 r9 `9 f% Im函数：

1. %file xyz.m5 d: a) y% t2 Z+ ?& P
   
2. function c=xyz(x,y,z)
   & F/ c+ r9 Q) b, F5 i& I' o- J
3. c=x+y+z;
   ) a5 T/ M9 l+ J( M9 S* ?1 B6 w8 N
4. end

复制代码

测试代码：

1. tic;\" j4 u# b2 e9 D
   
2. s=0;\" R  `6 [; w8 u
   
3. for x=0:1000
   9 j# m1 r/ `7 i- W& X\" t
4.   for y=0:1000 n\" F& W. j. y+ `! |\" E/ y
   
5.     for z=0:100
   $ t6 L$ H! k5 J8 {7 @
6.       s=s+xyz(x,y,z);
   & c) b' c6 E4 p9 s% {
7.     end
   ; o) o. {5 P5 a& x8 v7 \3 [
8.   end
   ; T  f% J9 ^. ]; D7 v2 W
9. end
   9 i( o\" L/ b% x: C; I) u
10. s: p3 D2 G; s$ B. @. c& s
    
11. toc
    0 R\" i; ?3 X! S, H/ S' ?. F
12. 
    ' H\" I0 b6 j7 k- B
13. s =& f) ~; L/ d( V4 G3 v
    
14.     6.126720600000000e+009
    1 [' ]6 B* H( R! N) M
15. 
    ! m, N! M4 ~( Y$ d1 g8 ?
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

( l" e- H8 P. c+ j5 ~  JForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & k: Z: @0 N7 J# |5 n6 y& z
2. mvar:
   3 ^5 [( _( x$ |1 Q0 }7 t& P* b! L
3. t=sys::clock(),
   8 O; F\" x7 A6 R6 B
4. s=0,. O. p& t2 r* g/ o5 F
   
5. x=0, while{x<=1000,: ]7 v$ l# _' m& c, y2 Z+ x7 Y! S
   
6.   y=0, while{y<=100,
   5 Q! f& e! o1 C
7.     z=0, while{z<=100,
   + ?' h+ ?2 B3 r\" k/ R) X
8.       s=s+f(x,y,z),: ?7 {1 B. F0 V; O1 t# G
   
9.       z++' e\" I- |' a- Q! s/ O. A
   
10.     },
    8 G- U, ?\" z& c- e\" V; M) F
11.     y++
    3 d: k6 n4 x8 J
12.   },
    , X7 l8 d& r8 p# U( o2 v
13.   x++
    % w8 S% J6 q2 ?3 {
14. },
    * m\" t. N& V: C' L. J$ Z
15. s;
    / r3 K, a, T1 i! J1 L8 c  Y
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
1 G( c5 h. A* ~0 h$ E" t7 v( u! |9 V
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

, w+ g0 r0 r0 A* T# ?


/ D: o# K+ {" i; B2 U
3 @% }9 K% ^9 X3 ]! k; L& Q7 M. x
7 ]$ T1 a# T: m/ y, b

! c2 O" X7 P  _! x5 p
7 _) A, `/ M2 I5 t! @  ~
1 Q6 Y- W: O" E5 j) r




/ Y& Y7 }: b8 G5 {