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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  O% p9 ], f) `
, B3 t% F( C; h! n3 L＝＝＝＝＝＝＝＝＝＝＝＝
: n: Y/ K# g. h0 o" O* g
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;( x/ N; C2 C1 @* D5 w
   
2. tic;
   / h' X; T; O+ M\" G6 t! X
3. s=0;5 j( y  L8 V  i/ m
   
4. for x=0:1000! z% W0 Y0 o# G( }* n
   
5.   for y=0:100% G  \# \9 E( |$ {- }
   
6.     for z=0:1005 `4 K$ z1 X* K2 B3 Q
   
7.       s=s+f(x,y,z);  ]* q1 A' X  O5 @0 M
   
8.     end% v# `; D2 [) s7 x
   
9.   end
   $ T$ m0 M  m+ a4 z9 a
10. end' }% R* _  U3 x& ?
    
11. s
    5 E' Z* q2 w; h  z& [5 n' E+ R
12. toc% X% A- T4 I( {  w\" c' I1 H* w
    
13. 1 k% S' D/ h) n: |
    
14. s =$ N/ C5 _, u& E7 I% x\" X
    
15.     6.126720600000000e+009
    , I2 x6 }0 |0 x0 ?9 e# ?
16. , ]0 G( W3 c/ u, z9 ~6 A
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
  C2 x% H* Y2 jm函数：

1. %file xyz.m
   ( ^. \  ?' N5 ?
2. function c=xyz(x,y,z)
   ( e6 ^$ c/ q' P- g) d5 J# V2 ^
3. c=x+y+z;& |# l; `/ R& |# v. Y, x
   
4. end

复制代码

测试代码：

1. tic;. y5 Y\" S; Z3 }# Y
   
2. s=0;( p. ]- q$ O- j0 q& o% D
   
3. for x=0:1000; A0 _# R5 ~% X) a
   
4.   for y=0:100: C3 |8 I4 S6 I) U' j* p\" K/ d
   
5.     for z=0:100' K. W% K# [7 V6 s( D+ D
   
6.       s=s+xyz(x,y,z);
   \" a# {  |' Z. {, e0 F: I1 w. A
7.     end
   ! b% t; p1 ?/ W3 I% q
8.   end\" x3 d. L0 [. p8 [! r# W
   
9. end
   / e6 t( x8 u4 B% @& a
10. s( B  E0 ]) b/ d9 K
    
11. toc% Y, W5 t% P. Q. q' E
    
12. 0 a  I+ Z6 f3 U' c\" }1 }  g: }7 l0 V
    
13. s =
    7 {9 a\" E0 Q4 K) `2 b; X! S
14.     6.126720600000000e+0097 \5 a: S7 |* d
    
15. % W! a5 y/ j- ]4 ^
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

4 o' B& H; ^5 w5 H; }- JForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;! q2 U) l+ g' F/ d) i
   
2. mvar:
   ) [2 T/ f( n$ s: H
3. t=sys::clock(),$ W8 F) E% x% g; O) J' N0 H
   
4. s=0,6 Y& k$ Y$ D: y
   
5. x=0, while{x<=1000,
   $ M3 x7 e$ Y6 @) l\" X/ K% P
6.   y=0, while{y<=100,
   ) E& o: @\" I3 z4 }
7.     z=0, while{z<=100,$ x) ]1 w$ R* t/ N$ u; H+ E
   
8.       s=s+f(x,y,z),8 N& i+ e8 R0 z' ]) q4 {/ m* a
   
9.       z++\" P. K- S# z0 e& ?: a7 F6 F
   
10.     },+ C; r3 f3 G6 L4 H
    
11.     y++
    + i; C- S+ X- E5 @% E( I- w
12.   },. a* t8 _+ K0 _5 t; e
    
13.   x++2 x+ ^# e% a0 B\" o: o
    
14. },6 k2 W; p* w' u# l$ U
    
15. s;
    + n/ {- U0 O$ ]
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
3 Q6 V' o7 ^* H% M/ h6 }1.109  秒
$ Y/ _3 G6 [' `; X1 f0 x
: x' t* R9 L( P+ j, f' o: o- n二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
3 u  g: a$ r- R* A

hxftntA

% j$ Q# @% n3 i8 ]. o  ]
% N& f* s% p$ Q$ s' w% U: p: i


: U8 Y: W2 @" c2 m4 y8 I' G1 @# g


( M1 n, r8 D, m/ w/ m- c8 L  [7 v
/ e# n6 M5 v. l& S/ z: u$ h' ~6 U
) j- [0 Y4 @) T& F

4 M7 P2 H4 I+ L