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

forcal

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

4 ?. K! @* s: k+ ]% q4 j＝＝＝＝＝＝＝＝＝＝＝＝
4 @) \, @; J3 u/ a6 B+ R: Q2 d
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ) F( {4 F* H$ V9 |' Z1 n
2. tic;
   ! Y- E! U! m# ^' t& h4 x
3. s=0;/ M% Z0 n\" x/ y7 R
   
4. for x=0:1000
   6 R  M  i* O5 X
5.   for y=0:100
   4 ?6 ~0 {! v( d2 r/ O5 L3 e7 `
6.     for z=0:100- z6 p! F\" N0 q* @* ~
   
7.       s=s+f(x,y,z);  ^4 X3 w# [& Z* Y0 I% O8 R( S
   
8.     end
   0 M8 e7 E$ ^( i3 C
9.   end
   6 q7 y4 V/ ^+ D
10. end
    * v/ M1 N- p2 r# z- M. K
11. s; @0 l  ~& j3 y  v+ q
    
12. toc
    % f8 J  d5 o: O4 A
13. 
    + f0 n1 D) E( u' X
14. s =5 o) q\" ]$ [$ O+ C
    
15.     6.126720600000000e+009
    5 j9 Y) ?, c6 G; w  W; F; J
16. 
    2 y8 N' V\" P' h  r
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m8 H# @- t* w9 n0 a
   
2. function c=xyz(x,y,z)
   6 F4 m2 i3 B6 q/ W5 S- J1 m
3. c=x+y+z;9 F6 A  ~4 ]8 S6 N2 K
   
4. end

复制代码

测试代码：

1. tic;  F  T' W8 q( _  M# T5 S\" M, m
   
2. s=0;- V6 ^2 m9 t( W: L- e
   
3. for x=0:1000! q+ b4 _2 a' d# {2 I) J0 B4 v( x
   
4.   for y=0:100
   * l, [# l2 v( y& A
5.     for z=0:100
   . \- @. q7 z\" A! j
6.       s=s+xyz(x,y,z);
   ( J& d2 \1 c\" P) H
7.     end
   : ]% L' b( e5 U7 \3 c$ r9 E
8.   end0 A) s! F* T8 z' a: R
   
9. end5 h( A0 h. G) o5 t+ |\" L% S
   
10. s  M- g5 x- j3 ]+ K' v
    
11. toc+ O. s- u\" g; e- n
    
12. 
    3 B\" L1 ?7 f, B( L, Y% j0 g
13. s =
    % Z* {( U0 E; X; q+ [( l
14.     6.126720600000000e+009
    , s4 j& U6 m\" ]
15. \" t\" I5 r  X7 A+ t$ N+ z3 Q  z+ S* a0 s
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
+ a2 B! e3 w- c; _, g4 r. u$ n3 J
( ~2 H. X3 H$ u, BForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;$ b1 I3 K! o5 O: y$ z2 w4 f: V% B
   
2. mvar:
   0 Y8 p, |/ n  L8 H7 F- M# }
3. t=sys::clock(),
   $ |; p' U0 B( L6 ?8 R
4. s=0,
   # Z* O( s. P( ]! V' a
5. x=0, while{x<=1000,
   2 k  P5 ^( N3 {# ~5 J/ t8 ^' ~
6.   y=0, while{y<=100,
   5 v( F\" [/ ?7 Q6 v# u4 T
7.     z=0, while{z<=100,9 W$ E, t0 m6 o, E6 S3 ?) i& f* `
   
8.       s=s+f(x,y,z),
   4 i! R. P! k0 F8 [3 v
9.       z++
   * C8 C4 p& h; W. x, {& ~
10.     },- e, D7 v- B+ M: r8 d; i
    
11.     y++
    / _! b8 U+ D% }7 o
12.   },( ~* g1 p% d# J# }9 S
    
13.   x++
      X+ B; g: s' d9 p: V0 Q& P$ v
14. },: L0 o% C! p/ x/ o2 s  I6 X- U\" u
    
15. s;
    7 D, A/ g' h$ H3 k  ?0 t
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
2 y: ^: b1 E: C, X# Z1.109  秒
7 v! O0 H, _4 t# }* s2 ?3 z
" N/ X$ O6 ?1 ~: e# }0 @二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

# V" d  ], u4 ?& ^3 U
1 X3 t4 z' D, p6 X2 B8 y
* e4 o0 c) R* @/ R9 L0 U$ b1 v1 _
/ C+ k- E; v2 c& H( a
6 D1 b% _, w$ \
& f; h1 Q* w. o( h6 K



) L. H1 E8 h4 o% W& \+ l
* C) e' M7 e9 o
6 u7 D" s) d9 h2 c- [! a


, z; V# v$ g! T* C