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

forcal

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

$ i, ?6 i# X% |6 C- E6 `: n8 b2 d＝＝＝＝＝＝＝＝＝＝＝＝
  @/ w" ?# L0 O
9 A+ S# B- s/ X( vMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;% @& I+ Y, `0 ]! T# Y
   
2. tic;
   5 u3 {) i3 a3 q1 E
3. s=0;' R, b. Z; q8 A0 @
   
4. for x=0:1000% d5 g0 z2 ~8 n\" P4 U
   
5.   for y=0:100
   7 O* c; d3 L0 K& z
6.     for z=0:1009 [9 p$ Q+ _8 _; O. O3 h. f8 N: n
   
7.       s=s+f(x,y,z);
   ( m1 j) O9 A* h, X5 G3 @
8.     end
   ; I: G( B8 N  U
9.   end% }! K7 j. Y* S& y
   
10. end$ j. R4 F# y, B* a% ]
    
11. s3 W) U% }$ S( ~4 B/ k# H$ E
    
12. toc
    8 ~; H. F\" v2 D) D5 E+ g
13. 
    : S( q1 |# `7 m\" u) i6 M
14. s =
    4 X( P8 A+ E6 Y0 r4 ^6 d\" J/ V+ Z
15.     6.126720600000000e+009
    3 z& i2 }7 H& _
16. 
    8 k  o( ~8 n% W& w5 Q& J
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
4 {- I1 f1 `+ g; X# {m函数：

1. %file xyz.m( H1 J& M% L( ^6 ^
   
2. function c=xyz(x,y,z)
   0 q! g4 E' l/ q
3. c=x+y+z;
   # `  X; I' {7 E, P2 U) u) M
4. end

复制代码

测试代码：

1. tic;
   4 K& L5 P, j, k: ?/ `
2. s=0;: i, H\" D& o( w\" a
   
3. for x=0:1000
   ) J& ~  D\" q% \' O; ?
4.   for y=0:100
     c. W# S% L) V5 P# d1 P* `
5.     for z=0:100' y5 @7 s: T& ^$ A4 P% R/ U  ~
   
6.       s=s+xyz(x,y,z);
   ( y! q9 a. ^0 l0 d
7.     end
   . h4 |7 L; y: D+ \# d' b. B6 W
8.   end
   : f0 P\" n, _% @/ V! D0 Q
9. end' c7 H, z3 x! b, c+ J. R  S' j0 c
   
10. s
    ; f) I4 ^* u' n
11. toc. Q* l: O# ?( d# J3 S* ~& Y
    
12. 
    ; b* ~. Q; m0 p! i) m! ~5 H9 Y6 b
13. s =
    9 R7 Y& A4 c% s5 s; G2 S1 e
14.     6.126720600000000e+009
    ! T6 ]9 k) n\" ?# d) g  t  r
15. 4 c6 D: c: V8 c
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

" e1 q' M, @7 O# A# T( h+ n' [Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   4 k+ B1 y, R6 E4 A' G
2. mvar:  z$ O5 d) a& w% [4 d4 l9 P$ V2 ]! L. R
   
3. t=sys::clock(),; T0 J( f2 ^$ G7 A- V1 ^
   
4. s=0,
   ' k5 g9 p4 c6 n\" ^9 m! o+ Q6 d
5. x=0, while{x<=1000,1 _; M( v1 ?7 F# h
   
6.   y=0, while{y<=100,
   6 l2 Y$ N0 c# Y' C: Q- \5 F: Y
7.     z=0, while{z<=100,
   : R7 z& p& p' L6 o
8.       s=s+f(x,y,z),3 g' Q2 W/ i3 x: |, I+ _( j
   
9.       z++0 }9 Y. C) v\" ]\" Y% ~) f; E' g
   
10.     },2 q; U6 o3 F8 C
    
11.     y++5 a, L% n6 X% D6 m& J
    
12.   },
    6 N5 K5 I1 r\" n. l
13.   x++
    4 ?9 J$ x- i% Q
14. },3 [4 a7 T1 g3 O\" d, r
    
15. s;- X. `) s( B3 y' B4 x
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
2 ~( T2 Y2 y: N6 Z$ Z1.109  秒

, H9 v0 I, L- f/ Y) |8 u; X二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

( A. |8 J- E# b; B1 |5 x  {


$ v5 i" m2 ~8 U" C1 Y

$ k5 ?: e1 h( Y: E9 Y
- U8 m) X' C5 A. J
: s3 D1 X, Y- z) Q9 Z
: b; w( ^# o( F1 a

9 q# A' l# o( t% ~' ?; e8 e- l

, ]7 Y2 K% N1 `

# C  n9 [  N, Y/ U  ?% t