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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
& c8 U$ y% ^! h" ]/ g
* S4 R' y( f$ T2 k: {: }& X. K＝＝＝＝＝＝＝＝＝＝＝＝
/ \7 ?' R" X% z4 Q) n6 s9 k
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   9 |- s6 K# n! L* K7 W
2. tic;
   \" [, Q8 ~\" o* k3 a
3. s=0;0 S/ ~% ]4 t9 z5 I  x. d5 |
   
4. for x=0:1000
   : U  e. q0 H+ M4 p& ]
5.   for y=0:100
   2 E' Q7 O3 z# e) w$ |& |. e
6.     for z=0:100
   * T$ P  e8 O  t( B  Z: y- R
7.       s=s+f(x,y,z);
   & _$ s9 n+ b. h
8.     end
   6 k* L+ T6 t3 Y+ q1 t; M
9.   end
   8 K- ^- G2 s/ R: t
10. end
    7 _5 Y9 {* A  i) X( K
11. s
    2 ?) `% T/ B+ j! s\" ~5 }) I2 m) O
12. toc
      q. S- ^\" L# g
13. 
    ) d' ?, v2 k6 }1 |1 n. C. t$ K$ c
14. s =
    4 |# c9 v5 e\" p; s0 @
15.     6.126720600000000e+0098 S( o) v* s\" H: T
    
16. 
    # U# n2 y* Q5 v7 O\" v1 ~. {
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   + v- ~6 l; p2 O: K
2. function c=xyz(x,y,z)* Z' i3 `$ e# ?; d
   
3. c=x+y+z;
   + E  l% G; i  ?0 O' m
4. end

复制代码

测试代码：

1. tic;
   0 O7 t! }* ^! I7 |8 J( W
2. s=0;
   # _2 c+ \0 O; ~6 w( k
3. for x=0:1000
   / h- l0 w% u2 B
4.   for y=0:100
   # J  U. X2 A2 N; b# p. Q
5.     for z=0:100; F/ c! c\" C\" z2 W6 @
   
6.       s=s+xyz(x,y,z);% \7 W. l3 k% e/ {7 O
   
7.     end# O0 i4 X7 B. R6 K
   
8.   end4 J2 l& v# v/ S, U& t5 K
   
9. end9 }6 F* w% E; |- y
   
10. s% l* u: s' f; `( W  c! V0 o( u! c# x
    
11. toc
    ! H  I\" H+ {5 E: L% r' R
12. & W  \; X: C; b5 o  ~* A- C7 ]
    
13. s =$ `9 R! l8 k+ y- R' F
    
14.     6.126720600000000e+009
    5 v* W) z3 P\" R6 \* ~8 v, E
15. 8 J) S$ }5 ?6 A4 L: @3 R, f6 @
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
% f! `  K8 P7 {/ n
  Y& F% c- Z# g4 M5 J' |% RForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   6 L* x1 m, U. V4 d, d/ G
2. mvar:' A8 P+ o- ?! p/ d9 y+ t3 c3 M; g7 k& j
   
3. t=sys::clock(),. F( x* _\" h$ n* ?6 @% D
   
4. s=0,
   / z\" t/ v  v$ x4 ?) X) k! M
5. x=0, while{x<=1000,
   ! u$ C0 M+ y\" R6 \
6.   y=0, while{y<=100,3 K3 x( a1 K( X0 a8 F
   
7.     z=0, while{z<=100,- a: a  B! j) }0 \+ U! |0 T  \' s, R4 R! s
   
8.       s=s+f(x,y,z),
   / a8 }4 M- F) x. O) T
9.       z++% Y& T3 N- B4 w5 F, Z* ], ]
   
10.     },: y+ K. Y! s6 Y5 N* _2 @' I
    
11.     y++0 _! A  z2 ]( }
    
12.   },
    3 F8 |9 ?$ B3 {! p
13.   x++
    $ a* j) r, i' s3 H% U( |2 d
14. },
    1 e: }\" s+ J! J
15. s;
    # l8 J  r$ k; `. m6 _# \0 V1 j  {
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
3 k5 [* R9 c  Q" t  K" L3 d7 X

hxftntA

% u. d5 O1 `& ]5 O$ r: h) Z, s1 @
# E. a' F' H, v2 V


3 i7 a) n1 Q( E) K2 Y0 b/ T3 I
% w( W; k+ ]5 {" k1 l

3 V+ S% H. K( D" K

4 ]0 w/ O3 N/ p# Z4 u7 ]( K1 T
/ J+ \% i$ p# t8 C5 t4 N

( Z! }! _& T6 r" n& Q
3 M  m3 a2 u$ A  z, z
0 u% g; D5 d: r# u. f6 d% F/ y