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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
' f& ?& ]7 ?( e* w: j5 k/ ~
＝＝＝＝＝＝＝＝＝＝＝＝
% ]) G9 r! m2 @8 J  p
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   5 C* y1 Y/ A. h3 L* H. r: Q. n
2. tic;7 K9 T9 z; ]8 H- d\" u. Z\" T- P. d
   
3. s=0;\" E' i* O  l7 g0 E) z/ ?& M\" _
   
4. for x=0:1000
   * h- e6 Y  L# s& {& w+ c* H# i2 C
5.   for y=0:100
   % `8 F& y4 E6 d, K; d
6.     for z=0:100; j) b# Q0 x. z, C
   
7.       s=s+f(x,y,z);
   - E8 n5 C6 p+ O9 l6 E, l) B2 u
8.     end
   & e) y$ i' Y( k+ g
9.   end; c, `' e& x5 S
   
10. end
    6 S9 t. [; \+ N6 u\" `1 \4 A
11. s\" G) d  X\" i/ B; q! X& x+ \
    
12. toc
    # O# `3 R; P0 f0 u: L( m/ N
13. 2 h/ b$ N- M5 X1 v3 a7 m\" O* d
    
14. s =
    4 Y\" b+ U2 F2 r7 W
15.     6.126720600000000e+0095 W: v9 Y3 W9 O! |. }# g
    
16. 8 w6 K) L) O- r8 y
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   + E/ A8 X0 c- v' j\" w! u; z
2. function c=xyz(x,y,z)
   # [6 \# Q: V1 F: \0 j
3. c=x+y+z;( U7 A# p2 F! t* O
   
4. end

复制代码

测试代码：

1. tic;
   \" b* s. p% v: u, u* G# J, _
2. s=0;
   $ p% ?6 S( ?: `: ?9 l4 O
3. for x=0:10009 U) ?+ j* E0 N4 d
   
4.   for y=0:100: S: S' L- P! }4 x
   
5.     for z=0:100
   3 r& U# ~/ d7 `  }
6.       s=s+xyz(x,y,z);# P5 Z1 h0 o  s
   
7.     end2 E( {: }7 d! i' f6 ^5 A
   
8.   end
   / K0 Y2 V% W- ]$ C% H7 y
9. end
     N* B  }/ A\" {\" F- P, u7 u
10. s# ?- f. O3 d( C- g
    
11. toc
    % V. T: A  }# t5 T/ _% A
12. ) Z9 F/ w* y- x% ?* c) Z6 M
    
13. s =2 P* ~& k2 i  C6 o7 U
    
14.     6.126720600000000e+009
    3 H( Q# p/ a) `- L( f
15. 1 w/ B$ h3 [) E2 Z  Z\" z: x
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
8 O! U; y/ ?* l- Y  s* E" _
+ [9 }+ J' r& R! ?5 \8 u* x; sForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;% B( b( o. T0 K
   
2. mvar:( q! l, Y: H/ ~' S4 Z! f' J
   
3. t=sys::clock(),
   $ S( P6 M; l+ Z- i2 ?. g: a, Z
4. s=0,, ^6 l7 M\" \, }/ V. ?
   
5. x=0, while{x<=1000,* p9 @8 h* |7 l
   
6.   y=0, while{y<=100,  K# g- v  g7 a- \1 h2 N
   
7.     z=0, while{z<=100,
     _; S, y3 o) k5 N. Q4 m\" p
8.       s=s+f(x,y,z),
   \" N/ w; d# F4 C2 \8 d
9.       z++
   ( w8 h5 D5 L% ^
10.     },
    0 I. D! l$ a; ]& W! S, y\" i
11.     y++
    ) E0 w1 q+ [  o$ P* d
12.   },. b  ]5 u6 P+ d, e\" O7 r
    
13.   x++
    , l5 Y$ b' Q* w
14. },) t- W\" I: s9 k0 }/ r0 C* p
    
15. s;6 t2 ~% q- e/ a- E
    
16. [sys::clock()-t]/1000;

复制代码

结果：
8 A5 q* i* j9 a9 k6126720600.
1.109  秒

: t3 y7 O5 u: u& B& f) x  j0 F二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
1 ], L: L8 R6 b7 l% m

hxftntA

- b3 G0 t5 _5 a
$ ~7 H) o/ K" j7 i! s
7 ?  e6 t8 b: }



6 }' }# _* c2 z* ~0 f! b
9 O2 f) l3 w, ~+ E: R) S! F9 k

4 o( ^) B2 l8 ]
" s3 ?1 w! w( ]1 o; B9 c