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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
, f6 J5 u! {: l
8 j" R& I1 c0 c4 X0 p1 a( YMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;1 I  p3 d* k\" m8 [9 q/ H8 U0 A
   
2. tic;0 l\" M2 ]4 D5 g( ?' i: e1 K
   
3. s=0;
   2 F/ o9 l: P' w; X
4. for x=0:1000: z% ]$ b, W) w: }
   
5.   for y=0:100( d, F. W8 Z, P# `
   
6.     for z=0:100/ d9 `$ r6 i9 g/ ~  J/ }) o
   
7.       s=s+f(x,y,z);
   . s2 U1 ?8 p' e5 K1 e
8.     end
   9 z% T' y, G  ^0 ^$ K5 x' S
9.   end
   ! K  q# C. u5 G* a1 p/ p
10. end
    ' e- d& L\" m) n8 Z' w: i; \
11. s
    1 X2 p3 u+ l) B3 U
12. toc
    ) J. x9 h\" t# K* K
13. ' Y7 |8 M& Z2 g* S3 e' p- x
    
14. s =
    7 j9 Q$ G; g8 x
15.     6.126720600000000e+0095 h. D) Y: D+ b1 u' y
    
16. ! W6 l1 e4 g- Z, y) ]
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
! p7 U, T  P* f  R9 j4 [9 a; g& c: z: Im函数：

1. %file xyz.m
   / G) T( K3 {( g5 Y& O
2. function c=xyz(x,y,z)
   ; C. {! g1 `9 D1 V- h\" m
3. c=x+y+z;
   7 v1 D, }\" W. ?/ G' o
4. end

复制代码

测试代码：

1. tic;) U* G' [/ j% N
   
2. s=0;6 ]! F( F/ ]% I- d! u* w
   
3. for x=0:10006 L) I6 d- {: X4 \# _
   
4.   for y=0:100
   3 E3 ?- T/ O6 ]6 o
5.     for z=0:1009 J% b5 H. ?/ {. l6 A
   
6.       s=s+xyz(x,y,z);
   7 N+ O: T+ K5 S6 x- b7 d& e) Y
7.     end% k/ W7 I# N\" g2 U2 r9 u$ F* R
   
8.   end! U5 T+ Z' R2 c4 G0 A9 K+ h
   
9. end
   9 i* y( z6 Y\" i( J! U) K9 S+ K7 `
10. s
    $ q+ g' E$ K- m\" u% w* T$ Y- r
11. toc4 w2 S& r\" [  V2 [0 V
    
12. $ E% d6 W+ k& Y  `( w, F
    
13. s =& \, x2 b- ~9 w9 U4 {! w
    
14.     6.126720600000000e+0093 L. N* f- S/ y\" `4 L# X
    
15. 
    ( ~4 [5 @3 g7 n, m& j
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

" h. e% P7 V- v. k; XForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;, M, }  A: z$ b9 N% E
   
2. mvar:/ s1 Y6 a6 w( a5 n% o4 n
   
3. t=sys::clock(),
   / f\" |( {+ c/ F7 e0 n
4. s=0,1 G1 v- }0 y  S; Y
   
5. x=0, while{x<=1000,+ q1 X5 R7 u+ Z3 Z: o
   
6.   y=0, while{y<=100,
     n. m0 L: H; e$ z5 X* U. _: M
7.     z=0, while{z<=100,# |( {/ k8 z: r3 C: l3 H
   
8.       s=s+f(x,y,z),1 Y/ I, y. H\" [  \4 ]( L
   
9.       z++
   1 C1 X7 Q/ g  f1 G8 b
10.     },( t3 L' Q+ o9 }' I2 e
    
11.     y++! g( P$ _! S% e5 V& J9 N& q
    
12.   },
    / x3 m, b6 K: L$ Q# {' E' p$ t
13.   x++7 c9 z5 Y; ?. i& Q0 L
    
14. },- r. m\" C; X( H2 L# m0 F
    
15. s;
    . X8 H& L4 Z6 M: Z7 n5 p, B% f  V; p
16. [sys::clock()-t]/1000;

复制代码

结果：
" G, b$ t9 Y) Q7 o) m6 @5 [, S6126720600.
1.109  秒
: q) F: U2 _' }+ G( p" T. n& i
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

5 z0 y- _- u, z! r$ s: ~


  p% ]9 d2 s  H, R) d- [
9 a* u. B* H: l) b
& d' |% I3 k& l1 s9 C
& H1 @; T2 G" l% ?( e; h$ I


+ z$ u4 ~$ k, H
0 s  f9 ^' B7 T
) ]/ n/ {) q+ X3 B
+ D/ ^6 f4 |7 o& ~
! w1 s4 X2 q& G/ {# O/ }

" }# Z. O# ^# V( N, L0 o4 A
6 x' H1 c# W8 U