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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
6 X; u1 G' g; U  ]! f3 T
. T: b. V, U4 h$ G8 K* _7 O) M" r＝＝＝＝＝＝＝＝＝＝＝＝

% v9 ]; v% I3 N5 S& GMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;  @/ W- g) k  @( S5 {$ B' h9 F
   
2. tic;
   + z3 D& X$ G$ _2 A0 O* o
3. s=0;
   / E1 j& J( c9 t+ U! }
4. for x=0:10008 C, y7 F5 \4 O# V
   
5.   for y=0:1001 q9 B2 d% P' I: b# j5 Y' V
   
6.     for z=0:100; c1 e) \# n) G4 Y
   
7.       s=s+f(x,y,z);2 {9 `3 r1 f\" ~/ I' S. ]
   
8.     end! G' O. {2 I7 A) W$ J
   
9.   end: y  p  A! W; u6 N: P' n
   
10. end
    : ?5 r, _# C( [# y# @
11. s  C' y) D+ R8 W! h5 t
    
12. toc
    & k% e4 A; T4 c9 y+ q6 }! _
13. 
    , d2 t7 @! ~\" G& M7 I' Y
14. s =6 ~4 Q  z, f7 d' t- d8 c, a
    
15.     6.126720600000000e+009
    7 P2 |3 i& H. m. K/ `; c$ f
16. - ^  u% @8 r6 o4 j) Q0 j
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
7 R& ~" N. E5 R. xm函数：

1. %file xyz.m
   $ v. v# A/ I- v, O$ a8 y  A8 b
2. function c=xyz(x,y,z)9 i\" \2 V( ^# [# _5 o9 ?
   
3. c=x+y+z;
   / N7 z: l# h9 i7 l1 E$ x; F
4. end

复制代码

测试代码：

1. tic;3 j, d1 R  C\" Z- K. |
   
2. s=0;
   % F3 ~! z2 K- o% n) ^
3. for x=0:1000. Q) L: M: Q3 g  }+ v
   
4.   for y=0:1008 L4 z& K: e/ T* j' {
   
5.     for z=0:100\" z\" e2 \  v) g/ G8 V$ x
   
6.       s=s+xyz(x,y,z);6 i3 N, ]0 D: P
   
7.     end8 A\" |+ S( P( `: ]$ O# b
   
8.   end$ t9 S- d9 k( p6 N3 W( Z( l, x! `
   
9. end
     P0 v( u; A- f. S; }- U
10. s7 p9 C5 f7 Q4 e5 Y, g
    
11. toc/ @/ w& M- v0 l- Z; g
    
12. 
    . u* |. R0 R# S6 [2 x( B: C
13. s =' w- @0 j3 H3 N$ F  P. @4 m
    
14.     6.126720600000000e+009
    $ ]7 W$ g. g: x( V
15. 
    1 b  y0 {  C2 \  O; v! \8 W
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;! G' u% U$ V) t; \: D  v8 h
   
2. mvar:6 a1 H( }' b5 ^5 q0 j/ ]; F
   
3. t=sys::clock(),1 f6 n/ a- g) l& f/ l% v# K
   
4. s=0,
   # I% b7 U- ~$ y9 }/ Z' F
5. x=0, while{x<=1000,
   - ~) b- L4 a+ b/ C
6.   y=0, while{y<=100,
   ( \\" ?) T1 j; ]4 r
7.     z=0, while{z<=100,, r) \5 Y1 P5 {4 N$ x! ?
   
8.       s=s+f(x,y,z),  D5 {/ G) K) B* q
   
9.       z++
   8 x9 w9 [7 b( M2 r; G7 _
10.     },
    4 C1 Q9 r# B3 O5 ^/ x
11.     y++' S1 s7 f9 E. B7 N6 r+ j\" g
    
12.   },! Q9 b# ~+ z& G) p, m
    
13.   x++
    # ?6 _6 s- ~; J' {2 G3 e# {
14. },
    ' ]7 e) E2 h2 h! K* N
15. s;$ q' J  `, U; O$ g5 k1 O- U9 e
    
16. [sys::clock()-t]/1000;

复制代码

结果：
/ i& m9 ?  z1 a" B0 ^6126720600.
1.109  秒
: D. V" I  O  e2 ]# C
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

5 J" p- D3 u1 ]& b7 D' Q
5 [* [! n% d. z' B; \( n' J5 U+ H
3 \& v% Q- r' ?  D




% u, V: ~- Y+ q3 R9 F' T- }: J
/ I2 R7 {8 g% f% J6 N, B3 j
+ t3 e& i8 `5 y# F; w, I. ?! _
: [/ @; s3 E1 l( L, B+ w+ a+ D
7 R) e, X' o& H% y/ {2 S