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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
# e8 y5 L% J% s' X& E/ k' b
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;4 n; B7 Y5 q3 n6 f
   
2. tic;1 K4 G& c5 T0 H; R/ C* ^# x! K
   
3. s=0;
   ! q: [# i) C, Z1 T# U# c* \0 C
4. for x=0:10002 M- t& _6 @- b7 w! b* u
   
5.   for y=0:100
   9 M$ L+ f7 K0 m\" D* |) k; q
6.     for z=0:1004 h6 M8 P9 t: g5 b9 K
   
7.       s=s+f(x,y,z);( F9 \' [) ~, c5 X4 X- _) K
   
8.     end
   ! ]3 Y5 y; x+ m7 t  n0 N1 o. p
9.   end. s, [\" k% |, p5 S\" U
   
10. end
    8 v2 k% t8 D2 d
11. s9 X9 s3 a1 B+ ]- O  M; M3 ?
    
12. toc# O  X+ K0 Z+ m5 k9 x+ I' R
    
13. \" `2 o( m+ d% p8 U2 V$ [3 `
    
14. s =
    . Q* f7 B* f* E9 f/ P
15.     6.126720600000000e+009
      R8 V/ H  P9 c5 r* e& U, p: Z: M& N
16. 0 W% D5 {2 s! G  \1 l! I  B
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
( i; h" J  K" E+ Dm函数：

1. %file xyz.m
     ?+ {' P  q- R& d- H5 t
2. function c=xyz(x,y,z)
   7 X8 i\" q8 C- ]' Z$ K
3. c=x+y+z;
   ; S\" Q  a! ^9 `5 I' P1 u7 j9 {
4. end

复制代码

测试代码：

1. tic;
   ; M4 b  {  B5 S& ^/ E5 h! Z
2. s=0;8 m8 d3 H: l. t( d- F
   
3. for x=0:1000% e- Y, y. c/ S+ n$ X
   
4.   for y=0:100. V5 ^0 M7 L6 g! `
   
5.     for z=0:100
   9 K( m; G) ?4 F/ Y! u; i& @: Z
6.       s=s+xyz(x,y,z);
   * b4 b  n$ [& U! H1 X3 x
7.     end+ o$ y: i/ H- v& B+ Y: }
   
8.   end) ?$ E* p, s+ Y
   
9. end
   4 H' g' C/ k* e3 M) v
10. s5 Q' C' F# z+ N' K
    
11. toc3 Y. h1 N2 p$ |\" ~) e; I% L7 S
    
12. 1 g# G' z7 H! `! b% L2 W
    
13. s =
    & b, X! T& K, V1 D
14.     6.126720600000000e+009- u9 P& W2 \7 U) Q' d
    
15. 7 B( B9 _; J! F4 g; K: f2 Q/ k/ {
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   1 ?2 n+ e; [: n0 X5 y
2. mvar:
   & z; J# X& o3 J
3. t=sys::clock(),
   . Q, c! ]\" m* m5 w) y6 S4 _1 _$ ]
4. s=0,4 y# G4 ^\" s\" P6 e0 K5 D( d/ n+ E, p$ x
   
5. x=0, while{x<=1000,\" M( f. ^* L6 K( Q4 s
   
6.   y=0, while{y<=100,
   + w! R& V) [8 m4 u+ d& i& I' @) O
7.     z=0, while{z<=100,
   9 c' l' X, r$ X\" d( u/ R$ t3 F; d5 {
8.       s=s+f(x,y,z),
   . ^' p; K- x% [& y: @
9.       z++# `# \. E: Z9 |& U  m/ E4 E( f
   
10.     },
    3 W$ @\" N4 J: T0 L: W4 x) j
11.     y++
    1 [, U) \3 l5 X! F7 x0 @
12.   },
    \" Q/ s& x/ b, G  x0 Z3 h  f
13.   x++3 d  h4 F' [- f/ `* B
    
14. },
      d5 u7 l. Q7 i: a
15. s;. B; k/ J- M8 @( q/ G1 K7 V\" p( t- B
    
16. [sys::clock()-t]/1000;

复制代码

结果：
3 h" p. h9 @: H8 h" ?* G6126720600.
- w3 `& Q1 m9 E& U1.109  秒
/ y: i# q" z  @* Q, N& ^, R* D
' ^5 }9 @% ^: Q+ v/ p/ q二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

& s, q0 x% I0 s1 u

$ H( B: i8 _1 a% m% W  G. H4 r


/ o. P8 {' i# U) k
! q& M* g% u: A! |2 V  V* R
9 @: l" Y$ O7 x8 Y) r7 B


/ n2 I- m' U8 X' M: F
+ _1 @% g# d8 V5 Q; C( x
: x8 Z% h$ l) `9 g# [+ M! i+ L5 ^
9 k7 }' a8 ?/ F* E* Q7 e
' Q2 j" C* ?4 @1 p! S9 H3 I1 ?
% r$ }3 s/ |8 v9 e3 H: f$ ~