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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
2 f" Z0 y' v6 q
  R+ ^, W. s$ C: w7 ?' F. j$ j2 j＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   \" ]( F5 f7 f$ Z( L& I
2. tic;; g+ E7 G& w4 }* z8 z
   
3. s=0;
   ; r; Z+ d; G$ C8 b. P, B
4. for x=0:1000' L9 J3 }9 \* x; ]. H; S! b2 h/ b
   
5.   for y=0:100
   7 [) R$ u, J! l8 v7 h
6.     for z=0:100
   $ G# @9 j: w1 V8 N. w
7.       s=s+f(x,y,z);
   ) Q1 d; Q7 w; r# ?4 y. ~
8.     end% b8 ?& {3 F) C: w\" R+ q
   
9.   end
   + S1 \/ t4 I0 m
10. end
    8 ]! ?2 A# h+ }7 h\" Z  K  z4 B
11. s
      g+ u: l, S  M
12. toc
    7 y' c$ y& @1 f
13. 
    ' {$ c3 S9 a+ s/ `/ y
14. s =- V/ _$ L  a! n# G2 @
    
15.     6.126720600000000e+009
    7 N# n1 k4 |, d8 m% }
16. 
    % s# C+ O: g* s; |# `$ w
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
5 n# L7 @) I) r) `, \+ m/ O7 _: Tm函数：

1. %file xyz.m
   . `; c3 Q( ?6 w* ~* W6 F8 m1 a
2. function c=xyz(x,y,z)5 t- O8 _& o4 B8 u
   
3. c=x+y+z;' X  m* q4 d/ M8 ^
   
4. end

复制代码

测试代码：

1. tic;' ?, E- H' R$ \8 L( ?  N
   
2. s=0;
   8 W) h4 M+ F) ?* B
3. for x=0:1000
   + A\" ^3 n# {1 b\" U
4.   for y=0:100
   3 W! R# d+ |0 I( n
5.     for z=0:100* y! X5 j9 Q9 R
   
6.       s=s+xyz(x,y,z);. B& ?$ D6 r* G$ h, p7 I
   
7.     end; l3 a( C2 e# v% \$ L
   
8.   end
   / V1 ]$ z4 \2 ~9 C, R3 c9 p
9. end
   7 ]\" O6 `! X& V9 q* D2 e6 K
10. s
    1 x% E  S% u- J7 }7 D
11. toc
    ( x9 S. O7 }- E: z
12. % q7 X& v3 ?) Y! E' v
    
13. s =
    \" F% t& J) S7 _' |7 t
14.     6.126720600000000e+009
    2 r1 `0 Q  Q. H0 N
15. 8 D2 F- k+ M( h3 R, A( W
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

: ]( Y; `! N  R8 S" ZForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;$ c. \; h9 @% A# [
   
2. mvar:
   4 o$ Q. A8 X: ]) x9 n! p
3. t=sys::clock(),
   - |; h; V: P\" C' n8 a; b
4. s=0,
   5 W. a% }9 _5 s
5. x=0, while{x<=1000,9 Z4 J( C  g2 f+ h) S8 R
   
6.   y=0, while{y<=100,! f7 k* ]  K& O
   
7.     z=0, while{z<=100,
   6 x. X4 o4 S9 E6 Y% M6 N
8.       s=s+f(x,y,z),  C& G& g. ?2 y' I+ _1 F( E- ?
   
9.       z++
   9 o( z, o2 G3 S; E$ @\" q7 a
10.     },. D; Z, a. x3 ^' z
    
11.     y++) [/ Q  s. k7 S+ ^$ O( G
    
12.   },1 a5 {; V7 S# j7 L
    
13.   x++
    2 z) o  Z, n1 Y! I& E
14. },
    $ }$ y2 g; h/ d3 Z& V5 `; x8 k
15. s;  s: o/ R( C3 i) ~- A
    
16. [sys::clock()-t]/1000;

复制代码

结果：
% h; B" N. O6 ]! Y( z9 c) F6126720600.
1.109  秒

2 w: Y0 u8 O7 i1 A. }& U8 |二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

2 q2 w  A2 [1 w2 D7 u& p2 a


1 w/ F0 v9 b! c# _0 ?# T" c( B
9 h% u, a. Y: U' T
( z7 t8 {% e2 j

6 N; V# t( H+ N3 t

; l5 R4 U" d$ l3 w" s
! D% ?( j/ q% ~( F9 C' \, o' R8 o

% L+ K8 J2 u9 f; _

1 j7 q2 y: `$ ?3 B- R