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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
! G, |3 D# U4 p8 W
! L' {3 L( K1 l# l8 G＝＝＝＝＝＝＝＝＝＝＝＝

/ d/ V3 `0 C8 P5 t) s$ k; oMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   & y3 v\" j+ \; L! o( E: S
2. tic;8 R( Q6 ~7 J6 _0 k& g' B# @0 C
   
3. s=0;/ F  v. {0 ?& V
   
4. for x=0:10005 a2 E# \3 N! z' x' P+ [
   
5.   for y=0:1002 \4 f* k' I4 y% A
   
6.     for z=0:100
   6 C% o, [+ M4 \' q6 G5 h
7.       s=s+f(x,y,z);  T\" ~1 l. _4 c/ |1 \* ^+ F9 E. n
   
8.     end) C\" v3 q* S  n; K$ d, U9 ~: ~
   
9.   end
   ! E, T% g7 h0 D  f* g3 s' A
10. end\" V) @( U% S. Z, B
    
11. s
    & f- }9 @6 L# S' Q! T! ?; V
12. toc
    2 q1 s' ]# G8 L: y( S/ L5 A3 K' @$ `
13. $ @+ c; w; c* Y
    
14. s =8 [2 E! p# I5 b9 k
    
15.     6.126720600000000e+009
    ' B- U1 B2 ~6 X- W6 L
16. \" h) X4 R0 p; O* O! ?% B
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
! g' C. m5 s6 Mm函数：

1. %file xyz.m* b+ t$ b- a$ i& g4 Z$ w
   
2. function c=xyz(x,y,z)
   $ G  i' [( @\" b0 W/ p' q$ L: ]
3. c=x+y+z;
   , I+ j& d- P6 X$ g& R
4. end

复制代码

测试代码：

1. tic;+ A( G: W+ h2 n( ]6 s% \
   
2. s=0;
   1 P* B5 G2 r, V% l. Y' }8 |2 k
3. for x=0:1000
   0 v  Z0 C; I/ U9 ?
4.   for y=0:1006 I% `3 z- m! ]5 }. Y
   
5.     for z=0:100$ f$ g+ B0 F. |4 m/ U0 l# @* |0 V
   
6.       s=s+xyz(x,y,z);
   8 H( H* f  i# _; ?- |
7.     end
   & m1 I1 w( @2 d
8.   end
     P) z* b; \% R7 [# p& r% S
9. end0 y2 ?- e& s/ \
   
10. s. w\" O2 K0 A+ W6 ?! J2 v
    
11. toc
    9 I8 Q6 i3 v, Y% B& o% s
12. 1 x, O2 ~9 d+ c6 |4 H- n
    
13. s =
    7 h3 |  W8 H- T5 `
14.     6.126720600000000e+009
    , R3 p1 f' o3 h9 x0 p2 Z5 u6 @
15. 4 a  _8 a& y* y\" n: z+ H; b
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
% g, U  r5 I0 M- X& s, E
# d8 V- ]8 @1 |0 YForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;) k* x/ D' a: ]/ S/ T0 |# y
   
2. mvar:3 C+ v( b/ ~/ v# {\" V* U  \7 b
   
3. t=sys::clock(),
   2 l) B0 b- g2 ~' o+ m
4. s=0,9 F- r7 X) _: J+ D
   
5. x=0, while{x<=1000,
   1 |' W( `/ j9 Q/ @/ {
6.   y=0, while{y<=100,
   0 F0 k3 V' J\" O1 i; N2 V9 v
7.     z=0, while{z<=100,
   . [$ s8 G: t8 x7 m* W
8.       s=s+f(x,y,z),\" i8 s& I3 w- {
   
9.       z++, f' u. i  S$ l: h! B
   
10.     },
    ' t. u& J0 x/ ?! R$ W+ v: V6 |# {
11.     y++
    % G2 ~! s/ C4 t4 j; ^1 g
12.   },
    ' R- Y. l* @7 k; o  L, ~. z, u
13.   x++: V0 M) S- S\" f1 y* r\" u+ z2 a
    
14. },
    4 U9 a5 U# e, y4 ]' ]5 ]' y
15. s;
    6 O1 j- w* b5 {. I1 I* I/ o
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
  w5 S0 w* q2 V) @4 Q7 A4 Y  v1.109  秒
% |1 S- [* y0 d" F. V$ s& o% V; N
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

  h; j# E, p1 h, v3 Q+ S4 \

2 q, I" T/ b# C' D) U6 p
, c2 K$ X1 k2 {+ D1 ^

% A5 P: h, v, U6 ^
9 k& G4 Y) _% Z3 u2 E% l
; j7 z* }4 G1 Q+ ?  _4 a3 P6 V
5 [: B/ N$ K' ?2 \
6 m% j0 z: p) @9 E' Z



: s$ n/ q- q" P3 K/ o" {$ u( f