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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
- H0 b$ k" w- F6 n4 [# K5 D
/ d, Z, O* V0 U+ M: d  p9 l6 A# Z, S＝＝＝＝＝＝＝＝＝＝＝＝
% E: S* T1 F' i$ {* Q: I' V
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;0 B' ]0 @- O2 {
   
2. tic;
   # y( u- y' v! `  }/ z2 m+ Z+ s
3. s=0;
   , J6 V' v  L0 i1 {6 O' T9 H
4. for x=0:1000
   7 N( _# y! F* N% ~4 t% Y6 n7 W
5.   for y=0:100
   $ n1 s- M2 ~! E
6.     for z=0:100* r. g2 f  [8 P0 P% S' u: [4 p8 H3 `* Z! V
   
7.       s=s+f(x,y,z);5 W* P+ ^\" ^( B8 [8 D+ l0 x
   
8.     end
   ! k/ {& k) Q\" h/ y9 U9 W
9.   end6 a* N, A/ Y' K! g) t9 X
   
10. end! E! B% L6 q/ u\" ?; @. }/ x
    
11. s2 i( C+ Q. ?1 k0 l8 r) F2 ^
    
12. toc+ L4 y3 e2 i+ _& h1 J- c\" J
    
13. ( F; \/ @  _' I% I
    
14. s =. i7 l- F2 `; h! i8 B9 u' f& m
    
15.     6.126720600000000e+009/ A$ w- c! [4 o
    
16. ! [9 C7 x$ U- ~$ A- J
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
% ^- F9 W4 L" A, {  V4 ^7 p9 {m函数：

1. %file xyz.m
   & d/ o8 |$ w& ]( u
2. function c=xyz(x,y,z)* s/ E6 |4 m( m& j% o  F+ Y0 ~; d! ^
   
3. c=x+y+z;
     J% u8 v5 `3 A1 C  m
4. end

复制代码

测试代码：

1. tic;
   9 x6 c: H7 g% K
2. s=0;; @. L) B: u9 b( J4 }' ]
   
3. for x=0:1000
   + v- o7 M8 X. F! k) `9 P8 N
4.   for y=0:100) b! D4 p# S3 A- v+ w$ N
   
5.     for z=0:100+ s$ ?8 F  _- H. E7 I
   
6.       s=s+xyz(x,y,z);% g6 x4 h( |+ J8 x9 R: c* j
   
7.     end% m. z: \; C/ ^4 P' f1 I
   
8.   end
   9 S2 w9 p/ p4 ^1 U. O$ h2 i1 Q1 C( S; O
9. end
   : S% I. j6 Z+ d. k' j& K
10. s\" F* W* L6 R4 ?2 G1 i\" @+ d
    
11. toc2 f# R- z; v+ G, j
    
12. 
    % O# t+ [, o. S) a, f
13. s =! Q, Q: l6 l, `+ A; \. m& k
    
14.     6.126720600000000e+009
    9 ~# }( R% Q, ^, I2 M\" x
15. - z2 I4 ?5 x9 E4 S! e
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
$ |" r# u3 F% v
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;\" c# k/ V6 E3 G- d! W! X
   
2. mvar:
   4 K8 w) C+ E' o8 ^' t( M) j3 _
3. t=sys::clock(),
   4 z' F8 s, |' G+ c+ r
4. s=0,. Z' _2 }& E: P. v6 e! K
   
5. x=0, while{x<=1000,* O( r% ]! j\" n
   
6.   y=0, while{y<=100,3 ]; J% X( K; N: \) [( R
   
7.     z=0, while{z<=100,2 Y+ V( l) e1 Y( Q& ~' E3 x) h. D% }6 Q
   
8.       s=s+f(x,y,z),/ o# m+ F\" G$ r) N! p( N
   
9.       z++
   7 U. P% i% N' q% H3 R- L
10.     },
    , C. E$ b% A4 [
11.     y++* h% Q+ |( I( J7 e
    
12.   },. _3 q' ~( b1 e: b( o
    
13.   x++, r0 D7 e+ ]2 q% B; M; `8 i
    
14. },+ v4 |% x( M3 A3 W# w5 m
    
15. s;# F0 g# z3 x) A/ J2 q  t
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
) r: }0 w, o) A: a- M$ P1.109  秒

# G+ ?# K' S. t+ q( }2 [" Y  \5 [二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

) z2 v* H! J* A# n; \
8 F. X+ A  p* _1 W. B3 C' ~0 W
7 G# k( m' P" ^* z4 M4 p
7 r& m* B# i) y" U. f

1 X. {0 J- u) r. |

  G7 X" P( m/ G9 X! K6 T
0 H' Y9 h3 r! \1 j& E1 C

- o4 m5 p3 \. c& `


) v; a  G7 D8 q" h
$ j! D1 Z% c& y0 ]+ t% ?
3 d/ Z( B6 a* [, L