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

forcal

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

# H' d% O0 p( A% p＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;. n4 w6 y7 c( _' f
   
2. tic;! n4 D8 o7 P1 q+ S
   
3. s=0;
   2 ~4 b8 ?9 F  e
4. for x=0:1000& T8 ]! O! t4 k4 c
   
5.   for y=0:1008 ?6 [- x( }\" ]4 k0 T
   
6.     for z=0:100! @: d6 {3 o3 Y) y( b! I. P
   
7.       s=s+f(x,y,z);
   $ X) u; r: B: m% G& R( j
8.     end; g! }) [9 h+ y/ q& t8 T/ c
   
9.   end$ Z3 `6 ~/ e  r  F' s% k
   
10. end
    : @3 ^$ }/ @& U  w8 C( E. i
11. s
    0 J$ C- L9 z$ `9 p
12. toc
    + N+ m  Y/ v8 S' ^/ s
13. 
      W1 _1 i2 d! U) ~) o
14. s =! {- _9 x; k& P( d\" c1 Z7 _
    
15.     6.126720600000000e+009  W) T2 n# {7 t3 e( l' J8 n
    
16. 
    - O! [& \* d1 g7 ^
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m6 |4 S: y& [& U; \- y# d/ r
   
2. function c=xyz(x,y,z)) p& J6 I$ ]2 n0 ^
   
3. c=x+y+z;6 t; D' g! U+ r7 @6 t! L
   
4. end

复制代码

测试代码：

1. tic;2 F& u( Q8 j, S# \
   
2. s=0;5 T- F: R9 }- p$ C  m\" B) @( F
   
3. for x=0:10009 F; Z5 Q4 M9 `
   
4.   for y=0:100% V% w+ j6 T4 s' I# ]- f
   
5.     for z=0:1000 P4 N4 C\" k, q' n
   
6.       s=s+xyz(x,y,z);
   & k+ v- e* u0 I
7.     end
   % @. Q) g) A  A( j) Y  ?$ ?. B
8.   end
   . U, |+ q4 K) B4 o
9. end+ r7 b- F* q- q/ s* {; H8 c- k
   
10. s
    % r/ F. h$ p! @! ]5 @3 {8 h1 }
11. toc
    / h) g8 [8 e% o' |1 B2 t) [' `
12. 1 a, ]\" y; b- d
    
13. s =; |( p3 G) H3 G
    
14.     6.126720600000000e+009
    / n1 E1 {  G! V
15. / ?6 T9 i& @1 G# G# Z
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
4 X  w  f# \  Q. A$ }+ I4 _
7 v1 s( D/ i$ r" _' ~6 }$ LForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   % m7 i/ C8 U9 s, d
2. mvar:
   0 e- K# e: O) l! E0 _
3. t=sys::clock(),
   / _3 ~* `- g# ]\" B6 g
4. s=0,
   0 e9 R% N6 \& E2 s8 Z. }
5. x=0, while{x<=1000,% q7 x& |. f- I4 Y6 P
   
6.   y=0, while{y<=100,; A& W' `# l; A  R3 V* \
   
7.     z=0, while{z<=100,
   \" c5 z; G3 H9 D0 A! ]0 u
8.       s=s+f(x,y,z),
   ! y7 a3 m( m; ^\" b* r' ~
9.       z++4 I# m2 A8 A\" q\" S1 d0 z
   
10.     },; H7 Y6 _\" i  W
    
11.     y++
    ( s$ }6 v9 {0 R: R
12.   },
    \" L8 d& X4 m+ S+ s* O: n6 c' T
13.   x++
    : I* y% h1 J$ c/ V6 l
14. },
    ( t/ M& R) @; r. o
15. s;; |5 P6 `( V# y6 ?
    
16. [sys::clock()-t]/1000;

复制代码

结果：
# J+ h! k' I2 s  s, k7 Z- R6126720600.
1.109  秒

- v3 k! I, F' N& L0 t: A  o0 T二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
! [, X" Q/ C# u

hxftntA

( t$ R! ^0 R5 N/ w8 H; e3 _
4 U/ ?9 n" u7 d+ Y" ]& ]
' i0 A6 ]9 K3 a2 a+ f
$ G% I, L8 T. W, e

' u* @4 z6 d. {, A5 n5 l6 M
, G% V" O: k5 k6 _! P' F; n

: o) T/ u9 W5 u! T) w2 j4 t% ^# B: N
: g+ [7 e! H) l


* s+ ^& c+ d' T% ^3 T
# s6 y! Y- a( u( k

$ c  I! H6 h8 w