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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
( W- s8 U% i/ a( T; n( ^0 U
- q) C4 L, T  ^8 n! R- \; f＝＝＝＝＝＝＝＝＝＝＝＝

% l/ P/ _) _& d0 {' EMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;6 f  g  b( u/ B2 k) x
   
2. tic;
   ( s. K- e+ G8 @5 O6 j
3. s=0;
   3 `  K4 `: R$ ~/ F: H3 P, E
4. for x=0:1000
   % s' r; D2 @8 H& f
5.   for y=0:100
   # H# _% f( k; Q$ Z+ m4 @
6.     for z=0:100- d1 |- d1 g/ w. ~0 j
   
7.       s=s+f(x,y,z);2 `5 N1 W0 s# u+ I\" @3 x\" R
   
8.     end5 K( |4 }( p) `1 S5 X
   
9.   end! y7 [1 g, \% ^6 q. Z3 x
   
10. end
    9 E- q\" n1 o, u* f. ]% _
11. s
    + z$ ~( O1 O, d
12. toc
      k+ z( B  P. T% |6 F7 z* l2 M
13. 0 `/ t\" z; R/ x. r0 U9 m  u
    
14. s =# v3 j+ p. L6 e1 u  b. k6 y
    
15.     6.126720600000000e+009
    3 ?2 Y4 ]2 m; q. B! n; U1 T' Z
16. 
    . M+ T4 U: H6 O5 X: B' i
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
# a8 _5 ^2 C* V# a# U, i+ k% r% Ym函数：

1. %file xyz.m
   2 @0 X& v& ^\" M9 o
2. function c=xyz(x,y,z)4 m# Y( a/ b6 P7 u  P3 n, Y) u
   
3. c=x+y+z;
   1 F) e( H& Z  g7 E0 i
4. end

复制代码

测试代码：

1. tic;
   \" M/ G  d& F0 l0 c
2. s=0;
   ' N  ?! M1 r' z+ }
3. for x=0:1000
   9 R- j  V9 Q  k& p& y
4.   for y=0:100
   8 u7 X4 R4 b1 Y% R7 g
5.     for z=0:100
   ' h6 c* A$ D  t/ [! m
6.       s=s+xyz(x,y,z);5 ~' }. H+ q7 H1 W& P
   
7.     end- s/ a* N3 U0 Y# I) `1 v) ?/ ?: g; i, J
   
8.   end+ r/ U\" r8 A1 ^8 S1 w( n
   
9. end
   + i- A5 l6 k3 @
10. s+ }( `. d+ I( U/ U$ }. N
    
11. toc. A' H# n  g( ]. `6 n
    
12. 4 y% d1 i/ K, l\" I5 s
    
13. s =' `& }& V8 I: P; D; {+ l/ x
    
14.     6.126720600000000e+009! C4 x+ ~: S9 u6 f
    
15. , w3 K, K2 c* Q+ g4 k) S# @! J1 G. R
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   0 K; N. d1 m% j: i9 E# b
2. mvar:
   6 e5 |- f% ], C  s& K5 A0 B5 l# {
3. t=sys::clock(),
   \" U, O0 @) s! p5 Z+ J
4. s=0,: z\" _& T$ k9 s) l
   
5. x=0, while{x<=1000,9 V6 m& K. Z\" K% M
   
6.   y=0, while{y<=100,
   + s9 ]9 t8 L. g/ X- M0 [6 ]
7.     z=0, while{z<=100,. F. r8 J6 ^3 g9 K2 B* U# P
   
8.       s=s+f(x,y,z),7 |& g' {% j9 C6 V
   
9.       z+++ j% L; O! V& }) r) e. M
   
10.     },
    ; }; |! O% w\" N/ W
11.     y++
    2 B1 B3 c7 Q* `8 f% ?: B
12.   },1 y* A: S0 P* l; S4 g\" K/ j8 D9 S4 Q
    
13.   x++
    - w' ~* E. P6 o7 B9 l9 S# ]7 V! M' U6 d& s
14. },
    9 @6 a7 H, t5 w1 w& Q0 t
15. s;* {+ _  [' y7 f' [9 w! d; B2 Y
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

( r% ^/ N4 C5 R! U  A" _二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

1 c: T" u" J! y


/ o. d; m7 K% D, l% ?% C
& g# B) c; X% h& o9 F) t4 s
( [6 \! C1 x& i% X; Y
, Y9 C. Z: B* y1 v% t0 F/ }% x



$ |9 C9 k* @* F0 c- c$ p/ ?( \4 s



# Z$ e8 m' `4 k$ c) ~
" n2 t/ n5 \6 W6 c* Y2 {' M