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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
) y' N" ]1 Q  u" R0 c# W7 Q
＝＝＝＝＝＝＝＝＝＝＝＝
' r0 s+ [8 H! P# \/ F
* F7 h; |7 T* T3 t' i4 M5 hMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;- O. K  S/ S5 u6 [, u' D9 l
   
2. tic;
   : X; ^. z) x\" ^  \, |
3. s=0;
   3 v+ `2 t4 b  J; x# Q
4. for x=0:1000
   9 B$ ?) E7 S& s9 f  \5 T* S
5.   for y=0:100; f8 _. m2 [0 y+ N
   
6.     for z=0:1003 r4 W/ O: L6 `- Z2 `
   
7.       s=s+f(x,y,z);
   4 z+ }( a) \: G; R/ D' A( M
8.     end
   4 Z. c: x4 B# j% W/ V
9.   end
   - z: K/ |3 O- E5 A6 F. @/ v) I
10. end( ]/ [, q, i5 |7 Z- ?
    
11. s
    8 T& c+ |# S) {+ F: v+ _2 [9 L
12. toc
    \" o$ H0 W7 l1 |; \6 o
13. 6 a; v( b- N) J2 e, E
    
14. s =
    : o1 C; q1 ?$ x( x\" L
15.     6.126720600000000e+009# @$ Q& h0 j& G, P5 g4 {7 r4 V
    
16. 4 `) Y! r7 y  ^, \4 y/ L
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m: H, n/ `( p5 p6 o% M
   
2. function c=xyz(x,y,z)
   8 _4 H% G( B% |
3. c=x+y+z;0 y* v; L0 h+ Q\" o$ R' {
   
4. end

复制代码

测试代码：

1. tic;
   * R! N* X* X/ a% Z
2. s=0;
   5 u% F4 C+ }\" ~* h2 Y
3. for x=0:1000
   ) n* Q: V  a; G: ]( f
4.   for y=0:1000 @2 t4 y\" }4 o4 v) N6 o
   
5.     for z=0:1007 t+ r# z$ t  q/ ?. B# I
   
6.       s=s+xyz(x,y,z);' y( w! ?, X) y3 C% d5 E! C
   
7.     end; }( C\" q+ Z, \( }# G
   
8.   end
   7 U* ~  [+ r- m% l) g+ K& Z
9. end
   # ?- D: m  w. W\" A& B8 \
10. s; f4 v\" V3 X9 ]. e
    
11. toc7 m5 |% ]7 N1 r; I
    
12. 7 i& ~6 h$ W6 `. d/ R3 l8 j, u
    
13. s =\" D. @: ~/ U5 Z0 U
    
14.     6.126720600000000e+009
    . I' v& K9 m) j& n
15. 0 _+ g4 ?\" ]# L) ?4 A6 q
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;( q$ h/ j+ P# r\" Y. @7 O' m9 M
   
2. mvar:
   * E  i& r1 n7 Z$ R- t4 Q
3. t=sys::clock(),
   4 G# U4 J6 @. t
4. s=0,
   # q6 E- R0 O2 H# s2 c2 h
5. x=0, while{x<=1000,  Y' S, N7 }* |- d: _
   
6.   y=0, while{y<=100,
   % p0 o+ m0 b5 n4 {
7.     z=0, while{z<=100,
   0 s, ]9 R  z, M- p1 y
8.       s=s+f(x,y,z),2 h( F+ I0 X$ J
   
9.       z++
   , P9 [4 e! O' w7 i- B- ~( c. t
10.     },
    6 s; [& F  p\" z' U  {
11.     y++
    8 T3 O\" A5 s0 Y- u
12.   },
    0 Y! S: y9 w; M2 _6 w  q4 o# m
13.   x++
    1 c  Z0 o+ Y4 E9 N3 k9 i8 b
14. },
    3 G( P: f3 u* v/ M: M* S
15. s;+ f- M' @6 w% ^8 N
    
16. [sys::clock()-t]/1000;

复制代码

结果：
7 G$ \3 a, ]# _$ P7 c4 |: G' j6126720600.
1.109  秒
! a( \5 g" `3 N/ f8 _" E
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
& b& [" E. t* a6 D& s. i

hxftntA

$ V# X* F' n( q$ `
) R2 r9 N5 ]& J' H1 v* |
$ ?9 T- g) f% d% C
) B/ D* e- |/ k' S* s6 p6 B4 T/ x, Z




( O& p3 P8 x! m( X- Y; \$ o



" j1 @4 ]# U1 \/ a

1 h* N( e; u( u
' ]# e/ ^( q7 a5 b4 Y& ?& a
# x8 B+ {3 }* B) i/ }* a- q