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

forcal

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

" Z: b% ^' v* z4 s% u3 u＝＝＝＝＝＝＝＝＝＝＝＝
' _! R5 F$ a# w8 \2 L% r6 R
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;\" Y; T# E  x! M/ c* _
   
2. tic;/ b0 s3 a+ j! g9 \! q- Y
   
3. s=0;
   : X; L( E( {/ ~5 ]/ J% I
4. for x=0:1000: p4 }2 |& n% M2 Z( R/ t, n. T
   
5.   for y=0:100
   ! m& X7 B  V\" ~7 u% M
6.     for z=0:100
   ) b9 b% \$ Y! }3 F8 p8 N
7.       s=s+f(x,y,z);# U0 K2 d, Y2 C& }
   
8.     end! h1 f\" y4 R) T* R
   
9.   end- h, y7 @6 q5 R+ I, }
   
10. end
    % U8 f! c! m1 o! n
11. s- z# a6 f! p# [0 A) p. y# K0 X( D
    
12. toc+ i; \1 W9 g6 t3 \
    
13. 
    ; K+ V6 h- o. O: [
14. s =: ?5 l! {\" k. R
    
15.     6.126720600000000e+0090 h. F$ K. n: M8 \) W) m* R4 ?
    
16. $ {7 B' Z% ~1 ~6 L# F
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m8 P8 X3 D' v# V) J& `% n1 N
   
2. function c=xyz(x,y,z)5 i8 A\" i% B; u4 V: R4 X* v
   
3. c=x+y+z;
   9 N: w\" X# `) E7 z
4. end

复制代码

测试代码：

1. tic;9 I. K  l2 x5 ]& y( L9 g9 w
   
2. s=0;
   5 P) F& @+ K/ G0 o
3. for x=0:10009 `8 k+ F- D  ^
   
4.   for y=0:100
   ' I1 \6 S6 Y- x
5.     for z=0:1001 E. k, l8 N- T# p) N7 H1 f
   
6.       s=s+xyz(x,y,z);
   7 s8 U8 I/ @4 I1 K0 F' Z' N7 C8 {- H
7.     end
   5 _: [2 Y/ W/ o( \  W2 c; P' @- }
8.   end
   + R& z) ^  \: V5 E; Z( W, z+ n
9. end
   8 |0 i) z- [* \
10. s9 \0 g  Q- {# [/ E
    
11. toc9 d8 A& b/ S3 e- @$ C- I, r
    
12. % A. y7 ?% t$ @; e# [4 j( r3 p; m
    
13. s =) E, d; I0 V/ ^  i
    
14.     6.126720600000000e+009: ~+ R2 |1 |9 |. p, T' N# \
    
15. 
    ! N1 W( [8 D: k$ k' d
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
( ~0 _& X5 [+ c$ v/ d
% A9 N( {' r8 \Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ; B( N) K! d0 ^5 f9 d\" O6 D& K
2. mvar:
   7 R8 U, x) X( E
3. t=sys::clock(),& ^7 u/ u, _) x, P$ ~2 C
   
4. s=0,) l; s4 k' m& ]' F
   
5. x=0, while{x<=1000,$ s+ z! O* Z- v+ T& ]8 p0 \
   
6.   y=0, while{y<=100,/ a  |\" l+ t, _/ H  @( W
   
7.     z=0, while{z<=100,
   ( n& Q$ I7 i. e% K, T& S+ L: `
8.       s=s+f(x,y,z),
   4 e: V$ [; \\" h  X' b2 k1 J
9.       z++( S* @' p+ l# c* N
   
10.     },
    & k) b* C% r: n& D- q, ~
11.     y++
    / A; `8 J0 ?+ {7 I- q4 T
12.   },
    : A& d! |& `- S
13.   x++' [5 f. k6 Y3 V) c3 Z, P6 ~+ a2 I5 I
    
14. },2 L1 V\" {0 @7 a# j! L4 h
    
15. s;
    , h1 K' h. X) L# ?
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
9 I& g  N5 q3 _& @/ i  m4 e$ ~( i# ?
6 `5 `: z  q. b9 L- N' q; Z二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

) q/ m4 [: w1 M$ D3 j9 ~/ `4 s






' n& L% _9 Z: s! w: n2 {

, u; o5 F# \  n5 n& D& ]7 L' I

0 ]1 Y8 L( W( W% M; c

) R! Y# F; z, Y0 Y# z

% W& k3 T. v& l/ I/ @9 I9 ]