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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
  n  x( T0 }" o- Q6 L
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ; g( r: |2 P& F7 ]% V4 _: n4 E
2. tic;! q% i0 ]% b1 W* Z
   
3. s=0;7 x( h% }7 y- |
   
4. for x=0:1000
   # a8 Q: f7 G  g( s+ y% r- k. \. E
5.   for y=0:100
   & d9 V/ ]* r; ~% _
6.     for z=0:100' \! C: V* V. b; C7 b4 g4 Q( F
   
7.       s=s+f(x,y,z);( q% A7 w6 M9 s! q3 u0 W  T! Y5 U& w
   
8.     end; m4 `7 W5 e6 \/ N
   
9.   end' q! |. R) Y& n# u4 ]/ q7 Y- b: _
   
10. end  l* q1 E+ P. P- u\" J
    
11. s1 u) f2 j# C; V! W) u
    
12. toc: ?- O5 B: Y1 W) y3 U2 @6 z
    
13. 5 W; A8 t4 O0 J: {: s
    
14. s =
      X1 Q* T) x( w  N5 l3 W
15.     6.126720600000000e+009
    ) A# V2 J/ w: r) v) Y
16. 
    * ^$ H# m) I: i4 M3 Q/ e6 g0 [
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
7 [! ^! w) i" x2 ^m函数：

1. %file xyz.m6 l  i* E( e. R' c0 r9 d4 J
   
2. function c=xyz(x,y,z)
   1 M) R' ]+ ?* F- {
3. c=x+y+z;0 e) s) Q8 W& Y, f  T! ?/ L) ]7 L+ P
   
4. end

复制代码

测试代码：

1. tic;
   9 X2 r% o! B; m& \
2. s=0;
   5 t  I* }; C7 y' m% @( u
3. for x=0:1000
   ! ]; R' ^( s% V+ Y1 F( K
4.   for y=0:100
   5 o. i# R8 u7 J! j9 y
5.     for z=0:100\" v\" }  o& w! l: @& C
   
6.       s=s+xyz(x,y,z);
   ' z3 m6 I! m9 l  W. [
7.     end4 k$ S, |4 F% [9 k$ D\" y\" {6 i! D, n5 q
   
8.   end2 S  P! S5 c( B8 o/ z
   
9. end5 ?, i) \2 d0 m1 u3 B/ a
   
10. s; I7 U# l% E6 e9 |- D1 K( S
    
11. toc
    % p& C+ J/ X) \: O& c% ^
12. 
      i. A, H1 l! d8 B
13. s =) {- t\" a  i2 }  x/ F
    
14.     6.126720600000000e+009  b: d/ x8 K- S) Q7 Z  L
    
15. # B7 _, t) p& `5 {) f. U+ J
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
1 d, A6 S, p! ^* X
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   6 O6 x$ r' g/ Z# _, r
2. mvar:; d5 s\" j: S* X! o$ e1 l
   
3. t=sys::clock(),) {  i\" q8 P+ S\" N
   
4. s=0,! \1 d: [) Y; K& A, r
   
5. x=0, while{x<=1000,
   7 U! ?( }; g1 |+ q$ o
6.   y=0, while{y<=100,
     @3 U\" z  }+ z
7.     z=0, while{z<=100,7 i5 f: m; g$ B
   
8.       s=s+f(x,y,z),
   ! h& {! d5 Z. {6 B4 W# @3 i
9.       z++6 W' \# u! P- p: {
   
10.     },* t$ _9 u5 `9 g6 u# i' z- x. i
    
11.     y++
    ! I& S/ }% A8 f  m% K2 h1 Q
12.   },# O, q5 q% E  M* D) k- J
    
13.   x++
    & ]; D7 O1 c/ T+ i\" V* y) T: Q' C
14. },
    ) T! |5 b9 D$ Z2 I6 P$ a1 Z. J
15. s;9 \* a5 ?$ Q' l\" O. {' Y
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
! y7 M# q1 N  Q3 Q, t1.109  秒
6 e6 C: H1 s. w5 h3 w( T' k
7 Y. A( U& P' H" y" z( m二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
% z5 ~8 d5 L7 X, @# n

hxftntA

/ Y  ]+ z' f6 F6 r* t, |
* V- t1 i, F* ?# G7 ?

* f+ w; n' c/ H
1 U) o  ^. i* x

- q' k2 R' D; O3 H4 S( U" t
( V" j- s) N9 V. Q# m% ]
2 b1 u1 z4 z6 Z$ Q, i

' `' o  B) g- f2 w  S; e3 K( A


; T# a" Q% ]' q! g* O

/ x* x' a) n: \, Y