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

forcal

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

5 c% c. p" H( h＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;/ a0 F8 A5 M: `3 \' j# Q/ h4 {
   
2. tic;! K  r- \5 h- p7 ~8 t9 c
   
3. s=0;
   7 B/ [, [1 F' o( p
4. for x=0:10000 i, j; q: z% C! S
   
5.   for y=0:1003 }' [6 q( Q3 A
   
6.     for z=0:100! K- Q# A+ M* X, ?* ?5 U
   
7.       s=s+f(x,y,z);- }* i* {# }! z3 E1 ]$ {5 M+ v% Q
   
8.     end
   + a% _: @. A2 r- H7 |/ {
9.   end, ^: S/ v+ @8 E4 o- ?) ]
   
10. end
    8 R( N' H! o: C7 F
11. s
    # g6 q/ S0 P) [3 R! I
12. toc
    9 F5 c3 k- q& Q: J/ t
13. 
    / p. s5 a( B9 k
14. s =
    / t5 W9 `; f( Z
15.     6.126720600000000e+009
    : D- M( M! t4 ?3 S% z9 N
16. 
    8 s, ^% o4 g/ P0 }( ^
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
" ~/ ^4 U  Y: K' |  ^9 k: A9 l+ B  ^4 em函数：

1. %file xyz.m4 d9 O. |+ y4 Q. n4 {
   
2. function c=xyz(x,y,z)2 O( y. J$ U8 K0 ~
   
3. c=x+y+z;
   2 H- I- M7 K. r4 v
4. end

复制代码

测试代码：

1. tic;# x& ]8 A; G( @
   
2. s=0;
   \" E! C) t/ `! N! J8 I7 S
3. for x=0:1000
   : Q  K; y+ S$ e- |  N/ _
4.   for y=0:1000 R8 t5 _$ z# p, l
   
5.     for z=0:100
   2 E. g, d8 D6 h. }) S- @. q
6.       s=s+xyz(x,y,z);4 d! e0 o. O3 \* r  O  e
   
7.     end
   ) ^& \0 D\" R7 P: l. @$ a
8.   end) V+ P\" o3 M! ]$ R
   
9. end8 D( C' k8 q3 c7 w0 }: }
   
10. s
    * j& e6 {! h. [  F' L' G. p* w
11. toc7 Q+ k\" m3 D5 d, B: E( y, Y* C! e5 l+ ~
    
12. 
    \" G2 \- Y6 [* T! ~
13. s =
    ' D# h2 U9 c  W! d\" ^
14.     6.126720600000000e+009
    6 m* K; e. A6 e8 ^- g) g
15. 6 ]) r  i! g! ~# `- w( Q
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

& d0 [, ^1 Y- w% {1 C( @2 ]Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;, y9 n0 p! D/ B5 n
   
2. mvar:0 g/ Y  J* e7 h  E/ |
   
3. t=sys::clock(),
   ! @: g5 w) ~! ~- _; ]
4. s=0,+ M$ H% z2 ~# r0 F* q2 U' z
   
5. x=0, while{x<=1000,
   9 [+ A  u0 z( [* ]* m
6.   y=0, while{y<=100,
   : F7 X9 R2 L& U/ z5 C- x9 O/ p
7.     z=0, while{z<=100,
   . ~\" q7 O7 d* t. x* O* g
8.       s=s+f(x,y,z),
   - @' y\" O) z6 v: ]  ~; [* f. s
9.       z++
   # ^8 r& l/ o/ x  z# @0 C0 m4 W! _7 U. d
10.     },  W7 J% ?- K- R/ D0 ^9 }' @
    
11.     y++2 b) X# z5 [  ~* r
    
12.   },
    \" `, Z' U; g% {, Q* \
13.   x++
    2 }( b  g; \( g$ v; v' d
14. },
    & u6 A+ m' y/ Q& E
15. s;7 @2 d8 E\" }\" m! }2 l0 r8 r
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
% Z) w2 w9 X1 a3 ?1.109  秒
, z5 l3 }& _5 p* s; V
7 r. j* G# Q! L: M6 Y二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

! {# k$ x1 d" k( f0 H7 B, B
9 \. h6 n% g- m0 s; ~4 p$ m


2 Z2 y/ _- e+ [! O2 r2 F7 A. S

, l3 _' R4 o( U3 l  ?, _, @
/ f( S8 A  O4 f* A& W. x
$ f3 k" w7 T% ^# U$ X' i4 F

/ ~6 S7 n! p$ e) q% V" z
# H2 k; N5 B- C. T: P" M

' _  a' r3 u" D5 R& n$ C

' f4 @7 ?7 t/ L- c5 s. V  ^
: @& n; S6 B  K, Q* \6 a$ Z  F  M