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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
& r  ?8 s0 c3 i; h" |
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;7 t! l9 H, D6 q  \
   
2. tic;3 |$ N  d3 j& S& b0 P( ~' u+ \
   
3. s=0;
   # |- G1 W( h  U9 _) ?) T5 ~) x  G; \9 ~# u
4. for x=0:1000
   9 u8 Q- V' M& C! R% U# D$ J, I8 P
5.   for y=0:100
   ' `\" U6 o/ ]( O$ E$ ^6 ~
6.     for z=0:100
   + a5 l) k) M# Z. t+ a, X
7.       s=s+f(x,y,z);- O- _+ J9 @; N
   
8.     end9 u% A# S$ m2 p6 u+ j! ?$ V' z
   
9.   end9 T0 e# P. a! |( Z+ T
   
10. end
    9 E+ O, B0 i& u! ?. d; a
11. s! _: W. \5 m  }9 T
    
12. toc
    + T) ~$ r0 ?% I7 j) r4 |( p$ \. ]
13. . u. _4 G2 B0 W9 L/ `6 Q
    
14. s =
    2 k* _* L  p% o/ V0 Y9 E% ]
15.     6.126720600000000e+0094 [* _; ^) \2 _/ t+ r# j
    
16. 
    8 `1 N2 D# w- d
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
$ _+ g! G6 K' h, x  y1 }! \m函数：

1. %file xyz.m
   + E3 n7 {# c0 b
2. function c=xyz(x,y,z)7 X5 B# J* z8 Q! t7 d- d\" Y$ D
   
3. c=x+y+z;# i! Z/ Q; o* L; `/ e
   
4. end

复制代码

测试代码：

1. tic;
   5 H- J. v/ W5 @\" @/ V2 G
2. s=0;
     p& t! E7 W' D
3. for x=0:1000/ w) m4 V' e0 E2 g3 J: H3 m
   
4.   for y=0:100& {4 `1 f' l( i- ?$ `% d2 X
   
5.     for z=0:100/ R' Q, z. T' j
   
6.       s=s+xyz(x,y,z);0 J, I6 {% B: h2 |4 f
   
7.     end
   ) H( f8 K# @$ G7 `6 E# D1 ^. E/ d
8.   end
   + Q\" a5 n8 M; H( S* n0 O) S
9. end
   8 ?& b\" b/ j+ Z; z
10. s- w6 X* C( D7 B2 b
    
11. toc\" W  i/ t; d6 M\" P) c
    
12. 7 M) Y  ^: v* t8 V. W5 x\" t
    
13. s =
    9 j! K\" j4 M5 Q  R5 u
14.     6.126720600000000e+009
    7 K% b2 r7 ~' u1 Q2 }0 }$ U
15. 
    % A. x7 v+ k. ^! _3 j% G
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

$ X' i  h) j; {( bForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;! F, _! L\" M, y
   
2. mvar:4 Z0 r# z. V. z, P0 C; x$ j4 q
   
3. t=sys::clock(),
   9 S% G- t$ I\" q% J
4. s=0,
   1 G% `8 I( C$ f
5. x=0, while{x<=1000,3 C$ A5 D( l6 [+ P$ k& X& E
   
6.   y=0, while{y<=100,
   , X. q% M/ ~8 c: N; Y5 I
7.     z=0, while{z<=100,
   3 s# q3 {9 [$ H% `1 R' ?  w
8.       s=s+f(x,y,z),
   3 O) g$ d  h& e& ~) X
9.       z++6 @# l. x4 Y6 E8 p, O
   
10.     },& K5 E, Y/ h9 f5 M' W+ B
    
11.     y++# u\" U6 K' O( o7 C/ S) Z
    
12.   },
    # s; ~7 H2 [$ c1 ]  D
13.   x++9 e! O% V! k: ]' o% {6 v
    
14. },; l/ C8 D: H- c4 H! i4 ^9 f
    
15. s;
    + S3 E; p( k: Z- Y7 t
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
# k' m- r, _8 A. d1.109  秒
6 h6 `6 }9 M8 o- [! F
4 V2 z# d5 j/ C: C) g9 k二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
7 u. _* S5 u) Y

hxftntA

6 |7 r! w3 F: v6 D




  v$ ]9 Y  \9 u/ g

. x, g/ G  N- ~



" [0 Q8 y% g4 q' W7 T& Z% j
0 x, F, l' }/ a0 O' h

. v  C! l# b& }; f; [, r