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

forcal

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

" {( Q# a3 E9 @＝＝＝＝＝＝＝＝＝＝＝＝
$ E! H! G( \5 X& S8 p
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;$ P& \1 z, q8 |# L+ k4 g
   
2. tic;+ g\" G* N. Q3 Q+ [
   
3. s=0;
   4 r. h+ }$ L) \2 O* Z. f
4. for x=0:1000
   0 N8 G# j: t  X1 X$ w
5.   for y=0:100  R1 e: z2 [7 K0 z4 J! i
   
6.     for z=0:100
   * d: H/ k6 t+ n\" w( ^4 B& C
7.       s=s+f(x,y,z);
   5 N( _  c* u% n  [- s
8.     end
   , |- ]! d$ d2 w/ F) k
9.   end+ C5 F\" b2 ]( @$ J3 n* a
   
10. end4 f0 C' @  y# m7 P1 B% n
    
11. s
    3 ]3 \) B6 j0 U) w6 Q* U
12. toc0 J) X2 ]4 w; z8 g% O: s
    
13. 
    6 u5 G) }, `9 K0 C3 k7 e
14. s =
    2 \\" w8 @2 [* ~% a& L# M( x
15.     6.126720600000000e+0091 a0 ^- ]6 C9 z' C
    
16. + v5 q- \2 W+ j\" z
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
6 o; ~/ ]9 i, fm函数：

1. %file xyz.m- Z- F5 b' d/ R5 B
   
2. function c=xyz(x,y,z)
   - }\" ~# ]# H5 e. k
3. c=x+y+z;
   0 k; b  d: u; P. l
4. end

复制代码

测试代码：

1. tic;
   9 u3 j+ k# y( m0 s: ]9 \5 L
2. s=0;9 n+ P; ^( R% t: }+ T# o
   
3. for x=0:1000) _4 y2 ?2 ~- K5 }( M
   
4.   for y=0:1008 Y+ \, N2 T- V. @) @
   
5.     for z=0:1006 o; j- [: c7 l& _9 ^& Y
   
6.       s=s+xyz(x,y,z);
   ) V+ E3 L  ?. c( e; G0 q
7.     end; T, o5 @* v/ i\" }4 ^7 d
   
8.   end
   , g7 }9 M3 F# o4 i3 v( f' q. m
9. end3 \5 G+ \4 s9 R+ z
   
10. s
    \" x: _' u+ a3 C1 Q' B3 I) H0 [5 T1 ^
11. toc
    / s% ~) I, c\" @; i0 ^\" R8 F
12. 
    , T' y0 }5 L0 U
13. s =
    6 ^; f$ Q& ^; N& l8 P2 A
14.     6.126720600000000e+0095 b+ X6 G0 ]; }/ b* E' t
    
15. $ n% U* p, i2 g1 `) z
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
. {3 A/ l6 E# E7 }  J+ _# q- w% i
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   / Y: a3 }' i$ p  h* u
2. mvar:
   5 ~$ R) s1 |0 |3 j1 Y
3. t=sys::clock(),3 o/ v5 U+ X! w
   
4. s=0,, ?; p9 v& z, f4 d) U7 ^3 d
   
5. x=0, while{x<=1000,6 M* c8 C& L; `9 H1 E+ N\" Z1 M
   
6.   y=0, while{y<=100,9 m( R  U- V# T9 L* J
   
7.     z=0, while{z<=100,* b7 z2 j* O) p
   
8.       s=s+f(x,y,z),9 l0 X' A2 I; P  t0 G, e1 Q( `& ~
   
9.       z++
   2 p2 j\" v& U0 n' i/ k$ k. W6 k
10.     },( O3 \. X( b# h% t3 ^# t
    
11.     y++( S7 a. j6 n. }& w# J
    
12.   },
    6 S! a6 u4 U1 |4 z
13.   x++# ~# @6 {\" u& v9 g, R
    
14. },1 e\" Z% n7 a. M3 D# w; p
    
15. s;; g! ]. r; |' u9 B! S7 l* ~7 P  \
    
16. [sys::clock()-t]/1000;

复制代码

结果：
- Q( b8 c6 H9 Z2 Z7 i: e; T1 d; y- p6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

2 b5 X& W+ N. Q* |. B, p5 K

) c0 v1 x# K6 W$ I# T
7 H) Y! G& r+ [$ S


+ [( q$ T: j7 ]8 Z; T. \, H
' ^# p2 g+ K6 |2 I* ?
) K. H9 H4 D- x5 x) b1 D7 P( Q% i

6 j- G. ~  D7 v
" z+ @/ {9 ^) {1 t9 p


( u. W4 B% u& \0 R( V5 `