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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
' G; d* q3 v9 O# w
7 b3 T3 a/ P; @9 j. h1 \) g＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   + {4 t( c/ V. z$ B$ v  [0 d
2. tic;- K7 ?; f# c1 E6 I
   
3. s=0;
   $ d3 a# G; O; `% N
4. for x=0:1000
   % S  H5 b9 \0 Y* N: S/ }
5.   for y=0:100
   : P/ t: R8 b7 c! V  \7 X8 |
6.     for z=0:100
   $ X9 P  C$ m- a) g- y& n0 _
7.       s=s+f(x,y,z);
   # Q: Z) T. ]5 S' S
8.     end3 V( H6 i3 C, K) z0 b( e5 f
   
9.   end
   \" c- H! N2 e0 `! x, U: w# h
10. end
    ( v' Q- ]# ?4 l5 n& U\" x
11. s7 P* u# B9 J7 M2 \
    
12. toc
    ' f: w3 S5 B) e- [/ l9 z; `# g4 }
13. 7 s. k8 @& z2 Q4 t
    
14. s =
    0 J# m5 N/ y2 c3 V8 w/ l6 w8 \
15.     6.126720600000000e+009( D' W' U; M3 k
    
16. / A/ @/ o4 p2 b  ]0 H
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m( L% J\" i/ R- h8 i
   
2. function c=xyz(x,y,z), @+ y: H6 |  p' ~: Y: Z6 g\" f
   
3. c=x+y+z;# R: e/ C  z& T4 W2 B\" T
   
4. end

复制代码

测试代码：

1. tic;9 c' @9 T# |\" H& q7 m5 T
   
2. s=0;
   ' U( u  B3 Z5 [, E1 B# T. J
3. for x=0:1000
   9 W5 C5 E+ D$ |, U! O# l2 \! g
4.   for y=0:100
   4 m# o. J9 [: W# U# [! `\" d- T% G; t
5.     for z=0:100
   $ \0 S( @6 Y' Q9 Z5 w
6.       s=s+xyz(x,y,z);; M% C9 k& v1 W7 @) ]
   
7.     end
   7 J: r5 z# \* i; K7 C* w
8.   end
   7 P0 X\" b\" ]* e' n4 R
9. end; E. Q% X- ]; z\" o! t0 |; \
   
10. s
    ) X  v% ]& Y9 s9 j
11. toc
    6 Z8 `6 K! ~# u9 u2 i7 \; `
12. ' W9 r4 h. Q3 b) o: H
    
13. s =
    7 }7 B6 E0 q$ V3 ~/ p: _* m1 M  X
14.     6.126720600000000e+009
    # u) b( ~$ {\" q2 u. p# K
15. - l4 }# k. y1 ^6 F# T3 [, e. z6 @
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

# D6 p* Y# J3 T# r  y: @9 |Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;' u% i. B$ M) l$ v+ a
   
2. mvar:
   6 J- k. l3 q# S' w. Z: y# ~2 D
3. t=sys::clock(),
   ; o; L. a\" Q' ^9 L6 |
4. s=0,
   \" n8 f! x( M+ [( |) f$ w5 g
5. x=0, while{x<=1000,
   ! A1 L4 z) ^\" A, S+ U5 q
6.   y=0, while{y<=100,
   # q7 `\" v8 g, z$ s\" C5 z- l
7.     z=0, while{z<=100,- I! p. `) b  i2 S5 x% O
   
8.       s=s+f(x,y,z),
   $ t% y- \  _! |7 D  T% f
9.       z++
   ) ~; f  _5 {: F
10.     },* E. Y  k7 [$ m% s
    
11.     y++1 y6 W2 V4 A' L' q1 x+ N) r6 ~
    
12.   },
    * o6 a! p0 i- m  r1 P2 U( \
13.   x++* P( |- u7 N. d- o0 X
    
14. },
      V2 \1 x4 Y7 {. D
15. s;1 }3 Z0 w6 E* @+ D( I) g
    
16. [sys::clock()-t]/1000;

复制代码

结果：
9 g+ Q8 F4 P' u9 d% v' r) a$ b6126720600.
1.109  秒
( t& I- c- Z" }2 ]
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
6 k, k1 x2 G5 e3 {! w

hxftntA

) S2 ^5 z5 |) P2 k
; y" o  i  ?$ K( u/ F

$ \6 X( i3 K1 V% q$ ~# d" h
% t; s/ y" m0 i: O
8 Q5 Z6 \0 [; W5 o' V4 V




! z) h5 W1 W. T5 k

8 C# g$ r6 f+ @3 U% ^6 j
7 m/ G, g% z: Q