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

forcal

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

% a, V% Q/ u0 r# T9 \3 X) a＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;8 L  [: w* w0 F7 N& B\" O: A
   
2. tic;
   : b5 H  o7 m/ N8 s3 a9 \. p
3. s=0;
   0 @! W. f3 [* Y$ p. ]
4. for x=0:1000, e* p2 V, A% G7 m2 V7 E7 N
   
5.   for y=0:100
   $ Q; X# I7 h+ \0 q. e0 x
6.     for z=0:1009 m4 U+ x. C8 j
   
7.       s=s+f(x,y,z);
   - P% M! y1 J( j$ g4 H1 u
8.     end0 L& m1 q# }( ~, P1 i2 u
   
9.   end* [# n9 b) T. d
   
10. end, F7 g4 Z8 m7 p6 N7 l# l
    
11. s5 S3 Q* a) Y! r( T9 n- c9 `7 A! k
    
12. toc
    \" f# Z( @\" k4 ?! ~
13. & L* V' P. q1 M. \
    
14. s =  l$ ]- B5 w; I% a
    
15.     6.126720600000000e+009& O# C% E) o- l+ k' h& p! H
    
16. 
    ) T/ V1 D! ~. e0 \) ]
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
4 Y1 ?" U; R6 M3 j# ]m函数：

1. %file xyz.m% p% u1 N: H- S4 i
   
2. function c=xyz(x,y,z)
   0 K* Y- E; E! B- s
3. c=x+y+z;
   9 N- x3 w6 C- F0 G% U, j9 D
4. end

复制代码

测试代码：

1. tic;( S! }2 {. Y$ ^\" d
   
2. s=0;* j1 e! J\" F! R  O( s. Y
   
3. for x=0:1000' W4 ?  `2 a2 s; f: Y  w/ \. m
   
4.   for y=0:100. i; ?5 M+ K2 D$ g1 p7 a
   
5.     for z=0:100
   4 A/ l9 e  E7 y8 W+ B
6.       s=s+xyz(x,y,z);9 x0 A6 Z1 E+ o, @0 x9 Z
   
7.     end
   & j! S0 R) P: b# p( B$ o' Y4 K# x
8.   end
   ! A. x7 A2 E0 d. I2 L
9. end
   ; Q& ]2 N9 V0 f4 i* z! n
10. s2 q\" Q! s( F1 O+ {# L' Z8 t$ B( E
    
11. toc
    ! \; B( ~9 Y- O  T1 d2 V( J
12. # x/ k. J. l* c3 \! b
    
13. s =
    / p/ o  K\" V\" U( }
14.     6.126720600000000e+009
    3 h# W1 J7 n% \/ _( Y
15. 4 P, y( u, D1 P
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
8 N  q4 @( U: H* N3 I- a9 c9 D1 ]) P* s  L
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;, }! C- D8 _, ?! m* l8 F5 k8 P
   
2. mvar:# s2 v\" h\" Y% _: B6 \; o4 I
   
3. t=sys::clock(),
   & V3 @, {+ F3 u0 e$ R\" L6 Z
4. s=0,1 k( ]: ?% T& p+ E8 K' y# V6 C9 S
   
5. x=0, while{x<=1000,
   + [9 ~# ~\" g( b3 G  \0 k& ?1 x5 M
6.   y=0, while{y<=100,
   - |+ e, N\" x2 ~) A( y* z0 a, h
7.     z=0, while{z<=100,
   6 m8 z+ f, ~' L/ s; @
8.       s=s+f(x,y,z),: ^\" ]4 n3 B( D) F! f\" P
   
9.       z++4 a3 @% y- w\" L\" m9 [( x4 z
   
10.     },
    + V8 y  p/ E* J8 _  }9 S\" u
11.     y++
    : ~; ~9 m. a* i4 B
12.   },4 `% X# j2 x. S
    
13.   x++2 |: `, S) N7 u* x
    
14. },
    $ ]! R; ]0 {$ S& G
15. s;
    : ^. n& A, k4 ]* A
16. [sys::clock()-t]/1000;

复制代码

结果：
+ {" `/ D( u# x9 ?6126720600.
. l. N4 C# D+ C1.109  秒
0 X+ [9 L- b$ E
3 Y) T( T6 N# G3 d3 g! {二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

0 M: i0 O& ]4 C" w) t/ d



' g- a! w% R8 q. t# T) V
' `# {" u+ G0 v+ d
% x5 u2 R/ {) a
( _' I" R/ c3 [- i2 }

3 |; A) `9 c% G


% d# _9 d- K6 A, L, H: ?  L" j
# B& x. [3 o' V; m2 N) \& I; ]% |3 `