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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
% E$ j5 y" b+ W" ]% \  n+ Y) H% Q
) j6 T8 Q: d& z* ?8 X- X＝＝＝＝＝＝＝＝＝＝＝＝
' h5 x; h8 q2 P. D" V
0 Q# t( q4 m! n9 lMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;9 d5 Y7 S\" m3 D7 J' u) B: Z
   
2. tic;$ M% ]4 u% y% G) i/ H# M
   
3. s=0;
   : H5 ^' K1 }& X5 ]4 @4 q
4. for x=0:1000& t- ]- f7 ^' L  `
   
5.   for y=0:1007 C) ^\" E* y; ~; u+ ~
   
6.     for z=0:100' l. x0 ]! K5 g) Q  l9 z. D. K
   
7.       s=s+f(x,y,z);( C: R7 i5 m+ E2 f% f
   
8.     end. b$ n! q\" K8 s! Q2 K8 x0 B; o
   
9.   end7 f\" a% ~: k$ m. r
   
10. end$ L9 @  t* g. }$ K# G
    
11. s) Q' E( p7 a9 F1 U# d6 N' t3 n
    
12. toc
    . k, L/ H3 R% g( x
13. 
    \" h\" C! e( n% l$ F9 w6 v
14. s =; I7 R5 H: j% k% w) G\" h
    
15.     6.126720600000000e+009- x+ f8 c. T9 u+ N5 o2 g) i4 \
    
16. 
    : [5 L* g( i$ ?6 Z
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
9 s" b3 S1 D; `- bm函数：

1. %file xyz.m
   3 r; z) ?0 w& ?' i
2. function c=xyz(x,y,z), _; y% h' d% D$ s( [2 N: ?\" H
   
3. c=x+y+z;0 @7 I8 K2 O+ W6 ]
   
4. end

复制代码

测试代码：

1. tic;% u( X  H% W( k! R4 `; q% K8 ?\" v! U, o/ v
   
2. s=0;; c2 I0 ^* M; h* l- z+ X/ D% D
   
3. for x=0:1000
   / G' c  i. f* J: K
4.   for y=0:100\" x2 A9 @* k\" B\" ^5 w& w
   
5.     for z=0:1006 g( ]9 w\" Q, W7 ?5 _) i
   
6.       s=s+xyz(x,y,z);\" n% g$ {% c/ d) `3 J) y# \
   
7.     end6 A+ x* j* {: W1 {! R1 o\" I0 }& v
   
8.   end
   ; q5 k1 s* K; f7 i3 L
9. end4 @: Z, W- ]4 H4 E7 L) Q
   
10. s
    / C. ]2 Q) V2 B5 l
11. toc
    . Z; D3 r) f5 f& w1 C$ G0 P
12. 
    * e* h, K  `; Q. y+ {
13. s =/ b, w, X9 Z8 }3 t
    
14.     6.126720600000000e+009
    ) @  b# `- G5 Y, U* p7 y4 L# G) Q
15. # C5 G4 i% Q& F0 B
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   2 y7 e5 @8 D- T2 R4 w2 `( b$ O
2. mvar:
   ! ~% E, b) S! `4 g, f
3. t=sys::clock(),$ V/ l: ]/ L( H( `$ q- m
   
4. s=0,) z\" w* J+ B. t
   
5. x=0, while{x<=1000,
   + J$ c\" e+ S2 ^2 k9 ~* Z! b
6.   y=0, while{y<=100,
   4 E. d* a' ~( y( Q
7.     z=0, while{z<=100,
   , E) i1 p% s$ G- \$ X& t) |
8.       s=s+f(x,y,z),: N' ~5 L8 Q  X7 h\" ?
   
9.       z++
   % f+ j( S\" r, j* I
10.     },; t7 g( ]7 l& H7 h\" {) N) H1 K7 d
    
11.     y++
    5 b# d/ N+ t2 y( U\" F6 r
12.   },
    # n) H' e: Y( P\" }' o; b
13.   x++: l8 g\" H$ E; J1 R3 D
    
14. },
    0 D( N' [\" ^: o: }+ J
15. s;
    # `0 x. o9 E4 V5 E
16. [sys::clock()-t]/1000;

复制代码

结果：
6 i" [9 l( x& i- E; v; H: ^6126720600.
5 i5 M/ }- [4 h$ ^, [( R1 T" P1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

/ s" g, p; i  W" Z( I3 i8 X( j/ d: v+ _



! J" ]' N- H2 v


  n) Y) h$ I- E0 R5 W6 a' x
  V$ I, N# I) O. x. L( p, G: J# g. ?



2 U' Q0 u  T2 b/ Z- \

& S+ c2 Z- p4 W8 W. O; [! ]
. Q; c9 b9 U7 g% ?/ h. B
) ~& T6 x1 @8 S2 C